Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Climate gradient”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Understanding and Predictability of Integrated Mountain Hydroclimate (Workshop Report)

Mountainous systems cover approximately 23% of Earth’s land and are distributed across all continents. They can capture and store atmospheric moisture that is then cycled through the terrestrial surface and subsurface system, released to downstream communities, and cycled back to the atmosphere. Mountain hydroclimate—characterized by steep gradients, geological, ecological, and biogeochemical diversity—is influenced by topographic forcing and elevated warming and susceptible to large subseasonal to multidecadal variability and rapid changes. Terrestrial hydrological and biogeochemical cycles also experience cascading effects from global warming impacts, such as multidecadal declines in mountain snowpack, longer growing seasons, and increased frequency and severity of extreme events like droughts and wildfires. However, little is known about the effects of these impacts and their feedbacks on climate systems and surface-subsurface compartments. Also unknown are the full implications of changing hydroclimate and extreme events on hydro biogeochemical cycles across atmosphere, terrestrial, and human systems in mountain regions and beyond. This knowledge gap is critical, given human reliance on mountain systems for stable water supply and quality. Mountain systems’ increasing vulnerability to climate change and human perturbations motivates the need to improve understanding of integrated mountain hydroclimate (IMHC) systems and their feedbacks and impacts on humans across scales. However, due to large heterogeneity and strong gradients, coupled natural-human processes in mountain regions present significant challenges for observations, modeling, predictions, and projections. Motivated by gaps in mountain hydroclimate understanding, observations, and modeling and the need for credible projections of future changes, the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program organized a virtual workshop on “Understanding and Predictability of Integrated Mountain Hydroclimate.” Sponsored by BER’s Earth and Environmental Systems Sciences Division (EESSD), the workshop aimed to inform and catalyze EESSD’s growing interests in enhancing predictive understanding of IMHC. Organizers structured the workshop to identify (1) knowledge gaps, (2) observational and modeling challenges, (3) short-term (1 to 3 years) and long-term (10 years and beyond) research opportunities, and (4) strategies for fostering collaboration and coordination. To address the outstanding challenges of IMHC, the workshop included two sessions organized by disciplinary, cross-disciplinary, and crosscutting science topics. The disciplinary and cross-disciplinary topics focused on essential IMHC elements: atmosphere, terrestrial, and human systems and their interactions. Breakout sessions on disciplinary and cross-disciplinary topics facilitated identification of crosscutting topics and central emerging themes. Session 1 focused on connecting existing DOE investments to accelerate progress related to scientific challenges in understanding mountain hydroclimate. In Session 2, participants further explored key Session 1 takeaways through the lens of multiagency collaborations and coordination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2024

Poleward Migration of the Latitude of Maximum Tropical Cyclone Intensity—Forced or Natural?

Abstract Past studies have shown a significant observed poleward trend in the latitude at which tropical cyclones reach their lifetime maximum intensity (LMI), especially in the northwest Pacific basin. Given the brevity of the historical record, it remains difficult to separate the forced trend from internal variability of the climate system. A recently developed tropical cyclone downscaling model is used to downscale the Community Earth System Model, version 2 (CESM2), preindustrial control simulation. It is found that the observed trend in the latitude at which tropical cyclones reach their LMI in the northwest Pacific is very unlikely to be caused by internal variability. The same downscaling model is then used to downscale CESM2 simulations under historical forcing. The resulting trend distribution shows a significant poleward migration of tropical cyclone LMI even after regressing out both natural variability and the part of the forced warming pattern that projects onto natural variability. The results indicate that the observed poleward migration of the latitude at which tropical cyclones reach their LMI in the northwest Pacific basin is likely to be, at least in part, forced. However, the magnitude of the projected poleward trend in climate models can be significantly modulated by the simulated spatial pattern of ocean warming. This highlights how discrepancies between models and observations, with regard to projected changes to the equatorial zonal sea surface temperature gradient under anthropogenic forcing, can lead to large uncertainties in projected changes to the LMI latitude of tropical cyclones. Significance Statement Observations in the northwest Pacific basin show that the latitude at which tropical cyclones are at their most intense has been trending northward in the recent half century. These changes are important since tropical cyclones could bring hazardous weather to coastal areas that are poorly equipped to handle them. Here, we show that natural variations in Earth’s climate are very unlikely to explain the observed poleward trend in the latitude that tropical cyclone reach their maximum intensity. We find that it is much more likely that the observed trend is forced by human-related emissions, though the spatial pattern of warming in response to greenhouse emissions can have significant impacts on the magnitude of the trend.

