Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “land surface temperature”

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 379 records · Page 21

An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast—Part 2: Circulation, Boundary Layer, and Clouds

The Western North Atlantic Ocean (WNAO) is a complex land-ocean-atmosphere system that experiences a wide range of atmospheric phenomena, which in turn drive unique aerosol transport pathways, cloud morphologies, and boundary layer variability. This work, Part 2 of a 2-part paper series, provides an overview of the atmospheric circulation, boundary layer variability, three-dimensional cloud structure, and precipitation over the WNAO; the companion paper (Part 1) focused on chemical characterization of aerosol, gas, and wet deposition. Seasonal changes in the circulation and sea surface temperature explain a clear transition in cloud morphologies from small shallow cumulus clouds, convective clouds, and tropical storms in summer, to stratus/stratocumulus and multi-layer cloud systems associated with winter storms. Synoptic variability in cloud fields is estimated using satellite-based weather states, and the role of postfrontal conditions (cold-air outbreaks) in the development of stratiform clouds is further analyzed. Precipitation is persistent over the ocean, with a regional peak over the Gulf Stream path, where offshore sea surface temperature gradients are large and surface fluxes reach a regional peak. Satellite data show a clear annual cycle in cloud droplet number concentration along the coast in winter, and minima in summer, suggesting a marked annual cycle in aerosol-cloud interactions. The performance of four climate models are assessed against satellite cloud retrievals. The models qualitatively reproduce the annual cycle in cloud cover and liquid water path, but with large discrepancies across models, especially in the extra-tropics. The paper concludes with a summary of outstanding issues and recommendations for future work.

54 ENVIRONMENTAL SCIENCES↗

CHELAX-BNF: Turbulent Parameters by airborne measurements

