Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Snow model”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydrologic Model Files (PRMS) for Historical Conditions in the East River Watershed, Colorado between 1987-2019

The data package contains model input and output file for the East River watershed (85 km2) located in the headwaters of the Colorado River. The model is the U.S. Geological Survey (USGS) Precipitation-Modeling Runoff System (PRMS) for the historical, or baseline, condition (1987-2019) to assess the efficiency of summer monsoon rains to generate streamflow. The finite difference grid resolution is 100 m with elevations resampled from the USGS National Elevation Dataset. LANDFIRE raster data is used to derive parameters of dominant cover type, summer and winter cover density, canopy interception characteristics for snow and rain and transmission coefficient for shortwave radiation. Climate forcing uses minimum and maximum daily temperature lapse rates defined by the two Snow Telemetry (SNOTEL) stations located in the vicinity of the East River adjusted for aspect. Schofield SNOTEL snowfall is spatially distributed using LiDAR derived snow depth observations from the Airborne Snow Observatory (ASO) flown 4 April 2016. Snow depths are converted to snow water equivalent based on ground surveys and density modeling. Rainfall is spatially distributed using the monthly Parameter-elevation Relationships on Independent Slopes Model 30-year monthly averages. Simulated solar radiation is calibrated to match weather station observations. Model verification of snow water equivalent (SWE) accumulation and ablation relies on the 2018 and 2019 ASO maps and weather station snow depth. The data package contains model input files in text format, model documentation in pfd and a readme.txt file to help users understand the file structure to run the program.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deep Learning Advances Arctic River Water Temperature Predictions

The accelerated warming in the Arctic poses serious risks to freshwater ecosystems by altering streamflow and river thermal regimes. However, limited research on Arctic River water temperatures exists due to data scarcity and the absence of robust methodologies, which often focus on large, major river basins. To address this, we leveraged the newly released, extensive AKTEMP data set and advanced machine learning techniques to develop a Long Short-Term Memory (LSTM) model. By incorporating ERA5-Land reanalysis data and integrating physical understanding into data-driven processes, our model advanced river water temperature predictions in ungauged, snow- and permafrost-affected basins in Alaska. Our model outperformed existing approaches in high-latitude regions, achieving a median Nash-Sutcliffe Efficiency of 0.95 and root mean squared error of 1.0°C. The LSTM model learned air temperature, soil temperature, solar radiation, and thermal radiation—factors associated with energy balance—were the most important drivers of river temperature dynamics. Soil moisture and snow water equivalent were highlighted as critical factors representing key processes such as thawing, melting, and groundwater contributions. Glaciers and permafrost were also identified as important covariates, particularly in seasonal river water temperature predictions. Our LSTM model successfully captured the complex relationships between hydrometeorological factors and river water temperatures across varying timescales and hydrological conditions. This scalable and transferable approach can be potentially applied across the Arctic, offering valuable insights for future conservation and management efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Experimental Whole-Ecosystem Warming Enables Novel Estimation of Snow Cover and Depth Sensitivities to Temperature, and Quantification of the Snow-Albedo Feedback Effect

Climate change is reducing the amount, duration, and extent of snow across high-latitude ecosystems. But, in landscapes where persistent winter snow cover develops, experimental platforms to specifically investigate interactions between warming and changes in snowpack, and impacts on ecosystem processes, have been lacking. We leveraged a whole-ecosystem warming experiment in a boreal peatland forest to quantify how snow duration, depth, and fractional cover vary with warming of up to +9°C. We found that every snow-related quantity we examined declined precipitously as the amount of warming increased. The importance of deep, continuous snow cover for moderating shallow soil temperature is highlighted by an increase in soil temperature variance and the frequency of short-duration freeze-thaw cycles in the warmer plots. We used a paired-plot approach to estimate the magnitude of the snow-albedo feedback effect, and demonstrate that albedo-driven warming linked to reduced snow cover varies between December (+0.4°C increase in maximum air temperature) and March (+1.2°C increase) because of differences in insolation. Overall, results show that even modest future warming will have profound impacts on northern winters and cold-season ecosystem processes. Plot-level data from this warming experiment, and emergent relationships between warming and quantities related to snow cover and duration, could be of enormous value for testing and improving the representation of snow processes in simulation models, especially under future climate scenarios that are outside of the range of historically observed variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling the Winter Heat Conduction Through the Sea Ice System During MOSAiC