Lin, Jonathan↗
DOE OSTI2023

An ecologically-constrained deep learning model for tropical leaf phenology monitoring using PlanetScope satellites

In tropical forests, leaf phenology signals leaf-on/off status and exhibits considerable variability across scales from a single tree-crown to the entire forest ecosystem. Such phenology signals importantly regulate large-scale biogeochemical cycles and regional climate. PlanetScope CubeSats data with a 3-m resolution and near-daily global coverage provide an unprecedented opportunity to monitor both fine- and ecosystem-scale phenology variability along large environmental gradients. However, a scalable method that accurately characterizes leaf phenology from PlanetScope with biophysically meaningful metrics remains lacking. Here we developed an index-guided, ecologically constrained autoencoder (IG-ECAE) method to automatically derive a deciduousness metric (i.e. percentage of tree canopies with leaf-off status within an image pixel) from PlanetScope. The IG-ECAE first estimated the reflectance spectra of leafy/leafless canopies based on their spectral indices characteristics, then used the derived reflectance spectra to guide an autoencoder deep learning method with additional ecological constraints to refine the reflectance spectra, and finally used linear spectral unmixing to estimate the relative abundance of leafless canopies (or deciduousness) per PlanetScope image pixel. We tested the IG-ECAE method at 16 tropical forest sites spanning multiple continents and a large precipitation gradient (1470-2819 mm year -1 ). Among these sites, we evaluated the PlanetScope-derived deciduousness with corresponding measures derived from WorldView-2 (n=9 sites) and local phenocams (n=9 sites). Our results show that PlanetScope-derived deciduousness agrees: 1) with that derived from WorldView-2 at the patch level (90m×90m) with r 2 =0.89 across all sites; and 2) with that derived from phenocams to quantify ecosystem-scale seasonality with r 2 ranging from 0.62-0.96. These results demonstrate the effectiveness and scalability of IG-ECAE in characterizing the wide variability in deciduousness across scales from pixels to forest ecosystems, and from a single date to the full annual cycle, indicating the potential for using high-resolution satellites to track the large-scale phenological patterns and response of tropical forests to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Controls on Soil Organic Matter Degradation and Subsequent Greenhouse Gas Emissions Across a Permafrost Thaw Gradient in Northern Sweden

Warming-induced permafrost thaw could enhance microbial decomposition of previously stored soil organic matter (SOM) to carbon dioxide (CO 2 ) and methane (CH 4 ), one of the most significant potential feedbacks from terrestrial ecosystems to the atmosphere in a changing climate. The environmental parameters regulating microbe-organic matter interactions and greenhouse gas (GHG) emissions in northern permafrost peatlands are however still largely unknown. The objective of this work is to understand controls on SOM degradation and its impact on porewater GHG concentrations across the Stordalen Mire, a thawing peat plateau in Northern Sweden. Here, we applied high-resolution mass spectrometry to characterize SOM molecular composition in peat soil samples from the active layers of a Sphagnum-dominated bog and rich fen sites in the Mire. Microbe-organic matter interactions and porewater GHG concentrations across the thaw gradient were controlled by aboveground vegetation and soil pH. An increasingly high abundance of reduced organic compounds experiencing greater humification rates due to enhanced microbial activity were observed with increasing thaw, in parallel with higher CH 4 and CO 2 porewater concentrations. Bog SOM however contained more Sphagnum-derived phenolics, simple carbohydrates, and organic- acids. The low degradation of bog SOM by microbial communities, the enhanced SOM transformation by potentially abiotic mechanisms, and the accumulation of simple carbohydrates in the bog sites could be attributed in part to the low pH conditions of the system associated with Sphagnum mosses. We show that Gibbs free energy of C half reactions based on C oxidation state for OM can be used as a quantifiable measure for OM decomposability and quality to enhance current biogeochemical models to predict C decomposition rates. We found a direct association between OM chemical diversity and δ 13 C-CH 4 in peat porewater; where higher substrate diversity was positively correlated with enriched δ 13 C-CH 4 in fen sites. Oxidized sulfur-containing compounds, produced by Sphagnum, were further hypothesized to control GHG emissions by acting as electron acceptors for a sulfate-reducing electron transport chain, inhibiting methanogenesis in peat bogs. These results suggest that warming-induced permafrost thaw might increase organic matter lability, in subset of sites that become wetlands, and shift biogeochemical processes toward faster decomposition with an increasing proportion of carbon released as CH 4 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impacts of climate change, population growth, and power sector decarbonization on urban building energy use

