skip to main content

Title: Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling
Abstract. Evapotranspiration (ET) is critical in linking global water, carbon andenergy cycles. However, direct measurement of global terrestrial ET is notfeasible. Here, we first reviewed the basic theory and state-of-the-artapproaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surfacemodels (LSMs). We then utilized 4 remote-sensing-based physical models,2 machine-learning algorithms and 14 LSMs to analyze the spatial andtemporal variations in global terrestrial ET. The results showed that theensemble means of annual global terrestrial ET estimated by these threecategories of approaches agreed well, with values ranging from 589.6 mm yr−1(6.56×104 km3 yr−1) to 617.1 mm yr−1(6.87×104 km3 yr−1). For the period from 1982 to 2011, boththe ensembles of remote-sensing-based physical models and machine-learningalgorithms suggested increasing trends in global terrestrial ET (0.62 mm yr−2 with a significance level of p<0.05 and 0.38 mm yr−2 with a significance level of p<0.05,respectively). In contrast, the ensemble mean of the LSMs showed nostatistically significant change (0.23 mm yr−2, p>0.05),although many of the individual LSMs reproduced an increasing trend.Nevertheless, all 20 models used in this study showed that anthropogenicEarth greening had a positive role in increasing terrestrial ET. Theconcurrent small interannual variability, i.e., relative stability, found inall estimates of global terrestrial ET, suggests that a potentialplanetary boundary exists in regulating global terrestrial ET, with the value more » of this boundary beingaround 600 mm yr−1. Uncertainties among approaches were identified inspecific regions, particularly in the Amazon Basin and arid/semiaridregions. Improvements in parameterizing water stress and canopy dynamics,the utilization of new available satellite retrievals and deep-learning methods,and model–data fusion will advance our predictive understanding of globalterrestrial ET. « less
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
1903722 1243232
Publication Date:
NSF-PAR ID:
10167126
Journal Name:
Hydrology and Earth System Sciences
Volume:
24
Issue:
3
Page Range or eLocation-ID:
1485 to 1509
ISSN:
1607-7938
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. The quantity and quality of river discharge in Arctic regions is influenced by many processes including climate, watershed attributes and, increasingly, hydrological cycle intensification and permafrost thaw. We used a hydrological model to quantify baseline conditions and investigate the changing character of hydrological elements for Arctic watersheds between Utqiagvik (formerly known as Barrow)) and just west of Mackenzie River over the period 1981–2010. A synthesis of measurements and model simulations shows that the region exports 31.9 km3 yr−1 of freshwater via river discharge, with 55.5 % (17.7 km3 yr−1) coming collectively from the Colville, Kuparuk, and Sagavanirktok rivers. The simulations point to significant (p<0.05) increases (134 %–212 % of average) in cold season discharge (CSD) for several large North Slope rivers including the Colville and Kuparuk, and for the region as a whole. A significant increase in the proportion of subsurface runoff to total runoff is noted for the region and for 24 of the 42 study basins, with the change most prevalent across the northern foothills of the Brooks Range. Relatively large increases in simulated active-layer thickness (ALT) suggest a physical connection between warming climate, permafrost degradation, and increasing subsurface flow to streams and rivers. A decline in terrestrial water storage (TWS) is attributed tomore »losses in soil ice that outweigh gains in soil liquid water storage. Over the 30-year period, the timing of peak spring (freshet) discharge shifts earlier by 4.5 d, though the time trend is only marginally (p=0.1) significant. These changing characteristics of Arctic rivers have important implications for water, carbon, and nutrient cycling in coastal environments.« less
  2. Global aridification is projected to intensify. Yet, our knowledge of its potential impacts on species ranges remains limited. Here, we investigate global aridity velocity and its overlap with three sectors (natural protected areas, agricultural areas, and urban areas) and terrestrial biodiversity in historical (1979 through 2016) and future periods (2050 through 2099), with and without considering vegetation physiological response to rising CO2. Both agricultural and urban areas showed a mean drying velocity in history, although the concurrent global aridity velocity was on average +0.05/+0.20 km/yr−1(no CO2effects/with CO2effects; “+” denoting wetting). Moreover, in drylands, the shifts of vegetation greenness isolines were found to be significantly coupled with the tracks of aridity velocity. In the future, the aridity velocity in natural protected areas is projected to change from wetting to drying across RCP (representative concentration pathway) 2.6, RCP6.0, and RCP8.5 scenarios. When accounting for spatial distribution of terrestrial taxa (including plants, mammals, birds, and amphibians), the global aridity velocity would be -0.15/-0.02 km/yr−1(“-” denoting drying; historical), -0.12/-0.15 km/yr−1(RCP2.6), -0.36/-0.10 km/yr−1(RCP6.0), and -0.75/-0.29 km/yr−1(RCP8.5), with amphibians particularly negatively impacted. Under all scenarios, aridity velocity shows much higher multidirectionality than temperature velocity, which is mainly poleward. These results suggest that aridification risks may significantlymore »influence the distribution of terrestrial species besides warming impacts and further impact the effectiveness of current protected areas in future, especially under RCP8.5, which best matches historical CO2emissions [C. R. Schwalmet al.,Proc. Natl. Acad. Sci. U.S.A.117, 19656–19657 (2020)].

    « less
  3. Abstract The physical processes of heat exchange between lakes and the surrounding atmosphere are important in simulating and predicting terrestrial surface energy balance. Latent and sensible heat fluxes are the dominant physical process controlling ice growth and decay on the lake surface, as well as having influence on regional climate. While one-dimensional lake models have been used in simulating environmental changes in ice dynamics and water temperature, understanding the seasonal to daily cycles of lake surface energy balance and its relationship to lake thermal properties, atmospheric conditions, and how those are represented in models is still an open area of research. We evaluated a pair of one-dimensional lake models, Freshwater Lake (FLake) and the General Lake Model (GLM), to compare modeled latent and sensible heat fluxes against observed data collected by an eddy covariance tower during a 1-yr period in 2017, using Lake Mendota in Madison, Wisconsin, as our study site. We hypothesized transitional periods of ice cover as a leading source of model uncertainty, and we instead found that the models failed to simulate accurate values for large positive heat fluxes that occurred from late August into late December. Our results ultimately showed that one-dimensional models are effective inmore »simulating sensible heat fluxes but are considerably less sensitive to latent heat fluxes than the observed relationships of latent heat flux to environmental drivers. These results can be used to focus future improvement of these lake models especially if they are to be used for surface boundary conditions in regional numerical weather models. Significance Statement While lakes consist of a small amount of Earth’s surface, they have a large impact on local climate and weather. A large amount of energy is stored in lakes during the spring and summer, and then removed from lakes before winter. The effect is particularly noticeable in high latitudes, when the seasonal temperature difference is larger. Modeling this lake energy exchange is important for weather models and measuring this energy exchange is challenging. Here we compare modeled and observed energy exchange, and we show there are large amounts of energy exchange happening in the fall, which models struggle to capture well. During periods of partial ice coverage in early winter, lake behavior can change rapidly.« less
  4. Abstract. At the leaf level, stomata control the exchange of water and carbon across the air–leaf interface. Stomatal conductance is typically modeledempirically, based on environmental conditions at the leaf surface. Recently developed stomatal optimization models show great skills at predictingcarbon and water fluxes at both the leaf and tree levels. However, how well the optimization models perform atlarger scales has not been extensively evaluated. Furthermore, stomatal models are often used with simple single-leaf representations of canopy radiative transfer (RT), such asbig-leaf models. Nevertheless, the single-leaf canopy RT schemes do not have the capability to model optical properties of the leaves nor the entirecanopy. As a result, they are unable to directly link canopy optical properties with light distribution within the canopy to remote sensing dataobserved from afar. Here, we incorporated one optimization-based and two empirical stomatal models with a comprehensive RT model in the landcomponent of a new Earth system model within CliMA, the Climate Modelling Alliance. The model allowed us to simultaneously simulate carbon and waterfluxes as well as leaf and canopy reflectance and fluorescence spectra. We tested our model by comparing our modeled carbon and water fluxes andsolar-induced chlorophyll fluorescence (SIF) to two flux tower observations (a gymnospermmore »forest and an angiosperm forest) and satellite SIFretrievals, respectively. All three stomatal models quantitatively predicted the carbon and water fluxes for both forests. The optimization model,in particular, showed increased skill in predicting the water flux given the lower error (ca. 14.2 % and 21.8 % improvement for thegymnosperm and angiosperm forests, respectively) and better 1:1 comparison (slope increases from ca. 0.34 to 0.91 for the gymnosperm forest andfrom ca. 0.38 to 0.62 for the angiosperm forest). Our model also predicted the SIF yield, quantitatively reproducing seasonal cycles for bothforests. We found that using stomatal optimization with a comprehensive RT model showed high accuracy in simulating land surface processes. Theever-increasing number of regional and global datasets of terrestrial plants, such as leaf area index and chlorophyll contents, will helpparameterize the land model and improve future Earth system modeling in general.« less
  5. BACKGROUND The availability of nitrogen (N) to plants and microbes has a major influence on the structure and function of ecosystems. Because N is an essential component of plant proteins, low N availability constrains the growth of plants and herbivores. To increase N availability, humans apply large amounts of fertilizer to agricultural systems. Losses from these systems, combined with atmospheric deposition of fossil fuel combustion products, introduce copious quantities of reactive N into ecosystems. The negative consequences of these anthropogenic N inputs—such as ecosystem eutrophication and reductions in terrestrial and aquatic biodiversity—are well documented. Yet although N availability is increasing in many locations, reactive N inputs are not evenly distributed globally. Furthermore, experiments and theory also suggest that global change factors such as elevated atmospheric CO 2 , rising temperatures, and altered precipitation and disturbance regimes can reduce the availability of N to plants and microbes in many terrestrial ecosystems. This can occur through increases in biotic demand for N or reductions in its supply to organisms. Reductions in N availability can be observed via several metrics, including lowered nitrogen concentrations ([N]) and isotope ratios (δ 15 N) in plant tissue, reduced rates of N mineralization, and reduced terrestrial Nmore »export to aquatic systems. However, a comprehensive synthesis of N availability metrics, outside of experimental settings and capable of revealing large-scale trends, has not yet been carried out. ADVANCES A growing body of observations confirms that N availability is declining in many nonagricultural ecosystems worldwide. Studies have demonstrated declining wood δ 15 N in forests across the continental US, declining foliar [N] in European forests, declining foliar [N] and δ 15 N in North American grasslands, and declining [N] in pollen from the US and southern Canada. This evidence is consistent with observed global-scale declines in foliar δ 15 N and [N] since 1980. Long-term monitoring of soil-based N availability indicators in unmanipulated systems is rare. However, forest studies in the northeast US have demonstrated decades-long decreases in soil N cycling and N exports to air and water, even in the face of elevated atmospheric N deposition. Collectively, these studies suggest a sustained decline in N availability across a range of terrestrial ecosystems, dating at least as far back as the early 20th century. Elevated atmospheric CO 2 levels are likely a main driver of declines in N availability. Terrestrial plants are now uniformly exposed to ~50% more of this essential resource than they were just 150 years ago, and experimentally exposing plants to elevated CO 2 often reduces foliar [N] as well as plant-available soil N. In addition, globally-rising temperatures may raise soil N supply in some systems but may also increase N losses and lead to lower foliar [N]. Changes in other ecosystem drivers—such as local climate patterns, N deposition rates, and disturbance regimes—individually affect smaller areas but may have important cumulative effects on global N availability. OUTLOOK Given the importance of N to ecosystem functioning, a decline in available N is likely to have far-reaching consequences. Reduced N availability likely constrains the response of plants to elevated CO 2 and the ability of ecosystems to sequester carbon. Because herbivore growth and reproduction scale with protein intake, declining foliar [N] may be contributing to widely reported declines in insect populations and may be negatively affecting the growth of grazing livestock and herbivorous wild mammals. Spatial and temporal patterns in N availability are not yet fully understood, particularly outside of Europe and North America. Developments in remote sensing, accompanied by additional historical reconstructions of N availability from tree rings, herbarium specimens, and sediments, will show how N availability trajectories vary among ecosystems. Such assessment and monitoring efforts need to be complemented by further experimental and theoretical investigations into the causes of declining N availability, its implications for global carbon sequestration, and how its effects propagate through food webs. Responses will need to involve reducing N demand via lowering atmospheric CO 2 concentrations, and/or increasing N supply. Successfully mitigating and adapting to declining N availability will require a broader understanding that this phenomenon is occurring alongside the more widely recognized issue of anthropogenic eutrophication. Intercalibration of isotopic records from leaves, tree rings, and lake sediments suggests that N availability in many terrestrial ecosystems has steadily declined since the beginning of the industrial era. Reductions in N availability may affect many aspects of ecosystem functioning, including carbon sequestration and herbivore nutrition. Shaded areas indicate 80% prediction intervals; marker size is proportional to the number of measurements in each annual mean. Isotope data: (tree ring) K. K. McLauchlan et al. , Sci. Rep. 7 , 7856 (2017); (lake sediment) G. W. Holtgrieve et al. , Science 334 , 1545–1548 (2011); (foliar) J. M. Craine et al. , Nat. Ecol. Evol. 2 , 1735–1744 (2018)« less