Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ecosystem modeling”

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.

At least 865 records · Page 48

Global changes in biogeochemical cycles in response to human activities

The main objective of our research was to characterize biogeochemical cycles at continental and global scales in both terrestrial and aquatic ecosystems. This characterization applied to both natural ecosystems and those disturbed by human activity. The primary elements of interest were carbon and nitrogen and the analysis sought to quantify standing stocks and dynamic cycling processes. The translocation of major nutrients from the terrestrial landscape to the atmosphere (via trace gases) and to fluvial systems (via leaching, erosional losses, and point source pollution) were of particular importance to this study. Our aim was to develop the first generation of Earth System Models. Our research was organized around the construction and testing of component biogeochemical models which treated terrestrial ecosystem processes, aquatic nutrient transport through drainage basins, and trace gas exchanges at the continental and global scale. A suite of three complementary models were defined within this construct. The models were organized to operate at a 1/2 degree latitude by longitude level of spatial resolution and to execute at a monthly time step. This discretization afforded us the opportunity to understand the dynamics of the biosphere down to subregional scales, while simultaneously placing these dynamics into a global context.

Moore, Berrien, III↗

Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?

Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily–monthly) timescales and reduces on longer (seasonal–annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. Finally, the relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.

54 ENVIRONMENTAL SCIENCES↗

Carbon Dynamics in Vegetation and Soils

The overall goals of CD-08 team in Phase I were to quantify the contributions of different components of the carbon cycle to overall ecosystem carbon balance in Amazonian tropical forests and to undertake process studies at a number of sites along the eastern LBA transect to understand how and why these fluxes vary with site, season, and year. We divided this work into a number of specific tasks: (1) determining the average rate (and variability) of tree growth over the past 3 decades; (2) determining age demographics of tree populations, using radiocarbon to determine tree age; (3) assessing the rate of production and decomposition of dead wood debris; (4) determining turnover rates for organic matter in soils and the mean age of C respired from soil using radiocarbon measurements; and (5) comparing our results with models and constructing models to predict the potential of tropical forests to function as sources or sinks of C. This report summarizes the considerable progress made towards our original goals, which have led to increased understanding of the potential for central Amazon forests to act as sources or sinks of carbon with altered productivity. The overall picture of tropical forest C dynamics emerging from our Phase I studies suggests that the fraction of gross primary production allocated to growth in these forests is only 25-30%, as opposed to the 50% assumed by many ecosystem models. Consequent slow tree growth rates mean greater mean tree age for a given diameter, as reflected in our measurements and models of tree age. Radiocarbon measurements in leaf and root litter suggest that carbon stays in living tree biomass for several years up to a decade before being added to soils, where decomposition is rapid. The time lags predicted from 14C, when coupled with climate variation on similar time scales, can lead to significant interannual variation in net ecosystem C exchange.

Trumbore, Susan↗

Canopy reflectance models illustrate varying NDVI responses to change in high latitude ecosystems

Multiyear trends in Normalized Difference Vegetation Index (NDVI) have been used as metrics of high latitude ecosystem change based on the assumption that NDVI change is associated with ecological change, generally as changes in green vegetation amount (green leaf area index [LAI] or plant cover). Further, no change in NDVI is often interpreted as no change in these variables. Three canopy reflectance models including linear mixture model, the SAIL (Scattering from Arbitrarily Inclined Leaves) model, and the GeoSail model were used to simulate scenarios representing high latitude landscape NDVI responses to changes in LAI and plant cover. The simulations showed inconsistent NDVI responses. Clear increases in NDVI are generally associated with increases in LAI and plant cover. At higher values of LAI, the change in NDVI per unit change in LAI decreases, with very little change in spruce forest NDVI where crown cover is >50% and at the tundra–taiga ecotone with transitions from shrub tundra to spruce woodland. These lower responses may bias the interpretation of greening/browning trends in boreal forests. Variations in water or snow coverage were shown to produce outsized nonbiological NDVI responses. Inconsistencies in NDVI responses exemplify the need for care in the interpretation of NDVI change as a metric of high latitude ecosystem change, and that landscape characteristics in terms of the type of cover and its characteristics, such as the initial plant cover, must be taken into account in evaluating the significance of any observed NDVI trends.

Karl F Huemmrich↗

Imaging spectroscopy algorithms for mapping canopy foliar chemical and morphological traits and their uncertainties

A major goal of remote sensing is the development of generalizable algorithms to repeatedly and accurately map ecosystem properties across space and time. Imaging spectroscopy has great potential to map vegetation traits that cannot be retrieved from broadband spectral data, but rarely have such methods been tested across broad regions. Here we illustrate a general approach for estimating key foliar chemical and morphological traits through space and time using NASA's Airborne Visible/Infrared Imaging Spectrometer (AVIRIS-Classic). We apply partial least squares regression (PLSR) to data from 237 field plots within 51 images acquired between 2008 and 2011. Using a series of 500 randomized 50/50 subsets of the original data, we generated spatially explicit maps of seven traits (leaf mass per area (M area ), percentage nitrogen, carbon, fiber, lignin, and cellulose, and isotopic nitrogen concentration, δ 15 N) as well as pixel-wise uncertainties in their estimates based on error propagation in the analytical methods. Both Marea and %N PLSR models had a R 2 > 0.85. Root mean square errors (RMSEs) for both variables were less than 9% of the range of data. Fiber and lignin were predicted with R 2 > 0.65 and carbon and cellulose with R 2 > 0.45. Although R 2 of %C and cellulose were lower than Marea and %N, the measured variability of these constituents (especially %C) was also lower, and their RMSE values were beneath 12% of the range in overall variability. Model performance for δ 15 N was the lowest (R 2 = 0.48, RMSE = 0.95‰), but within 15% of the observed range. The resulting maps of chemical and morphological traits, together with their overall uncertainties, represent a first-of-its-kind approach for examining the spatiotemporal patterns of forest functioning and nutrient cycling across a broad range of temperate and sub-boreal ecosystems. These results offer an alternative to categorical maps of functional or physiognomic types by providing non-discrete maps (i.e., on a continuum) of traits that define those functional types. A key contribution of this work is the ability to assign retrieval uncertainties by pixel, a requirement to enable assimilation of these data products into ecosystem modeling frameworks to constrain carbon and nutrient cycling projections.

54 ENVIRONMENTAL SCIENCES↗

Functionally Assembled Terrestrial Ecosystem Simulator (FATES) for Hurricane Disturbance and Recovery

Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post-hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES) to study the processes and the key factors regulating post-hurricane forest recovery. We implemented hurricane-induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post-hurricane forest recovery, and quantified the sensitivity of key model parameters to the post-hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre-hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post-hurricane forest recovery in simulations. This research improves understanding of the ELM-FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities.