Climate, technologies, and socio-economic changes will influence future building energy use in cities. However, current low-resolution regional and state-level analyses are insufficient to reliably assist city-level decision-making. Here we estimate mid-century hourly building energy consumption in 277 U.S. urban areas using a bottom-up approach. The projected future climate change results in heterogeneous changes in energy use intensity (EUI) among urban areas, particularly under higher warming scenarios, with on average 10.1–37.7% increases in the frequency of peak building electricity EUI but over 110% increases in some cities. For each 1 °C of warming, the mean city-scale space-conditioning EUI experiences an average increase/decrease of ~14%/ ~ 10% for space cooling/heating. Heterogeneous city-scale building source energy use changes are primarily driven by population and power sector changes, on average ranging from –9% to 40% with consistent south–north gradients under different scenarios. Across the scenarios considered here, the changes in city-scale building source energy use, when averaged over all urban areas, are as follows: –2.5% to –2.0% due to climate change, 7.3% to 52.2% due to population growth, and –17.1% to –8.9% due to power sector decarbonization. Our findings underscore the necessity of considering intercity heterogeneity when developing sustainable and resilient urban energy systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Measurements of TRACER pre-convective conditions and mesoscale circulations using small unmanned aircraft systems (sUAS)

Improved comprehension of the physical processes governing convective cloud formation and lifecycle are of critical importance for understanding and predicting future climate states. The influence of these clouds on the planetary energy budget, including on precipitation, is significant. Things are particularly complex in coastal regimes, where gradients in aerosol particle properties, localized circulations such as sea breezes, and large population centers are found. To date, numerical models struggle to accurately represent these critical clouds and are therefore challenged to provide a realistic view on the planetary energy budget. Through the proposed research, we deployed two uncrewed aircraft systems (UAS) equipped with a variety of instruments alongside sensors deployed by the US Department of Energy Atmospheric Radiation Measurement (ARM) program for the TRACER (Tracking Aerosol Convection Interactions Experiment) field campaign. The two small UAS platforms consisted of a CU RAAVEN fixed-wing airplane and an OU CopterSonde system, with the copter collecting frequent vertical profiles of thermodynamic and kinematic variables such as temperature, pressure, wind and humidity. At the same time, the fixed-wing captured horizontal gradients of these quantities and aerosol size distribution. These systems were deployed south of the Houston metro area, in an area that is impacted by the Gulf of Mexico sea breeze on a daily basis. These observations offer enhanced and complementary perspectives to those provided by the DOE ARM Mobile Facility (AMF) which is was deployed in southeast Houston, and an ancillary site in a more rural location west of the urban Houston area. Quality-controlled versions of the UAS data were collected and posted on the DOE ARM data archive after the conclusion of the campaign where they are accessible by the research community and general public. The UAS perspective offers revolutionary insight into key spatial and temporal effects that have not been evaluated previously.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Recovery of Forest Structure Following Large-Scale Windthrows in the Northwestern Amazon

