skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 10 until 2:00 AM ET on Friday, October 11 due to maintenance. We apologize for the inconvenience.


Title: A differentiable, physics-informed ecosystem modeling and learning framework for large-scale inverse problems: demonstration with photosynthesis simulations

Abstract. Photosynthesis plays an important role in carbon,nitrogen, and water cycles. Ecosystem models for photosynthesis arecharacterized by many parameters that are obtained from limited in situmeasurements and applied to the same plant types. Previous site-by-sitecalibration approaches could not leverage big data and faced issues likeoverfitting or parameter non-uniqueness. Here we developed an end-to-endprogrammatically differentiable (meaning gradients of outputs to variablesused in the model can be obtained efficiently and accurately) version of thephotosynthesis process representation within the Functionally AssembledTerrestrial Ecosystem Simulator (FATES) model. As a genre ofphysics-informed machine learning (ML), differentiable models couplephysics-based formulations to neural networks (NNs) that learn parameterizations(and potentially processes) from observations, here photosynthesis rates. Wefirst demonstrated that the framework was able to correctly recover multiple assumedparameter values concurrently using synthetic training data. Then, using areal-world dataset consisting of many different plant functional types (PFTs), welearned parameters that performed substantially better and greatly reducedbiases compared to literature values. Further, the framework allowed us togain insights at a large scale. Our results showed that the carboxylationrate at 25 ∘C (Vc,max25) was more impactful than a factorrepresenting water limitation, although tuning both was helpful inaddressing biases with the default values. This framework could potentiallyenable substantial improvement in our capability to learn parameters andreduce biases for ecosystem modeling at large scales.

 
more » « less
Award ID(s):
2103942
NSF-PAR ID:
10527055
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
European Geophysical Union
Date Published:
Journal Name:
Biogeosciences
Volume:
20
Issue:
13
ISSN:
1726-4189
Page Range / eLocation ID:
2671 to 2692
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Incorporating microbial processes into soil biogeochemical models has received growing interest. However, determining the parameters that govern microbially driven biogeochemical processes typically requires case‐specific model calibration in various soil and ecosystem types. Here each case refers to an independent and individual experimental unit subjected to repeated measurements. Using the Microbial‐ENzyme Decomposition model, this study aimed to test whether a common set of microbially‐relevant parameters (i.e., generalized parameters) could be obtained across multiple cases based on a two‐year incubation experiment in which soil samples of four distinct soil series (i.e., Coland, Kesswick, Westmoreland, and Etowah) collected from forest and grassland were subjected to cellulose or no cellulose amendment. Results showed that a common set of parameters controlling microbial growth and maintenance as well as extracellular enzyme production and turnover could be generalized at the soil series level but not land cover type. This indicates that microbial model developments need to prioritize soil series type over plant functional types when implemented across various sites. This study also suggests that, in addition to heterotrophic respiration and microbial biomass data, extracellular enzyme data sets are needed to achieve reliable microbial‐relevant parameters for large‐scale soil model projections.

     
    more » « less
  2. Primary production and photoacclimation models are two important classes of physiological models that find applications in remote sensing of pools and fluxes of carbon associated with phytoplankton in the ocean. They are also key components of ecosystem models designed to study biogeochemical cycles in the ocean. So far, these two classes of models have evolved in parallel, somewhat independently of each other. Here we examine how they are coupled to each other through the intermediary of the photosynthesis–irradiance parameters. We extend the photoacclimation model to accommodate the spectral effects of light penetration in the ocean and the spectral sensitivity of the initial slope of the photosynthesis–irradiance curve, making the photoacclimation model fully compatible with spectrally resolved models of photosynthesis in the ocean. The photoacclimation model contains a parameterθ<#comment/>m, which is the maximum chlorophyll-to-carbon ratio that phytoplankton can attain when available light tends to zero. We explore how size-class-dependent values ofθ<#comment/>mcould be inferred from field data on chlorophyll and carbon content in phytoplankton, and show that the results are generally consistent with lower bounds estimated from satellite-based primary production calculations. This was accomplished using empirical models linking phytoplankton carbon and chlorophyll concentration, and the range of values obtained in culture measurements. We study the equivalence between different classes of primary production models at the functional level, and show that the availability of a chlorophyll-to-carbon ratio facilitates the translation between these classes. We discuss the importance of the better assignment of parameters in primary production models as an important avenue to reduce model uncertainties and to improve the usefulness of satellite-based primary production calculations in climate research.

     
    more » « less
  3. Abstract

    Vegetation water content (VWC) plays a key role in transpiration, plant mortality, and wildfire risk. Although land surface models now often contain plant hydraulics schemes, there are few direct VWC measurements to constrain these models at global scale. One proposed solution to this data gap is passive microwave remote sensing, which is sensitive to temporal changes in VWC. Here, we test that approach by using synthetic microwave observations to constrain VWC and surface soil moisture within the Climate Modeling Alliance Land model. We further investigate the possible utility of sub‐daily observations of VWC, which could be obtained through a satellite in geostationary orbit or combinations of multiple satellites. These high‐temporal‐resolution observations could allow for improved determination of ecosystem parameters, carbon and water fluxes, and subsurface hydraulics, relative to the currently available twice‐daily sun‐synchronous observational patterns. We find that incorporating observations at four different times in the diurnal cycle (such as could be available from two sun‐synchronous satellites) provides a significantly better constraint on water and carbon fluxes than twice‐daily observations do. For example, the root mean square error of projected evapotranspiration and gross primary productivity during drought periods was reduced by approximately 40%, when using four‐times‐daily relative to twice‐daily observations. Adding hourly observations of the entire diurnal cycle did not further improve the inferred parameters and fluxes. Our comparison of observational strategies may be informative in the design of future satellite missions to study plant hydraulics, as well as when using existing remotely sensed data to study vegetation water stress response.

     
    more » « less
  4. As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change. 
    more » « less
  5. Climate models predict that water limited regions around the world will become drier and warmer in the near future, including southwestern North America. We developed a large-scale experimental system that allows testing of the ecosystem impacts of precipitation changes. Four treatments were applied to 1600 m2 plots (40 m × 40 m), each with three replicates in a piñon pine (Pinus edulis) and juniper (Juniper monosperma) ecosystem. These species have extensive root systems, requiring large-scale manipulation to effectively alter soil water availability.  Treatments consisted of: 1) irrigation plots that receive supplemental water additions, 2) drought plots that receive 55% of ambient rainfall, 3) cover-control plots that receive ambient precipitation, but allow determination of treatment infrastructure artifacts, and 4) ambient control plots. Our drought structures effectively reduced soil water potential and volumetric water content compared to the ambient, cover-control, and water addition plots. Drought and cover control plots experienced an average increase in maximum soil and air temperature at ground level of 1-4° C during the growing season compared to ambient plots, and concurrent short-term diurnal increases in maximum air temperature were also observed directly above and below plastic structures. Our drought and irrigation treatments significantly influenced tree predawn water potential, sap-flow, and net photosynthesis, with drought treatment trees exhibiting significant decreases in physiological function compared to ambient and irrigated trees.  Supplemental irrigation resulted in a significant increase in both plant water potential and xylem sap-flow compared to trees in the other treatments. This experimental design effectively allows manipulation of plant water stress at the ecosystem scale, permits a wide range of drought conditions, and provides prolonged drought conditions comparable to historical droughts in the past – drought events for which wide-spread mortality in both these species was observed.  Obviously, one of the important areas of interest in this experiment was the effects of elevated (greater-than-average) and decreased (less-than-average) precipitation levels on soil moisture.  The volumetric water content of the soil was monitored across all twelve plots, all four treatment types, and all three cover types.  The record created through these monitoring activities not only noted the initial “wetting-up” of the soil after a precipitation event but also tracked the “drying-down” of the soil after the event.  The water content of the soil and its associated storage capacity could then provide a frame of reference in which changes in the physiological properties of our two target tree species, such as water potential and sapflow rate, could be interpreted.  
    more » « less