Abstract Models struggle to accurately simulate observed sea ice thickness changes, which could be partially due to inadequate representation of thermodynamic processes. We analyzed co‐located winter observations of the Arctic sea ice from the Multidisciplinary Drifting Observatory for the Study of the Arctic Climate for evaluating and improving thermodynamic processes in sea ice models, aiming to enable more accurate predictions of the warming climate system. We model the sea ice and snow heat conduction for observed transects forced by realistic boundary conditions to understand the impact of the non‐resolved meter‐scale snow and sea ice thickness heterogeneity on horizontal heat conduction. Neglecting horizontal processes causes underestimating the conductive heat flux of 10% or more. Furthermore, comparing model results to independent temperature observations reveals a ∼5 K surface temperature overestimation over ice thinner than 1 m, attributed to shortcomings in parameterizing surface turbulent and radiative fluxes rather than the conduction. Assessing the model deficiencies and parameterizing these unresolved processes is required for improved sea ice representation.

Geology↗
DOE OSTI2025

Observations of high-time-resolution and size-resolved aerosol chemical composition and microphysics in the central Arctic: implications for climate-relevant particle properties

Abstract. Aerosols play a critical role in the Arctic's radiative balance, influencing solar radiation and cloud formation. Limited observations in the central Arctic leave gaps in understanding aerosol dynamics year-round, affecting model predictions of climate-relevant aerosol properties. Here, we present the first annual high-time-resolution observations of submicron aerosol chemical composition in the central Arctic during the Arctic Ocean 2018 (AO2018) and the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expeditions. Seasonal variations in the aerosol mass concentrations and chemical composition in the central Arctic were found to be driven by typical Arctic seasonal regimes and resemble those of pan-Arctic land-based stations. Organic aerosols dominated the pristine summer, while anthropogenic sulfate prevailed in autumn and spring under haze conditions. Ammonium, which impacts aerosol acidity, was consistently less abundant, relative to sulfate, in the central Arctic compared to lower latitudes of the Arctic. Cyclonic (storm) activity was found to have a significant influence on aerosol variability by enhancing emissions from local sources and the transport of remote aerosol. Local wind-generated particles contributed up to 80 % (20 %) of the cloud condensation nuclei population in autumn (spring). While the analysis presented herein provides the current central Arctic aerosol baseline, which will serve to improve climate model predictions in the region, it also underscores the importance of integrating short-timescale processes, such as seasonal wind-driven aerosol sources from blowing snow and open leads/ocean in model simulations. This is particularly important, given the decline in mid-latitude anthropogenic emissions and the increase in local ones.

Heutte, Benjamin (ORCID:0000000173727460)↗
DOE OSTI2021

Permafrost Promotes Shallow Groundwater Flow and Warmer Headwater Streams

Abstract The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from 11 headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7% of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5% moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data from: "Towards CONUS-Wide ML-Augmented Conceptually-Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics"

This data package was generated to support the manuscript “Towards CONUS-Wide Machine Learning-Augmented Conceptually Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics.” It provides input files, model outputs, plotting data, scripts, notebooks, and documentation used to develop, evaluate, and reproduce Mass-Conserving Perceptron (MCP)-based hydrologic modeling experiments across 513 selected Catchment Attributes and Meteorology for Large-sample Studies in the United States (CAMELS-US) basins. The files are organized by modeling component and analysis purpose, including rainfall–runoff experiments, snow module experiments, coupled hydrologic-snow experiments, Long Short-Term Memory (LSTM) benchmark results, model skill metrics, initialization and epoch records, cell-state normalization files, Akaike Information Criterion (AIC)-based model comparison files, and data used to generate manuscript figures. Tabular files can be opened using standard spreadsheet software or Python/R data-analysis tools. Python scripts, Jupyter notebooks, and selected MATLAB scripts are included for model execution, postprocessing, plotting, and statistical analysis. Quality assurance and quality control were conducted through the source-data selection and modeling workflow. Meteorological forcing, streamflow, and static catchment attributes were derived from the CAMELS-US dataset, and snow water equivalent data were derived from the University of Arizona (UA) Snow Water Equivalent dataset. Selected basins and time periods were screened during the associated research workflow to avoid missing observations or poor-quality cases. Static geospatial features were processed primarily using Quantum Geographic Information System (QGIS) and Geospatial Data Abstraction Library (GDAL) workflows. Additional details are provided in the associated manuscript and documentation.

ESS-DIVE CSV File Formatting Guidelines Reporting ↗
DOE OSTI2020

What Polarimetric Weather Radars Offer to Cloud Modelers: Forward Radar Operators and Microphysical/Thermodynamic Retrievals