The dynamics of forest recovery after windthrows (i.e., broken or uprooted trees by wind) are poorly understood in tropical forests. The Northwestern Amazon (NWA) is characterized by a higher occurrence of windthrows, greater rainfall, and higher annual tree mortality rates (~2%) than the Central Amazon (CA). We combined forest inventory data from three sites in the Iquitos region of Peru, with recovery periods spanning 2, 12, and 22 years following windthrow events. Study sites and sampling areas were selected by assessing the windthrow severity using remote sensing. At each site, we recorded all trees with a diameter at breast height (DBH) ≥ 10 cm along transects, capturing the range of windthrow severity from old-growth to highly disturbed (mortality > 60%) forest. Across all damage classes, tree density and basal area recovered to >90% of the old-growth values after 20 years. Aboveground biomass (AGB) in old-growth forest was 380 (±156) Mg ha -1 . In extremely disturbed areas, AGB was still reduced to 163 (±68) Mg ha -1 after 2 years and 323 (± 139) Mg ha -1 after 12 years. This recovery rate is ~50% faster than that reported for Central Amazon forests. The faster recovery of forest structure in our study region may be a function of its higher productivity and adaptability to more frequent and severe windthrows. These varying rates of recovery highlight the importance of extreme wind and rainfall on shaping gradients of forest structure in the Amazon, and the different vulnerabilities of these forests to natural disturbances whose severity and frequency are being altered by climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Permafrost microbial communities and functional genes are structured by latitudinal and soil geochemical gradients

Permafrost underlies approximately one quarter of Northern Hemisphere terrestrial surfaces and contains 25–50% of the global soil carbon (C) pool. Permafrost soils and the C stocks within are vulnerable to ongoing and future projected climate warming. The biogeography of microbial communities inhabiting permafrost has not been examined beyond a small number of sites focused on local-scale variation. Permafrost is different from other soils. Perennially frozen conditions in permafrost dictate that microbial communities do not turn over quickly, thus possibly providing strong linkages to past environments. Thus, the factors structuring the composition and function of microbial communities may differ from patterns observed in other terrestrial environments. Here, we analyzed 133 permafrost metagenomes from North America, Europe, and Asia. Permafrost biodiversity and taxonomic distribution varied in relation to pH, latitude and soil depth. The distribution of genes differed by latitude, soil depth, age, and pH. Genes that were the most highly variable across all sites were associated with energy metabolism and C-assimilation. Specifically, methanogenesis, fermentation, nitrate reduction, and replenishment of citric acid cycle intermediates. This suggests that adaptations to energy acquisition and substrate availability are among some of the strongest selective pressures shaping permafrost microbial communities. The spatial variation in metabolic potential has primed communities for specific biogeochemical processes as soils thaw due to climate change, which could cause regional- to global- scale variation in C and nitrogen processing and greenhouse gas emissions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A research agenda for nonvascular photoautotrophs under climate change

Nonvascular photoautotrophs (NVP), including bryophytes, lichens, terrestrial algae, and cyanobacteria, are increasingly recognized as being essential to ecosystem functioning in many regions of the world. Current research suggests that climate change may pose a substantial threat to NVP, but the extent to which this will affect the associated ecosystem functions and services is highly uncertain. Here, we propose a research agenda to address this urgent question, focusing on physiological and ecological processes that link NVP to ecosystem functions while also taking into account the substantial taxonomic diversity across multiple ecosystem types. Accordingly, we developed a new categorization scheme, based on microclimatic gradients, which simplifies the high physiological and morphological diversity of NVP and world-wide distribution with respect to several broad habitat types. We found that habitat-specific ecosystem functions of NVP will likely be substantially affected by climate change, and more quantitative process understanding is required on (1) potential for acclimation, (2) response to elevated CO 2 , (3) role of the microbiome, and (4) feedback to (micro)climate. We suggest an integrative approach of innovative, multimethod laboratory and field experiments and ecophysiological modelling, for which sustained scientific collaboration on NVP research will be essential.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Energy‐Water Asynchrony Principally Determines Water Available for Runoff From Snowmelt in Continental Montane Forests

