Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cold climate”

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.

918 records · Page 51

Sun, the Earth, and Near-Earth Space: A Guide to the Sun-Earth System

In a world of warmth and light and living things we soon forget that we are surrounded by a vast universe that is cold and dark and deadly dangerous, just beyond our door. On a starry night, when we look out into the darkness that lies around us, the view can be misleading in yet another way: for the brightness and sheer number of stars, and their chance groupings into familiar constellations, make them seem much nearer to each other, and to us, that in truth they are. And every one of them--each twinkling, like a diamond in the sky--is a white-hot sun, much like our own. The nearest stars in our own galaxy--the Milky Way-- are more than a million times further away from us than our star, the Sun. We could make a telephone call to the Moon and expect to wait but a few seconds between pieces of a conversation, or but a few hours in calling any planet in our solar system.

Sun↗

Postglacial migration shaped the genomic diversity and global distribution of the wild ancestor of lager-brewing hybrids

The wild, cold-adapted parent of hybrid lager-brewing yeasts, Saccharomyces eubayanus, has a complex and understudied natural history. The exploration of this diversity can be used both to develop new brewing applications and to enlighten our understanding of the dynamics of yeast evolution in the wild. Here, we integrate whole genome sequence and phenotypic data of 200 S. eubayanus strains, the largest collection known to date. S. eubayanus has a multilayered population structure, consisting of two major populations that are further structured into six subpopulations. Four of these subpopulations are found exclusively in the Patagonian region of South America; one is found predominantly in Patagonia and sparsely in Oceania and North America; and one is specific to the Holarctic ecozone. Plant host associations differed between subpopulations and between S. eubayanus and its sister species, Saccharomyces uvarum. S. eubayanus is most abundant and genetically diverse in northern Patagonia, where some locations harbor more genetic diversity than is found outside of South America, suggesting that northern Patagonia east of the Andes was a glacial refugium for this species. All but one subpopulation shows isolation-by-distance, and gene flow between subpopulations is low. However, there are strong signals of ancient and recent outcrossing, including two admixed lineages, one that is sympatric with and one that is mostly isolated from its parental populations. Using our extensive biogeographical data, we build a robust model that predicts all known and a handful of additional regions of the globe that are climatically suitable for S. eubayanus, including Europe where host accessibility and competitive exclusion by other Saccharomyces species may explain its continued elusiveness. We conclude that this industrially relevant species has rich natural diversity with many factors contributing to its complex distribution and natural history.

59 BASIC BIOLOGICAL SCIENCES↗

Environment-sensitivity functions for gross primary productivity in light use efficiency models

The sensitivity of photosynthesis to environmental changes is essential for understanding carbon cycle responses to global climate change and for the development of modeling approaches that explains its spatial and temporal variability. We collected a large variety of published sensitivity functions of gross primary productivity (GPP) to different forcing variables to assess the response of GPP to environmental factors. These include the responses of GPP to temperature; vapor pressure deficit, some of which include the response to atmospheric CO 2 concentrations; soil water availability (W); light intensity; and cloudiness. These functions were combined in a full factorial light use efficiency (LUE) model structure, leading to a collection of 5600 distinct LUE models. Each model was optimized against daily GPP and evapotranspiration fluxes from 196 FLUXNET sites and ranked across sites based on a bootstrap approach. The GPP sensitivity to each environmental factor, including CO 2 fertilization, was shown to be significant, and that none of the previously published model structures performed as well as the best model selected. From daily and weekly to monthly scales, the best model's median Nash-Sutcliffe model efficiency across sites was 0.73, 0.79 and 0.82, respectively, but poorer at annual scales (0.23), emphasizing the common limitation of current models in describing the interannual variability of GPP. Although the best global model did not match the local best model at each site, the selection was robust across ecosystem types. The contribution of light saturation and cloudiness to GPP was observed across all biomes (from 23% to 43%). Temperature and W dominates GPP and LUE but responses of GPP to temperature and W are lagged in cold and arid ecosystems, respectively. The findings of this study provide a foundation towards more robust LUE-based estimates of global GPP and may provide a benchmark for other empirical GPP products.

