Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sea level”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Performance of reanalysis and mesoscale models off the coast of Hawai'i

The eastern Hawai'i coast in the United States is characterized by considerable wind resource fuelled by persistent trade winds, making it an important area for energy research. The need is strong for reanalyses and higher-resolution regional simulations where observations have been historically limited, such as Hawai'i's offshore environments. However, studies using offshore observations in other parts of the world have shown that significant errors can occur in reanalyses and wind datasets, which can lead to inaccurate estimates of wind energy generation, payback periods, and extreme weather risks at project locations. The degree of such errors is influenced by a number of factors, including spatial resolution and the handling of processes within the planetary boundary layer (PBL). In this work, we provide a wind resource characterization from year-long lidar buoy measurements off the eastern coast of O'ahu, Hawai'i, an environment previously unobserved at the rotor level, and use the characterization to evaluate the performance of two simulation datasets. The O'ahu deployment location is meteorologically unique and less complex than land-based wind resource characterizations, being strongly characterized by trade winds with minimal land–atmosphere interaction influences. Despite the unique and fairly consistent meteorological conditions, we hypothesize that distinct simulation datasets will exhibit diverse ranges of errors similar to those that have been seen for other offshore locations. We find the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) to strongly underestimate observed wind speeds at the O'ahu location (bias = −1.54 m s −1 at a height of 140 m above sea level), while a regional Weather Research and Forecasting Model (WRF) simulation produced by the University of Hawai'i (UH-WRF) provides a significantly smaller wind speed bias (−0.25 m s −1 ), highlighting the value of running regional, higher-resolution simulations. The large bias noted for ERA5 is driven by significant underestimation of fast wind speeds (>9 m s −1 ), which the study site is largely characterized by, along with discontinuities in the ERA5 diurnal cycle. We also speculate that the relative sparsity of observations for data assimilation in this remote part of the world could influence the performance of ERA5 and that challenges with characterizing island effects could impact the performance of both datasets.

17 WIND ENERGY↗
DOE OSTI2022

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Improved Understanding of Multicentury Greenland Ice Sheet Response to Strong Warming in the Coupled CESM2‐CISM2 With Regional Grid Refinement

The simulation of ice sheet‐climate interactions, such as surface mass balance fluxes, is sensitive to model grid resolution. Here we simulate the multi‐century evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO 2 increases by 1% yr -1 until quadrupling the pre‐industrial level and then is held fixed). A variable‐resolution (VR) grid with 1/4° regional refinement over the broader Arctic and 1° resolution elsewhere is applied to the atmosphere and land components, and the results are compared with conventional 1° lat‐lon grid simulations to investigate the impact of grid refinement. Compared with the 1° runs, the VR run features a slower rate of surface melt, especially over the western and northern GrIS, where the ice surface slopes gently toward the periphery. This difference pattern originates primarily from higher snow albedo and, thus, weaker albedo feedback in the VR run. The VR grid better captures the CISM ice sheet topography by reducing elevation discrepancies between CAM and CISM and is, therefore, less reliant on the downscaling algorithm, which is known to underestimate albedo gradients. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% less than that of the 1° runs, respectively.

Earth System Model↗
DOE OSTI2023

Decoupling of Arctic variability from the North Pacific in a warmer climate

The sea level pressure variability in the North Pacific modulates the climate of the Arctic and surrounding continents, substantially impacting ecosystems and indigenous communities. Our analysis based on data from the CESM2 Large Ensemble and different Model Intercomparison Project datasets reveals that the interannual variability of Arctic surface air temperature (SAT) gradually decouples from the contemporaneous atmospheric conditions over the North Pacific as external forcing increases in intensity in the future. Future projections show that the North Pacific-Arctic relationship during the fall season consistently weakens in magnitude until the end of this century, and in the 22nd and 23rd centuries, the relationship is negligible throughout the year. We show that under increased greenhouse gas emissions, the regional heat fluxes extensively control the Arctic temperature variability, and the strength of the projected North Pacific-Arctic relationship is strongly dependent on the Arctic sea ice extent. Our results suggest that under future warming, a strong coupling of Arctic SAT with the underlying ocean and a weakening of the meridional pressure gradient driven by an enhanced rate of sea ice retreat will weaken the interannual footprint of North Pacific variability on Arctic SAT. Therefore, we propose that the alarming rate of sea-ice decline over recent decades and projected in the near future could accelerate the rate of decoupling. Further, we suggest that mitigation strategies for the Arctic should focus on regional mechanisms operating on interannual and seasonal timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Wintertime extreme warming events in the high Arctic: characteristics, drivers, trends, and the role of atmospheric rivers