The original data were collected on board the ARM Aerial Facility ArcticShark uncrewed aerial system (UAS; https://www.arm.gov/capabilities/observatories/aaf/uas ) during the “Characterizing HEterogeneous Land-Atmosphere eXchanges at BNF” field campaign (CHEAX-BNF; https://arm.gov/research/campaigns/aaf2025CHELAX-BNF ). The ARM Aerial Facility ArcticShark UAS was based at the public-use airport of Posey Field, AL (FAA LID: 1M4, 34.28027778° N, 87.60055556° W, 283m MSL) from May 28 through June 23, 2025. The ArcticShark UAS performed 5 flights, including 4 research flights over the BNF Main Site (ARM Mobile Facility 3, https://arm.gov/capabilities/observatories/amf ) and Supplemental Facilities to measure atmospheric state, turbulence, surface IR temperature and imagery, aerosol number concentration, and aerosol size distribution. The current data set presents a collection of turbulent parameters in the atmospheric boundary layer or lower free troposphere based on airborne measurement throughout the field campaign. The primary instruments used to create the current data set were the Aircraft Integrated Meteorological Measurement System (AIMMS-30) and the fine-wire thermocouple probe.

Aircraft Integrated Meteorological Measurement Sys↗

NASA Dragonfly Mission MMRTG

The Radioisotope Power Systems (RPS) group at the Idaho National Laboratory (INL) is supporting NASA’s New Frontier class mission, Dragonfly, to explore the Saturn moon Titan. Titan has been observed by both Voyager spacecrafts in 1979 and 1980 and during over 100 close flyby encounters by the Cassini orbiter and subsequent landing of the Huygens probe which landed on Titan in 2005. The probe lasted about 2.5 hours without an RPS due to the thick atmosphere and cold climate but gave scientists valuable insight into the geology of Titan. The surface is composed of vast sand dunes and liquid methane with an atmosphere composed primarily of nitrogen that is 60% more dense as the Earth’s with 1/7th of Earth’s gravity. These conditions allow for exploration using a dual-quadcopter. Due to the extreme temperatures (-290°F or -180°C), changes to the Multi-Mission Radioisotope Thermoelectric Generator (MMRTG) are required. This generator, designated flight unit 4 or MMRTG-F4, is being built by a team lead by INL and includes the prime contractor Aerojet Rocketdyne and their supplier Teledyne Energy System. This poster outlines the mission, the unique challenges and updates the build progress.

42 ENGINEERING↗

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc↗

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

Shrub Expansion Can Counteract Carbon Losses From Warming Tundra

Arctic warming is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and tree-line expansion and densification. However, predicting the carbon trajectories of the changing Arctic is challenging due to interacting feedbacks between vegetation composition and structure, and surface characteristics. We conduct a sensitivity analysis of the current-date to 2100 projected surface energy fluxes, soil carbon pools, and CO 2 fluxes to different shrub expansion rates under future emission scenarios (intermediate—RCP4.5, and high—RCP8.5) using the Arctic-focused configuration of E3SM Land Model (ELM). We focus on Trail Valley Creek (TVC), an upland tundra site in the western Canadian Arctic, which is experiencing shrub densification and expansion. We find that shrub expansion did not significantly alter the modeled surface energy and water budgets. However, the carbon balance was sensitive to shrub expansion, which drove higher rates of carbon sequestration as a consequence of higher shrubification rates. Thus, at low shrub expansion rates, the site would become a carbon source, especially under RCP8.5, due to higher temperatures, which deepen the active layer and enhance soil respiration. At higher shrub expansion rates, TVC would become a net CO 2 sink under both Representative Concentration Pathway scenarios due to higher shrub productivity outweighing temperature-driven respiration increase. Our simulations highlight the effect of shrub expansion on Arctic ecosystem carbon fluxes and stocks. We predict that at TVC, shrubification rate would interact with climate change intensity to determine whether the site would become a carbon sink or source under projected future climate.

Yazbeck, Theresia [The Ohio State Univ., Columbus,↗

Blodgett 13C–labeled litter incubation 2016-2019

The dataset is from 13C-labelled (stable isotope of carbon) root-litter in-situ field incubation experiment based on the whole-soil warming experiment at the Blodgett Forest Research Station, CA, USA. The files are in both ".csv" and ".xlsx" versions, and can be opened in "maCOS numbers", and "Microsoft Excel". The files includes several sheets with all the data published in the paper: Sun, B., Zosso, C., Wiesenberg, G. L. B., Pegoraro, E., Torn, M. S., and Schmidt, M. W. I.: Warming accelerates the decomposition of root-derived hydrolysable lipids in a temperate forest and is depth- and compound class-dependent, SOIL, 11, 1077–1093, https://doi.org/10.5194/soil-11-1077-2025, 2025. This dataset includes bulk soil carbon, nitrogen, delta 13C values, the normalized concentration (to organic carbon) of hydrolysable lipids identified, the absolute concentration (normalized to bulk soil) of hydrolysable lipids, hydrolysable lipids recovery, and weighted 13C-excess of bulk soil carbon, and weighted 13C-excess of each compound class in hydrolysable lipids. These data aim to answer two research questions: 1) How will warming affect the decomposition of 13C-labelled root-litter at different depth? 2) Will the decomposition of root-derived hydrolysable lipids under warming differ among different compound classes? The experiment sites located on the foothills of the Sierra Nevada near Georgetown, CA (120°3904000W; 38°5404300 N) at 1370m above see level. The Blodgett Forest is a mixed-coniferous forest. The site has a Mediterranean climate with a mean annual air temperature of 12.5 °C and a mean annual precipitation of 1774mm.

54 ENVIRONMENTAL SCIENCES↗

A Decadal Hybrid GCM Simulation Using Deep‐Learning‐Based Cloud and Convection Parameterization Generalized to a Warm Climate

A critical challenge for machine‐learning (ML) parameterization in global climate models (GCMs) is to achieve stable, accurate simulations under climates not seen during training. Previous studies have demonstrated promising offline performance and year‐long online stability in aquaplanet simulations but have encountered difficulties in real geography and under climate warming. Here we report that a GCM with real geography configuration using neural‐network‐based cloud and convection parameterization, trained exclusively with present‐day climate data, successfully performs a stable, decade‐long simulation of a warm climate with +4 K sea surface temperature (SST). The neural network (NN) is based on Han et al. (2023, https://doi.org/10.1029/2022ms003508 ) with additional inputs. The simulation captures the global precipitation distribution, surface temperatures, vertical atmospheric structures, and extreme precipitation very well, closely matching simulations from both the superparameterized CAM (SPCAM) and the conventional CAM5 in the warm climate without accuracy degradation compared to those in the baseline climate. Moreover, it produces a climate response to +4 K SST in atmospheric thermodynamic states and circulations similar to those from SPCAM and CAM5. Prognostic ablation tests on NN input variables show that the NN without convective memory as input suffers from numerical instability, and the NN without considering radiative variables and land fraction as input, or with reduced training samples produce less accurate results. To our knowledge, this is the first time an ML parameterization successfully achieves online extrapolation to a warm climate without using additional warm‐climate data for training. It demonstrates the potential of ML‐driven parameterizations for credible long‐term climate projections.

Atmosphere model↗

Warming Response of Deep Soil Carbon (LDRD Final Report)

The overarching objective of this LDRD project was to determine the vulnerability of deep soil organic carbon (SOC) to warming in the Sierra Nevada region. Deep soils (>30 cm) store more than 70% of global SOC, and increased SOC decomposition and CO 2 emissions caused by warming are potentially large climate change feedbacks. According to the Intergovernmental Panel on Climate Change, temperatures are expected to increase by 4°C by the year 2100, warming the land and underlying soil, and making understanding of how warming will influence deep SOC storage and persistence critical to projecting the land carbon sink. However, uncertainty remains in our process-level understanding and ability to quantify how projected warming will impact the stability of carbon in deep soils. We investigated warming effects on deep (up to 16 m) SOC stability across climate, vegetation, and soil mineralogy gradients in California. We sampled soils from the surface to bedrock (down to 16 meters) at four sites representing vastly different ecosystems across the Northern and Southern Sierra Nevada mountains. This study quantified the response of SOC concentration, distribution, and vulnerability to warming using a soil incubation experiment to warm the whole soil profile and radiocarbon and stable isotopes to assess which pools are vulnerable to loss under warming. Our results refine our understanding of the terrestrial carbon cycle by revealing the vulnerability of deep SOC to future changes in climate and how minerology may influence that vulnerability.

58 GEOSCIENCES↗

Implications of carbon content for the processing, stability, and mechanical properties of cast and wrought Ni-based superalloy Nimonic 105

As the temperature and pressure requirements in land-based turbines for power generation increase, the shift towards more advanced alloys, e.g., Ni-based superalloys, requires improved understanding of their long-term stability and mechanical properties. Additionally, the larger size of the components also presents fabrication and cost-reduction problems. In this study, we investigated Nimonic 105 as a potential rotor material at two carbon content levels – the commonly used maximum of the alloy specification and a lower carbon content variant. Reducing the carbon content presented fabrication challenges, resulting in surface crack formation during hot working. Additionally, a lower M 6 C carbide fraction was present after fabrication/solutioning, leading to a significantly larger grain size. After the standard aging treatment, similar γ' populations formed in the two variants, however the M 6 C carbides in the high-carbon variant already started transforming to intergranular M 23 C 6 . Although we did not observe a noticeable effect on the tensile properties of the two variants, the lower carbon content alloy possessed a superior creep property – mainly due to the larger grain size as the microstructural stability was inferior. Further, upon prolonged thermal exposure, the M 6 C carbides in the low-carbon variant further decomposed to intergranular and intragranular σ and μ precipitates, as opposed to mainly M 23 C 6 as in the high-carbon variant.

36 MATERIALS SCIENCE↗

Five Years of Dissolved Oxygen, Temperature, Salinity, Depth, Weather Data from a Transitioning Wetland at Beaver Creek, Washington, USA

Groundwater dissolved oxygen (DO) variability in coastal system remains poorly understood despite its importance for biogeochemical cycling and ecosystem modeling. Here we investigate the temporal variability in groundwater DO and its hydro-climatic drivers across hourly to seasonal timescales in a transitioning wetland at Beaver Creek, Washington, USA. The site is transitioning from a freshwater forest to a brackish tidal wetland following removal of a barrier in 2014 that prevented tides from accessing the freshwater creek. By utilizing novel optical dissolved oxygen instrumentation (Opti O2, LLC) we obtained continuous, high-frequency (5-minute), in-situ measurements of DO from the flood-plain from June 26th, 2019 through September 30th, 2024. This 63 month dataset is comprised of groundwater dissolved oxygen, temperature, water level and salinity timeseries from the floodplain. This dataset also includes rainfall, air pressure, air temperature, and solar radiation data collected with a co-located Campbell ClimaVUE50 weather sensor. All data is contained within a single csv (2019-06-26 to 2024-09-30 Beaver Creek DO, saln, BGS, temp, weather.csv) that can easily be viewed either using software such as Excel or using any text editor.

54 ENVIRONMENTAL SCIENCES↗

Observational benchmarks inform representation of soil organic carbon dynamics in land surface models

Abstract. Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Machine learning models can help identify dominant environmental controllers and establish their functional relationships with SOC stocks. The resulting knowledge can be integrated into ESMs to reduce uncertainty and improve predictions of SOC dynamics over space and time. In this study, we used a large number of SOC field observations (n=54 000), geospatial datasets of environmental factors (n=46), and two machine learning approaches (namely random forest, RF, and generalized additive modeling, GAM) to (1) identify dominant environmental controllers of global and biome-specific SOC stocks, (2) derive functional relationships between environmental controllers and SOC stocks, and (3) compare the identified environmental controllers and predictive relationships with those in models used in Phase 6 of the Coupled Model Intercomparison Project (CMIP6). Our results showed that the diurnal temperature, drought index, cation exchange capacity, and precipitation were important observed environmental predictors of global SOC stocks. While the RF model identified 14 environmental factors that describe climatic, vegetation, and edaphic conditions as important predictors of global SOC stocks (R2=0.61, RMSE = 0.46 kg m−2), current ESMs oversimplify the relationships between environmental factors and SOC, with precipitation, temperature, and net primary productivity explaining > 96 % of the variability in ESM-modeled SOC stocks. Further, our study revealed notable disparities among the functional relationships between environmental factors and SOC stocks simulated by ESMs compared with observed relationships. To improve SOC representations in ESMs, it is imperative to incorporate additional environmental controls, such as the cation exchange capacity, and refine the functional relationships to align more closely with observations.

54 ENVIRONMENTAL SCIENCES↗

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗

The underappreciated importance of solar radiation in constraining spring phenology of temperate ecosystems in the Northern and Eastern United States

Spring phenology of temperate ecosystems is highly sensitive to climate change, generating various impacts on many important terrestrial surface biophysical processes. Although various prognostic models relying on environmental variables of temperature and photoperiod have been developed for spring phenology, comprehensive ecosystem-scale evaluations over large landscapes and long-time periods remain lacking. Further, environmental variables other than temperature and photoperiod might also importantly constrain spring phenology modelling but remain under-investigation. To address these issues, we leveraged around 20-years datasets of environmental variables (from Daymet and GLDAS products) and the spring phenology metric (i.e., the greenup date) respectively derived from MODIS and PhenoCams across 108 sites in the Northern and Eastern United States. We firstly cross-compared MODIS-derived greenup date with official PhenoCams product with high accuracy (R 2 = 0.70). Then, we evaluated the three prognostic models (i.e., Growing Degree Date (GDD), Sequential (SEQ) and optimality-based (OPT)) with MODIS-derived spring phenology, assessed the model residuals and their associations with soil moisture, rainfall, and solar radiation, and revised the two photoperiod-relevant models (SEQ, OPT) by replacing the daylength variable with solar radiation, which was found to contribute the most to model residuals. We found that 1) all models demonstrated good capability in characterizing spring phenology, with OPT performing the best (RMSE = 8.04 ± 5.05 days), followed by SEQ (RMSE = 10.57 ± 7.77 days) and GDD (RMSE = 10.84 ± 8.42 days), 2) all models displayed high model residuals showing tight correlation with solar radiation (r = 0.45–0.75), and 3) the revised models that included solar radiation significantly performed better with an RMSE reduction by 22.08%. Such results are likely because solar radiation better constrains early growing season plant photosynthesis than photoperiod, supporting the hypothesis of spring phenology as an adaptive strategy to maximize photosynthetic carbon gain (approximated by solar radiation) while minimizing frost damage risk (captured by temperature). Collectively, our study here reveals the underappreciated importance of solar radiation in constraining spring phenology of temperate ecosystems, and suggests ways to improve spring phenology modelling and other phenology-related ecological processes.