54 ENVIRONMENTAL SCIENCES↗

Using Remotely Sensed Data and Hydrologic Models to Evaluate the Effects of Climate Change on Shallow Aquatic Ecosystems in the Mobile Bay, AL Estuary

Coastal systems in the northern Gulf of Mexico, including the Mobile Bay, AL estuary, are subject to increasing pressure from a variety of activities including climate change. Climate changes have a direct effect on the discharge of rivers that drain into Mobile Bay and adjacent coastal water bodies. The outflows change water quality (temperature, salinity, and sediment concentrations) in the shallow aquatic areas and affect ecosystem functioning. Mobile Bay is a vital ecosystem that provides habitat for many species of fauna and flora. Historically, submerged aquatic vegetation (SAV) and seagrasses were found in this area of the northern Gulf of Mexico; however the extent of vegetation has significantly decreased over the last 60 years. The objectives of this research are to determine: how climate changes affect runoff and water quality in the estuary and how these changes will affect habitat suitability for SAV and seagrasses. Our approach is to use watershed and hydrodynamic modeling to evaluate the impact of climate change on shallow water aquatic ecosystems in Mobile Bay and adjacent coastal areas. Remotely sensed Landsat data were used for current land cover land use (LCLU) model input and the data provided by Intergovernmental Panel on Climate Change (IPCC) of the future changes in temperature, precipitation, and sea level rise were used to create the climate scenarios for the 2025 and 2050 model simulations. Project results are being shared with Gulf coast stakeholders through the Gulf of Mexico Data Atlas to benefit coastal policy and climate change adaptation strategies.

Estes, M. G.↗

Observations and modeling reveal that heatwaves reduce photosynthesis, plant carbon reserves, and net carbon uptake