Abstract. An extreme warming event near the North Pole, with 2 m temperature rising above 0 °C, was observed in late December 2015. This specific event has been attributed to cyclones and their associated moisture intrusions. However, little is known about the characteristics and drivers of similar events in the historical record. Here, using data from European Centre for Medium-Range Weather Forecasts Reanalysis, version 5 (ERA5), we study these winter extreme warming events with 2 m temperature over a grid point above 0 °C over the high Arctic (poleward of 80° N) that occurred during 1980–2021. In ERA5, such wintertime extreme warming events can only be found over the Atlantic sector. They occur rarely over many grid points, with a total absence during some winters. Furthermore, even when occurring, they tend to be short-lived, with the majority of the events lasting for less than a day. By examining their surface energy budget, we found that these events transition with increasing latitude from a regime dominated by turbulent heat flux into the one dominated by downward longwave radiation. Positive sea level pressure anomalies which resemble blocking over northern Eurasia are identified as a key ingredient in driving these events, as they can effectively deflect the eastward propagating cyclones poleward, leading to intense moisture and heat intrusions into the high Arctic. Using an atmospheric river (AR) detection algorithm, the roles of ARs in contributing to the occurrence of these extreme warming events defined at the grid-point scale are explicitly quantified. The importance of ARs in inducing these events increases with latitude. Poleward of about 83° N, 100 % of these events occurred under AR conditions, corroborating that ARs were essential in contributing to the occurrence of these events. Over the past 4 decades, both the frequency, duration, and magnitude of these events have been increasing significantly. As the Arctic continues to warm, these events are likely to increase in both frequency, duration, and magnitude, with great implications for the local sea ice, hydrological cycle, and ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Offshore low-level jet observations and model representation using lidar buoy data off the California coast

Abstract. Low-level jets (LLJs) occur under a variety of atmospheric conditions and influence the available wind resource for wind energy projects. In 2020, lidar-mounted buoys owned by the US Department of Energy (DOE) were deployed off the California coast in two wind energy lease areas administered by the Bureau of Ocean Energy Management: Humboldt and Morro Bay. The wind profile observations from the lidars and collocated near-surface meteorological stations (4–240 m) provide valuable year-long analyses of offshore LLJ characteristics at heights relevant to wind turbines. At Humboldt, LLJs were associated with flow reversals and north-northeasterly winds, directions that are more aligned with terrain influences than the predominant northerly flow. At Morro Bay, coastal LLJs were observed primarily during northerly flow as opposed to the predominant north-northwesterly flow. LLJs were observed more frequently in colder seasons within the lowest 250 m a.s.l. (above sea level), in contrast with the summertime occurrence of the higher-altitude California coastal jet influenced by the North Pacific High, which typically occurs at heights of 300–400 m. The lidar buoy observations also validate LLJ representation in atmospheric models that estimate potential energy yield of offshore wind farms. The European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) was unsuccessful at identifying all observed LLJs at both buoy locations within the lowest 200 m. An extension of the National Renewable Energy Laboratory (NREL) 20-year wind resource dataset for the Outer Continental Shelf off the coast of California (CA20-Ext) yielded marginally greater captures of observed LLJs using the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary layer (PBL) scheme than the 2023 National Offshore Wind dataset (NOW-23), which uses the Yonsei University (YSU) scheme. However, CA20-Ext also produced the most LLJ false alarms, which are instances when a model identified an LLJ but no LLJ was observed. CA20-Ext and NOW-23 exhibited a tendency to overestimate the duration of LLJ events and underestimate LLJ core heights.

