Exploring Optimal Complexity for Water Stress Representation in Terrestrial Carbon Models: A Hybrid‐Machine Learning Model Approach
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A simple realistic biosphere model for calculating the transfer of energy, mass and momentum between the atmosphere and the vegetated surface of the earth has been developed for use in atmospheric general circulation models. The vegetation in each terrestrial model grid is represented by an upper level, representing the perennial canopy of trees and shrubs, and a lower level, representing the annual cover of grasses and other heraceous species. The vegetation morphology and the physical and physiological properties of the vegetation layers determine such properties as: the reflection, transmission, absorption and emission of direct and diffuse radiation; the infiltration, drainage, and storage of the residual rainfall in the soil; and the control over the stomatal functioning. The model, with prescribed vegetation parameters and soil interactive soil moisture, can be used for prediction of the atmospheric circulation and precipitaion fields for short periods of up to a few weeks.
The ability to accurately predict ecosystem drought response and recovery is necessary to produce reliable forecasts of land carbon uptake and future climate. Using a suite of models from the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), we assessed modeled net primary productivity (NPP) response to, and recovery from, drought events against a benchmark derived from tree ring observations between 1948 and 2008 across forested regions of the US and Europe. We find short lag times (0–6 months) between climate anomalies and modeled NPP response. Although models accurately simulate the direction of drought legacy effects (i.e. NPP decreases), projected effects are approximately four times shorter and four times weaker than observations suggest. This discrepancy between observed and simulated vegetation recovery from drought reveals a potential critical model deficiency. Since productivity is a crucial component of the land carbon balance, models that underestimate drought recovery time could overestimate predictions of future land carbon sink strength and, consequently, underestimate forecasts of atmospheric CO2.
A generalized carbon flow model for the major terrestrial ecosystems of the world is reported. The model is a simplification of the Century model and the Forest-Biogeochemical model. Topics covered include plant production, decomposition and nutrient cycling, biomes, the utility of the carbon flow model for predicting carbon dynamics under global change, and possible applications to state-and-transition models and environmentally driven global vegetation models.
The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.
In the northern hemisphere, terrestrial ecosystems transition from net sources of CO2 to the atmosphere in winter to net ecosystem carbon sinks during spring. The timing (or phase) of this transition, determined by the balance between ecosystem respiration (RECO) and primary production, is key to estimating the amplitude of the terrestrial carbon sink. We diagnose an apparent phase bias in the RECO and net ecosystem exchange (NEE) seasonal cycles estimated by the Terrestrial Carbon Flux (TCF) model framework and investigate its link to soil respiration mechanisms. Satellite observations of vegetation canopy conditions, surface meteorology, and soil moisture from the NASA SMAP Level 4 Soil Moisture product are used to model a daily carbon budget for a global network of eddy covariance flux towers. Proposed modifications to TCF include: the inhibition of foliar respiration in the light (the Kok effect); a seasonally varying litterfall phenology; an O2 diffusion limitation on heterotrophic respiration (RH); and a vertically resolved soil decomposition model. We find that RECO phase bias can result from bias in RECO magnitude and that mechanisms which reduce northern spring RECO, like substrate and O2 diffusion limitations, can mitigate the phase bias. A vertically resolved soil decomposition model mitigates this bias by temporally segmenting and lagging RH. Applying these model enhancements at Continuous Soil Respiration (COSORE) sites verifies their improvement of RECO and NEE skill compared to in situ observations (up to ∆RMSE = −0.76 g C m−2 d −1 35 ). Ultimately, these mechanisms can improve prior estimates of NEE for atmospheric inversion studies.
Terrestrial respiration is one of the most poorly understood processes in the global carbon cycle, making respiration predictions uncertain. However, expanding observations and machine learning approaches have led to a proliferation of estimates. We compiled total ecosystem and heterotrophic respiration estimates derived from top-down atmospheric inversions and bottom-up upscaling of ecosystem observations and compared them with dynamic vegetation models (DGVM) simulations over the 1980-2020 period. Our analysis revealed a convergence in mean annual global total ecosystem respiration estimates between top-down 97.1 (± SD 6.8) PgC yr-1 and bottom-up 98.5 (+/-13.4) PgC yr-1, which were both significantly lower than the ensemble mean from DGVMs estimates 133.7 (±4.7) PgC yr-1. We also found similar temporal trends between top-down estimates with a mean of 0.075 (±0.05) PgC yr-2, and bottom-up estimates of 0.05 (±0.05) PgC yr-2, which were 5 to 7 times smaller than the ensemble mean trend of 0.34 PgC yr-2 simulated by DGVMs. Global heterotrophic respiration showed much less agreement, ranging from top-down estimates of 42.7 (±4.0) PgC yr-1 to bottom-up estimates of 51.5 (±4.0) PgC yr-1 and a significantly larger ensemble model mean estimate of 60.8 PgC yr-1 (±1.9). The temporal trends in observation-based bottom-up estimates of heterotrophic respiration of 0.03 PgC yr-2 were five times lower than the model ensemble mean 0.15 PgC yr-2. Large regional disagreements in heterotrophic respiration estimates and simulations were evident in tropical and boreal latitudes. Therefore, improved regional and heterotrophic respiration estimates are necessary to reduce uncertainties regarding the future vulnerability of soil carbon.
Recent models for accretion of terrestrial bodies involve metal-silicate equilibrium as the metallic core formed during growth. Most elements considered are either refractory or well studied elements for which effects of pressure, temperature, oxygen fugacity, and metallic liquid composition are well known. There are a large number of elements that are both siderophile and volatile, whose fate in such models is unknown, largely due to a lack of data at comparable conditions and com-positions (FeNi core with light elements such as S, C, Si, and O). We have focused on Ge, In, As, Sb and determined the effect of Si and C on metal-silicate partitioning, and developed a thermo-dynamic model that allows application of these new data to a wide range of planetary bodies. New experiments: We have previously carried out experiments with FeSi metallic liquid at C-saturated conditions at 1600 and 1800 C [4]. In a new series of experiments we investigate the effect of Si in carbon-free systems at 1600 C for comparison. Experiments were carried out at 1 GPa in MgO capsules using the same basaltic starting composition as in previous studies. The MgO capsule reacts with the silicate melt to form more MgO-rich liquids that have 22-26 wt% MgO. Experimental met-als and silicates were analyzed using a combination of electron microprobe analysis and laser ablation ICP-MS. Results: The new results can be interpreted by considering Ge as an example, in the simple exchange equilibrium Fe + GeO = FeO + Ge, where the equilibrium constant Kd can be examined as a function of Si content of the metal. The slope of lnKd vs. (1-XSi) for this new series allows derivation of the epsilon interaction parameter for each of these four elements and Si (both C-saturated and C-free).All four elements have positive epsilon values, indicating that Si causes a decrease in the partition coefficients; values are 6.6, 6.5, 27.8 and 25.2 for In, Ge, As, and Sb, respectively, at 1 GPa and 1600 C. As an example of how large the effect of Si can be, these epsilon values correspond to activity coefficients (gamma) for As of 0.01 when XSi = 0, and up to gamma = 23 when XSi = 0.2. Combining these new results with previous determinations [5,6] of epsilon parameters for S and C for these elements allows us calculate activity of Ge, In, As, and Sb in Fe-Ni-Si-S-C-O metallic liquids. We apply this new model to sever-al terrestrial bodies such as Earth (Si-rich core), Mars (S-rich core), Moon (S-, C-, and Si-poor core), and Vesta, and examine the resulting core and mantle concentrations of these elements. Mantle concentrations of these four elements are well explained for Earth and Mars in models that call for mid-mantle equilibration between Si-bearing and S-bearing FeNi cores, respectively. Modeling results for the Moon and Vesta will also be presented.
Inland waters serve as important hydrological connections between the terrestrial landscape and oceans but are often overlooked in global carbon (C) budgets and Earth System Models. Terrestrially derived C entering inland waters from the watershed can be transported to oceans but over 83% is either buried in sediments or emitted to the atmosphere before reaching oceans. Anthropogenic pressures such as climate and landscape changes are altering the magnitude of these C fluxes in inland waters. Here, we synthesize the most recent estimates of C fluxes and the differential contributions across inland waterbody types (rivers, streams, lakes, reservoirs, and ponds), including recent measurements that incorporate improved sampling methods, small waterbodies, and dried areas. Across all inland waters, we report a global C emission estimate of 4.40 Pg C/year (95% confidence interval: 3.95–4.85 Pg C/year), representing a 13% increase from the most recent estimate. In this work, we also review the mechanisms by which the most globally widespread anthropogenically driven climate and landscape changes influence inland water C fluxes. The majority of these drivers are expected to influence terrestrial C inputs to inland waters due to alterations in terrestrial C quality and quantity, hydrological pathways, and biogeochemical processing. We recommend four research priorities for the future study of anthropogenic alterations to inland water C fluxes: (1) before-and-after measurements of C fluxes associated with climate change events and landscape changes, (2) better quantification of C input from land, (3) improved assessment of spatial coverage and contributions of small inland waterbodies to C fluxes, and (4) integration of dried and drawdown areas to global C flux estimates. Improved measurements of inland water C fluxes and quantification of uncertainty in these estimates will be vital to understanding both terrestrial C losses and the “moving target” of inland water C emissions in response to rapid and complex anthropogenic pressures.
Thermal evolution models based on subsolidus whole mantle convection which indicate that the surface heat flows of the earth and the moon do not necessarily provide good measures of the total amounts of radioactives in these bodies have been constructed. These models assume an initially hot state, but with a wide variety of choices for the parameters characterizing the rheology and convective vigor. All models are constrained to be consistent with present-day surface heat fluxes, and many of the terrestrial models are consistent with the mantle viscosities indicated by postglacial rebound. In the lunar models, heat generation is typically only 70-80% of the surface heat flow, even with allowance for the strong near-surface enhancement of radioactives. Despite the simplicity of these models, the persistence of a significant difference between heat generation and heat output indicates that this difference is real and should be incorporated in geochemical modeling of planets.
The First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) was conducted in a 15 km by 15 km research area located 8 km south of Manhattan, Kansas. The site consists primarily of native tallgrass prairie mixed with gallery oak forests and croplands. The objectives of FIFE are to better understand the role of biology in controlling the interactions between the land and the atmosphere, and to determine the value of remotely sensed data for estimating climatological parameters. The goals of FIFE are twofold: the upscale integration of models, and algorithm development for satellite remote sensing. The specific objectives of the field campaigns carried out in 1987 and 1989 were the simultaneous acquisition of satellite, atmospheric, and surface data; and the understanding of the processes controlling surface energy and mass exchange. Collected data were used to study the dynamics of various ecosystem processes (photosynthesis, evaporation and transpiration, autotrophic and heterotrophic respiration, etc.). Modelling terrestrial ecosystems at scales larger than that of a homogeneous plot led to the development of simple, generalized models of biogeochemical cycles that can be accurately applied to different biomes through the use of remotely sensed data. A model was developed called BIOME-BGC (for BioGeochemical Cycles) from a coniferous forest ecosystem model, FOREST-BGC, where a biome is considered a combination of a life forms in a specified climate. A predominately C4-photosynthetic grassland is probably the most different from a coniferous forest possible, hence the FIFE site was an excellent study area for testing BIOME-BGC. The transition from an essentially one-dimensional calculation to three-dimensional, landscape scale simulations requires the introduction of such factors as meteorology, climatology, and geomorphology. By using remotely sensed geographic information data for important model inputs, process-based ecosystem simulations at a variety of scales are possible. The second objective of this study is concerned with determining the accuracy of the estimated fluxes from BIOME-BGC, when extrapolated spatially over the entire 15-km by 15-km FIFE site. To accomplish this objective, a topographically distributed map of soil depth at the FIFE site was developed. These spatially-distributed fluxes were then tested with data from aircraft by eddy-flux correlation obtained during the FIFE experiment.
This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domainscale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospherichydrological process within a land surface model for use in HMA.
This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domain-scale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospheric-hydrological process within a land surface model for use in HMA.