54 ENVIRONMENTAL SCIENCES↗

Cirrus Susceptibility to Changes in Ice Nuclei: Physical Processes, Model Uncertainties, and Measurement Needs

One of the proposed concepts for mitigating the warming effect of increasing greenhouse gases is seeding cirrus cloud with ice nuclei (IN) in order to reduce the lifetime and coverage of cold cirrus that have a net warming impact on the earth's surface. Global model simulations of the net impact of changing upper tropospheric IN have given widely disparate results, partly as a result of poor understanding of ice nucleation processes in the current atmosphere, and partly as a result of poor representation of these processes in global models. Here, we present detailed process-model simulations of tropical tropopause layer (TTL) transport and cirrus formation with ice nuclei properties based on recent laboratory nucleation experiments and field measurements of aerosol composition. The model is used to assess the sensitivity of TTL cirrus occurrence frequency and microphysical properties to the abundance and efficacy of ice nuclei. The simulated cloud properties compared with recent high-altitude aircraft measurements of TTL cirrus and ice supersaturation. We find that abundant effective IN (either from glassy organic aerosols or crystalline ammonium sulfate with concentrations greater than about 100/L) prevent the occurrences of large ice concentration and large ice supersaturations, both of which are clearly indicated by the in situ observations. We find that concentrations of effective ice nuclei larger than about 50/L can drive significant changes in cirrus microphysical properties and occurrence frequency. However, the cloud occurrence frequency can either increase or decrease, depending on the efficacy and abundance of IN added to the TTL. We suggest that our lack of information about ice nuclei properties in the current atmosphere, as well as uncertainties in ice nucleation processes and their representations in global models, preclude meaningful estimates of climate impacts associated with addition of ice nuclei in the upper troposphere. We will briefly discuss the key field measurements needed to constrain ice nucleation processes.

Jensen, Eric↗

Evaluation of Aqua MODIS Thermal Emissive Bands Stability through Radiative Transfer Modeling

Moderate Resolution Imaging Spectroradiometer (MODIS) on Aqua has been inoperation providing continuous global observations for science research and applications since2002. The long-term stability of thermal emissive bands (TEBs) of Aqua MODIS was monitoredthrough inter-comparisons with measurements by hyperspectral or multi-spectral infraredsensors or through vicarious monitoring over cold targets such as Dome-C and deep convectiveclouds. The radiative transfer modeling (RTM)-based simulation models are developed to per-form the long-term monitoring of the stability of Aqua MODIS TEBs. Long-term EuropeanCentre for Medium-Range Weather Forecasts global atmospheric reanalysis data are used asinputs to the RTM simulation. By confining ocean in low to middle latitudes as the area ofinterest, the long-term stabilities of Aqua MODIS TEBs are monitored through observationminus background (O-B) brightness temperature (BT) bias between MODIS measurements andsimulations from the RTM. In general, the RTM-based O-B analysis shows that the long-termstability of Aqua MODIS TEBs are maintained well. It is shown that Aqua MODIS surfacebands are all radiometrically stable with yearly BT bias drifts <0.005 K∕yearfor B20, B22,and B23 and∼0.01 K∕yearfor B31 and B32. The carbon dioxide (CO2) absorption channelsof Aqua MODIS, e.g., B33 to B36, are stable with a BT bias drift <0.005 K∕year. It is also foundthat after accounting for the long-term growth of greenhouse gasN2O, the O-B bias trends forB24 and B25 are quite stable with yearly BT bias drifts around 0.0002 and0.0055 K∕year,respectively. The stabilities of the moisture-sensitive channels B27 to B29 and the ozone channelB30 of Aqua MODIS are also evaluated and the remnant small variations in the O-B bias trend-ing of these channels are further discussed. The RTM-based inter-comparison of MODIS TEBmeasurements with global atmospheric and climate re-analysis data provides validation of thelong-term radiometric stability of Aqua MODIS TEBs. The analysis in this paper also demon-strates that the comparison of well-calibrated MODIS TEB measurements with RTM simulationhelps quantify the long-term impacts of global greenhouse gas concentration growth on the BTdecreases in the MODIS greenhouse gas-sensitive channels.