17 WIND ENERGY↗
DOE OSTI2025

Back Trajectories during SAIL

This data set contains air mass back trajectories generated during the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign using the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT, Stein et al. 2015, Rolph et al. 2017). For each hour from January 1, 2021 to June 30, 2023, a 96-hour back trajectory was initiated at the SAIL sampling site, starting 100 m above ground level. The calculations used the GDAS meteorological data set with model vertical velocity. The files within this dataset show location and meteorological parameters of the back trajectory analysis. There are 21 columns within each data file. The first 2 columns pertain to back trajectory parameters: the trajectory number (which is always 1) and the meteorological grid number. The next 9 columns pertain to spatiotemporal information: the year, month, day, hour (0-23), and minute of the trajectory, the previous time of the trajectory (in negative hours), and the airmass location (latitude, longitude, and height in m AGL). The rest of the parameters within each file pertains to surface and meteorological characteristics of the airmass: pressure (in hPa), potential temperature (in K), temperature (in K), precipitation rate (in mm/hr), mixing depth (in m), relative humidity (in % relative to liquid water), specific humidity (in g/kg), H2O mixing ratio (in g/kg), underlying terrain height (in meters above sea level), and incoming radiation (in W/m^2).

back trajectory↗
DOE OSTI2024

Implementing Superresolution of Nonstationary Tides with Wavelets: An Introduction to CWT_Multi

Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.

Lobo, Matthew↗
DOE OSTI2022

The role of föhn winds in eastern Antarctic Peninsula rapid ice shelf collapse

Ice shelf collapse reduces buttressing and enables grounded glaciers to contribute more rapidly to sea-level rise in a warming climate. The abrupt collapses of the Larsen A (1995) and B (2002) ice shelves on the Antarctic Peninsula (AP) occurred, at least for Larsen B, when long-period ocean swells damaged the calving front and the ice shelf was inundated with melt lakes that led to large-scale hydrofracture cascades. During collapse, field and satellite observations indicate föhn winds were present on both ice shelves. Here we use a regional climate model and machine learning analyses to evaluate the contributory roles of föhn winds and associated melt events prior to and during the collapses for ice shelves on the AP. Föhn winds caused about 25 % ± 3 % of the total annual melt in just 9 d on Larsen A prior to and during collapse and were present during the Larsen B collapse, which helped form extensive melt lakes. At the same time, the off-coast wind direction created by föhn winds helped melt and physically push sea ice away from the ice shelf calving fronts that allowed long-period ocean swells to reach and damage the front, which has been theorized to have ultimately triggered collapse. Collapsed ice shelves experienced enhanced surface melt driven by föhn winds over a large spatial extent and near the calving front, whereas SCAR inlet and the Larsen C ice shelves are affected less by föhn-wind-induced melt and do not experience large-scale melt ponds. These results suggest SCAR inlet and the Larsen C ice shelves may be less likely to experience rapid collapse due to föhn-driven melt so long as surface temperatures and föhn occurrence remain within historical bounds.

58 GEOSCIENCES↗
DOE OSTI2023

Cloud condensation nuclei activity of internally mixed particle populations at a remote marine free troposphere site in the North Atlantic Ocean