The utilization of polarimetric weather radars for optimizing cloud models is a next frontier of research. It is widely understood that inadequacies in microphysical parameterization schemes in numerical weather prediction (NWP) models is a primary cause of forecast uncertainties. Due to its ability to distinguish between hydrometeors with different microphysical habits and to identify “polarimetric fingerprints” of various microphysical processes, polarimetric radar emerges as a primary source of needed information. There are two approaches to leverage this information for NWP models: (1) radar microphysical and thermodynamic retrievals and (2) forward radar operators for converting the model outputs into the fields of polarimetric radar variables. In this paper, we will provide an overview of both. Polarimetric measurements can be combined with cloud models of varying complexity, including ones with bulk and spectral bin microphysics, as well as simplified Lagrangian models focused on a particular microphysical process. Combining polarimetric measurements with cloud modeling can reveal the impact of important microphysical agents such as aerosols or supercooled cloud water invisible to the radar on cloud and precipitation formation. Some pertinent results obtained from models with spectral bin microphysics, including the Hebrew University cloud model (HUCM) and 1D models of melting hail and snow coupled with the NSSL forward radar operator, are illustrated in the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗
DOE OSTI2020

Model of Operation-and-Maintenance Costs for Photovoltaic Systems

This article presents a method for calculating costs associated with operation and maintenance (O&M) of photovoltaic (PV) systems. It compiles details regarding the cost and frequency of multiple O&M services to estimate annual O&M costs ($\$$/year) for each year of an analysis period, the net present value ($\$$) of life cycle costs accumulated over the analysis period, and the reserve account amount ($\$$). Here we show that this method is an improvement over the previous averaged or levelized per-unit ($\$$/kW/year) valuations for estimating PV O&M costs, because it allows a detailed selection of services to perform based on system size, market served (e.g., residential, commercial, or utility), type and configuration of system components (e.g., micro-, string or central inverter) and site and environmental conditions (e.g., snow, pollen, bird populations). This model also distinguishes costs that vary from year to year and increase at different rates over time because of heuristic failure distributions (e.g. Weibull or Lognormal distribution) based on actuarial data for many of the services. This cost model was created by the PV O&M Working Group of researchers and industry, sponsored by DOE Solar Energy Technologies Office, and has been published in an on-line version hosted by SunSpec Alliance at apsuite.sunspec.org. A spreadsheet version is included with this paper as supplementary material.

14 SOLAR ENERGY↗
DOE OSTI2023

Evolution of global snow drought characteristics from 1850 to 2100

Abstract Seasonal snow is an integral part of the global water supply and storage system. Snow droughts impact ecological, agricultural, and urban systems by altering the amount and timing of meltwater delivery. These droughts are characterized by a lack of on-the-ground snow (snow water equivalent, SWE) that can be caused by low total precipitation (dry drought) or low proportion of precipitation falling as snowfall (warm drought), often combined with an early melt. The standardized SWE index (SWEI) ranks the current status of SWE for a given location compared to a baseline condition and identifies the existence, but not the cause, of snow drought. In this work, we use estimates of SWE, temperature, and precipitation from nine coupled model intercomparison project phase 6 (CMIP6) models to quantify the frequency, severity, and type of snow droughts globally for historical and future scenarios. Compared to a historical baseline (1850–1900) total snow drought frequency more than doubles under socio-economic pathway (SSP)2-4.5 and SSP5-8.5; all of the increase in snow drought frequency comes from an increase in warm droughts. The probability distribution of future SWEI in major snowy basins around the world are likely to be centered on more negative values, which corresponds to more severe drought and, with only moderate changes in distribution spread, more frequent drought. CMIP6 simulations pinpoint snow drought as an emerging global threat to water resources and highlight the need to explore higher resolution future models that better capture complex mountain topography, wildland fires, and snow-forest interactions.

CMIP6↗
DOE OSTI2021

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗
DOE OSTI2022

Overview of the MOSAiC expedition: Snow and sea ice