MODIS↗

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗

Trends and meteorological drivers of extreme daily reservoir evaporation events in the western United States

Extreme daily evaporation from reservoir surfaces can lead to significant short-term water losses, affecting water quality, water supply, and reservoir operation strategies. Historical trends in daily reservoir evaporation events have eluded the scientific and operational communities, largely due to a lack of long-term, consistent data record. This study quantifies trends in extreme daily reservoir evaporation events at 165 major reservoirs located in the western U.S. Here, we use the place-based energy balance and aerodynamic Daily Lake Evaporation Model (DLEM) driven by multiple meteorological data products (RTMA, gridMET, Daymet) to estimate daily evaporation rates at these reservoirs from 1981 to 2022. The results—while are based on different meteorological forcing datasets—consistently indicate that the California, Lower Colorado, and Rio Grande hydrologic regions are more prone to higher daily evaporation extremes. Compared to the rest of western U.S, these three regions also experience a more pronounced increasing trend in the annual maximum daily evaporation rate, at about 0.3 mm day -1 decade -1 during 1981-2022. The results show that heat and dryness are the main drivers to the increasing trend of extreme evaporation, while extreme wind speed is the primary contributor to exceptionally high daily evaporation events across all regions. This phenomenon is particularly prominent in the arid Lower Colorado region, but less significant in the cold and humid Pacific Northwest region. We also find that the correlation between extreme wind speed and extreme evaporation degrades as the time scale increases from daily, to monthly and seasonal. Our findings have strong implications for the pattern and distribution of extreme evaporation events at the western U.S. reservoirs, and illustrate how various drivers influence extreme evaporation across different time scales.

13 HYDRO ENERGY↗

Changes in atmospheric circulation amplify extreme snowfall fueled by Arctic sea ice loss over high-latitude land

Arctic sea-ice retreat has been linked to increased winter precipitation and heavy snowfall over land, likely due to a combination of enhanced evaporation from ice-free Arctic marginal seas (AMS) and changes in atmospheric circulation. However, their relative roles and contributions remain uncertain. Here, we show that a greater proportion of AMS evaporative moisture reached high-latitude land during the cold seasons from 1980–1989 to 2012–2021. Atmospheric circulation changes added an additional 13 % increase in the AMS moisture contribution, accounting for 11 % of the total increase in AMS-sourced land precipitation. Notably, 46 % of the increase in AMS-sourced extreme snowfall is attributed to circulation-driven landward moisture transport, representing an 84 % increase beyond the effect of enhanced AMS evaporation alone. Further analysis indicates that both the rise in Arctic moisture and the atmospheric circulation shifts are primarily driven by anthropogenic forcing. These findings highlight how atmospheric circulation changes amplify extreme snowfall fueled by AMS evaporation, underscoring the synergistic effects of Arctic sea ice loss and circulation change on high-latitude winter precipitation.

Anthropogenic forcing↗

Sensing Water Vapon via Spacecraft Radio Occultation Observations

