Search

Societal Impact Statement Nitrogen, phosphorus, and potassium are essential elements for plant growth and are the primary nutrients in commercial fertilizers. However, the extent to which these nutrients individually limit plant growth is still unclear, as is the influence of their interactions. Using a meta‐analysis approach, plant growth was found to be co‐limited by nitrogen, phosphorus, and potassium. Furthermore, these findings demonstrate that multiple nutrient additions are primarily additive, with the exception of synergistic effects observed with NP additions on plant growth. These findings are important as they provide a timely and valuable reference for assessing the ecosystem‐scale consequences of increasing nutrient deposition caused by human activities. Summary Numerous field studies and a few meta‐analyses have confirmed that terrestrial primary production is individually or jointly limited by nitrogen (N) and phosphorus (P). Meanwhile, a recent synthesis emphasized that potassium (K) limitation in terrestrial ecosystems is more widespread than previously thought. However, it remains unclear whether N, P, and K co‐limit terrestrial production and how their interactions influence plant growth. Therefore, we conducted a globally comprehensive synthesis of a total 2877 observations to explore the individual and interactive effects of N, P, and K additions on plant growth. We found that N, P, K, NP, NK, PK, and NPK additions significantly increased aboveground biomass production (AGP) by 30.8%, 15.3%, 12.3%, 68.4%, 36.5%, 21.1%, and 78.7%, respectively. The combined addition of multiple nutrients resulted in a larger increase in AGP than a single addition of nutrients. These findings suggested that terrestrial production is widely co‐limited by N, P, and K. However, no significant responses of belowground biomass production (BGP) to nutrient additions were observed. Except for synergistic interactions between N and P additions, the responses of AGP to other pairs of nutrient additions were predominantly additive. This meta‐analysis highlights the importance of multiple‐nutrient limitation in regulating terrestrial primary production while pointing out that non‐additive effects (including synergistic and antagonistic) of multiple‐nutrient additions on terrestrial production are relatively rare.

Abstract Large datasets of carbon dioxide, energy, and water fluxes were measured with the eddy-covariance (EC) technique, such as FLUXNET2015. These datasets are widely used to validate remote-sensing products and benchmark models. One of the major challenges in utilizing EC-flux data is determining the spatial extent to which measurements taken at individual EC towers reflect model-grid or remote sensing pixels. To minimize the potential biases caused by the footprint-to-target area mismatch, it is important to use flux datasets with awareness of the footprint. This study analyze the spatial representativeness of global EC measurements based on the open-source FLUXNET2015 data, using the published flux footprint model (SAFE-f). The calculated annual cumulative footprint climatology (ACFC) was overlaid on land cover and vegetation index maps to create a spatial representativeness dataset of global flux towers. The dataset includes the following components: (1) the ACFC contour (ACFCC) data and areas representing 50%, 60%, 70%, and 80% ACFCC of each site, (2) the proportion of each land cover type weighted by the 80% ACFC (ACFCW), (3) the semivariogram calculated using Normalized Difference Vegetation Index (NDVI) considering the 80% ACFCW, and (4) the sensor location bias (SLB) between the 80% ACFCW and designated areas (e.g. 80% ACFCC and window sizes) proxied by NDVI. Finally, we conducted a comprehensive evaluation of the representativeness of each site from three aspects: (1) the underlying surface cover, (2) the semivariogram, and (3) the SLB between 80% ACFCW and 80% ACFCC, and categorized them into 3 levels. The goal of creating this dataset is to provide data quality guidance for international researchers to effectively utilize the FLUXNET2015 dataset in the future.

Accurate evapotranspiration (ET) estimation is crucial for understanding ecosystem dynamics and managing water resources. Existing methodologies, including traditional techniques like the Penman–Monteith model, remote sensing approaches utilizing Solar-Induced Fluorescence (SIF), and machine learning algorithms, have demonstrated varying levels of effectiveness in ET estimation. However, these methods often face significant challenges, such as reliance on empirical coefficients, inadequate representation of canopy dynamics, and limitations due to cloud cover and sensor constraints. These issues can lead to inaccuracies in capturing ET’s spatial and temporal variability, highlighting the need for improved estimation techniques. This study introduces a novel approach to enhance ET estimation by integrating SIF partitioning with Photosynthetically Active Radiation (PAR) and leaf area index (LAI) data, utilizing the TL-LUE model (Two-Leaf Light Use Efficiency). Partitioning SIF data into sunlit and shaded components allows for a more detailed representation of the canopy’s functional dynamics, significantly improving ET modelling. Our analysis reveals significant advancements in ET modelling through SIF partitioning. At Xiaotangshan Station, the correlation between modelled ET and SIFsu is 0.71, while the correlation between modelled ET and SIFsh is 0.65. The overall correlation (R2) between the modelled ET and the combined SIF partitioning (SIF(P)) is 0.69, indicating a strong positive relationship at Xiaotangshan Station. The correlations between SIFsh and SIFsu with modelled ET show notable patterns, with R2 values of 0.89 and 0.88 at Heihe Daman, respectively. These findings highlight the effectiveness of SIF partitioning in capturing canopy dynamics and its impact on ET estimation. Comparing modelled ET with observed ET and the Penman–Monteith model (PM model) demonstrates substantial improvements. R2 values for modelled ET against observed ET were 0.68, 0.76, and 0.88 across HuaiLai, Shangqiu, and Yunxiao Stations. Modelled ET correlations to the PM model were 0.75, 0.73, and 0.90, respectively, at three stations. These results underscore the model’s capability to enhance ET estimations by integrating physiological and remote sensing data. This innovative SIF-partitioning approach offers a more nuanced perspective on canopy photosynthesis, providing a more accurate and comprehensive method for understanding and managing ecosystem water dynamics across diverse environments.