Heatwaves threaten ecosystem carbon balances, yet the mechanisms driving short-term carbon flux responses remain poorly understood. Here, integrating high-frequency eddy covariance (EC) data from 140 global flux tower sites (872 site-years) with detailed process-based modeling, we examine ecosystem responses during and immediately after heatwaves. We show that heatwaves caused a −40% (range [−29%, −128%]) reduction in net ecosystem productivity (NEP) compared to pre-heatwave values, with this reduction persisting over the following two weeks (−38% range [+3%, −154%]). We attributed NEP decreases to photosynthesis decreases more than to ecosystem respiration (RE) increases. Forest sites had greater NEP decreases during heatwaves than non-forest sites, but remained carbon sinks afterwards, indicating resilience. Our modeling analysis of extreme heatwaves at selected EC sites shows that decreased photosynthesis, increased maintenance respiration, and decreased plant non-structural carbon reserves during heatwaves drive carbon cycle changes that persist for weeks. Consistent with phenocam observations, we modeled a reduction in leaf area index caused by reduced non-structural carbon reserves, leading to early leaf senescence and longer-term impacts. Ongoing increases in heatwaves are therefore likely to reduce NEP across a range of ecosystems, exacerbating carbon cycle feedback.

carbon cycle↗

Impacts of benchmarking choices on inferred model skill of the Arctic–Boreal terrestrial carbon cycle

Abstract Land surface models require continuous validation against observations to improve and reduce simulation uncertainty. However, inferred model performance can be heavily influenced by subjective choices made in the selection and application of observational data products. A key area often misrepresented by models is the Arctic–Boreal region, which is a potential tipping point region in Earth’s climate system due to large permafrost carbon stocks that are vulnerable to release with climate warming. We use the International Land Model Benchmarking (ILAMB) framework to evaluate how the model skill of TRENDY-v9 models varies based on the choice of observational-based benchmark and how benchmarks are applied in model evaluation. This analysis uses global datasets integrated into ILAMB and new, regionally-specific observational products from the Arctic–Boreal Vulnerability Experiment. Our results cover the overall time period of 1979–2019 and show that model scores can vary substantially depending on the data product applied, with higher model scores indicating better model performance against observations. The lowest model scores occur when benchmarked against regional, compared to global, datasets. We also evaluate observed and modeled functional relationships between ecosystem respiration and air temperature and between gross primary production and precipitation. Here, we find that the magnitude and shape of the responses are strongly impacted by the choice of observational dataset and the approach used to construct the functional relationship benchmark. These results suggest that model evaluation studies could conclude a false sense of model skill if only using a single benchmark data product or if not applying regional data products when performing a regional model analysis. Collectively, our findings highlight the influence of benchmarking choices on model evaluation and point to the need for benchmarking guidelines when assessing model skill.

Poe, Jeralyn (ORCID:0000000318495278)↗

Drivers of Phytoplankton Variability in and Near the Pearl River Estuary, South China Sea During Typhoon Hato (2017): A Numerical Study

The Pearl River estuary-coastal system in the Northern South China Sea is frequently affected by tropical cyclones (TCs) every year. Previous in-situ and remote sensing studies have found typhoons can enhance phytoplankton biomass and induce blooms in this region. However, the mechanistic links between phytoplankton blooms and typhoons have not been well elucidated due to the interplay of multiple processes along the land-ocean-atmosphere interfaces. Unravelling these interactions will have to rely on an integrated modelling system. In this work, we constructed a realistic, 3-dimensional, land-ocean-atmosphere modelling system with the marine ecosystem and sediment components for the China Great Bay Area. By using the integrated modelling system, we quantitatively investigated phytoplankton response to hydrological conditions variations under Typhoon Hato (2017), a strong typhoon case. Passive tracer experiments showed that with high river discharge induced by heavy rainfall, the residence time of Lingding Bay is as short as 15-day, less than half of that under the climatological discharge. The increase in freshwater pulse washed out the phytoplankton biomass within Lingding Bay. While for the offshore region, the source and sink terms analysis showed that the increase of phytoplankton biomass in the first week was because of the uplift of nutrient-rich subsurface water, while in the second week was because of the seaward propagated nearshore high phytoplankton biomass water. While riverine nutrients support phytoplankton growth in the third week, a large part of phytoplankton biomass was lost to zooplankton grazing, showing the system shifted from the bottom-up control to the top-down control.