The radio occultation technique has been used to characterize planetary atmospheres since the 1960's spanning atmospheric pressures from 16 microbars to several bars. In 1988, the use of GPS signals to make occultation observations of Earth's atmosphere was realized by Tom Yunck and Gunnar Lindal at JPL. In the GPS to low-Earth-orbiter limb- viewing occultation geometry, Fresnel diffraction yield a unique combination of high vertical resolution of 100 m to 1 km at long wavelengths (approx. 20 cm) insensitive to particulate scattering which allows routine limb sounding from the lower mesosphere through the troposphere. A single orbiting GPS/GLONASS receiver can observe - 1000 to 1400 daily occultations providing as many daily, high vertical resolution soundings as the present global radiosonde network, but with far more evenly distributed, global coverage. The occultations yield profiles of refractivity as a function of height. In the cold, dry conditions of the upper troposphere and above (T less than 240 K), profiles of density, pressure (geopotential), and temperature can be derived. Given additional temperature information, water vapor can be derived in the midddle and lower troposphere with a unique combination of vertical resolution, global distribution and insensitivity to clouds and precipitation to an accuracy of approx. 0.2 g/kg. At low latitudes, moisture profiles will be accurate to 1-5% within the convective boundary layer and better than 20% below 6 to 7 km. Accuracies of climatological averages should be approx. 0. 1 g/kg limited by the biases in the temperature estimates. To use refractivity to constrain water vapor, knowledge of temperature is required. The simplest approach is to use the temperature field from an analysis such as the 6 hour ECMWF global analysis interpolated to the locations of each occultation. A better approach is to combine the temperature and moisture fields from such an analysis with the occultation refractivity in a weighting scheme based on the errors in each data field. A ID variational combinational approach has been developed at the UKMO. We win present results from both approaches from GPS/MET data taken in June and July 1995 and compare them with the ECMWF global 6 hour moisture analyses which are derived largely from TOVS and radiosonde data. Overall, the atmosphere below the 500 mb level appears somewhat drier in general than the ECNIWF humidity field. A 2-D (latitude vs. height) climatological snapshot derived from a 2-week span of GPS/MET data will be compared to the humidity climatology of Peixoto and Oort derived from radiosonde data from 1963-1973. Differences between the GPS results and Peixoto and Oort may be the signature of a climate trend over the past 30 years.

Kursinski, E. Robert↗

Precursors to Extreme Wintertime Cold in the Midwest United States

ABSTRACT Cold air outbreaks (CAOs) are extreme weather events that affect millions of people every winter, including those in the US Midwest. Per our criteria (a sufficiently cold event over a wide area and a long period), we identified 43 wintertime CAOs in the Midwest region using data from 1980 to 2021 (approximately one per year). These occurred in the Midwest during three of five North American weather regimes. Two regimes' CAOs are related to a weak stratospheric polar vortex, consistent with previous research on CAOs in the neighbouring Great Plains region. A different regime, characterised by anomalous ridging along the west coast of North America, is the most common for Midwest CAOs, unlike in regions further to the west. Unlike the other two regimes, these so‐called West Coast Ridge CAOs have no clear stratospheric connection. West Coast Ridge CAOs are instead linked to tropospheric processes on synoptic timescales, usually preceded by a weaker than average mid‐tropospheric height gradient in the western Pacific about 10 days prior to CAO onset. Our results demonstrate a particular challenge for predicting some extreme events in the Midwest, with important implications for early warning systems.