54 ENVIRONMENTAL SCIENCES↗

A permafrost implementation in the simple carbon–climate model Hector v.2.3pf

Abstract. Permafrost currently stores more than a fourth of global soil carbon. A warming climate makes this carbon increasingly vulnerable to decomposition and release into the atmosphere in the form of greenhouse gases. The resulting climate feedback can be estimated using land surface models, but the high complexity and computational cost of these models make it challenging to use them for estimating uncertainty, exploring novel scenarios, and coupling with other models. We have added a representation of permafrost to the simple, open-source global carbon–climate model Hector, calibrated to be consistent with both historical data and 21st century Earth system model projections of permafrost thaw. We include permafrost as a separate land carbon pool that becomes available for decomposition into both methane (CH4) and carbon dioxide (CO2) once thawed; the thaw rate is controlled by region-specific air temperature increases from a preindustrial baseline. We found that by 2100 thawed permafrost carbon emissions increased Hector’s atmospheric CO2 concentration by 5 %–7 % and the atmospheric CH4 concentration by 7 %–12 %, depending on the future scenario, resulting in 0.2–0.25 ∘C of additional warming over the 21st century. The fraction of thawed permafrost carbon available for decomposition was the most significant parameter controlling the end-of-century temperature change in the model, explaining around 70 % of the temperature variance, and was distantly followed by the initial stock of permafrost carbon, which contributed to about 10 % of the temperature variance. The addition of permafrost in Hector provides a basis for the exploration of a suite of science questions, as Hector can be cheaply run over a wide range of parameter values to explore uncertainty and can be easily coupled with integrated assessment and other human system models to explore the economic consequences of warming from this feedback.

