Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “winter”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A full year of aerosol size distribution data from the central Arctic under an extreme positive Arctic Oscillation: insights from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition

Abstract. The Arctic environment is rapidly changing due to accelerated warming in the region. The warming trend is driving a decline in sea ice extent, which thereby enhances feedback loops in the surface energy budget in the Arctic. Arctic aerosols play an important role in the radiative balance and hence the climate response in the region, yet direct observations of aerosols over the Arctic Ocean are limited. In this study, we investigate the annual cycle in the aerosol particle number size distribution (PNSD), particle number concentration (PNC), and black carbon (BC) mass concentration in the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. This is the first continuous, year-long data set of aerosol PNSD ever collected over the sea ice in the central Arctic Ocean. We use a k-means cluster analysis, FLEXPART simulations, and inverse modeling to evaluate seasonal patterns and the influence of different source regions on the Arctic aerosol population. Furthermore, we compare the aerosol observations to land-based sites across the Arctic, using both long-term measurements and observations during the year of the MOSAiC expedition (2019–2020), to investigate interannual variability and to give context to the aerosol characteristics from within the central Arctic. Our analysis identifies that, overall, the central Arctic exhibits typical seasonal patterns of aerosols, including anthropogenic influence from Arctic haze in winter and secondary aerosol processes in summer. The seasonal pattern corresponds to the global radiation, surface air temperature, and timing of sea ice melting/freezing, which drive changes in transport patterns and secondary aerosol processes. In winter, the Norilsk region in Russia/Siberia was the dominant source of Arctic haze signals in the PNSD and BC observations, which contributed to higher accumulation-mode PNC and BC mass concentrations in the central Arctic than at land-based observatories. We also show that the wintertime Arctic Oscillation (AO) phenomenon, which was reported to achieve a record-breaking positive phase during January–March 2020, explains the unusual timing and magnitude of Arctic haze across the Arctic region compared to longer-term observations. In summer, the aerosol PNCs of the nucleation and Aitken modes are enhanced; however, concentrations were notably lower in the central Arctic over the ice pack than at land-based sites further south. The analysis presented herein provides a current snapshot of Arctic aerosol processes in an environment that is characterized by rapid changes, which will be crucial for improving climate model predictions, understanding linkages between different environmental processes, and investigating the impacts of climate change in future Arctic aerosol studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Long-range transported continental aerosol in the eastern North Atlantic: three multiday event regimes influence cloud condensation nuclei