54 ENVIRONMENTAL SCIENCES↗

A Multi-Scale Sampling Strategy for Detecting Physiologically Significant Signals in AVIRIS Imagery

Models of photosynthetic production at ecosystem and global scales require multiple input parameters specifying physical and physiological surface features. While certain physical parameters (e.g., absorbed photosynthetically active radiation) can be derived from current satellite sensors, other physiologically relevant measures (e.g., vegetation type, water status, carboxylation capacity, or photosynthetic light-use efficiency), are not generally directly available from current satellite sensors at the appropriate geographic scale. Consequently, many model parameters must be assumed or derived from independent sources, often at an inappropriate scale. An abundance of ecophysiological studies at the leaf and canopy scales suggests strong physiological control of vegetation-atmosphere CO2 and water vapor fluxes, particularly in evergreen vegetation subjected to diurnal or seasonal stresses. For example hot, dry conditions can lead to stomatal closure, and associated "downregulation" of photosynthetic biochemical processes, a phenomenon often manifested as a "midday photosynthetic depression". A recent study with the revised simple biosphere (SiB2) model demonstrated that photosynthetic downregulation can significantly impact global climate. However, at the global scale, the exact significance of downregulation remains unclear, largely because appropriate physiological measures are generally unavailable at this scale. Clearly, there is a need to develop reliable ways of extracting physiologically relevant information from remote sensing. Narrow-band spectrometers offer many opportunities for deriving physiological parameters needed for ecosystem and global scale photosynthetic models. Experimental studies on the ground at the leaf- to stand-scale have indicated that several narrow-band features can be used to detect plant physiological status. One physiological signal is caused by xanthophyll cycle pigment activity, and is often expressed as the Photochemical Reflectance Index (PRI). Because the xanthophyll cycle pigments are photoregulatory pigments closely linked to photosynthetic function, this index can be used to derive relative photosynthetic rates. An additional signal with physiological significance is the 970 nm water absorption band, which provides a measure of liquid water content. This feature has been quantified both using a simple 2-band ratio (900/970 nm, here referred to as the "Water Band Index" or WBI;), and using the "continuum removal" method. Current atmospheric correction methods for AVIRIS imagery also obtain quantitative expressions of surface liquid water absorption based on the 970 nm water band and may be comparable to ground-based estimates of water content using this feature. However, physiological interpretations of both the PRI and the WBI are best understood at the leaf and canopy scales, where complications of atmospheric interference and complex stand and landscape features can be minimized, and where experimental manipulations can be readily applied. Currently it is not known whether these physiological indices can be used to derive meaningful physiological information from AVIRIS imagery. In addition to the problem of atmospheric interference, another challenge is that any simple physiological index can be confounded by multiple factors unrelated to physiology, and this problem can become more severe at progressively larger spatial scales. For example, previous work has suggested that both the PRI and the WBI, are strongly correlated with other optical measures of canopy structure (e.g., the Normalized Difference Vegetation Index or green vegetation fraction), indicating a confounding effect of structure on physiological signals at the larger, landscape scale. Furthermore, the normal operating mode of most imaging spectrometers does not allow simultaneous, ground truthing at a level of detail needed for physiological sampling. Additionally, manipulative experiments of physiology are difficult to apply at a geographic scale suitable for comparison with remote imagery, which often works at spatial scales that are several orders of magnitude larger than those typically used for physiological studies. These limitations require the consideration of alternative approaches to validating physiological information derived from AVIRIS data. In this report, we present a multi-scale sampling approach to detecting physiologically significant signals in narrow-band spectra. This approach explores the multi-dimensional data space provided by narrow-band spectrometry, and combines AVIRIS imagery at a large scale, with ground spectral sampling at an intermediate scale, and detailed ecophysiological measurements at a fine scale, to examine seasonally and spatially changing relationships between multiple structural and physiological variables. Examples of this approach are provided by simultaneous sampling of the Normalized Difference Vegetation Index (NDVI), an index of fractional PAR interception and green vegetation cover, the Water Band Index (WBI, an index of liquid water absorption, and the Photochemical Reflectance Index (PRI, an index of xanthophyll cycle pigment activity and photosynthetic light-use efficiency. By directly linking changing optical properties sampled on the ground with measurable physiological states, we hope to develop a basis for interpreting similar signals in AVIRIS imagery.

Gamon, John A.↗

Allometric relationships and trade‐offs in 11 common M editerranean‐climate grasses

Abstract Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant–environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co‐occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process‐based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above‐ and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade‐offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast‐growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf‐to‐stem ratio. However, contrary to our hypothesis, there were no differences in root‐to‐leaf ratio between perennial and annual or C 4 and C 3 plants. Our study provides species‐specific and functional‐type‐specific allometry equations for both above‐ and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade‐offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

54 ENVIRONMENTAL SCIENCES↗

Evaluations of CMIP5 simulations over cropland

Cropland is one of major sources of carbon lost to the atmosphere and directly contributes to the emissions of greenhouse gases. There is, however, large potential for cropland to reduce its carbon flux to the atmosphere and sequester soil carbon through conservative agriculture management including no-tillage, perennial and/or deep root crops, irrigation, and organic fertilization etc. But these estimations on carbon emissions and sequestrations over cropland under future climate changes and variability remain largest uncertain among all other terrestrial biomes. Global climate and earth system models are an effective tool to study the cropland responses and feedbacks to present and future climate, yet most models in the latest couple model intercomparsion project phase 5 (CMIP5), generally treat cropland similarly as grassland with tuned parameters and do not account for realistic crop phenology, physiology, and management. In this study, we will evaluate the limitations and deficiencies of the CMIP5 models without process-based crop growth models over cropland by comparing their simulations against FLUXNET observations at eight cropland sites. The results show that: (1) the observed and simulated annual cycles generally are not consistent in either phase or amplitude; (2) the MPI and IPSL model families have better skills in the annual cycles of gross primary product (GPP), net ecosystem production (NEP), and terrestrial ecosystem respiration (TER) than other models at the corn/soybean and cereal sites respectively; (3) none of the CMIP5 models successfully simulate the observed two-peak pattern in the annual cycles of sensible heat fluxes at the corn/soybean sites; (4) the simulated GPPs and NEPs of the CESM1, BCC model families and NorESM1-M are much smaller than the observations for entire year; (5) model members from same model family normally simulate similar annual cycles both in phase and magnitude, but the model members from the CESM1 model family with different atmospheric models have different annual cycles; (6) the biases both in phases and magnitudes of annual cycles for the biogeochemical variables (GPP, NEP and TER) are generally larger than those biogeophysical variables (sensible and latent heat fluxes). Because of the above limitations and deficiencies of CMIP5 simulations over cropland, it is essential to incorporate process-based crop growth models into ESMs to improve the model physics and performance over cropland.

AGRICULTURE↗

Historical Carbon Dioxide Emissions Caused by Land-Use Changes are Possibly Larger than Assumed

The terrestrial biosphere absorbs about 20% of fossil-fuel CO2 emissions. The overall magnitude of this sink is constrained by the difference between emissions, the rate of increase in atmospheric CO2 concentrations, and the ocean sink. However, the land sink is actually composed of two largely counteracting fluxes that are poorly quantified: fluxes from land-use change andCO2 uptake by terrestrial ecosystems. Dynamic global vegetation model simulations suggest that CO2 emissions from land-use change have been substantially underestimated because processes such as tree harvesting and land clearing from shifting cultivation have not been considered. As the overall terrestrial sink is constrained, a larger net flux as a result of land-use change implies that terrestrial uptake of CO2 is also larger, and that terrestrial ecosystems might have greater potential to sequester carbon in the future. Consequently, reforestation projects and efforts to avoid further deforestation could represent important mitigation pathways, with co-benefits for biodiversity. It is unclear whether a larger land carbon sink can be reconciled with our current understanding of terrestrial carbon cycling. Our possible underestimation of the historical residual terrestrial carbon sink adds further uncertainty to our capacity to predict the future of terrestrial carbon uptake and losses.

Arneth, A.↗

Sensitivity of Global Terrestrial Gross Primary Production to Hydrologic States Simulated by the Community Land Model Using Two Runoff Parameterizations

The terrestrial water and carbon cycles interact strongly at various spatio-temporal scales. To elucidate how hydrologic processes may influence carbon cycle processes, differences in terrestrial carbon cycle simulations induced by structural differences in two runoff generation schemes were investigated using the Community Land Model 4 (CLM4). Simulations were performed with runoff generation using the default TOPMODEL-based and the Variable Infiltration Capacity (VIC) model approaches under the same experimental protocol. The comparisons showed that differences in the simulated gross primary production (GPP) are mainly attributed to differences in the simulated leaf area index (LAI) rather than soil moisture availability. More specifically, differences in runoff simulations can influence LAI through changes in soil moisture, soil temperature, and their seasonality that affect the onset of the growing season and the subsequent dynamic feedbacks between terrestrial water, energy, and carbon cycles. As a result of a relative difference of 36% in global mean total runoff between the two models and subsequent changes in soil moisture, soil temperature, and LAI, the simulated global mean GPP differs by 20.4%. However, the relative difference in the global mean net ecosystem exchange between the two models is small (2.1%) due to competing effects on total mean ecosystem respiration and other fluxes, although large regional differences can still be found. Our study highlights the significant interactions among the water, energy, and carbon cycles and the need for reducing uncertainty in the hydrologic parameterization of land surface models to better constrain carbon cycle modeling.

carbon↗

Seasonal and Interannual Patterns and Controls of Hydrological Fluxes in an Amazon Floodplain Lake With a Surface-Subsurface Process Model

Floodplain lakes represent important aquatic ecosystems, and field-based estimates of their water budgets are difficult to obtain, especially over multiple years. We examine the hydrological fluxes for an Amazon floodplain lake connected to the Solimões River using a process-based hydrologic model. Water exchanges between the river and lake agree well with field estimates, including the timing of different hydrological phases. However, beyond available field data, modeling results show that the seven simulated years all differed from each other. These interannual differences were caused by the interplay between phases when water levels were rising with river-water flowing into the lake (RWRI), versus rising with lake-water flowing out to the river (RWLO). This exchange determines the river-water content in the lake ( C L ). Maximum C L occurred before river levels peaked because local catchment contributions can be sufficient to push lake-water out to the river, even as river levels rise. Numerical experiments show that the seasonal distribution of local rainfall, local catchment size, and interannual variability in both climate and river stage can contribute to differing dynamics of C L in a floodplain lake. Their impacts vary among phases: river-rise dominates the RWRI, whereas local hydrological processes dominate the RWLO and receding-water phases. Intermediate-to-long-term rainfall accumulation controls C L during the RWLO phase, whereas annual precipitation accumulation is important for C L during low water. Finally, our model generalizes beyond limited available field studies and offers potential to better understand floodplain lakes in other areas and how regional versus local changes in climate may affect their hydrological dynamics.

54 ENVIRONMENTAL SCIENCES↗

Lower soil moisture and deep soil temperatures in thermokarst features increase old soil carbon loss after 10 years of experimental permafrost warming

Almost half of the global terrestrial soil carbon (C) is stored in the northern circumpolar permafrost region, where air temperatures are increasing two times faster than the global average. As climate warms, permafrost thaws and soil organic matter becomes vulnerable to greater microbial decomposition. Long-term soil warming of ice-rich permafrost can result in thermokarst formation that creates variability in environmental conditions. Consequently, plant and microbial proportional contributions to ecosystem respiration may change in response to long-term soil warming. Natural abundance δ 13 C and Δ 14 C of aboveground and belowground plant material, and of young and old soil respiration were used to inform a mixing model to partition the contribution of each source to ecosystem respiration fluxes. We employed a hierarchical Bayesian approach that incorporated gross primary productivity and environmental drivers to constrain source contributions. We found that long-term experimental permafrost warming introduced a soil hydrology component that interacted with temperature to affect old soil C respiration. Old soil C loss was suppressed in plots with warmer deep soil temperatures because they tended to be wetter. When soil volumetric water content significantly decreased in 2018 relative to 2016 and 2017, the dominant respiration sources shifted from plant aboveground and young soil respiration to old soil respiration. The proportion of ecosystem respiration from old soil C accounted for up to 39% of ecosystem respiration and represented a 30-fold increase compared to the wet-year average. Our findings show that thermokarst formation may act to moderate microbial decomposition of old soil C when soil is highly saturated. However, when soil moisture decreases, a higher proportion of old soil C is vulnerable to decomposition and can become a large flux to the atmosphere. Furthermore, as permafrost systems continue to change with climate, we must understand the thresholds that may propel these systems from a C sink to a source.

permafrost↗