Here, this study reports results from research conducted at the Observatory of Mount Pico (OMP), 2225 m above mean sea level on Pico Island in the Azores archipelago in June and July 2017. We investigated the chemical composition, mixing state, and cloud condensation nuclei (CCN) activities of long-range transported free tropospheric (FT) particles. FLEXible PARTicle Lagrangian particle dispersion model (FLEXPART) simulations reveal that most air masses that arrived at the OMP during the sampling period originated in North America and were highly aged (average plume age > 10 days). We probed size-resolved chemical composition, mixing state, and hygroscopicity parameter ($\mathcal{κ}$) of individual particles using computer-controlled scanning electron microscopy with an energy-dispersive X-ray spectrometer (CCSEM-EDX). Based on the estimated individual particle mass from elemental composition, we calculated the mixing state index, $χ$. During our study, FT particle populations were internally mixed ($χ$ of samples are between 53% and 87%), owing to the long atmospheric aging time. We used data from a miniature Cloud Condensation Nucleus Counter (miniCCNC) to derive the hygroscopicity parameter, $\mathcal{κ}$ CCNC . Combining $\mathcal{κ}$ CCNC and FLEXPART, we found that air masses recirculated above the North Atlantic Ocean with lower mean altitude had higher $\mathcal{κ}$ CCNC due to the higher contribution of sea salt particles. We used CCSEM-EDX and phase state measurements to predict single-particle $\mathcal{κ}$ ($\mathcal{κ}$ CCSEM-EDX ) values, which overlap with the lower range of $\mathcal{κ}$ CCNC measured below 0.15 % SS. Therefore, CCSEM-EDX measurements can be useful in predicting the lower bound of $\mathcal{κ}$, which can be used in climate models to predict CCN activities, especially in remote locations where online CCN measurements are unavailable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Human‐induced salinity changes impact marine organisms and ecosystems

Abstract Climate change is fundamentally altering marine and coastal ecosystems on a global scale. While the effects of ocean warming and acidification on ecology and ecosystem functions and services are being comprehensively researched, less attention is directed toward understanding the impacts of human‐driven ocean salinity changes. The global water cycle operates through water fluxes expressed as precipitation, evaporation, and freshwater runoff from land. Changes to these in turn modulate ocean salinity and shape the marine and coastal environment by affecting ocean currents, stratification, oxygen saturation, and sea level rise. Besides the direct impact on ocean physical processes, salinity changes impact ocean biological functions with the ecophysiological consequences are being poorly understood. This is surprising as salinity changes may impact diversity, ecosystem and habitat structure loss, and community shifts including trophic cascades. Climate model future projections (of end of the century salinity changes) indicate magnitudes that lead to modification of open ocean plankton community structure and habitat suitability of coral reef communities. Such salinity changes are also capable of affecting the diversity and metabolic capacity of coastal microorganisms and impairing the photosynthetic capacity of (coastal and open ocean) phytoplankton, macroalgae, and seagrass, with downstream ramifications on global biogeochemical cycling. The scarcity of comprehensive salinity data in dynamic coastal regions warrants additional attention. Such datasets are crucial to quantify salinity‐based ecosystem function relationships and project such changes that ultimately link into carbon sequestration and freshwater as well as food availability to human populations around the globe. It is critical to integrate vigorous high‐quality salinity data with interacting key environmental parameters (e.g., temperature, nutrients, oxygen) for a comprehensive understanding of anthropogenically induced marine changes and its impact on human health and the global economy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ecohydrology of wetland plant communities along an estuarine to tidal river gradient

Abiotic filters that interact with wetland plant communities along tidal–fluvial gradients are highly dynamic, and understanding their quantitative thresholds and relationships to interspecific competition is important during an era of sea-level rise and watershed hydrologic change. Yet, landscape-scale studies of major coastal rivers from the river mouth to the head of tide, such as this study, remain rare. Here, we develop a new predictive framework for estuarine–tidal river research and management using a river-specific low-water datum and the wetland inundation indicator SEVg, the growing-season sum exceedance value of hourly surface-water depth. The distribution and variability of the wetland species pool (n = 203) on the 234 river kilometer (rkm) lower Columbia River and estuary floodplain are described for the first time. 4,940 quadrats at 50 marshes were surveyed (2005–2016). Throughout the estuarine–tidal river system, SEVg was well suited to describe the wetland inundation regime and its variability based on the combination of longitudinal river position and elevation. SEVg increased significantly landward. Two primary wetland inundation regimes were identified: the seaward-tidal, usually greatest during the winter months, and landward-fluvial, greatest during the growing season. Nearest the ocean, salinity is the abiotic factor limiting species richness and non-native species. Farther upriver, the daily, seasonal, and interannual variability of the wetting and drying cycle encourage disturbance-tolerant species and non-natives and limit the number of hydrophytes and total vegetative cover. Hence, the average between-year similarity of site-scale areal cover significantly decreased landward. Hierarchical cluster analysis indicated five vegetative groups and five ecohydrologic zones between rkm 0 and 234 were discriminated with 76 significant species–zone associations. All zones had unique indicator species. Species with high indicator values were Carex lyngbyei throughout the estuarine zones, and Eleocharis palustris, Sagittaria latifolia, and the invasive non-native Phalaris arundinacea in the upper estuarine and lower, middle, and upper tidal river zones (IV > 0.90). Competition from C. lyngbyei nearest the ocean and P. arundinacea in the tidal river was associated with reduced species richness when total cover was >65%. This framework of filters informs the design and prediction of future wetland plant communities on coastal river floodplains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Whole plant traits of coastal dune vegetation and implications for interactions with dune dynamics

