Developing a FAIR digital ecosystem for DOE materials modeling
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
With leading energy experts, high-performance computing, and powerful analysis and modeling capabilities, the National Renewable Energy Laboratory (NREL) empowers aviation stakeholders with tools and data to project future energy demand and explore supply options. As a result, NREL can help electric utilities and airport operations identify aviation electrification opportunities and prepare for future electricity needs.
The satellite-driven carbon flux estimation has been playing important role to estimate continental-scale carbon budget. One of the biggest recent advances in the satellite-driven carbon flux modeling is utilization of high-frequent geostationary satellites to estimate diurnal cycle in carbon fluxes. The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. The new generation of geostationary satellite sensors provide frequent observations, often less than every 10 minutes. Here, we use GOES Advanced Baseline Imager (ABI) data to estimate hourly NEE over the conterminous US. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and solar radiation were derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Ameriflux data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.
The interaction between terrestrial carbon (C) and nitrogen (N) cycles has been incorporated into more and more land surface models. However, the scheme of C-N coupling differs greatly among models, and how these diverse representations of C-N interactions will affect C-cycle modeling remains unclear. In this study, we explored how the simulated ecosystem C storage capacity in the terrestrial ecosystem (TECO) model varied with three different commonly used schemes of C-N coupling. The three schemes (SM1, SM2, and SM3) have been used in three different coupled C-N models (i.e., TECO-CN, CLM 4.5, and O-CN, respectively). They differ mainly in the stoichiometry of C and N in vegetation and soils, plant N uptake strategies, down regulation of photosynthesis, and the pathways of N import. We incorporated the three C-N coupling schemes into the C-only version of the TECO model and evaluated their impacts on the C cycle with a traceability framework. Our results showed that all three of the C-N schemes caused significant reductions in steady-state C storage capacity compared with the C-only version with magnitudes of −23%, −30%, and −54% for SM1, SM2, and SM3, respectively. This reduced C storage capacity was mainly derived from the combined effects of decreases in net primary productivity (NPP; −29%, −15%, and −45%) and changes in mean C residence time (MRT; 9%, −17%, and −17%) for SM1, SM2, and SM3, respectively. The differences in NPP are mainly attributed to the different assumptions on plant N uptake, plant tissue C : N ratio, down regulation of photosynthesis, and biological N fixation. In comparison, the alternative representations of the plant vs. microbe competition strategy and the plant N uptake, combined with the flexible C : N ratio in vegetation and soils, led to a notable spread in MRT. These results highlight the fact that the diverse assumptions on N processes represented by different C-N coupled models could cause additional uncertainty for land surface models. Understanding their difference can help us improve the capability of models to predict future biogeochemical cycles of terrestrial ecosystems.
Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbonfluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle,we must also be able to understand and model its effects on dryland vegetation. Dynamic Global VegetationModels (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduceseasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to bothtotal productivity and competitive interactions with other plants, errors in how process-based models representphenology hinder our ability to predict climate change impacts. This may be particularly problematic in drylandecosystems where many species have developed a complex phenology in response to seasonal variability in bothmoisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure ofexisting phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) andgross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized modelparameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, wemodified the model to include a semi-deciduous phenology type and added flexibility to the representation ofgrass phenology. While optimizing parameters reduced model bias, the largest gains in model performance wereassociated with the development of our new representation of phenology. This modified model was able to bettercapture seasonal patterns of both leaf area index (R2=0.75) and gross primary productivity (R2=0.84), thoughits ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new modelalso resulted in a more realistic outcome of modeled competition between grass and shrubs. These findingsdemonstrate the importance of improving how DGVMs represent phenology in order to accurately forecastclimate change impacts in dryland ecosystems.
Primary goal: Provide resource managers with remote sensing products that support ecosystem forecasting models requiring salinity and inundation data. Work supports the habitat-switching modules in the Coastal Louisiana Ecosystem Assessment and Restoration (CLEAR) model, which provides scientific evaluation for restoration management (Visser et al., 2008). Ongoing work to validate flooding with radar (NWRC/USGS) and enhance persistence estimates through "fusion" of MODIS and Landsat time series (ROSES A.28 Gulf of Mexico). Additional work will also investigate relationship between saltwater dielectric constant and radar returns (Radarsat) (ROSES A.28 Gulf of Mexico).
An important potential use of ocean color chlorophyll data is to determine other important properties of the marine biosphere, such as primary productivity, new production, and particulate fluxes at spatial scales larger and temporal scales longer than those possible with ground-based observations. Such determinations will likely progress from relatively simple empirical correlations to algorithms that are actually predictive models of ecosystem dynamics. As an example, this paper demonstrates how an empirical correlation between nitrate concentration and new production can be understood by a simple productivity model. Several models are then constructed to examine the functional relationship between total production and surface chlorophyll. The empirical correlation is substantially different than the analogous relation in the model. Understanding the relationship between surface chlorophyll and productivity on a global scale will probably require families of models for various marine ecosystems.
Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbon fluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle, we must also be able to understand and model its effects on dryland vegetation. Dynamic Global Vegetation Models (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduce seasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to both total productivity and competitive interactions with other plants, errors in how process-based models represent phenology hinder our ability to predict climate change impacts. This may be particularly problematic in dryland ecosystems where many species have developed a complex phenology in response to seasonal variability in both moisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure of existing phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) and gross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized model parameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, we modified the model to include a semi-deciduous phenology type and added flexibility to the representation of grass phenology. While optimizing parameters reduced model bias, the largest gains in model performance were associated with the development of our new representation of phenology. This modified model was able to better capture seasonal patterns of both leaf area index (R(exp 2) = 0.75) and gross primary productivity (R(exp 2) = 0.84), though its ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new model also resulted in a more realistic outcome of modeled competition between grass and shrubs. These findings demonstrate the importance of improving how DGVMs represent phenology in order to accurately forecast climate change impacts in dryland ecosystems.
Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.
Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.
Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.
Lidar instruments provide scientists with the unique opportunity to characterize the 3D structure of forest ecosystems. This information allows us to estimate properties such as wood volume, biomass density, stocking density, canopy cover, and leaf area. Structural information also can be used as drivers for photosynthesis and ecosystem demography models to predict forest growth and carbon sequestration. All lidars use time-in-flight measurements to compute accurate ranging measurements; however, there is a wide range of instruments and data types that are currently available, and instrument technology continues to advance at a rapid pace. This seminar will present new technologies that are in use and under development at NASA for airborne and space-based missions. Opportunities for instrument and data fusion will also be discussed, as Dr. Cook is the PI for G-LiHT, Goddard's LiDAR, Hyperspectral, and Thermal airborne imager. Lastly, this talk will introduce radiative transfer models that can simulate interactions between laser light and forest canopies. Developing modeling capabilities is important for providing continuity between observations made with different lidars, and to assist the design of new instruments. Dr. Bruce Cook is a research scientist in NASA's Biospheric Sciences Laboratory at Goddard Space Flight Center, and has more than 25 years of experience conducting research on ecosystem processes, soil biogeochemistry, and exchange of carbon, water vapor and energy between the terrestrial biosphere and atmosphere. His research interests include the combined use of lidar, hyperspectral, and thermal data for characterizing ecosystem form and function. He is Deputy Project Scientist for the Landsat Data Continuity Mission (LDCM); Project Manager for NASA s Carbon Monitoring System (CMS) pilot project for local-scale forest biomass; and PI of Goddard's LiDAR, Hyperspectral, and Thermal (G-LiHT) airborne imager.
Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.
Carbon (C) allocation and nonstructural carbon (NSC) dynamics play essential roles in plant growth and survival under stress and disturbance. However, quantitative understanding of these processes remains limited. Here, in this work, we propose a framework where we connect commonly measured carbon cycle components (eddy covariance fluxes of canopy CO 2 exchange, soil CO 2 efflux, and allometry-based biomass and net primary production) by a simple mass balance model to derive ecosystem-level NSC dynamics (NSC i ), C translocation (dC i ), and the biomass production efficiency (BPE i ) in above- and belowground plant (i = agp and bgp) compartments. We applied this framework to two long-term monitored loblolly pine (Pinus taeda) plantations of different ages in North Carolina and characterized the variations of NSC and allocation in years under normal and drought conditions. The results indicated that the young stand did not have net NSC flux at the annual scale, whereas the mature stand stored a near-constant proportion of new assimilates as NSC every year under normal conditions, which was comparable in magnitude to new structural growth. Roots consumed NSC in drought and stored a significant amount of NSC post drought. The above- and belowground dC i and BPE i varied more from year to year in the young stand and approached a relatively stable pattern in the mature stand. The belowground BPE bgp differed the most between the young and mature stands and was most responsive to drought. With the internal C dynamics quantified, this framework may also improve biomass production estimation, which reveals the variations resulting from droughts. Overall, these quantified ecosystem-scale dynamics were consistent with existing evidence from tree-based manipulative experiments and measurements and demonstrated that combining the continuous fluxes as proposed here can provide additional information about plant internal C dynamics. Given that it is based on broadly available flux data, the proposed framework is promising to improve the allocation algorithms in ecosystem C cycle models and offers new insights into observed variability in soil–plant–climate interactions.
State of the art Engineering and Science codes have grown in complexity dramatically over the last two decades. As a consequence application teams have adopted more sophisticated development strategies, leveraging third party libraries, deploying comprehensive testing and using advanced debugging and profiling tools. In todays environment of diverse hardware platforms, these applications also desire performance portability - avoiding the need to duplicate work for various platforms - which makes it necessary that these tools and libraries also work across the various systems. The Kokkos EcoSystem provides that portable software stack. Based on the Kokkos Core Programming Model, the EcoSystem provides math libraries, interoperability capabilities with Python and Fortran, and Tools for analysing, debugging, and optimizing applications. In this paper we will provide an overview of the components, discuss some specific use cases, and highlight how co-designing these components enables a more developer friendly experience.
The third primary production algorithm round robin (PPARR3) compares output from 24 models that estimate depth-integrated primary production from satellite measurements of ocean color, as well as seven general circulation models (GCMs) coupled with ecosystem or biogeochemical models. Here we compare the global primary production fields corresponding to eight months of 1998 and 1999 as estimated from common input fields of photosynthetically-available radiation (PAR), sea-surface temperature (SST), mixed-layer depth, and chlorophyll concentration. We also quantify the sensitivity of the ocean-color-based models to perturbations in their input variables. The pair-wise correlation between ocean-color models was used to cluster them into groups or related output, which reflect the regions and environmental conditions under which they respond differently. The groups do not follow model complexity with regards to wavelength or depth dependence, though they are related to the manner in which temperature is used to parameterize photosynthesis. Global average PP varies by a factor of two between models. The models diverged the most for the Southern Ocean, SST under 10 C, and chlorophyll concentration exceeding 1mg Chlm-3. Based on the conditions under which the model results diverge most, we conclude that current ocean-color-based models are challenged by high-nutrient low-chlorophyll conditions, and extreme temperatures or chlorophyll concentrations. The GCM-based models predict comparable primary production to those based on ocean color: they estimate higher values in the Southern Ocean, at low SST, and in the equatorial band, while they estimate lower values in eutrophic regions (probably because the area of high chlorophyll concentrations is smaller in the GCMs). Further progress in primary production modeling requires improved understanding of the effect of temperature on photosynthesis and better parameterization of the maximum photosynthetic rate.
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.
Surface conditions and their impacts on surface mass balance (SMB) play important roles in ice sheet dynamics and ice sheet interactions with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) – an ecosystem of models and reanalyses – represent ice sheet and glacier SMB components, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. We find the successful representation of ice sheet surface mass balance in MERRA-2, and similar systems, is partly due to compensating biases in surface energy budget. In the atmospheric system, subsequent changes to cloud microphysics and in the longwave radiative transfer model are found to have reduced biases in surface net longwave radiation fluxes while exacerbating surface net shortwave radiation biases, producing erroneously high near-surface temperatures and surface melt in summer months. These issues can be exacerbated by poor experiment initialization and the use of two-moment cloud microphysics within the GEOS ocean-atmosphere coupled system. Here, we document the spatial and temporal extent and causes of these biases in the ice sheet surface energy budget across GEOS systems and implement a range of model improvements to mitigate these issues. We examine the combined impact on these changes on regional energy budget and SMB using both free running and replay experiments (to simulate the impact in reanalyses). Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet SMB in GMAO’s future reanalyses and forecasting systems.