54 ENVIRONMENTAL SCIENCES↗

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗

Plant Physiology Increases the Magnitude and Spread of the Transient Climate Response to CO 2 in CMIP6 Earth System Models

Increasing concentrations of CO 2 in the atmosphere influence climate both through CO 2 ’s role as a greenhouse gas and through its impact on plants. Plants respond to atmospheric CO 2 concentrations in several ways that can alter surface energy and water fluxes and thus surface climate, including changes in stomatal conductance, water use, and canopy leaf area. These plant physiological responses are already embedded in most Earth system models, and a robust literature demonstrates that they can affect global-scale temperature. However, the physiological contribution to transient warming has yet to be assessed systematically in Earth system models. Here this gap is addressed using carbon cycle simulations from phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP) to isolate the radiative and physiological contributions to the transient climate response (TCR), which is defined as the change in globally averaged near-surface air temperature during the 20-yr window centered on the time of CO 2 doubling relative to preindustrial CO 2 concentrations. In CMIP6 models, the physiological effect contributes 0.12°C ( σ : 0.09°C; range: 0.02°–0.29°C) of warming to the TCR, corresponding to 6.1% of the full TCR ( σ : 3.8%; range: 1.4%–13.9%). Moreover, variation in the physiological contribution to the TCR across models contributes disproportionately more to the intermodel spread of TCR estimates than it does to the mean. The largest contribution of plant physiology to CO 2 -forced warming—and the intermodel spread in warming—occurs over land, especially in forested regions.

54 ENVIRONMENTAL SCIENCES↗