Abstract Coastal dunes are important protective features against sea level rise and coastal storms. Interactions between dune plant aboveground structures and sediment trapping that allow for dune building and maintenance are well established. More recently, studies documenting belowground biomass for promoting erosion resistance in dominant dune species have been conducted, yet a knowledge gap remains regarding species‐specific characterization of whole plants, specifically with respect to roots, rhizomes, and belowground stems. Our objective was to quantify above‐ and belowground traits of four dominant dune grasses to document the potential for species‐specific effects on dune growth, maintenance, and erosion resistance. We examined above‐ and belowground traits among four prominent dune grasses of the Atlantic and Gulf Coasts of North America: Ammophila breviligulata , Panicum amarum , Spartina patens , and Uniola paniculata . Whole plant samples of each species were collected from the foredune at the US Army Engineer Research and Development Center's Field Research Facility in Duck, North Carolina, USA, and quantified for several above‐ and belowground traits (e.g., stem height, rhizome number and length, root surface area by diameter class, root tensile strength, and mycorrhizal percent infection). Belowground factors known to impact important dune processes, such as rhizome length, mycorrhizal percent infection, and root traits, differed substantially among species. When visualized in multivariate space, all species significantly differed in suites of above‐ and belowground traits. When considering belowground only, Ammophila and Spartina were similar, despite differences in biomass allocation. Species separated along axes related to mycorrhizal association, biomass allocation, and root construction. The four co‐occurring dune grass species were dissimilar in suites of plant traits. Belowground trait differences were driven by those describing root construction, biomass allocation, and mycorrhizal infection. Dissimilarity in above‐ and belowground suites of traits may demonstrate different approaches for surviving the dune environment. Incorporating belowground traits into modeling will enhance predictions of dune response to climate change through interactions between vegetation and dune dynamics that facilitate coastal resistance and resilience.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Reciprocal transplants reveal asymmetric local adaptation of Himalayan Rhododendron approaching elevational range limit

Abstract As plant species expand their upper limits of distribution under current warming, some retain both traditional climate space and biotic environment while others encounter novel conditions. The latter is the case for Rhododendron campanulatum , a woody shrub that grows both above and below treeline at our study site in the Eastern Himalayas where a very conspicuous, stable treeline was defined by a nearly contiguous canopy of tall Abies spectabilis trees, many of which are over a century old. Prior work showed that treeline had remained static in this region while R. campanulatum expanded its elevational range limit. We tested local adaptation of R. campanulatum by performing reciprocal transplants between the species' current elevational range limit (4023 m above sea level [asl]) and just above treeline (3876 m asl). Contrary to expectation, the coldest temperatures of late winter and early mid‐spring were experienced by plants at the lower elevation: R. campanulatum at species' limit (upper site) were covered by snow for a longer period (40 more days) and escaped the coldest temperatures suffered by conspecifics at treeline (lower site). The harsher spring conditions at treeline likely explain why leaves were smaller at treeline (15.3 cm 2 ) than at species limit (21.3 cm 2 ). Contrary to results from equivalent studies in other regions, survival was reduced more by downslope than by upslope movement, again potentially due to extreme cold temperatures observed at treeline in spring. Upslope transplantation had no effect on mortality, but mortality of species limit saplings transplanted downslope was three times higher than that of residents at both sites. A general expectation is that locals should survive better than foreign transplants, but survival of locals and immigrants at our species limit site was identical. However, those species limit saplings that survived the transplant to treeline grew faster than both locals at treeline and the transplants at species limit. Overall, we found asymmetric adaptation: Compared with treeline saplings, those at species limit (147 m above treeline) were more tolerant of extremes in the growing season but less tolerant of extremes in winter and early mid‐spring, displaying local adaptation in a more complex manner than simply home advantage, and complicating predictions about impacts of future regional climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Stream discharge and water quality data for East Fork Poplar Creek beginning 2012