Ryan, James M. [Department of Earth and Atmospheri↗

Amplified Mesoscale and Submesoscale Variability and Increased Concentration of Precipitation under Global Warming over Western North America

Abstract Cold-season precipitation statistics in simulations from the storm-resolving WRF Model at 6-km and 1-h resolution over western North America are analyzed. Pseudo–global warming future simulations for the 2041–80 period, constrained by GCMs under the RCP8.5 scenario, are compared to the 1981–2020 historical simulation. The analysis focuses on the dynamical properties of precipitation time series at subdaily scales and on the morphology of storms. The statistical distribution of precipitation intensities in each pixel of the simulation domain is characterized through nonparametric statistical indicators: frequency of wet hours, mean wet-hour precipitation intensity, and Gini coefficient as a measure of the temporal concentration of the precipitation volume. Additionally, the temporal and spatial Fourier power spectra of precipitation time series and precipitation fields are analyzed. The half-power period (HPP) and half-power wavelength (HPW) are defined as spectral measures of the characteristic scales of precipitation’s temporal and spatial patterns. The results show statistically significant increases in the mean wet-hour precipitation intensity and in the Gini coefficient in 99% of the pixels, indicating that the seasonal precipitation volume becomes more concentrated within a smaller number of hours with higher precipitation intensity. The statistics of change in the frequency of wet hours are more contrasted across the simulation domain. The changes are also reflected in the power spectra, which show the spatial and temporal variability increasing proportionally more with finer spatial and temporal scales and the HPW and HPP decreasing. These projected changes are expected to have consequences, not only in terms of hydrologic impacts but also in terms of the predictability of precipitation patterns. Significance Statement The precipitation characteristics of winter storms over the western United States and southwestern Canada are analyzed in future climate simulations for the 2041–80 period. As compared to present-day climate, the most intense parts of the storms are projected to produce a higher rainfall volume, with increased concentration over smaller areas and shorter time intervals. The propensity of rainfall intensity to vary rapidly over time will be enhanced in the future according to the simulations. These model predictions imply an increased risk of rapid flooding in small basins. They also suggest that predicting several hours ahead the time and location at which a storm will produce maximum rainfall may become more challenging in the future.

Climate change↗

Microbial Ecology of Permafrost Soils: Populations, Processes, and Perspectives

Permafrost microbial research has flourished in the past decades, due in part to improvements in sampling and molecular techniques, but also the increased focus on the permafrost greenhouse gas feedback to climate change and other ecological processes in high latitude and alpine permafrost soils. Permafrost microorganisms are adapted to these extreme environments and remain active at low temperatures and when resources are limited. They are also an important component of global elemental cycles as they regulate organic matter turnover and greenhouse gas production, particularly as permafrost thaws. Here we review the permafrost microbiology literature coupled with an exploration of its historical aspects, with a particular focus on a new understanding advanced by molecular biology techniques. We further identify knowledge gaps and ways forward to improve our understanding of microbial contributions to ecosystem biogeochemistry of permafrost-affected systems.

54 ENVIRONMENTAL SCIENCES↗

Robust but weak winter atmospheric circulation response to future Arctic sea ice loss

The possibility that Arctic sea ice loss weakens mid-latitude westerlies, promoting more severe cold winters, has sparked more than a decade of scientific debate, with apparent support from observations but inconclusive modelling evidence. Here we show that sixteen models contributing to the Polar Amplification Model Intercomparison Project simulate a weakening of mid-latitude westerlies in response to projected Arctic sea ice loss. We develop an emergent constraint based on eddy feedback, which is 1.2 to 3 times too weak in the models, suggesting that the real-world weakening lies towards the higher end of the model simulations. Still, the modelled response to Arctic sea ice loss is weak: the North Atlantic Oscillation response is similar in magnitude and offsets the projected response to increased greenhouse gases, but would only account for around 10% of variations in individual years. We further find that relationships between Arctic sea ice and atmospheric circulation have weakened recently in observations and are no longer inconsistent with those in models.

54 ENVIRONMENTAL SCIENCES↗

Vertical Moist Thermodynamic Structure and Spatial-Temporal Evolution of the MJO in AIRS Observations

The atmospheric moisture and temperature profiles from the Atmospheric Infrared Sounder (AIRS)/Advanced Microwave Sounding Unit on the NASA Aqua mission, in combination with the precipitation from the Tropical Rainfall Measuring Mission (TRMM), are employed to study the vertical moist thermodynamic structure and spatial-temporal evolution of the Madden-Julian oscillation (MJO). The AIRS data indicate that, in the Indian Ocean and western Pacific, the temperature anomaly exhibits a trimodal vertical structure: a warm (cold) anomaly in the free troposphere (800-250 hPa) and a cold (warm) anomaly near the tropopause (above 250 hPa) and in the lower troposphere (below 800 hPa) associated with enhanced (suppressed) convection. The AIRS moisture anomaly also shows markedly different vertical structures as a function of longitude and the strength of convection anomaly. Most significantly, the AIRS data demonstrate that, over the Indian Ocean and western Pacific, the enhanced (suppressed) convection is generally preceded in both time and space by a low-level warm and moist (cold and dry) anomaly and followed by a low-level cold and dry (warm and moist) anomaly. The MJO vertical moist thermodynamic structure from the AIRS data is in general agreement, particularly in the free troposphere, with previous studies based on global reanalysis and limited radiosonde data. However, major differences in the lower-troposphere moisture and temperature structure between the AIRS observations and the NCEP reanalysis are found over the Indian and Pacific Oceans, where there are very few conventional data to constrain the reanalysis. Specifically, the anomalous lower-troposphere temperature structure is much less well defined in NCEP than in AIRS for the western Pacific, and even has the opposite sign anomalies compared to AIRS relative to the wet/dry phase of the MJO in the Indian Ocean. Moreover, there are well-defined eastward-tilting variations of moisture with height in AIRS over the central and eastern Pacific that are less well defined, and in some cases absent, in NCEP. In addition, the correlation between MJO-related mid-tropospheric water vapor anomalies and TRMM precipitation anomalies is considerably more robust in AIRS than in NCEP, especially over the Indian Ocean. Overall, the AIRS results are quite consistent with those predicted by the frictional Kelvin-Rossby wave/conditional instability of the second kind (CISK) theory for the MJO.

troposhere↗

A regional nuclear conflict would compromise global food security

A limited nuclear war between India and Pakistan could ignite fires large enough to emit more than 5 Tg of soot into the stratosphere. Climate model simulations have shown severe resulting climate perturbations with declines in global mean temperature by 1.8 °C and precipitation by 8%, for at least 5 y. Here we evaluate impacts for the global food system. Six harmonized state-of-the-art crop models show that global caloric production from maize, wheat, rice, and soybean falls by 13 (±1)%, 11 (±8)%, 3 (±5)%, and 17 (±2)% over 5 y. Total single-year losses of 12 (±4)% quadruple the largest observed historical anomaly and exceed impacts caused by historic droughts and volcanic eruptions. Colder temperatures drive losses more than changes in precipitation and solar radiation, leading to strongest impacts in temperate regions poleward of 30°N, including the United States, Europe, and China for 10 to 15 y. Integrated food trade network analyses show that domestic reserves and global trade can largely buffer the production anomaly in the first year. Persistent multiyear losses, however, would constrain domestic food availability and propagate to the Global South, especially to food-insecure countries. By year 5, maize and wheat availability would decrease by 13% globally and by more than 20% in 71 countries with a cumulative population of 1.3 billion people. In view of increasing instability in South Asia, this study shows that a regional conflict using <1% of the worldwide nuclear arsenal could have adverse consequences for global food security unmatched in modern history.

54 ENVIRONMENTAL SCIENCES↗

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗

Impact of the Arctic Oscillation from March on summertime sea ice

Current understanding of the cold season Arctic oscillation (AO) impact on the summertime sea ice is revisited in this study by analyzing the role from each month. Earlier studies examined the prolonged AO impact using a smooth average over 1–2 seasons (e.g. December–March, December–April, March–May), ignoring large month-to-month AO variability. This study finds that the March AO is most influential on the summertime sea ice loss. First, the March AO is most highly negative-correlated with the AO in summer. Secondly, surface energy budget, sea level pressure, and low-tropospheric circulation exhibit that their time-lagged responses to the positive (negative) phase of the March AO grow with time, transitioning to the patterns associated with the negative (positive) phase of the AO that induces sea ice decrease (increase) in summer. Time evolution of the surface energy budget explains the growth of the sea ice concentration anomaly in summer, and a warming-to-cooling transition in October. The regional difference in sea ice anomaly distribution can be also explained by circulation and surface energy budget patterns. The sea ice concentration along the pan-Arctic including the Laptev, East Siberian, Chukchi, and Beaufort Sea decreases (increases) in summer in response to the positive (negative) phase of the March AO, while the sea ice to the northeast of Greenland increases (decreases). This sea ice response is better represented by the March AO than by the seasonally averaged winter AO, suggesting that the March AO can play more significant role. This study also finds that the sea ice decrease in response to the positive AO is distinctively smaller in the 20th century than in the 21st century, along with the opposite sea ice response over the Canada Basin due to circulation difference between the two periods.

Arctic Oscillation↗

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