Land surface models (LSMs) have prominent advantages for exploring the best agricultural practices in terms of both economic and environmental benefits with regard to different climate scenarios. However, their applications to optimizing fertilization and irrigation have not been well discussed because of their relatively underdeveloped crop modules. We used a CLM5-Crop LSM to optimize fertilization and irrigation schedules that follow actual agricultural practices for the cultivation of maize and wheat, as well as to explore the most economic and environmental-friendly inputs of nitrogen fertilizer and irrigation (FI), in the North China Plain (NCP), which is a typical intensive farming area. The model used the indicators of crop yield, farm gross margin (FGM), nitrogen use efficiency (NUE), water use efficiency (WUE), and soil nitrogen leaching. The results showed that the total optimal FI inputs of FGM were the highest (230 ± 75.8 kg N ha−1 and 20 ± 44.7 mm for maize; 137.5 ± 25 kg N ha−1 and 362.5 ± 47.9 mm for wheat), followed by the FIs of yield, NUE, WUE, and soil nitrogen leaching. After multi-objective optimization, the optimal FIs were 230 ± 75.8 kg N ha−1 and 20 ± 44.7 mm for maize, and 137.5 ± 25 kg N ha−1 and 387.5 ± 85.4 mm for wheat. By comparing our model-based diagnostic results with the actual inputs of FIs in the NCP, we found excessive usage of nitrogen fertilizer and irrigation during the current cultivation period of maize and wheat. The scientific collocation of fertilizer and water resources should be seriously considered for economic and environmental benefits. Overall, the optimized inputs of the FIs were in reasonable ranges, as postulated by previous studies. This result hints at the potential applications of LSMs for guiding sustainable agricultural development.

Evaluating the potential impacts of land use and land cover change (LULCC) and climate change on air pollution is crucial to unravel the driving forces and mechanisms behind changes in air quality. A multi-faceted approach was adopted, including a land change model (LCM) and Mann–Kendall (MK) test, to evaluate the transition of land cover type, changes in climate, and atmospheric pollutants during 2004–2021 in Pakistan. Moreover, a multiscale geographically weighted regression (MGWR) model and a mathematical model were used to assess the potential contribution of LULCC and climate dynamics to atmospheric pollution. It was revealed that during 2004, croplands covered an area of 9.72 × 104 mile2, accounting for 38% of the total area. However, the area of the croplands increased to 10.1 × 104 mile2, accounting for 40% of the total area in 2021. The MK test showed that the north and west–south regions significantly experienced air pollution, with the increasing trend for nitrogen dioxide (NO2) and sulfur dioxide (SO2) being 0.89× 1015 molecules/cm2 per year and 0.54 DU/year, respectively. For climate variability, mean precipitation (Precp) and mean surface pressure (SP) showed a prominent increasing trend, with a maximum value of 1 mm/year and 0.01 Kpa/year, respectively. The mean temperature maximum (Tmax) showed an increasing and decreasing trend, with the highest value of 0.28°C/year and 0.08°C/year, respectively. In the context of contribution, the conversion of cropland to grasslands increased the trend for SO2 concentrations. The highest increasing trend of 1.5 DU for ozone (O3) was found due to conversion of grasslands to shrublands. Additionally, regional climate played a significant role in making air pollution stagnant across the country. Precp and wind speed (WS) contributed significantly in escalating NO2 concentrations in Pakistan, while Precp contributed most (0.004 DU) to increasing SO2 concentrations. For O3, the most influential climate factor was Precp. These results on a long-term temporal scale demonstrated how maintaining climate variability through comprehensive land use management can help improve ambient air quality in Pakistan.

Droughts have an adverse influence on agriculture, the environment, water supplies, and the global economy. The drought risk was computed using an integrated prospective approach: drought hazard, exposure, and vulnerability based on biophysical and socio-economic conditions over Karachi, Pakistan during 2000–2019. Drought hazard map (DHM) was created using annual Palmer drought severity Index (PDSI). Drought exposure map (DEM) was derived using population density and gross domestic product (GDP), as well as land surface temperature (LST), Normal difference vegetation index (NDVI), Night light images (NTL), land use land cover (LULC), and Distance to water were used for drought vulnerability map (DVM). An estimation of drought Risk (EDR) was derived by integrating layers of DHM, DEM, and DVM. Results showed that Central, South, and East regions of Karachi were at high risk, whereas the North East and North were less affected by the drought. The estimated average drought hazard (EDH) was 0.84, with minimum (maximum) value of 0.68 (1). Similarly, the average estimated drought exposure (estimated drought vulnerability) for EDE (EDV) was 0.27 (0.42), with the maximum value of 0.55 (0.84) and the minimum value of 0 (0). The drought risk assessment map (DRAM) shows that the average risk values is 0.18 while highest value is 0.36.

