Monitoring of carbon-water fluxes at Eurasian meteorological stations using random forest and remote sensing.

Simulating the carbon-water fluxes at more widely distributed meteorological stations based on the sparsely and unevenly distributed eddy covariance flux stations is needed to accurately understand the carbon-water cycle of terrestrial ecosystems. We established a new framework consisting of machine learning, determination coefficient (R2), Euclidean distance, and remote sensing (RS), to simulate the daily net ecosystem carbon dioxide exchange (NEE) and water flux (WF) of the Eurasian meteorological stations using a random forest model or/and RS. The daily NEE and WF datasets with RS-based information (NEE-RS and WF-RS) for 3774 and 4427 meteorological stations during 2002-2020 were produced, respectively. And the daily NEE and WF datasets without RS-based information (NEE-WRS and WF-WRS) for 4667 and 6763 meteorological stations during 1983-2018 were generated, respectively. For each meteorological station, the carbon-water fluxes meet accuracy requirements and have quasi-observational properties. These four carbon-water flux datasets have great potential to improve the assessments of the ecosystem carbon-water dynamics.

Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown.

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

Uncovering the critical soil moisture thresholds of plant water stress for European ecosystems.

Understanding the critical soil moisture (SM) threshold (θcrit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θcrit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θcrit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θcrit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θcrit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θcrit . We further found that soil texture dominates the spatial variability of θcrit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θcrit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models.

Altered energy partitioning across terrestrial ecosystems in the European drought year 2018.

Drought and heat events, such as the 2018 European drought, interact with the exchange of energy between the land surface and the atmosphere, potentially affecting albedo, sensible and latent heat fluxes, as well as CO2 exchange. Each of these quantities may aggravate or mitigate the drought, heat, their side effects on productivity, water scarcity and global warming. We used measurements of 56 eddy covariance sites across Europe to examine the response of fluxes to extreme drought prevailing most of the year 2018 and how the response differed across various ecosystem types (forests, grasslands, croplands and peatlands). Each component of the surface radiation and energy balance observed in 2018 was compared to available data per site during a reference period 2004-2017. Based on anomalies in precipitation and reference evapotranspiration, we classified 46 sites as drought affected. These received on average 9% more solar radiation and released 32% more sensible heat to the atmosphere compared to the mean of the reference period. In general, drought decreased net CO2 uptake by 17.8%, but did not significantly change net evapotranspiration. The response of these fluxes differed characteristically between ecosystems; in particular, the general increase in the evaporative index was strongest in peatlands and weakest in croplands. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data.

The FLUXNET2015 dataset provides ecosystem-scale data on CO2, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Climate change impact assessment on the hydrology of a large river basin in Ethiopia using a local-scale climate modelling approach.

Local-scale climate change adaptation is receiving more attention to reduce the adverse effects of climate change. The process of developing adaptation measures at local-scale (e.g., river basins) requires high-quality climate information with higher resolution. Climate projections are available at a coarser spatial resolution from Global Climate Models (GCMs) and require spatial downscaling and bias correction to drive hydrological models. We used the hybrid multiple linear regression and stochastic weather generator model (Statistical Down-Scaling Model, SDSM) to develop a location-based climate projection, equivalent to future station data, from GCMs. Meteorological data from 24 ground stations and the most accurate satellite and reanalysis products identified for the region, such as Climate Hazards Group InfraRed Precipitation with Station Data were used. The Soil Water Assessment Tool (SWAT) was used to assess the impacts of the projected climate on hydrology. Both SDSM and SWAT were calibrated and validated using the observed climate and streamflow data, respectively. Climate projection based on SDSM, in one of the large and agricultural intensive basins in Ethiopia (i.e., Awash), show high variability in precipitation but an increase in maximum (Tmax) and minimum (Tmin) temperature, which agrees with global warming. On average, the projection shows an increase in annual precipitation (>10%), Tmax (>0.4 °C), Tmin (>0.2 °C) and streamflow (>34%) in the 2020s (2011-2040), 2050s (2041-2070), and 2080s (2071-2100) under RCP2.6-RCP8.5. Although no significant trend in precipitation is found, streamflow during March-May and June-September is projected to increase throughout the 21 century by an average of more than 1.1% and 24%, respectively. However, streamflow is projected to decrease during January-February and October-November by more than 6%. Overall, considering the projected warming and changes in seasonal flow, local-scale adaptation measures to limit the impact on agriculture, water and energy sectors are required.

Long-term trends in rainfall and temperature using high-resolution climate datasets in East Africa.

Detecting changes in climate is a prerequisite for a better understanding of the climate and developing adaptation and mitigation measures at a regional and local scale. In this study long-term trends in rainfall and maximum and minimum temperature (T-max and T-min) were analysed on seasonal and annual time scales for East Africa. High resolution gridded rainfall (1981-2016) and temperature (1979-2010) data from international databases like the Climate Hazards Group are used. Long-term seasonal trend analysis shows a non-significant (except for small areas), decreasing (increasing) trend in rainfall in eastern (western) parts of Ethiopia and Kenya and a decreasing trend in large parts of Tanzania during the long rainy season. On the other hand, a non-significant increasing trend in large parts of the region is observed during the short rain season. With regard to annual trends, results largely confirm seasonal analyses: only a few significant trends in rainfall, but significant increasing trends in T-max (up to 1.9 °C) and T-min (up to 1.2 °C) for virtually the whole region. Our results demonstrate the need and added value of analysing climate trends based on data with high spatial resolution allowing sustainable adaptation measures at local scales.

Impacts of projected change in climate on water balance in basins of East Africa.

In East Africa, climate change and variability have shown a strong impact on sectors such as agriculture, energy, and water. To allow mitigation and adaptation of the possible impacts of the projected change in climate, this study applies a Statistical Downscaling Model (SDSM) to generate a high-resolution climate projection, equivalent to future station data, to drive impact assessment models in selected, agricultural intensive, basins of Ethiopia (EthShed), Kenya (KenShed), and Tanzania (TanShed). Observed and large-scale climate variables (predictors) are obtained from the national meteorological agency of Ethiopia and international databases. BROOK90, a physical-based hydrological model, is used to assess the impacts of the projected change in precipitation and maximum and minimum temperature (T-max, and T-min) on the water balance. Based on SDSM, the results show an increase in precipitation, relative to the baseline period (1961-1990), in EthShed (14% - 50%) and KenShed (15% - 86%) and a decrease in TanShed (1.3% - 6.3%) in the 20s (2011-2040), 50s (2041-2070), and 80s (2071-2100) under the three Representative Concentration Pathways (RCPs; RCP2.6, RCP4.5, and RCP8.5). T-max (anomalies up to 3.7 °C) and T-min (anomalies up to 2.76 °C) will be warmer than the baseline period throughout the 21 century in all three basins. In line with the projected change in precipitation and temperature, an increase (decrease) in seasonal and annual streamflow, soil-water, and evaporation in EthShed and KenShed (TanShed) is projected in the 20s, 50s, and 80s. In general, sustainable adaptation measures are required to be developed in a site-specific manner, considering the projected increase in temperature and evaporation in all three basins and a decrease in soil-water and streamflow in TanShed.

Statistically downscaled climate dataset for East Africa.

For many regions of the world, current climate change projections are only available at coarser spatial resolution from Global Climate Models (GCMs) that cannot directly be used in impact assessment and adaptation studies at regional and local scale. Impact assessment studies require high-resolution climate data to drive impact assessment models. To overcome this data challenge, we produced a station based climate projection (precipitation and maximum and minimum temperature) for Ethiopia, Kenya, and Tanzania using observed daily data from 211 stations obtained from the National Meteorological Agency of Ethiopia and international databases. Moreover, 26 large-scale climate variables derived from the National Centers for Environmental Prediction reanalysis data (1961-2005) and second generation Canadian Earth System Model (CanESM2, 1961-2100) are used. Statistical Down-Scaling Model (SDSM) is used to produce the required high-resolution climate projection by developing a statistical relationship between the large- and local-scale climate variables. The predictors are analysed more than 16458 times and we provided 20 ensembles for the current (1961-2005) and future (2006-2100, under RCP2.6, RCP4.5, and RCP8.5) climate.

The uncertain climate footprint of wetlands under human pressure.

Significant climate risks are associated with a positive carbon-temperature feedback in northern latitude carbon-rich ecosystems, making an accurate analysis of human impacts on the net greenhouse gas balance of wetlands a priority. Here, we provide a coherent assessment of the climate footprint of a network of wetland sites based on simultaneous and quasi-continuous ecosystem observations of CO2 and CH4 fluxes. Experimental areas are located both in natural and in managed wetlands and cover a wide range of climatic regions, ecosystem types, and management practices. Based on direct observations we predict that sustained CH4 emissions in natural ecosystems are in the long term (i.e., several centuries) typically offset by CO2 uptake, although with large spatiotemporal variability. Using a space-for-time analogy across ecological and climatic gradients, we represent the chronosequence from natural to managed conditions to quantify the "cost" of CH4 emissions for the benefit of net carbon sequestration. With a sustained pulse-response radiative forcing model, we found a significant increase in atmospheric forcing due to land management, in particular for wetland converted to cropland. Our results quantify the role of human activities on the climate footprint of northern wetlands and call for development of active mitigation strategies for managed wetlands and new guidelines of the Intergovernmental Panel on Climate Change (IPCC) accounting for both sustained CH4 emissions and cumulative CO2 exchange.

Above-ground woody carbon sequestration measured from tree rings is coherent with net ecosystem productivity at five eddy-covariance sites.

• Attempts to combine biometric and eddy-covariance (EC) quantifications of carbon allocation to different storage pools in forests have been inconsistent and variably successful in the past. • We assessed above-ground biomass changes at five long-term EC forest stations based on tree-ring width and wood density measurements, together with multiple allometric models. Measurements were validated with site-specific biomass estimates and compared with the sum of monthly CO2 fluxes between 1997 and 2009. • Biometric measurements and seasonal net ecosystem productivity (NEP) proved largely compatible and suggested that carbon sequestered between January and July is mainly used for volume increase, whereas that taken up between August and September supports a combination of cell wall thickening and storage. The inter-annual variability in above-ground woody carbon uptake was significantly linked with wood production at the sites, ranging between 110 and 370 g C m(-2) yr(-1) , thereby accounting for 10-25% of gross primary productivity (GPP), 15-32% of terrestrial ecosystem respiration (TER) and 25-80% of NEP. • The observed seasonal partitioning of carbon used to support different wood formation processes refines our knowledge on the dynamics and magnitude of carbon allocation in forests across the major European climatic zones. It may thus contribute, for example, to improved vegetation model parameterization and provides an enhanced framework to link tree-ring parameters with EC measurements.

Variability and distribution of spatial evapotranspiration in semi arid Inner Mongolian grasslands from 2002 to 2011.

Grasslands in Inner Mongolia are important for livestock farming while ecosystem functioning and water consumption are dominated by evapotranspiration (ET). In this paper we studied the spatiotemporal distribution and variability of ET and its components in Inner Mongolian grasslands over a period of 10 years, from 2002 to 2011. ET was modelled pixel-wise for more than 3000 1 km(2) pixels with the physically-based hydrological model BROOK90. The model was parameterised from eddy-covariance measurements and daily input was generated from MODIS leaf area index and surface temperatures. Modelled ET was also compared with the ET provided by the MODIS MOD16 ET data. The study showed ET to be highly variable in both time and space in Inner Mongolian grasslands. The mean coefficient of variation of 8-day ET in the study area varied between 25% and 40% and was up to 75% for individual pixels indicating a high innerannual variability of ET. Generally, ET equals or exceeds P during the vegetation period, but high precipitation in 2003 clearly exceeded ET in this year indicating a recharge of soil moisture and groundwater. Despite the high interannual and innerannual variations of spatial ET, the study also showed the existence of an intrinsic long-term spatial pattern of ET distribution, which can be explained partly by altitude and longitude (R(2) = 0.49). In conclusion, the results of this research suggest the development of dynamic and productive rangeland management systems according to the inherent variability of rainfall, productivity and ET in order to restore and protect Inner Mongolian grasslands.

Thermal optimality of net ecosystem exchange of carbon dioxide and underlying mechanisms.

• It is well established that individual organisms can acclimate and adapt to temperature to optimize their functioning. However, thermal optimization of ecosystems, as an assemblage of organisms, has not been examined at broad spatial and temporal scales. • Here, we compiled data from 169 globally distributed sites of eddy covariance and quantified the temperature response functions of net ecosystem exchange (NEE), an ecosystem-level property, to determine whether NEE shows thermal optimality and to explore the underlying mechanisms. • We found that the temperature response of NEE followed a peak curve, with the optimum temperature (corresponding to the maximum magnitude of NEE) being positively correlated with annual mean temperature over years and across sites. Shifts of the optimum temperature of NEE were mostly a result of temperature acclimation of gross primary productivity (upward shift of optimum temperature) rather than changes in the temperature sensitivity of ecosystem respiration. • Ecosystem-level thermal optimality is a newly revealed ecosystem property, presumably reflecting associated evolutionary adaptation of organisms within ecosystems, and has the potential to significantly regulate ecosystem-climate change feedbacks. The thermal optimality of NEE has implications for understanding fundamental properties of ecosystems in changing environments and benchmarking global models.

Surface characteristics of grasslands in Inner Mongolia as detected by micrometeorological measurements.

A roving tower concept was used to compare a semi-arid grassland site in Inner Mongolia (China), which was fenced in 1979 and ungrazed thereafter (UG79) with differently grazed semi-arid steppe ecosystems. The study was conducted during three consecutive years characterised by contrasting precipitation. The different grazing intensities included continuously and moderately grazed (CG), winter grazed (WG), and heavily grazed (HG). Here, we compare the energy fluxes and surface parameters that characterise the differently managed plots. The main focus is on sensible heat flux (H), available energy (AE), surface temperature (T ( s )), and surface albedo (alpha). Systematic errors were excluded by a side-to-side intercomparison of the instruments, and systematic climatic differences were minimised by the close distance between the fixed and the roving eddy covariance tower. Statistically, AE and T ( s ) were always significantly different between two simultaneously measured grazing intensities. Whereas AE was higher at UG79 in all years (mean difference of about 19Wm(-2)), T ( s ) was typically lower at UG79 (mean differences of 0.4 degrees C to about 2 degrees C). The exception was the end of the vegetation period in 2004 when T ( s ) was 0.6 degrees C higher at UG79 compared to CG. At UG79 alpha was typically significantly lower, and H was typically significantly higher. Consequently, latent heat fluxes (both as energy balance residual and directly measured) do not differ much between the different grazing intensities. It is concluded, that (1) the roving tower concept is able to detect differences due to grazing, (2) differences between the sites can be attributed to real surface differences, and (3) differences due to grazing intensities are small compared to interannual differences in surface fluxes.