Abstract. The eastern North Atlantic (ENA) is a region dominated by pristine marine environment and subtropical marine boundary layer clouds. Under unperturbed atmospheric conditions, the regional aerosol regime in the ENA varies seasonally due to different seasonal surface-ocean biogenic emissions, removal processes, and meteorological regimes. However, during periods when the marine boundary layer aerosol in the ENA is impacted by particles transported from continental sources, aerosol properties within the marine boundary layer change significantly, affecting the concentration of cloud condensation nuclei (CCN). Here, we investigate the impact of long-range transported continental aerosol on the regional aerosol regime in the ENA using data collected at the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility on Graciosa Island in 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. We develop an algorithm that integrates number concentrations of particles with optical particle dry diameter (Dp) between 100 and 1000 nm, single scattering albedo, and black carbon concentration to identify multiday events (with duration >24 consecutive hours) of long-range continental aerosol transport in the ENA. In 2017, we detected nine multiday events of long-range transported particles that correspond to ∼ 7.5 % of the year. For each event, we perform HYSPLIT 10 d backward trajectories analysis, and we evaluate CALIPSO aerosol products to assess, respectively, the origins and compositions of aerosol particles arriving at the ENA site. Subsequently, we group the events into three categories, (1) mixture of dust and marine aerosols, (2) mixture of marine and polluted continental aerosols from industrialized areas, and (3) biomass burning aerosol from North America and Canada, and we evaluate their influence on aerosol population and cloud condensation nuclei in terms of potential activation fraction and concentrations at supersaturation of 0.1 % and 0.2 %. The arrival of plumes dominated by the mixture of dust and marine aerosol in the ENA in the winter caused significant increases in baseline Ntot. Simultaneously, the baseline particle size modes and CCN potential activation fraction remained almost unvaried, while cloud condensation nuclei concentrations increased proportionally to Ntot. Events dominated by a mixture of marine and polluted continental aerosols in spring, fall, and winter led to a statistically significant increase in baseline Ntot, a shift towards larger particular sizes, a higher CCN potential activation fractions, and cloud condensation nuclei concentrations of >170 % and up to 240 % higher than during baseline regime. Finally, the transported aerosol plumes characterized by elevated concentration of biomass burning aerosol from continental wildfires detected in the summertime did not statistically contribute to increase baseline aerosol particle concentrations in the ENA. However, particle diameters were larger than under baseline conditions, and CCN potential activation fractions were >75 % higher. Consequentially, cloud concentration nuclei concentrations increased by ∼ 115 % during the period affected by the biomass burning events. Our results suggest that, through the year, multiday events of long-range continental aerosol transport periodically affect the ENA and represent a significant source of CCN in the marine boundary layer. Based on our analysis, in 2017, the multiday aerosol plume transport dominated by a mixture of dust and marine aerosol, a mixture of marine and polluted continental aerosols, and biomass burning aerosols caused increases in the NCCN baseline regime of, respectively, 6.6 %, 8 %, and 7.4 % at SS 0.1 % (and, respectively, 6.5 %, 8.2 %, and 7.3 % at SS 0.2 %) in the ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Physicochemical characterization of free troposphere and marine boundary layer ice-nucleating particles collected by aircraft in the eastern North Atlantic

Abstract. Atmospheric ice nucleation impacts the hydrological cycle and climate by modifying the radiative properties of clouds. To improve our predictive understanding of ice formation, ambient ice-nucleating particles (INPs) need to be collected and characterized. Measurements of INPs at lower latitudes in a remote marine region are scarce. The Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign, in the region of the Azores islands, provided the opportunity to collect particles in the marine boundary layer (MBL) and free troposphere (FT) by aircraft during the campaign's summer and winter intensive operation period. The particle population in samples collected was examined by scanning transmission X-ray microscopy with near-edge X-ray absorption fine structure spectroscopy. The identified INPs were analyzed by scanning electron microscopy with energy-dispersive X-ray analysis. We observed differences in the particle population characteristics in terms of particle diversity, mixing state, and organic volume fraction between seasons, mostly due to dry intrusion events during winter, as well as between the sampling locations of the MBL and FT. These differences are also reflected in the temperature and humidity conditions under which water uptake, immersion freezing (IMF), and deposition ice nucleation (DIN) proceed. Identified INPs reflect typical particle types within the particle population on the samples and include sea salt, sea salt with sulfates, and mineral dust, all associated with organic matter, as well as carbonaceous particles. IMF and DIN kinetics are analyzed with respect to heterogeneous ice nucleation rate coefficients, Jhet, and ice nucleation active site density, ns, as a function of the water criterion Δaw. DIN is also analyzed in terms of contact angles following classical nucleation theory. Derived MBL IMF kinetics agree with previous ACE-ENA ground-site INP measurements. FT particle samples show greater ice nucleation propensity compared to MBL particle samples. This study emphasizes that the types of INPs can vary seasonally and with altitude depending on sampling location, thereby showing different ice nucleation propensities, which is crucial information when representing mixed-phase cloud and cirrus cloud microphysics in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Constraints on simulated past Arctic amplification and lapse rate feedback from observations