Accurate global monitoring of carbon dioxide (CO2) is essential for understanding climate change and informing policy decisions. This study compares column-averaged dry-air mole fractions of CO2 (XCO2) between ACOS_L2_Lite_FP V9r for Japan’s Greenhouse Gases Observing Satellite (GOSAT), OCO-2_L2_Lite_FP V10r for the USA’s Orbiting Carbon Observatory-2 (OCO-2), and IAPCAS V2.0 for China’s Carbon Dioxide Observation Satellite (TANSAT) collectively referred to as GOT, with data from the Total Carbon Column Observing Network (TCCON). Our findings are as follows: (1) Significant data quantity differences exist between OCO-2 and the other satellites, with OCO-2 boasting a data volume 100 times greater. GOT shows the highest data volume between 30–45°N and 20–30°S, but data availability is notably lower near the equator. (2) XCO2 from GOT exhibits similar seasonal variations, with lower concentrations during June, July, and August (JJA) (402.72–403.74 ppm) and higher concentrations during December, January, and February (DJF) (405.74–407.14 ppm). XCO2 levels are higher in the Northern Hemisphere during March, April, and May (MAM) and DJF, while slightly lower during JJA and September, October, and November (SON). (3) The differences in XCO2 (ΔXCO2) reveal that ΔXCO2 between OCO-2 and TANSAT are minor (−0.47 ± 0.28 ppm), whereas the most significant difference is observed between GOSAT and TANSAT (−1.13 ± 0.15 ppm). Minimal differences are seen in SON (with the biggest difference between GOSAT and TANSAT: −0.84 ± 0.12 ppm), while notable differences occur in DJF (with the biggest difference between GOSAT and TANSAT: −1.43 ± 0.17 ppm). Regarding latitudinal variations, distinctions between OCO-2 and TANSAT are most pronounced in JJA and SON. (4) Compared to TCCON, XCO2 from GOT exhibits relatively high determination coefficients (R2 > 0.8), with GOSAT having the highest root mean square error (RMSE = 1.226 ppm, 2. Furthermore, it offers insights that can inform the analysis of differences in the inversion of carbon sources and sinks within assimilation systems when incorporating XCO2 data from satellite observations.

Abstract Global terrestrial vegetation dynamics have been rapidly altered by climate change. A widespread vegetation greenness over a large part of the planet from the 1980s to early this century has been reported, whereas weakening of CO2 fertilization effects and increasing climate extremes and the adverse impact of increasing rate of warming and severity of drought on vegetation growth were also reported. Earth system models project that the land carbon sink will decrease in size in response to an increase in warming during this century. How global vegetation is changing during this century in response to global warming and water availability across spatial and temporal scales remains uncertain. Our understanding of the widespread vegetation greening or browning processes and identifying the biogeochemical mechanisms remain incomplete. Here we use multiple long‐term satellite leaf area index (LAI) records to investigate vegetation growth trends from 1982 to 2018. We find that the widespread increase of growing‐season integrated LAI (greening) since 1980s was reversed (p‐value 45°N). Overall, the reversal of greening trend since 2000 weakened the negative feedback of carbon sequestration on the climatic system and should be considered in the strategies for climate warming mitigation and adaptation. Our findings of the diversity of processes that drive browning across bioclimatic‐zones and ecosystems and of how those driving processes are changing would enhance our ability to project global future vegetation change and its climatic and abiotic consequences.

To reduce air pollution, China adopted rigorous control mechanisms and announced the Air Pollution Prevention and Control Action Plan (APPCAP) in 2013. Here, using OMI satellite, the NASA Socioeconomic Data and Application Center (SEDAC), and Fifth ECMWF (ERA5) data at a 0.25° × 0.25° resolution, we explored changes in NO2, PM, SO2, and O3 and climatology over China in response to the Action Plan between 2004 and 2021. This study attempts to investigate the long term trend analysis of air pollution and climatic variations during two scenarios before (2004–2013) and after (2013–2021) APPCAP. We investigated the climatic effects of air pollution in China before and after APPCAP adoption using geographically weighted regression (GWR) and differential models to assess the contribution of air pollution. The spatial representation analysis demonstrated how air pollution affected climatic factors before and after the APPCAP. Several important findings were derived: (1) the APPCAP significantly influenced air pollution reduction in China post-scenario (2013–2021); (2) the Mann Kendall test investigated that all pollutants showed an increasing trend pre-APPCAP, while they showed a decreasing trend, except for O3, post-APPCAP; (3) for climatic factors, the MK test showed an increasing trend of precipitation and mean minimum air temperature tmin post-APPCAP; (4) innovative trend analysis (ITA) showed a reduction in NO2, SO2, and PM, although O3 showed no trend post-APPCAP; and (5) pre-scenario, NO2 contributed to an increase in the mean maximum air temperature (tmax) by 0.62 °C, PM contributed to raising tmin by 0.41 °C, while O3 reduced the tmax(tmin) by 0.15 °C (0.05 °C). PM increased tmax and precipitation with a magnitude 0.38 °C (7.38 mm), and NO2 contributed to increasing tmin by (0.35 °C), respectively, post-scenario. In particular, post-scenario led to an increase in tmin and precipitation across China. The results and discussion presented in this study can be beneficial for policymakers in China to establish long-term management plans for air pollution and climatological changes.

This study investigates the relationship of urban thermal environment (UTE) with various influential factors as well as ecological conditions. The relation between LST and land use and land cover (LULC) changes was explored in terms of remote-sensing (RS) based indices; heat effect contribution index (HECI), Urban thermal field variance index (UTFVI), Surface urban heat island intensity (SUHII), Normal Difference Built-up Index (NDBI), and Normal Difference Vegetation Index (NDVI). LULC maps were classified using the unsupervised classification technique and made error matrix to determine the accuracy. Results revealed that the vegetated area in Faisalabad decreased by 230 km2 due to an expansion in the urban area of 124-320 km2 during the period 1992-2014. An average LST in the rural buffers is increasing rapidly as compare to urban buffer and varied over the eight years with a range of 0.68-2.57 (°C). After 2007, SUHII's linear trend was negative because rural temperatures were still rising. Based on HECI, we found that urban expansion mainly led to increase in LST. UTFVI has shown poor ecological conditions in all urban buffers. In addition, there is a positive correlation between LST and NDBI, while NDVI indicates a negative correlation with LST.

Fine particulate matter (PM2.5) has a serious impact on human health. Forecasting PM2.5 levels and analyzing the pollution sources of PM2.5 are of great significance. In this study, the Lagrangian particle dispersion (LPD) model was developed by combining the FLEXPART model and the Bayesian inventory optimization method. The LPD model has the capacity for real-time forecasting and determination of pollution sources of PM2.5, which refers to the contribution ratio and spatial distribution of each type of pollution (industry, power, residential, and transportation). In this study, we applied the LPD model to the Beijing-Tianjin-Hebei (BTH) region to optimize the a priori PM2.5 emission inventory estimates during 15–20 March 2018. The results show that (1) the a priori estimates have a certain degree of overestimation compared with the a posteriori flux of PM2.5 for most areas of BTH; (2) after optimization, the correlation coefficient (R) between the forecasted and observed PM2.5 concentration increased by an average of approximately 10%, the root mean square error (RMSE) decreased by 30%, and the IOA (index of agreement) index increased by 16% at four observation sites (Aotizhongxin_Beijing, Beichenkejiyuanqu_Tianjin, Dahuoquan_Xintai, and Renmingongyuan_Zhangjiakou); and (3) the main sources of pollution at the four sites mainly originated from industrial and residential emissions, while power factory and transportation pollution accounted for only a small proportion. The concentration of PM2.5 forecasts and pollution sources in each type of analysis can be used as corresponding reference information for environmental governance and protection of public health.

Data from Landsat-8 and Sentinel-2A/2B are often combined for terrestrial monitoring because of their similar spectral bands. The bidirectional reflectance distribution function (BRDF) effect has been observed in both Landsat-8 and Sentinel-2A/2B reflectance data. However, there is currently no definition of solar zenith angle (θsz) that is suitable for the normalization of the BRDF-adjusted reflectance from the three sensors’ combined data. This paper describes the use of four machine learning (ML) models to predict a global θsz that is suitable for the normalization of bidirectional reflectance from the combined data in 2018. The observed θsz collected globally, and the three locations in the Democratic Republic of Congo (26.622°E, 0.356°N), Texas in the USA (99.406°W 30.751°N), and Finland (25.194°E, 61.653°N), are chosen to compare the performance of the ML models. At a global scale, the ML models of Support Vector Regression (SVR), Multi-Layer Perception (MLP), and Gaussian Process Regression (GPR) exhibit comparably good performance to that of polynomial regression, considering center latitude as the input to predict the global θsz. GPR achieves the best overall performance considering the center latitude and acquisition time as inputs, with a root mean square error (RMSE) of 1.390°, a mean absolute error (MAE) of 0.689°, and a coefficient of determination (R2) of 0.994. SVR shows an RMSE of 1.396°, an MAE of 0.638°, and an R2 of 0.994, following GPR. For a specific location, the SVR and GPR models have higher accuracy than the polynomial regression, with GPR exhibiting the best performance, when center latitude and acquisition time are considered as inputs. GPR is recommended for predicting the global θsz using the three sensors’ combined data.

Understanding the response of vegetation and ecosystem resilience to climate variability and drought conditions is essential for ecosystem planning and management. In this study, we assessed the vegetation changes and ecosystem resilience in the Horn of Africa (HOA) since 2000 and detected their drivers based mainly on analysis of the Moderate Resolution Imaging Spectroradiometer (MODIS) products. We found that the annual and seasonal trends of NDVI (Normalized Difference Vegetation Index) generally increased during the last two decades over the Horn of Africa particularly in western parts of Ethiopia and Kenya. The weakest annual and seasonal NDVI trends were observed over the grassland cover and tropical arid agroecological zones. The NDVI variation negatively correlated with Land Surface Temperature (LST) and positively correlated with precipitation at a significant level (p < 0.05) account for 683,197 km2 and 533,385 km2 area, respectively. The ecosystem Water Use Efficiency (eWUE) showed overall increasing trends with larger values for the grassland biome. The precipitation had the most significant effect on eWUE variation compared to LST and annual SPEI (Standardized Evapotranspiration Index). There were about 54.9% of HOA resilient to drought disturbance, whereas 32.6% was completely not-resilient. The ecosystems in the humid agroecological zones, the cropland, and wetland were slightly not-resilient to severe drought conditions in the region. This study provides useful information for policy makers regarding ecosystem and dryland management in the context of climate change at both national and regional levels.

Coastal wetlands provide essential ecosystem services and are closely related to human welfare. However, they can experience substantial degradation, especially in regions in which there is intense human activity. To control these increasingly severe problems and to develop corresponding management policies in coastal wetlands, it is critical to accurately map coastal wetlands. Although remote sensing is the most efficient way to monitor coastal wetlands at a regional scale, it traditionally involves a large amount of work, high cost, and low spatial resolution when mapping coastal wetlands at a large scale. In this study, we developed a workflow for rapidly mapping coastal wetlands at a 10 m spatial resolution, based on the recently emergent Google Earth Engine platform, using a machine learning algorithm, open-access Synthetic Aperture Radar (SAR) and optical images from the Sentinel satellites, and two terrain indices. We then generated a coastal wetland map of the Bohai Rim (BRCW10) based on the workflow. It has a producer accuracy of 82.7%, according to validation using 150 wetland samples. The BRCW10 data reflected finer information when compared to wetland maps derived from two sets of global high-spatial-resolution land cover data, due to the fusion of multiple data sources. The study highlights the benefits of simultaneously merging SAR and optical remote sensing images when mapping coastal wetlands.

The combination of Landsat-8, Landsat-9, Sentinel-2A and Sentinel-2B data provides a new perspective in remote sensing application for terrestrial monitoring. Jointly, these four sensors together offer global 10–30-m multi-spectral data coverage at a higher temporal revisit frequency. In this study, combinations of four sensors were used to examine the revisit interval by modelled orbit swath information. To investigate different factors that could influence data availability, an analysis was carried out for one year based on daytime surface observations of Landsat-8 and Sentinel-2A -2B. We found that (i) the global median average of revisit intervals for the combination of four sensors was 2.3 days; (ii) the global mean average number of surface observations was 141.4 for the combination of Landsat-8 and Sentinel-2A -2B; (iii) the global mean average cloud-weighted number of observations for the three sensors combined was 81.9. Three different locations were selected to compare with the cloud-weighted number of observations, and the results show an appropriate accuracy. The utility of combining four sensors together and the implication for terrestrial monitoring are discussed.

Coal mining, as one of the key drivers of land degradation worldwide, caused land subsidence problems. In this study, we conducted experimental research to explore the reclaimed mine soil (RMS) water dynamics and its sources in relation to reclaimed land use types using stable water isotopes in the Longdong coal mining area with high groundwater level in east China. We collected water samples seven times in 2017 from all of these water bodies (precipitation, surface waters (river water and water from subsidence pits (WSP)), groundwater and soil water). Our main findings are three fold: (1) the values of slope and intercept of the local meteoric water line of Craig (LMWL) of precipitation for the study area are higher than the global meteoric water line of Craig (GMWL) because of the humid monsoon climate zoon, and the values of δD and δ18O of surface waters and soil water and groundwater deviated from LMWL to some extent with a range of 5−30%, and the D and 18O of precipitation and the surface waters have higher seasonal variation than groundwater; (2) the values of δD and δ18O of RMS for the whole soil profile (0−100 cm) are lower than that of precipitation and have obvious seasonal variations and great fluctuation in the topsoil (0−30/40 cm) and decrease at depth (30/40−70 cm) and stable in deep soil layers (below 70 cm deep); (3) the RMS with forest and crop enhanced water infiltration capacity and soil water mixing strength compared with the waste RMS, so establishment of forest and crops should be encouraged in the RMS; (4) the main sources of topsoil (0−30 cm for crop and 0−40 cm for forest) of RMS are precipitation through infiltration, the main supply for deep soil water (below 70 cm deep) is groundwater, and the soil water for the middle deep soil layers (30/40−70 cm) is mainly from mixing sources of precipitation, groundwater, and river water through pant root water absorbing and groundwater upshifting.

Based on an improved risk assessment index system, this paper constructs an eco-environmental disaster risk assessment model during coalbed methane industrialization development in China by using intuitionistic fuzzy sets to describe the uncertain risk information, and the transformed interval value of Mamdani intuitionistic fuzzy neural networks model. Then, the validity of the model is verified by simulation tests. Furthermore, the assessment results are compared with those obtained by fuzzy neural networks model. Results show that suggested model has multidimensional nonlinearity and global approximation characteristics. By the procedure of “fuzzification - fuzzy-rules - defuzzifier”, the output conversion from uncertainty quantitative indicators to accurate risk assessment values can be effectively realized. Compared with the fuzzy neural networks model, the suggested model has better accuracy and stability. The risk assessment value calculated by the suggested model fairly matches the expected one. The study supplies a decision support for routine supervision and risk precaution and management on one hand, and enriches the theoretical research of the eco-environmental risk assessment of coalbed methane industrialization development on the other hand.

Global land cover is an important parameter of the land surface and has been derived by various researchers based on remote sensing images. Each land cover product has its own disadvantages and limitations. Data fusion technology is becoming a notable method to fully integrate existing land cover information. In this paper, we developed a method to generate a synergetic global land cover map (synGLC) based on Bayes theorem. A state probability vector was defined to precisely and quantitatively describe the land cover classification of every pixel and reduce the errors caused by legends harmonization and spatial resampling. Simple axiomatic approaches were used to generate the prior land cover map, in which pixels with high consistency were regarded to be correct and then used as benchmark to obtain posterior land cover map. Validation results show that our hybrid land cover map (synGLC, the dataset is available on request) has the best overall performance compared with the existing global land cover products. Closed shrub-lands and permanent wetlands have the highest uncertainty in our fused land cover map. This novel method can be extensively applied to fusion of land cover maps with different legends, spatial resolutions or geographic ranges.

Accurately modeling the land surface phenology based on satellite data is very important to the study of vegetation ecological dynamics and the related ecosystem process. In this study, we developed a Sigmoid curve (S-curve) function by integrating an asymmetric Gaussian function and a logistic function to fit the leaf area index (LAI) curve. We applied the resulting asymptotic lines and the curvature extrema to derive the vegetation phenophases of germination, green-up, maturity, senescence, defoliation and dormancy. The new proposed S-curve function has been tested in a specific area (Shangdong Province, China), characterized by a specific pattern in leaf area index (LAI) time course due to the dominant presence of crops. The function has not yet received any global testing. The identified phenophases were validated against measurement stations in Shandong Province. (i) From the site-scale comparison, we find that the detected phenophases using the S-curve (SC) algorithm are more consistent with the observations than using the logistic (LC) algorithm and the asymmetric Gaussian (AG) algorithm, especially for the germination and dormancy. The phenological recognition rates (PRRs) of the SC algorithm are obviously higher than those of two other algorithms. The S-curve function fits the LAI curve much better than the logistic function and asymmetric Gaussian function; (ii) The retrieval results of the SC algorithm are reliable and in close proximity to the green-up observed data whether using the AVHRR LAI or the improved MODIS LAI. Three inversion algorithms shows the retrieval results based on AVHRR LAI are all later than based on improved MODIS LAI. The bias statistics reveal that the retrieval results based on the AVHRR LAI datasets are more reasonable than based on the improved MODIS LAI datasets. Overall, the S-curve algorithm has the advantage of deriving vegetation phenophases across time and space as compared to the LC algorithm and the AG algorithm. With the SC algorithm, the vegetation phenophases can be extracted more effectively.

How well parameterization will improve gross primary production (GPP) estimation using the MODerate-resolution Imaging Spectroradiometer (MODIS) algorithm has been rarely investigated. We adjusted the parameters in the algorithm for 21 selected eddy-covariance flux towers which represented nine typical plant functional types (PFTs). We then compared these estimates of the MOD17A2 product, by the MODIS algorithm with default parameters in the Biome Property Look-Up Table, and by a two-leaf Farquhar model. The results indicate that optimizing the maximum light use efficiency (εmax) in the algorithm would improve GPP estimation, especially for deciduous vegetation, though it could not compensate the underestimation during summer caused by the one-leaf upscaling strategy. Adding the soil water factor to the algorithm would not significantly affect performance, but it could make the adjusted εmax more robust for sites with the same PFT and among different PFTs. Even with adjusted parameters, both one-leaf and two-leaf models would not capture seasonally photosynthetic dynamics, thereby we suggest that further improvement in GPP estimaiton is required by taking into consideration seasonal variations of the key parameters and variables.

Understanding how the dynamics of vegetation growth respond to climate change at different temporal and spatial scales is critical to projecting future ecosystem dynamics and the adaptation of ecosystems to global change. In this study, we investigated vegetated growth dynamics (annual productivity, seasonality and the minimum amount of vegetated cover) in China and their relations with climatic factors during 1982–2011, using the updated Global Inventory Modeling and Mapping Studies (GIMMS) third generation global satellite Advanced Very High Resolution Radiometer (AVHRR) Normalized Difference Vegetation Index (NDVI) dataset and climate data acquired from the National Centers for Environmental Prediction (NCEP). Major findings are as follows: (1) annual mean NDVI over China significantly increased by about 0.0006 per year from 1982 to 2011; (2) of the vegetated area in China, over 33% experienced a significant positive trend in vegetation growth, mostly located in central and southern China; about 21% experienced a significant positive trend in growth seasonality, most of which occurred in northern China (>35°N); (3) changes in vegetation growth dynamics were significantly correlated with air temperature and precipitation (p < 0.001) at a region scale; (4) at the country scale, changes in NDVI was significantly and positively correlated with annual air temperature (r = 0.52, p < 0.01) and not associated with annual precipitation (p > 0.1); (5) of the vegetated area, about 24% showed significant correlations between annual mean NDVI and air temperature (93% positive and remainder negative), and 12% showed significant correlations of annual mean NDVI with annual precipitation (65% positive and 35% negative). The spatiotemporal variations in vegetation growth dynamics were controlled primarily by temperature and secondly by precipitation. Vegetation growth was also affected by human activities; and (6) monthly NDVI was significantly correlated with the preceding month’s temperature and precipitation in western, central and northern China. The effects of a climate lag of more than two months in southern China may be caused mainly by the abundance of precipitation. These findings suggest that continuing efforts to monitor vegetation changes (in situ and satellite observations) over time and at broad scales are greatly needed, and are critical for the management of ecosystems and adapting to global climatic changes. It is likewise difficult to predict well future vegetation growth without linking these observations to mechanistic terrestrial ecosystem processes models that integrate all the satellite and in situ observations.

Global retrieval of solar-induced chlorophyll fluorescence (SIF) using remote sensing by means of satellites has been developed rapidly in recent years. Exploring how SIF could improve the characterization of photosynthesis and its role in the land surface carbon cycle has gradually become a very important and active area. However, compared with other gross primary production (GPP) models, the robustness of the parameterization of the SIF model under different circumstances has rarely been investigated. In this study, we examined and compared the effects of temporal aggregation and meteorological conditions on the stability of model parameters for the SIF model ( ε / S I F yield ), the one-leaf light-use efficiency (SL-LUE) model ( ε max ), and the two-leaf LUE (TL-LUE) model ( ε msu and ε msh ). The three models were parameterized based on a maize–wheat rotation eddy-covariance flux tower data in Yucheng, Shandong Province, China by using the Metropolis–Hasting algorithm. The results showed that the values of the ε / S I F yield and ε max were similarly robust and considerably more stable than ε msu and ε msh for all temporal aggregation levels. Under different meteorological conditions, all the parameters showed a certain degree of fluctuation and were most affected at the mid-day scale, followed by the monthly scale and finally at the daily scale. Nonetheless, the averaged coefficient of variation ( C V ) of ε / S I F yield was relatively small (15.0%) and was obviously lower than ε max ( C V = 27.0%), ε msu ( C V = 43.2%), and ε msh ( C V = 53.1%). Furthermore, the SIF model’s performance for estimating GPP was better than that of the SL-LUE model and was comparable to that of the TL-LUE model. This study indicates that, compared with the LUE-based models, the SIF-based model without climate-dependence is a good predictor of GPP and its parameter is more likely to converge for different temporal aggregation levels and under varying environmental restrictions in croplands. We suggest that more flux tower data should be used for further validation of parameter convergence in other vegetation types.

A fully coupled simulation of ecophysiological, hydrological and biochemical processes is significant for better understanding the individual and interactional impact of sophisticated land surface processes under future disturbances from nature and human beings. In this study, we spatially explicitly modelled evapotranspiration (ET) and photosynthesis (GPP) using a distributed hydrological model, Dynamic Land Model DLM-Ecohydro, over the Upper Heihe watershed for the years of 2013 and 2014. After considering the lateral water movements, the model fairly captured the variations in ET (R2 = 0.82, RMSE = 1.66 mm/day for 2013; R2 = 0.83, RMSE = 1.53 mm/day for 2014) and GPP (R2 = 0.71, RMSE = 5.25 gC/m2/day for 2013; R2 = 0.81, RMSE = 3.38 gC/m2/day for 2014) compared with the measurements from the Arou monitoring station. Vegetation transpiration accounted for total ET of around 65% and 64% in 2013 and 2014, respectively. A large spatial variability was found in these two indicators (14.30−885.36 mm/year for annual ET and 0−2174 gC/m2/day for annual GPP) over the watershed. Soil texture and vegetation functional types were the major factors affecting ET and GPP spatial variability, respectively. The study manifested a coupled water−carbon mechanism through the strong linear relationship between the variations in ET and GPP and the control of hydrological processes on the carbon cycle at the watershed scale. Although the model had a reasonable performance during most parts of the growing seasons, the lack of a soil freezing−thawing scenario caused inevitable discrepancies for the simulation of soil water and heat transfer mechanisms, hence inaccurately estimating the biophysiological processes in the transition period of winter to spring, which should be further improved especially for alpine regions.

The global carbon cycle research requires precise and sufficient observations of the column-averaged dry-air mole fraction of CO 2 (XCO 2 ) in addition to conventional surface mole fraction observations. In addition, assessing the consistency of multi-satellite data are crucial for joint utilization to better infer information about CO 2 sources and sinks. In this work, we evaluate the consistency of long-term XCO 2 retrievals from the Greenhouse Gases Observing Satellite (GOSAT), Orbiting Carbon Observatory 2 (OCO-2) in comparison with Total Carbon Column Observing Network (TCCON) and the 3D model of CO 2 mole fractions data from CarbonTracker 2017 (CT2017). We create a consistent joint dataset and compare it with the long-term model data to assess their abilities to characterize the carbon cycle climate. The results show that, although slight increasing differences are found between the GOSAT and TCCON XCO 2 in the northern temperate latitudes, the GOSAT and OCO-2 XCO 2 retrievals agree well in general, with a mean bias ± standard deviation of differences of 0.21 ± 1.3 ppm. The differences are almost within ±2 ppm and are independent of time, indicating that they are well calibrated. The differences between OCO-2 and CT2017 XCO 2 are much larger than those between GOSAT and CT XCO 2 , which can be attributed to the significantly different spatial representatives of OCO-2 and the CT-transport model 5 (TM5). The time series of the combined OCO-2/GOSAT dataset and the modeled XCO 2 agree well, and both can characterize significantly increasing atmospheric CO 2 under the impact of a large El Niño during 2015 and 2016. The trend calculated from the dataset using the seasonal Kendall (S-K) method indicates that atmospheric CO 2 is increasing by 2−2.6 ppm per year.

There is a growing concern over change in vegetation dynamics and drought patterns with the increasing climate variability and warming trends in Africa, particularly in the semiarid regions of East Africa. Here, several geospatial techniques and datasets were used to analyze the spatio-temporal vegetation dynamics in response to climate (precipitation and temperature) and drought in Eritrea from 2000 to 2017. A pixel-based trend analysis was performed, and a Pearson correlation coefficient was computed between vegetation indices and climate variables. In addition, vegetation condition index (VCI) and standard precipitation index (SPI) classifications were used to assess drought patterns in the country. The results demonstrated that there was a decreasing NDVI (Normalized Difference Vegetation Index) slope at both annual and seasonal time scales. In the study area, 57.1% of the pixels showed a decreasing annual NDVI trend, while the significance was higher in South-Western Eritrea. In most of the agro-ecological zones, the shrublands and croplands showed decreasing NDVI trends. About 87.16% of the study area had a positive correlation between growing season NDVI and precipitation (39.34%, p < 0.05). The Gash Barka region of the country showed the strongest and most significant correlations between NDVI and precipitation values. The specific drought assessments based on VCI and SPI summarized that Eritrea had been exposed to recurrent droughts of moderate to extreme conditions during the last 18 years. Based on the correlation analysis and drought patterns, this study confirms that low precipitation was mainly attributed to the slowly declining vegetation trends and increased drought conditions in the semi-arid region. Therefore, immediate action is needed to minimize the negative impact of climate variability and increasing aridity in vegetation and ecosystem services.

Exposure to fine particulate matter (PM2.5) is associated with adverse health impacts on the population. Satellite observations and machine learning algorithms have been applied to improve the accuracy of the prediction of PM2.5 concentrations. In this study, we developed a PM2.5 retrieval approach using machine-learning methods, based on aerosol products from the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the NASA Earth Observation System (EOS) Terra and Aqua polar-orbiting satellites, near-ground meteorological variables from the NASA Goddard Earth Observing System (GEOS), and ground-based PM2.5 observation data. Four models, which are orthogonal regression (OR), regression tree (Rpart), random forests (RF), and support vector machine (SVM), were tested and compared in the Beijing⁻Tianjin⁻Hebei (BTH) region of China in 2015. Aerosol products derived from the Terra and Aqua satellite sensors were also compared. The 10-repeat 5-fold cross-validation (10 × 5 CV) method was subsequently used to evaluate the performance of the different aerosol products and the four models. The results show that the performance of the Aqua dataset was better than that of the Terra dataset, and that the RF algorithm has the best predictive performance (Terra: R = 0.77, RMSE = 43.51 μg/m3; Aqua: R = 0.85, RMSE = 33.90 μg/m3). This study shows promise for predicting the spatiotemporal distribution of PM2.5 using the RF model and Aqua aerosol product with the assistance of PM2.5 site data.

High spatiotemporal resolution satellite imagery is useful for natural resource management and monitoring for land-use and land-cover change and ecosystem dynamics. However, acquisitions from a single satellite can be limited, due to trade-offs in either spatial or temporal resolution. The spatial and temporal adaptive reflectance fusion model (STARFM) and the enhanced STARFM (ESTARFM) were developed to produce new images with high spatial and high temporal resolution using images from multiple sources. Nonetheless, there were some shortcomings in these models, especially for the procedure of searching spectrally similar neighbor pixels in the models. In order to improve these models’ capacity and accuracy, we developed a modified version of ESTARFM (mESTARFM) and tested the performance of two approaches (ESTARFM and mESTARFM) in three study areas located in Canada and China at different time intervals. The results show that mESTARFM improved the accuracy of the simulated reflectance at fine resolution to some extent.

In this study, we evaluated estimates and predictions of the PM2.5 (fine particulate matter) concentrations and emissions in Xuzhou, China, using a coupled Lagrangian particle dispersion modeling system (FLEXPART-WRF). A Bayesian inversion method was used in FLEXPART-WRF to improve the emission calculation and mixing ratio estimation for PM2.5. We first examined the inversion modeling performance by comparing the model predictions with PM2.5 concentration observations from four stations in Xuzhou. The linear correlation analysis between the predicted PM2.5 concentrations and the observations shows that our inversion forecast system is much better than the system before calibration (with correlation coefficients of R = 0.639 vs. 0.459, respectively, and root mean square errors of RMSE = 7.407 vs. 9.805 µg/m3, respectively). We also estimated the monthly average emission flux in Xuzhou to be 4188.26 Mg/month, which is much higher (by ~10.12%) than the emission flux predicted by the multiscale emission inventory data (MEIC) (3803.5 Mg/month). In addition, the monthly average emission flux shows obvious seasonal variation, with the lowest PM2.5 flux in summer and the highest flux in winter. This pattern is mainly due to the additional heating fuels used in the cold season, resulting in many fine particulates in the atmosphere. Although the inversion and forecast results were improved to some extent, the inversion system can be improved further, e.g., by increasing the number of observation values and improving the accuracy of the a priori emission values. Further research and analysis are recommended to help improve the forecast precision of real-time PM2.5 concentrations and the corresponding monthly emission fluxes.

Coupled terrestrial carbon (C), nitrogen (N) and hydrological processes play a crucial role in the climate system, providing both positive and negative feedbacks to climate change. In this review we summarize published research results to gain an increased understanding of the dynamics between vegetation and atmosphere processes. A variety of methods, including monitoring (e.g., eddy covariance flux tower, remote sensing, etc.) and modeling (i.e., ecosystem, hydrology and atmospheric inversion modeling) the terrestrial carbon and water budgeting, are evaluated and compared. We highlight two major research areas where additional research could be focused: (i) Conceptually, the hydrological and biogeochemical processes are closely linked, however, the coupling processes between terrestrial C, N and hydrological processes are far from well understood; and (ii) there are significant uncertainties in estimates of the components of the C balance, especially at landscape and regional scales. To address these two questions, a synthetic research framework is needed which includes both bottom-up and top-down approaches integrating scalable (footprint and ecosystem) models and a spatially nested hierarchy of observations which include multispectral remote sensing, inventories, existing regional clusters of eddy-covariance flux towers and CO2 mixing ratio towers and chambers.

Background The global network of eddy-covariance (EC) flux-towers has improved the understanding of the terrestrial carbon (C) cycle, however, the network has a relatively limited spatial extent compared to forest inventory data and plots. Developing methods to use inventory-based and EC flux measurements together with modeling approaches is necessary evaluate forest C dynamics across broad spatial extents. Methods Changes in C stock change (ΔC) were computed based on repeated measurements of forest inventory plots and compared with separate measurements of cumulative net ecosystem productivity (ΣNEP) over four years (2003 – 2006) for Douglas-fir (Pseudotsuga menziesii var menziesii) dominated regeneration (HDF00), juvenile (HDF88 and HDF90) and near-rotation (DF49) aged stands (6, 18, 20, 57 years old in 2006, respectively) in coastal British Columbia. ΔC was determined from forest inventory plot data alone, and in a hybrid approach using inventory data along with litter fall data and published decay equations to determine the change in detrital pools. These ΔC-based estimates were then compared with ΣNEP measured at an eddy-covariance flux-tower (EC-flux) and modelled by the Carbon Budget Model - Canadian Forest Sector (CBM-CFS3) using historic forest inventory and forest disturbance data. Footprint analysis was used with remote sensing, soils and topography data to evaluate how well the inventory plots represented the range of stand conditions within the area of the flux-tower footprint and to spatially scale the plot data to the area of the EC-flux and model based estimates. Results The closest convergence among methods was for the juvenile stands while the largest divergences were for the regenerating clearcut, followed by the near-rotation stand. At the regenerating clearcut, footprint weighting of CBM-CFS3 ΣNEP increased convergence with EC flux ΣNEP, but not for ΔC. While spatial scaling and footprint weighting did not increase convergence for ΔC, they did provide confidence that the sample plots represented site conditions as measured by the EC tower. Conclusions Methods to use inventory and EC flux measurements together with modeling approaches are necessary to understand forest C dynamics across broad spatial extents. Each approach has advantages and limitations that need to be considered for investigations at varying spatial and temporal scales.

Climate change inevitably leads to changes in hydrothermal circulation. However, thermal-hydrologic exchanging caused by land cover change has also undergone ineligible changes. Therefore, studying the comprehensive effects of climate and land cover changes on land surface water and heat exchanges enables us to well understand the formation mechanism of regional climate and predict climate change with fewer uncertainties. This study investigated the land surface thermal-hydrologic exchange across southern China for the next 40 years using a land surface model (ecosystem-atmosphere simulation scheme (EASS)). Our findings are summarized as follows. (i) Spatiotemporal variation patterns of sensible heat flux (H) and evapotranspiration (ET) under the land cover scenarios (A2a or B2a) and climate change scenario (A1B) are unanimous. (ii) Both H and ET take on a single peak pattern, and the peak occurs in June or July. (iii) Based on the regional interannual variability analysis, H displays a downward trend (10%) and ET presents an increasing trend (15%). (iv) The annual average H and ET would, respectively, increase and decrease by about 10% when woodland converts to the cultivated land. Through this study, we recognize that land surface water and heat exchanges are affected greatly by the future climate change as well as land cover change.

The potential use of TerraSAR-X and Radarsat-2 data for soil moisture retrieval over bare agricultural areas was investigated using both empirical and semi-empirical approaches. For the empirical approach, the Support Vector Regression (SVR) model was used with two cases: (1) using only one C-band or X-band image; and (2) using a pair of C-band and X-band images jointly. For the semi-empirical approach, the modified Dubois model based on C-band and X-band SAR data was developed to estimate soil moisture content. The experiments were implemented over two bare agricultural areas, and in-situ measurements were carried out to assess the methods. The results showed that the TerraSAR-X and Radarsat-2 are suitable remote sensing tools for the estimation of surface soil moisture, with an accuracy of about 3 vol % (root mean square error, RMSE) over bare agricultural areas. Compared with the results obtained by Radarsat-2 data, TerraSAR-X data gives a slight improvement in estimating soil moisture. The accuracy of the soil moisture estimation was improved further when the two bands SAR data were used (RMSE of about 2.2 vol %) instead of only one. Moreover, the modified Dubois model showed comparable accuracy to the empirical model independent of the surface roughness.