Year-round observations of the physical snow and ice properties and processes that govern the ice pack evolution and its interaction with the atmosphere and the ocean were conducted during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition of the research vessel Polarstern in the Arctic Ocean from October 2019 to September 2020. This work was embedded into the interdisciplinary design of the 5 MOSAiC teams, studying the atmosphere, the sea ice, the ocean, the ecosystem, and biogeochemical processes. The overall aim of the snow and sea ice observations during MOSAiC was to characterize the physical properties of the snow and ice cover comprehensively in the central Arctic over an entire annual cycle. This objective was achieved by detailed observations of physical properties and of energy and mass balance of snow and ice. By studying snow and sea ice dynamics over nested spatial scales from centimeters to tens of kilometers, the variability across scales can be considered. On-ice observations of in situ and remote sensing properties of the different surface types over all seasons will help to improve numerical process and climate models and to establish and validate novel satellite remote sensing methods; the linkages to accompanying airborne measurements, satellite observations, and results of numerical models are discussed. We found large spatial variabilities of snow metamorphism and thermal regimes impacting sea ice growth. We conclude that the highly variable snow cover needs to be considered in more detail (in observations, remote sensing, and models) to better understand snow-related feedback processes. The ice pack revealed rapid transformations and motions along the drift in all seasons. The number of coupled ice–ocean interface processes observed in detail are expected to guide upcoming research with respect to the changing Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

LSTM-Based Data Integration to Improve Snow Water Equivalent Prediction and Diagnose Error Sources

Accurate prediction of snow water equivalent (SWE) can be valuable for water resource managers. Recently, deep learning methods such as long short-term memory (LSTM) have exhibited high accuracy in simulating hydrologic variables and can integrate lagged observations to improve prediction, but their benefits were not clear for SWE simulations. Here we tested an LSTM network with data integration (DI) for SWE in the western United States to integrate 30-day-lagged or 7-day-lagged observations of either SWE or satellite-observed snow cover fraction (SCF) to improve future predictions. SCF proved beneficial only for shallow-snow sites during snowmelt, while lagged SWE integration significantly improved prediction accuracy for both shallow- and deep-snow sites. The median Nash–Sutcliffe model efficiency coefficient (NSE) in temporal testing improved from 0.92 to 0.97 with 30-day-lagged SWE integration, and root-mean-square error (RMSE) and the difference between estimated and observed peak SWE values d max were reduced by 41% and 57%, respectively. DI effectively mitigated accumulated model and forcing errors that would otherwise be persistent. Moreover, by applying DI to different observations (30-day-lagged, 7-day-lagged), we revealed the spatial distribution of errors with different persistent lengths. For example, integrating 30-day-lagged SWE was ineffective for ephemeral snow sites in the southwestern United States, but significantly reduced monthly-scale biases for regions with stable seasonal snowpack such as high-elevation sites in California. These biases are likely attributable to large interannual variability in snowfall or site-specific snow redistribution patterns that can accumulate to impactful levels over time for nonephemeral sites. These results set up benchmark levels and provide guidance for future model improvement strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Atmospheric particle abundance and sea salt aerosol observations in the springtime Arctic: a focus on blowing snow and leads

Abstract. Sea salt aerosols play an important role in the radiation budget and atmospheric composition over the Arctic, where the climate is rapidly changing. Previous observational studies have shown that Arctic sea ice leads are an important source of sea salt aerosols, and modeling efforts have also proposed blowing snow sublimation as a source. In this study, size-resolved atmospheric particle number concentrations and chemical composition were measured at the Arctic coastal tundra site of Utqiaġvik, Alaska, during spring (3 April–7 May 2016). Blowing snow conditions were observed during 25 % of the 5-week study period and were overpredicted by a commonly used blowing snow parameterization based solely on wind speed and temperature. Throughout the study, open leads were present locally. During periods when blowing snow was observed, significant increases in the number concentrations of 0.01–0.06 µm particles (factor of 6, on average) and 0.06–0.3 µm particles (67 %, on average) and a significant decrease (82 %, on average) in 1–4 µm particles were observed compared to low wind speed periods. These size distribution changes were likely caused by the generation of ultrafine particles from leads and/or blowing snow, with scavenging of supermicron particles by blowing snow. At elevated wind speeds, both submicron and supermicron sodium and chloride mass concentrations were enhanced, consistent with wind-dependent local sea salt aerosol production. At moderate wind speeds below the threshold for blowing snow as well as during observed blowing snow, individual sea spray aerosol particles were measured. These individual salt particles were enriched in calcium relative to sodium in seawater due to the binding of this divalent cation with organic matter in the sea surface microlayer and subsequent enrichment during seawater bubble bursting. The chemical composition of the surface snowpack also showed contributions from sea spray aerosol deposition. Overall, these results show the contribution of sea spray aerosol production from leads on both aerosols and the surface snowpack. Therefore, if blowing snow sublimation contributed to the observed sea salt aerosol, the snow being sublimated would have been impacted by sea spray aerosol deposition rather than upward brine migration through the snowpack. Sea spray aerosol production from leads is expected to increase, with thinning and fracturing of sea ice in the rapidly warming Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

Meteorology & Atmospheric Sciences↗