The Arctic has warmed more rapidly than the global mean during the past few decades. The lapse rate feedback (LRF) has been identified as being a large contributor to the Arctic amplification (AA) of climate change. This particular feedback arises from the vertically non-uniform warming of the troposphere, which in the Arctic emerges as strong near-surface and muted free-tropospheric warming. Stable stratification and meridional energy transport are two characteristic processes that are evoked as causes for this vertical warming structure. Our aim is to constrain these governing processes by making use of detailed observations in combination with the large climate model ensemble of the sixth Coupled Model Intercomparison Project (CMIP6). We build on the result that CMIP6 models show a large spread in AA and Arctic LRF, which are positively correlated for the historical period of 1951–2014. Thus, we present process-oriented constraints by linking characteristics of the current climate to historical climate simulations. In particular, we compare a large consortium of present-day observations to co-located model data from subsets that show a weak and strong simulated AA and Arctic LRF in the past. Our analyses suggest that the vertical temperature structure of the Arctic boundary layer is more realistically depicted in climate models with weak (w) AA and Arctic LRF (CMIP6/w) in the past. In particular, CMIP6/w models show stronger inversions in the present climate for boreal autumn and winter and over sea ice, which is more consistent with the observations. These results are based on observations from the year-long Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic, long-term measurements at the Utqiagvik site in Alaska, USA, and dropsonde temperature profiling from aircraft campaigns in the Fram Strait. In addition, the atmospheric energy transport from lower latitudes that can further mediate the warming structure in the free troposphere is more realistically represented by CMIP6/w models. In particular, CMIP6/w models systemically simulate a weaker Arctic atmospheric energy transport convergence in the present climate for boreal autumn and winter, which is more consistent with fifth generation reanalysis of the European Centre for Medium-Range Weather Forecasts (ERA5). We further show a positive relationship between the magnitude of the present-day transport convergence and the strength of past AA. With respect to the Arctic LRF, we find links between the changes in transport pathways that drive vertical warming structures and local differences in the LRF. This highlights the mediating influence of advection on the Arctic LRF and motivates deeper studies to explicitly link spatial patterns of Arctic feedbacks to changes in the large-scale circulation.

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

The correlation between Arctic sea ice, cloud phase and radiation using A-Train satellites

Abstract. Climate warming has a stronger impact on Arctic climate and sea ice cover (SIC) decline than previously thought. Better understanding and characterization of the relationship between sea ice and clouds and the implications for surface radiation is key to improving our confidence in Arctic climate projections. Here we analyze the relationship between sea ice, cloud phase and surface radiation over the Arctic, defined as north of 60° N, using active- and passive-sensor satellite observations from three different datasets. We find that all datasets agree on the climatology of and seasonal variability in total and liquid-bearing (liquid and mixed-phase) cloud covers. Similarly, our results show a robust relationship between decreased SIC and increased liquid-bearing clouds in the lowest levels (below 3 km) for all seasons (strongest in winter) but summer, while increased SIC and ice clouds are positively correlated in two of the three datasets. A refined map correlation analysis indicates that the relationship between SIC and liquid-bearing clouds can change sign over the Bering, Barents and Laptev seas, likely because of intrusions of warm air from low latitudes during winter and spring. Finally, the increase in liquid clouds resulting from decreasing SIC is associated with enhanced radiative cooling at the surface. Our findings indicate that the newly formed liquid clouds reflect more shortwave (SW) radiation back to space compared to the surface, generating a cooling effect of the surface, while their downward longwave (LW) radiation is similar to the upward LW surface emission, which has a negligible radiative impact on the surface. This overall cooling effect should contribute to dampening future Arctic surface warming as SIC continues to decline.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Application of PRIM for understanding patterns in carbon dioxide model-observation differences

Reducing uncertainties in regional carbon balances requires a better understanding of CO 2 transport in synoptic weather systems. Here, we apply the Patient Rule Induction Method (PRIM), a data-mining method to identify high-density regions for a target-class within an input parameter space, to airborne observations of potential temperature, wind speed, water vapor mixing ratio, and CO 2 dry mol fraction gathered during the Atmospheric Carbon and Transport (ACT)-America Summer 2016 and Winter 2017 campaigns. ACT observations were targeted at expert-designated cases of fair weather and near-frontal warm and cold sector air at atmospheric boundary-layer, lower-, and higher free tropospheric levels (ABL, LFT, and HFT, respectively). We investigate atmospheric characteristics of these pre-defined cases and associated CO 2 model-observation-differences in the mesoscale WRF-Chem model. PRIM results separate winter- and summertime observations as well as observations from ABL, LFT, and HFT with enrichment factors of 4.0–20.5 inside the PRIM box compared to the entire dataset but cannot distinguish between near-frontal warm and cold sector observations in the higher free troposphere. Analyzing of the parameter space constrained by PRIM, we find that large magnitude model observation differences preferentially associated with times when atmospheric conditions are less typical. This association suggests that PRIM could provide a useful tool for isolating atmospheric conditions with large-magnitude and non-Gaussian CO 2 -residuals for targeted transport model evaluation and to potentially improve inversion results during synoptically active periods.