ABSTRACT Changes in the volume, rate, and timing of the snowmelt water pulse have profound implications for seasonal soil moisture, evapotranspiration (ET), groundwater recharge, and downstream water availability, especially in the context of climate change. Here, we present an empirical analysis of water available for runoff using five eddy covariance towers located in continental montane forests across a regional gradient of snow depth, precipitation seasonality, and aridity. We specifically investigated how energy‐water asynchrony (i.e., snowmelt timing relative to atmospheric demand), surface water input intensity (rain and snowmelt), and observed winter ET (winter AET) impact multiple water balance metrics that determine water available for runoff (WAfR). Overall, we found that WAfR had the strongest relationship with energy‐water asynchrony (adjusted r 2 = 0.52) and that winter AET was correlated to total water year evapotranspiration but not to other water balance metrics. Stepwise regression analysis demonstrated that none of the tested mechanisms were strongly related to the Budyko‐type runoff anomaly (highest adjusted r 2 = 0.21). We, therefore, conclude that WAfR from continental montane forests is most sensitive to the degree of energy‐water asynchrony that occurs. The results of this empirical study identify the physical mechanisms driving variability of WAfR in continental montane forests and are thus broadly relevant to the hydrologic management and modelling communities.

Webb, Ryan William↗
DOE OSTI2023

Using long‐term data from a whole ecosystem warming experiment to identify best spring and autumn phenology models

Abstract Predicting vegetation phenology in response to changing environmental factors is key in understanding feedbacks between the biosphere and the climate system. Experimental approaches extending the temperature range beyond historic climate variability provide a unique opportunity to identify model structures that are best suited to predicting phenological changes under future climate scenarios. Here, we model spring and autumn phenological transition dates obtained from digital repeat photography in a boreal Picea ‐ Sphagnum bog in response to a gradient of whole ecosystem warming manipulations of up to +9°C, using five years of observational data. In spring, seven equally best‐performing models for Larix utilized the accumulation of growing degree days as a common driver for temperature forcing. For Picea , the best two models were sequential models requiring winter chilling before spring forcing temperature is accumulated. In shrub, parallel models with chilling and forcing requirements occurring simultaneously were identified as the best models. Autumn models were substantially improved when a CO 2 parameter was included. Overall, the combination of experimental manipulations and multiple years of observations combined with variation in weather provided the framework to rule out a large number of candidate models and to identify best spring and autumn models for each plant functional type.

Schädel, Christina↗
DOE OSTI2025

Resource assessment of ocean thermal energy conversion in Puerto Rico and U.S. Virgin Islands

Island communities often struggle to establish and maintain traditional electric grids and are therefore heavily reliant on costly imported fossil fuels. In the case of Puerto Rico, these challenges are enhanced by extreme weather and other natural hazards that threaten the local electricity generation and transmission infrastructure. Ocean thermal energy conversion (OTEC) could play an important role in establishing a more resilient electrical grid in the region. Here, in this study, a detailed analysis is conducted to characterize the ocean thermal resource and power potential of OTEC in Puerto Rico based on a 14-year dataset of modeled ocean temperature. The assessment considers seasonal and interannual variability in the region's thermal resource and examines the operational limitations associated with minimal thermal gradients required to run a typical OTEC heat engine. Notably, the local thermal resource is found to be sensitive to El Niño-Southern Oscillation (ENSO) climate patterns, with La Niña conditions linked to greater OTEC power availability. Seven areas of opportunity are identified based on their resource potential and proximity to existing electrical distribution lines, including two that could benefit the nearby U.S. Virgin Islands. The greatest OTEC power potential is observed to the south of the main island of Puerto Rico in the Caribbean Sea with an estimated capacity of 138 MW for a plant pumping cold water from a depth of 1,000 m, or the equivalent amount of electricity required to power 219,000 households.

OTEC↗
DOE OSTI2021

Carbon Fluxes and Microbial Activities From Boreal Peatlands Experiencing Permafrost Thaw