This data note describes 15-min discharge and in situ water quality data at two locations along East Fork Poplar Creek in east Tennessee, USA. Data records include temperature, gauge height, water surface elevation above mean sea level, and volumetric discharge. Water quality measurements include temperature, specific conductance, pH, dissolved oxygen (percent saturation and concentration), turbidity, and less extensive fDOM data at one site. The data records begin in 2012 at one site and 2015 at the second site; monitoring at both sites is ongoing (as of 2021). The goal of this data collection is to improve understanding of watershed functions, hydrologic dynamics, and material flux. The data will contribute to site conceptual and numerical models, exposure and risk evaluation, remediation selection and design, and performance monitoring. The data are publicly available and can be accessed via unique url or DOI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Seasonal drivers of dissolved oxygen across a tidal creek–marsh interface revealed by machine learning

Abstract Dissolved oxygen (DO) is a key biogeochemical control in coastal systems, and its concentration and drivers vary markedly through time and space. This makes it difficult to accurately represent coastal DO and associated biogeochemical processes in models, limiting our ability to predict how these systems will respond to global change. We obtained high‐frequency (5‐min) in situ measurements of DO collected at three locations across the interface of a tidal creek and coastal marsh in the Pacific Northwest, USA. Random Forest machine learning models quantified the importance of three categories of environmental drivers (Aquatic, Climatic, and Terrestrial) of DO variability across the creek–marsh interface. We selected two 4‐month datasets representing Summer and Winter seasonal periods to test two hypotheses on the dominant drivers of DO at the coastal interface. We found that the Terrestrial driver—characterized by long periods of anaerobic conditions and episodic pulses in DO after floods—was most important during the Winter, whereas the Aquatic driver—characterized by variability over tidal, diel, and lunar cycles—was most important during the Summer. We explored how future climate change scenarios could alter the drivers of DO variability using a cumulative sums driver–response framework. Our results suggest that under climate change, Aquatic and Climatic drivers may increase in importance during the Summer, potentially linked to changing metabolic regimes and sea level, with Terrestrial driver importance potentially increasing during the Winter. Our approach highlights useful methods for understanding the spatiotemporal complexity of oxygen across coastal interfaces and quantifying the relative importance of distinct environmental drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Alternating salt and freshwater floods of coastal soils impact soil structure, hydraulic properties, and oxygen dynamics

Coastal soils are increasingly impacted by hydrologic intensification in the form of rising sea level and flooding from storm surges and precipitation. Alternating saltwater (SW) and freshwater (FW) exposure has the potential to disperse colloids, which can lead to disintegration of soil structure, clogging of pore spaces, and reduction in ecologically and biogeochemically important functions like infiltration or gas exchange. To investigate how hydrologic intensification affects soil structure and oxygen dynamics, we conducted a series of laboratory-based flood simulations. Intact soil cores of the A and B horizons (22 total) from the toe slope of an upland coastal forest along the western shore of Chesapeake Bay. We subjected the cores to 24 h of saturation with either FW or alternating brackish SW and FW, with a 24-h draining event in between flood events. Oxygen diffusion into soil during draining was reduced by up to 30% for soils flooded alternating SW-FW compared to soils flooded only with FW. We attributed the reduced oxygen diffusion to clogging of smaller pores by colloids, with colloid redistribution observed as an increase in specific surface area down the core profile of up to 59%. After three SW-FW floods (six floods total), there were significant changes in pore size distribution, significant redistribution of colloids, and the A horizon became sodic. We concluded that a small number of SW flooding events can induce a measurable change in soil physical properties that directly impacts the biogeochemical dynamics.

Rod, Kenton A. [Pacific Northwest National Laborat↗