Gerken, Tobias [James Madison Univ., Harrisonburg,↗
DOE OSTI2023

Evaluation of four ground-based retrievals of cloud droplet number concentration in marine stratocumulus with aircraft in situ measurements

Abstract. Cloud droplet number concentration (Nd) is crucial for understanding aerosol–cloud interactions (ACI) and associated radiative effects. We present evaluations of four ground-based Nd retrievals based on comprehensive datasets from the Atmospheric Radiation Measurement (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign. The Nd retrieval methods use ARM ENA observatory ground-based remote sensing observations from a micropulse lidar, Raman lidar, cloud radar, and the ARM NDROP (Droplet Number Concentration) value-added product (VAP), all of which also retrieve cloud effective radius (re). The retrievals are compared against aircraft measurements from the fast cloud droplet probe (FCDP) and the cloud and aerosol spectrometer (CAS) obtained from low-level marine boundary layer clouds on 12 flight days during summer and winter seasons. Additionally, the in situ measurements are used to validate the assumptions and characterizations used in the retrieval algorithms. Statistical comparisons of the probability distribution function (PDF) of the Nd and cloud re retrievals with aircraft measurements demonstrate that these retrievals align well with in situ measurements for overcast clouds, but they may substantially differ for broken clouds or clouds with low liquid water path (LWP). The retrievals are applied to 4 years of ground-based remote sensing measurements of overcast marine boundary layer clouds at the ARM ENA observatory to find that Nd (re) values exhibit seasonal variations, with higher (lower) values during the summer season and lower (higher) values during the winter season. The ensemble of various retrievals using different measurements and retrieval algorithms such as those in this paper can help to quantify Nd retrieval uncertainties and identify reliable Nd retrieval scenarios. Of the retrieval methods, we recommend using the micropulse lidar-based method. This method has good agreement with in situ measurements, less sensitivity to issues arising from precipitation and low cloud LWP and/or optical depth, and broad applicability by functioning for both daytime and nighttime conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of Four Ground-based Retrievals of Cloud Droplet Number Concentration in Marine Stratocumulus with Aircraft In Situ Measurements

Cloud droplet number concentration (N d ) is crucial for understanding aerosol-cloud interactions (ACI) and associated radiative effects. We present evaluations of four ground-based N d retrievals based on comprehensive datasets from the Atmospheric Radiation Measurements (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign. The N d retrieval methods use ARM ENA observatory ground-based remote sensing observations from a Micropulse lidar, Raman lidar, cloud radar, and the ARM NDROP Value-added Product (VAP), all of which also retrieve cloud effective radius (r e ). The retrievals are compared against aircraft measurements from the Fast-Cloud Droplet Probe (FCDP) and the Cloud and Aerosol Spectrometer (CAS) obtained from low-level marine boundary layer clouds on 12 flight days during summer and winter seasons. Additionally, the in situ measurements are used to validate the assumptions and characterizations used in the retrieval algorithms. Statistical comparisons of the probability distribution function (PDF) of the N d and cloud r e retrievals with aircraft measurements demonstrate that these retrievals align well with in situ measurements for overcast clouds, but they may substantially differ for broken clouds or clouds with low liquid water path (LWP). The retrievals are applied to four years of ground-based remote sensing measurements of overcast marine boundary layer clouds at the ARM ENA observatory to find that N d (r e ) values exhibit seasonal variations, with higher (lower) values during the summer season and lower (higher) values during the winter season. The ensemble of various retrievals using different measurements and retrieval algorithms such as those in this paper can help to quantify N d retrieval uncertainties and identify reliable N d retrieval scenarios. Of the retrieval methods, we recommend using the using the Micropulse lidar-based method given its good agreement with in situ measurements, it has less sensitivity to issues arising from precipitation and low cloud LWP/optical depth, and it has broad applicability by functioning for both day and nighttime conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating fine-resolution, regional outputs of a variable resolution global climate model

Climate models have been used to study water resources and regional hydrologic responses to climate change, but climate model outputs must be downscaled to provide relevant regional data. However, the accuracy of this regional data is limited by uncertainties across and within downscaling methods, uncertainty across global outputs, and discontinuities at downscaled boundaries. A new alternative to traditional downscaling is a variable resolution model that incorporates fine-resolution regions directly into a coarse-resolution, global climate simulation in order to capture contiguous dynamics across resolution boundaries. In this study, we used the Variable-Resolution Community Earth System Model (VR-CESM) to generate one-eighth degree (14 km) fine-resolution outputs for the western U.S. and eastern China from 1970-2006. We focus our evaluation on precipitaiton, temperature, snow pack, solar radiation, and wind. We compare the model outputs with remote-sensing-based precipitation data, and both reanalysis and gridded weather station data for precipitation and temperature. VR-CESM generally has a cold bias in winter and a warm bias in summer in the western U.S., which compensate each other to reduce the annual bias. In eastern China, however, the sign of temperature biases are more consistent throughout the year with cold biases in the higher mountains and warm biases throughout most of the rest of the region. Precipitation biases are dependent upon reference data, and show slight overestimation in high mountain regions in both the U.S. and China with respect to gridded weather station data. Simulated snow cover in the western U.S. is reasonable compared to remote sensing data, but snow cover and snow water equivalent have larger biases when compared to reanalysis data. In eastern China there are widespread snow cover biases compared to remote sensing data. VR-CESM underestimates downward shortwave radiation to a greater degree in summer than in winter, and underestimates surface layer windspeed over mountains to a greater degree than in other areas. Comparison between VR-CESM and a coarser simulation (1-degree Beijing Climate Center model) shows reduced precipitation biases in the mountainous regions with finer resolution, indicating the value of variable-resolution modeling for reigonal studies.

Di Vittorio, Alan↗
DOE OSTI2022

Atmospheric river representation in the Energy Exascale Earth System Model (E3SM) version 1.0

Abstract. The Energy Exascale Earth System Model (E3SM) project is an ongoing, state-of-the-science Earth system modeling, simulation, and prediction project developed by the US Department of Energy (DOE). With an emphasis on supporting the DOE's energy mission, understanding and quantifying how well the model simulates water cycle processes is of particular importance. Here, we evaluate E3SM version 1.0 (v1.0) for its ability to represent atmospheric rivers (ARs), which play significant roles in water vapor transport and precipitation. The characteristics and precipitation associated with global ARs in E3SM at standard resolution (1∘ × 1∘) are compared to the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA2). Global patterns of AR frequencies in E3SM show high degrees of correlation (≥0.97) with MERRA2 and low mean absolute errors (MAEs; <1 %) annually, seasonally, and across different ensemble members. However, some large-scale condition biases exist, leading to AR biases – most significant of which are the double intertropical convergence zone (ITCZ), a stronger and/or equatorward-shifted subtropical jet during boreal and austral winters, and enhanced Northern Hemisphere westerlies during summer. By comparing atmosphere-only and fully coupled simulations, we attribute the sources of the biases to the atmospheric component or to a coupling response. Using relationships revealed in Dong et al. (2021), we provide evidence showing the stronger North Pacific jet in winter and the enhanced Northern Hemisphere westerlies during summer, associated with E3SM's double ITCZ and related weaker Atlantic meridional overturning circulation (AMOC), respectively, which are significant sources of the AR biases found in the coupled simulations.

58 GEOSCIENCES↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
DOE OSTI2022

Inclusion of a cold hardening scheme to represent frost tolerance is essential to model realistic plant hydraulics in the Arctic–boreal zone in CLM5.0-FATES-Hydro

As temperatures decrease in autumn, vegetation of temperate and boreal ecosystems increases its tolerance to freezing. This process, known as hardening, results in a set of physiological changes at the molecular level that initiate modifications of cell membrane composition and the synthesis of anti-freeze proteins. Together with the freezing of extracellular water, anti-freeze proteins reduce plant water potentials and xylem conductivity. To represent the responses of vegetation to climate change, land surface schemes increasingly employ “hydrodynamic” models that represent the explicit fluxes of water from soil and through plants. The functioning of such schemes under frozen soil conditions, however, is poorly understood. Nonetheless, hydraulic processes are of major importance in the dynamics of these systems, which can suffer from, e.g., winter “frost drought” events. In this study, we implement a scheme that represents hardening into CLM5.0-FATES-Hydro. FATES-Hydro is a plant hydrodynamics module in FATES, a cohort model of vegetation physiology, growth, and dynamics hosted in CLM5.0. We find that, in frozen systems, it is necessary to introduce reductions in plant water loss associated with hardening to prevent winter desiccation. This work makes it possible to use CLM5.0-FATES-Hydro to model realistic impacts from frost droughts on vegetation growth and photosynthesis, leading to more reliable projections of how northern ecosystems respond to climate change.

58 GEOSCIENCES↗
DOE OSTI2022

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of E3SM land model snow simulations over the western United States

Abstract. Seasonal snow has crucial impacts on climate, ecosystems, and humans, but it is vulnerable to global warming. The land component (ELM) of the Energy Exascale Earth System Model (E3SM) mechanistically simulates snow processes from accumulation, canopy interception, compaction, and snow aging to melt. Although high-quality field measurements, remote sensing snow products, and data assimilation products with high spatio-temporal resolution are available, there has been no systematic evaluation of the snow properties and phenology in ELM. This study comprehensively evaluates ELM snow simulations over the western United States at 0.125∘ resolution during 2001–2019 using the Snow Telemetry (SNOTEL) in situ networks, MODIS remote sensing products (i.e., MCD43 surface albedo product), the spatially and temporally complete (STC) snow-covered area and grain size (MODSCAG) and MODIS dust and radiative forcing in snow (MODDRFS) products (STC-MODSCAG/STC-MODDRFS), and the snow property inversion from remote sensing (SPIReS) product and two data assimilation products of snow water equivalent and snow depth – i.e., University of Arizona (UA) and SNOw Data Assimilation System (SNODAS). Overall the ELM simulations are consistent with the benchmarking datasets and reproduce the spatio-temporal patterns, interannual variability, and elevation gradients for different snow properties including snow cover fraction (fsno), surface albedo (αsur) over snow cover regions, snow water equivalent (SWE), and snow depth (Dsno). However, there are large biases of fsno with dense forest cover and αsur in the Rocky Mountains and Sierra Nevada in winter, compared to the MODIS products. There are large discrepancies of snow albedo, snow grain size, and light-absorbing particle-induced snow albedo reduction between ELM and the MODIS products, attributed to uncertainties in the aerosol forcing data, snow aging processes in ELM, and remote sensing retrievals. Against UA and SNODAS, ELM has a mean bias of −20.7 mm (−35.9 %) and −20.4 mm (−35.5 %), respectively, for spring, and −13.8 mm (−27.8 %) and −10.2 mm (−22.2 %), respectively, for winter. ELM shows a relatively high correlation with SNOTEL SWE, with mean correlation coefficients of 0.69 but negative mean biases of −122.7 mm. Compared to the snow phenology of STC-MODSCAG and SPIReS, ELM shows delayed snow accumulation onset dates by 17.3 and 12.4 d, earlier snow end dates by 35.5 and 26.8 d, and shorter snow durations by 52.9 and 39.5 d, respectively. This study underscores the need for diagnosing model biases and improving ELM representations of snow properties and snow phenology in mountainous areas for more credible simulation and future projection of mountain snowpack.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Variations in boundary layer stability across Antarctica: a comparison between coastal and interior sites

The range of boundary layer stability profiles, from the surface to 500 m a.g.l. (above ground level), present in radiosonde observations from two continental-interior (South Pole Station and Dome Concordia Station) and three coastal (McMurdo Station, Georg von Neumayer Station III, and Syowa Station) Antarctic sites, is examined using the self-organizing maps (SOMs) neural network algorithm. A wide range of potential temperature profiles is revealed, from shallow boundary layers with strong near-surface stability to deeper boundary layers with weaker or near-neutral stability, as well as profiles with weaker near-surface stability and enhanced stability aloft, above the boundary layer. Boundary layer regimes were defined based on the range of profiles revealed by the SOM analysis; 20 boundary layer regimes were identified to account for differences in stability near the surface as well as above the boundary layer. Strong, very strong, or extremely strong stability, with vertical potential temperature gradients of 5 to in excess of 30 K per 100 m, occurred more than 80 % of the time at South Pole and Dome Concordia in the winter. Weaker stability was found in the winter at the coastal sites, with moderate and strong stability (vertical potential temperature gradients of 1.75 to 15 K per 100 m) occurring 70 % to 85 % of the time. Even in the summer, moderate and strong stability is found across all five sites, either immediately near the surface or aloft, just above the boundary layer. While the mean boundary layer height at the continental-interior sites was found to be approximately 50 m, the mean boundary layer height at the coastal sites was deeper, around 110 m. Further, a commonly described two-stability-regime system in the Arctic associated with clear or cloudy conditions was applied to the 20 boundary layer regimes identified in this study to understand if the two-regime behavior is also observed in the Antarctic. It was found that moderate and strong stability occur more often with clear- than cloudy-sky conditions, but weaker stability regimes occur almost equally for clear and cloudy conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DOE Cold Climate Heat Pump Challenge: Development, Metrics, and Early Field Observations

Space heating in residential buildings is a major contributor of Greenhouse Gas (GHG) emissions in the United States. New, advanced electric heat pumps are poised to provide a low carbon alternative to traditional fossil-based heating, especially in colder climates. Widespread deployment of cold climate heat pumps could help address the significant portion of building emissions and primary energy used in American households, but these gains will require broader acceptance from consumers and decision makers. As part of the “Energy, Emissions, and Equity” (E3) initiative, the U.S. Department of Energy launched the Cold Climate Heat Pump Challenge (CCHP) in 2021 to accelerate the deployment of the next generation of air source heat pumps. The Challenge is currently focused on residential, centrally ducted, electric heat pumps, with a nominal cooling capacity greater than or equal to 24,000 Btu/h (2 tons) and less than or equal to 65,000 Btu/h (5 tons). The Challenge specifications represent a best-in-class heat pump product that provides high-efficiency heating performance in cold climates, employs environmentally friendly low-Global Warming Potential (GWP) refrigerants, and is designed to be grid interactive. Spearheaded by DOE in partnership with the US Environmental Protection Agency (EPA) and Natural Resources Canada (NRCan), the CCHP Challenge brings together numerous major heat pump manufacturers and key stakeholders including utilities and state agencies across the country. In addition to performance testing in a laboratory environment, a key component of this research is evaluating the in-field performance of the prototype heat pumps developed as part of the Challenge. This paper will discuss the development of the CCHP Challenge, key performance specifications, and the energy and non-energy metrics that will be evaluated through the field study. The types of data that are being collected from the units installed in homes across North America as part of the field validation will be described, including details about the sensors and data acquisition system. Data cleaning and analysis methodologies will be discussed, along with early observations from the winter 2022-2023, spring and summer performance periods, and an assessment of non-energy metrics through pre- and post-installation homeowner surveys. Challenges uncovered and lessons learned throughout the process will also be discussed. Finally, the paper will discuss next steps for the second winter assessment period, areas needing additional research, and potential applications of the data collected and analyzed through this work.

Mendon, Vrushali V.↗
DOE OSTI2024

BENEFIT with Northeastern University: HVAC Hardware-in-the-Loop Experimental Testing of a Heat Pump and Air Conditioner

This dataset includes HVAC Hardware-in-the-Loop (HIL) experimental results for a single stage, SEER 16, HSPF 9.5, 3-ton single-speed air source heat pump with 15 kW of backup auxiliary heating tested in both cooling and heating mode, and a two stage, SEER 21, 2-ton central air conditioner tested in cooling mode for a set of outdoor temperatures and indoor setpoint temperatures. In addition to these tests, experimental tests focused on the operation of auxiliary heating for the heat pump for winter condition were also conducted. The laboratory experiments for transient testing of the heat pump and air conditioner were conducted using the two HIL systems in the Systems Performance Laboratory (SPL) at NREL’s Energy Systems Integration Facility (ESIF). Further information on laboratory design and capabilities of the SPL along with the architecture of HVAC HIL system can be found in: Sparn, B. F. 2018. Laboratory Resources and Techniques to Evaluate Smart Home Technology (No. NREL/CP-5500-71696). National Renewable Energy Laboratory (NREL), Golden, CO (United States). https://www.nrel.gov/docs/fy18osti/71696.pdf and the experimental setup and validation of HVAC HIL platform can be found in: Ramaraj, S. and Sparn, B. 2022. Validation of HVAC Hardware-In-the-Loop Simulation for Advanced Control Strategies in Smart Homes (No. NREL/CP-5500-82562). National Renewable Energy Lab (NREL), Golden, CO (United States). https://www.nrel.gov/docs/fy22osti/82562.pdf. These experimental results can be used to validate how we currently model the cycling behavior of heat pumps and air conditioners. Additionally, many demand response programs implement heat pump and air conditioner control by changing the thermostat set point – these data may also be used to verify our models for heat pump and air conditioner demand response control are implemented correctly. The Test_Matrix file describes all the indoor and outdoor test conditions for heat pump and air conditioner and the file names of data sets include information about the test conditions. A wide range of outdoor air temperatures were chosen to accommodate summer and winter conditions. In addition to operating the HVAC equipment with different outdoor temperatures, we also operate the system with different indoor temperature set points to represent different grid signals or different operating conditions. For cooling conditions, the baseline set point is 72°F. To represent Load Up signals, the setpoint is changed to 68°F. The Load Shed set point is 76°F. For heating conditions, the baseline set point was assumed to be 68°F. The Load add set point is 72°F and the Load shed set point is 64°F. The starting indoor temperature for cooling conditions was set ~2°F above the indoor setpoint temperature so that the equipment turned on quickly. Similarly, the initial indoor temperature was set ~2°F lower than setpoint for heating mode tests to ensure that heating began quickly. The return air temperature was assumed to be equal to the indoor setpoint temperature in all cases. The experimental data are sampled at 1-second intervals. The data from ecobee thermostat at 5-minute interval are resampled and added to the corresponding file. The content of each data set is as follows: • T_Return (C): Measured return air temperature [C] • T_Return_SP (C): Return air temperature setpoint from E+ model, sent to HIL [C] • T_Supply (C): Measured supply air temperature at evaporator outlet [C] • T_Outdoor (C): Measured outdoor air temperature [C] • T_Outdoor_SP (C): Outdoor air temperature setpoint from weather file, sent to HIL [C] • T_Indoor (C): Measured indoor air temperature [C] • T_Indoor_SP (C): Indoor air temperature setpoint from E+ model, sent to HIL [C] • Outdoor Unit Power (W): Measured power of the outdoor unit [W] • Indoor Unit Power (W): Measured power of the indoor unit [W] • Evaporator Airflow Rate (CFM): Measured evaporator or indoor unit airflow rate sent to E+ model [CFM] • Cooling/Heating Capacity (kW): Calculated cooling/heating capacity sent to E+ model [kW] • T_SP_Thermostat (C): Thermostat cooling/heating setpoint temperature [C] • T_Indoor_Thermostat (C): Thermostat indoor air temperature [C]

24 POWER TRANSMISSION AND DISTRIBUTION↗