Permafrost thaw is occurring in Alaska which may result in a positive feedback to climate warming, due to the release of greenhouse gases such as CO 2 and CH 4 from soils. Here we examined greenhouse gas production along a gradient of “time since thaw,” hypothesizing that fluxes and microbial activities would be highest soon after thaw, and then decline. We observed highest rates of microbial activities, particularly methanogenesis, soon after thaw, coinciding with less decomposed organic matter and higher concentrations of dissolved carbon and nitrogen in soil, possibly of permafrost origin. However, field fluxes were higher in the young thaw site, compared to the older sites, in winter and not summer, a phenomenon that is currently not well understood.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Active populations and growth of soil microorganisms are framed by mean annual precipitation in three California annual grasslands

Climate influences soil microbial composition and function, but the relative importance of a site's historic climate versus its more immediate environmental conditions is unclear. Using quantitative stable isotope probing (qSIP), we characterized actively growing soil microbial communities and soil properties in three California annual grasslands that span a rainfall gradient and have developed on similar parent material. The soils were assayed in the wet winter season, when environmental conditions are most similar across sites. Since growing populations might be expected to be most responsive to contemporary environmental conditions, we hypothesized that the structure of growing microbial communities would be more similar across the gradient than that of total communities (i.e., including non-growing populations). In addition, we hypothesized that population growth rates would be slowest in the driest site, reflecting a legacy effect of low soil moisture on microbial growth. Soils along the rainfall gradient differed in pH, texture, and cation exchange capacity, but not in total C, C:N or dominant minerals. The radiocarbon (14C) age of soil C (reflecting turnover time) increased with mean annual precipitation but soil respiration was uniformly modern, reflecting microbial reliance on recent C inputs across the sites. The structure of both total and growing microbial communities differed across sites. Across major microbial phyla, including the Actinobacteria, Acidobacteria, Bacteroidetes, Gemmatimonadetes and Proteobacteria, bacterial growth rates were consistently lower in the site with the lowest mean annual precipitation. Taxa that were growing at the dry site alone grew more slowly than taxa that grew at multiple sites. These results reflect the influence of climate history and point to the role of environmental filtering at the driest site in shaping its slower growing microbial community, possibly reflecting adaptation to repeated exposure to water stress. Lastly, across taxa, the growth rate of a taxon at one site was correlated with its growth rate in the other sites. Furthermore, this growth rate coherence is likely a consequence of genetically determined physiological traits and is consistent with the idea that evolutionary history constrains growth rate.

Environmental filtering↗
DOE OSTI2024

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗
DOE OSTI2021

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Process-based climate change assessment for European winds using EURO-CORDEX and global models

Wind energy is an important pillar of decarbonization strategies and potentially vulnerable to climate change. Existing wind climate change assessments rely on climate models but a systematic investigation of the global-to-regional climate modeling chain is missing. In this study, I highlight key limitations, namely (a) the differing representation of land use change in global and regional climate models which compromises comparability, and (b) the consistency of large-scale features along the global-to-regional climate modeling chain. To this end, I analyze the large European Coordinated Downscaling Experiment (EURO-CORDEX) ensemble (rcp85: N = 49; rcp45: N= 18; rcp26: N = 22) along with the driving global models (rcp85: N = 7; rcp45: N = 5; rcp26: N = 7), finding evidence that climate change reduces mean wind speeds by up to -0.8 ms -1 (offshore) and -0.3 ms -1 (onshore). I provide physical explanations for these changes by identifying two key drivers. First, onshore wind speeds drop in the driving global models in regions and scenarios with strong land use change but show no drop in EURO-CORDEX where land use is held constant. Second, offshore wind reductions follow decreases in the equator-to-pole temperature gradient remarkably well with correlations reaching around 0.9 in resource-rich European countries like Ireland, the United Kingdom and Norway, implying that arctic amplification is a severe risk for European offshore wind energy. My results suggest that earlier conclusions of negligible climate change impacts on wind energy might be premature if either land use changes strongly or polar amplification is at or above the range sampled in global climate models.

54 ENVIRONMENTAL SCIENCES↗