Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snow contribution”

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 73 records · Page 4

Extensive Liquid Meltwater Storage in Firn Within the Greenland Ice Sheet

The accelerating loss of mass from the Greenland ice sheet is a major contribution to current sea level rise. Increased melt water runoff is responsible for half of Greenlands mass loss increase. Surface melt has been increasing in extent and intensity, setting a record for surface area melt and runoff in 2012. The mechanisms and timescales involved in allowing surface melt water to reach the ocean where it can contribute to sea level rise are poorly understood. The potential capacity to store this water in liquid or frozen form in the firn (multi-year snow layer) is significant, and could delay its sea-level contribution. Here we describe direct observation of water within a perennial firn aquifer persisting throughout the winter in the southern ice sheet,where snow accumulation and melt rates are high. This represents a previously unknown storagemode for water within the ice sheet. Ice cores, groundairborne radar and a regional climatemodel are used to estimate aquifer area (70 plue or minus 10 x 10(exp 3) square kilometers ) and water table depth (5-50 m). The perennial firn aquifer represents a new glacier facies to be considered 29 in future ice sheet mass 30 and energy budget calculations.

meltwater↗

New Developments for Physically-based Falling Snow Retrievals over Land in Preparation for GPM

The NASA Global Precipitation Measurement mission (GPM) concept centers on deploying a Core spacecraft carrying a dual-frequency precipitation radar and a microwave radiometric imager with channels from 10 to 183 GHz to serve as a precipitation physics observatory and a calibration reference to unify a constellation of dedicated and operational passive microwave sensors. Because of the extended orbit of the Core (plus or minus 65 deg) and the enhanced dual frequency radar and high frequency radiometer, GPM will be able to sense falling snow precipitation and light rain over land. Accordingly, GPM has partnered with the Canadian CloudSat/CALIPSO Validation Project (C3VP) to obtain observations to provide one of several important ground-based validation data sets around which the falling snow models and retrieval algorithms can be further developed and tested. In this work we compare and correlate the long time series (Nov.'06 - March '07) measurements of precipitation rate from parsivels to the passive (89, 150, 183 plus or minus 1, plus or minus 3, plus or minus 7 GHz) observations of NOAA's AMSU-B radiometer. We separate the comparisons into categories of no precipitation, liquid rain and falling snow precipitation. We found that there are similar TBs (especially at 89 and 150 GHz) for cases with falling snow and for non-precipitating cases. The comparisons indicate that surface emissivity contributions to the satellite observed TB over land can add uncertainty in detecting and estimating falling snow. The newest results show that by computing brightness temperatures based on CARE radiosonde data and a rough estimate of surface emissivity show that the cloud ice scattering signal in the AMSU-B data is detected. That is the differences in computed TB and AMSU-B TB for precipitating and non-precipitating cases are unique such that the precipitating and non-precipitating cases can be identified. These results require that the radiosonde releases are within an hour of the AMSU-B data. Forest fraction, snow cover, and measured emissivities were combined to calculate the surface emissivities.

Jackson, Gail S.↗

GPM Pre-Launch Algorithm Development for Physically-Based Falling Snow Retrievals

In this work we compare and correlate the long time series (Nov.-March) n~easurements of precipitation rate from the Parsivels and 2DVD to the passive (89, 150, 183+/-1, +/-3, +/-7 GHz) observations of NOAA's AMSU-B radiometer. There are approximately 5-8 AMSU-B overpass views of the CARE site a day. We separate the comparisons into categories of no precipitation, liquid rain and falling snow precipitation. Scatterplots between the Parsivel snowfall rates and AMSU-B brightness temperatures (TBs) did not show an exploitable relationship for retrievals. We further compared and contrasted brightness temperatures to other surface measurements such as temperature and relative humidity with equally unsatisfying results. We found that there are similar TBs (especially at 89 and 150 GHz) for cases with falling snow and for non-precipitating cases. The comparisons indicate that surface emissivity contributions to the satellite observed TB over land can add uncertainty in detecting and estimating falling snow. The newest results show that the cloud icc scattering signal in the AMSU-B data call be detected by computing clear air TBs based on CARE radiosonde data and a rough estimate of surface emissivity. That is the differences in computed TI3 and AMSU-B TB for precipitating and nonprecipitating cases are unique such that the precipitating versus ~lon-precipitating cases can be identified. These results require that the radiosonde releases are within an hour of the AMSU-B data and allow for three surface types: no snow on the ground, less than 5 cm snow on the ground, and greater than 5 cm on the ground (as given by ground station data). Forest fraction and measured emissivities were combined to calculate the surface emissivities. The above work and future work to incorporate knowledge about falling snow retrievals into the framework of the expected GPM Bayesian retrievals will be described during this presentation.

Gail Skofronick-Jackson↗

Comparisons of Bin and Bulk Microphysics Schemes in Simulations of Topographic Winter Precipitation with Radar and Radiometer Measurements

The Hebrew University Cloud Model (HUCM) bin scheme and the Thompson bulk scheme in the Weather Research and Forecasting (WRF) model are compared to assess biases often found in simulated brightness temperature and radar reflectivity. Compared to our preceding study that evaluated several bulk schemes in the WRF model, the current study obtains a reduction of the bias from excessive microwave scattering by precipitation ice for both HUCM bin and the Thompson bulk microphysics schemes for a topographic winter precipitation event associated with an atmospheric river. The Thompson particle size distributions (PSDs) and snow particle density assumption are implemented into the Goddard Satellite Data Simulator Unit (G‐SDSU) and have produced improvements. Despite the greater sophistication of the bin scheme in representing cloud and precipitation processes, the simulation with the Thompson bulk scheme is generally in better agreement with observations for this winter event. The explicitly resolved hydrometeor PSDs in HUCM enable analysis of mass spectra variations in response to changes in microphysics assumptions. Two HUCM sensitivity runs tested the enhancement of snow particle breakup and the influence of ice nuclei (IN) concentration. Higher IN concentration resulted in increased snow mass and broadened the spectrum toward small‐size particles. Modified snow mass spectra and resultant changes in graupel contributed to modifications in scattering and reflectivity simulations. The article demonstrates the bin scheme's capability to provide a new means to improve our understanding of uncertainties in mesoscale weather models and radiative transfer models.

microwave remote sensing↗

North Atlantic Simulations in Coordinated Ocean-Ice Reference Experiments Phase II (CORE-II). Part I: Mean States

Simulation characteristics from eighteen global ocean-sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Oceanice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort. atmospheric forcing atmospheric temperatures

ocean model comparisons↗

Seasonal Variability in the Mechanisms Behind the 2020 Siberian Heatwaves

Record-breaking heatwaves and wildfires immersed Siberia during the boreal spring of 2020 following an anomalously warm winter. Springtime heatwaves are becoming more common in the region, with statistically significant trends in the frequency, magnitude, and duration of heatwave events over the past four decades. Mechanisms by which the heatwaves occur and contributing factors differ by season. Winter heatwave frequency is correlated with the atmospheric circulation, particularly the Arctic Oscillation, while the frequency of heatwaves during the spring months is highly correlated with aspects of the land surface including snow cover, albedo, and latent heat flux. Idealized AMIP-style experiments are used to quantify the contribution of suppressed Arctic sea ice and snow cover over Siberia on the atmospheric circulation, surface energy budget, and surface air temperature in Siberia during the winter and spring of 2020. Sea ice concentration contributed to the strength of the stratospheric polar vortex and Arctic Oscillation during the winter months, thereby influencing the tropospheric circulation and surface air temperature over Siberia. Warm temperatures across the region resulted in an earlier-than-usual recession of the winter snowpack. The exposed land surface contributed to up to 20% of the temperature anomaly during the spring through the albedo feedback and changes in the ratio of the latent and sensible heat fluxes. This, in combination with favorable atmospheric circulation patterns, resulted in record-breaking heatwaves in Siberia in the spring of 2020.

Allison B. Marquardt Collow↗

A cleaner snow future mitigates Northern Hemisphere snowpack loss from warming

Abstract Light-absorbing particles (LAP) deposited on seasonal snowpack can result in snow darkening, earlier snowmelt, and regional climate change. However, their future evolution and contributions to snowpack change relative to global warming remain unclear. Here, using Earth System Model simulations, we project significantly reduced black carbon deposition by 2081-2100, which reduces the December-May average LAP-induced radiative forcing in snow over the Northern Hemisphere from 1.3 Wm −2 during 1995-2014 to 0.65 (SSP126) and 0.49 (SSP585) Wm −2 . We quantify separately the contributions of climate change and LAP evolution on future snowpack and demonstrate that projected LAP changes in snow over the Tibetan Plateau will alleviate future snowpack loss due to climate change by 52.1 ± 8.0% and 8.0 ± 1.1% at the end of the century for the two scenarios, mainly due to reduced black carbon contamination. Our findings highlight a cleaner snow future and its benefits for future water supply from snowmelt especially under the sustainable development pathway of SSP126.

54 ENVIRONMENTAL SCIENCES↗

An AeroCom Assessment of Black Carbon in Arctic Snow and Sea Ice

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng/g for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng/g for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng/g. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model-measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60-90degN) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07-0.25) W/sq m and 0.18 (0.06-0.28) W/sq m in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W/sq m for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.

carbon↗

Impact of Seasonal Snow-Cover Change on the Observed and Simulated State of the Atmospheric Boundary Layer in a High-Altitude Mountain Valley

The structure and evolution of the atmospheric boundary layer (ABL) under clear-sky fair weather conditions over mountainous terrain is dominated by the diurnal cycle of the surface energy balance and thus strongly depends on surface snow cover. We use data from three passive ground-based infrared spectrometers deployed in the East River Valley in Colorado's Rocky Mountains to investigate the response of the thermal ABL structure to changes in surface energy balance during the seasonal transition from low to high snow cover. Temperature profiles were retrieved from the infrared radiances using the optimal estimation physical retrieval Tropospheric Remotely Observed Profiling via Optimal Estimation. A nocturnal surface inversion formed in the valley during clear-sky days, which was subsequently mixed out during daytime with the development of a convective boundary layer when snow cover was low. Over high snow cover, a very shallow convective boundary layer formed, above which the inversion persisted through the daytime hours. We compare these observations to NOAA's operational High-Resolution-Rapid-Refresh model and find large warm biases on clear-sky days resulting from the model's inability to form strong nocturnal inversions and to maintain the stable stratification in the valley during daytime when there was snow on the ground. We suggest several factors contributing to the large model errors. These are (a) the inability of the model to represent well-developed thermally driven flows likely due to the too coarse horizontal grid spacing (3 km), (b) too much convective mixing during daytime, and (c) too strong vertical coupling between the valley atmosphere and the free troposphere.

54 ENVIRONMENTAL SCIENCES↗

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, there has become an increasing interest in the use of millimeter-wave (mm-wave) and sub-millimeter-wave (submm-wave) radiometer observations to investigate the properties of ice particles in clouds. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. One methodology used since the 1960's to discern the relationship between the physical state observed and the brightness temperature (TB) is through the temperature weighting function profile. In this research, the temperature weighting function concept is exploited to analyze the sensitivity of various characteristics of the cloud profile, such as relative humidity, ice water path, liquid water path, and surface emissivity. In our numerical analysis, we compute the contribution (in Kelvin) from each of these cloud and surface characteristics, so that the sum of these various parts equals the computed TB. Furthermore, the percentage contribution from each of these characteristics is assessed. There is some intermingling/contamination of the contributions from various components due to the integrated nature of passive observations and the absorption and scattering between the vertical layers, but all in all the knowledge gained is useful. This investigation probes the sensitivity over several cloud classifications, such as cirrus, blizzards, light snow, anvil clouds, and heavy rain. The focus is on mm-wave and submm-wave frequencies, however discussions of the effects of cloud variations to frequencies as low as 10 GHz and up to 874 GHz will also be presented. The results show that nearly 60% of the TB value at 89 GHz comes from the earth's surface for even the heaviest blizzard snow rates. On the other hand, a significant percentage of the TB value comes from the snow in the cloud for 166, and 183 plus or minus 7 GHz for the heavy and medium snow rates. For submm-wave channels, there is no contribution from the surface because these channels cannot probe through clouds, nor normal water vapor amounts in clear air regions. This work is extremely valuable in physically-based retrieval algorithm development research.

Gail Skofronick-Jackson↗

Simulation of Seasonal Snow Microwave TB Using Coupled Multi-Layered Snow Evolution and Microwave Emission Models

The accurate quantification of SWE has important societal benefits, including improving domestic and agricultural water planning, flood forecasting and electric power generation. However, passive-microwave SWE algorithms suffer from variations in TB due to snow metamorphism, difficult to distinguish from those due to SWE variations. Coupled snow evolution-emission models are able to predict snow metamorphism, allowing us to account for emissivity changes. They can also be used to identify weaknesses in the snow evolution model. Moreover, thoroughly evaluating coupled models is a contribution toward the assimilation of TB, which leads to a significant increase in the accuracy of SWE estimates.

microwave↗

Warming and snow loss increase reliance on old groundwater in a Colorado River headwater

Atmospheric warming is reducing snowpack, with uncertain effects on mountainous streamflow, a crucial water resource. Despite limited historical observations of groundwater–streamflow interactions above 2,500 m, new measurements in the Upper Colorado River headwaters indicate declining groundwater storage that is dated decades to millennia old. Here we use integrated hydrologic modelling spanning water years 2015–2021 to determine whether the loss of old-age groundwater buffers streamflow during low-snow years and whether that loss is exacerbated with warming. Results show that old-groundwater contributions to streams remain relatively steady through time, unlike the more variable contributions from young groundwater. Numerical experiments of increased surface air temperatures (+2.5 °C and +4 °C) increase rain–snow fractions and evapotranspiration and decrease runoff ratio by 2–3% per degree Celsius increase. As streamflow declines with warming, the age of groundwater supporting it gets older, in part owing to intermediate-aged (1–3 year) groundwater declining twice as fast. Simulations show that water table depths at higher elevations (>3,700 m) decline disproportionately and fail to recover even during wet years. These findings suggest altered groundwater–streamflow interactions with warming and snow loss, with implications for water resources.

Siirila-Woodburn, Erica R. [Lawrence Berkeley Nati↗

Snowex 2017 Community Snow Depth Measurements: A Quality-Controlled, Georeferenced Product

Snow depth was one of the core ground measurements required to validate remotely-sensed data collected during SnowEx Year 1, which occurred in Colorado. The use of a single, common protocol was fundamental to produce a community reference dataset of high quality. Most of the nearly 100 Grand Mesa and Senator Beck Basin SnowEx ground crew participants contributed to this crucial dataset during 6-25 February 2017. Snow depths were measured along ~300 m transects, whose locations were determined according to a random-stratified approach using snowfall and tree-density gradients. Two-person teams used snowmobiles, skis, or snowshoes to travel to staked transect locations and to conduct measurements. Depths were measured with a 1-cm incremented probe every 3 meters along transects. In shallow areas of Grand Mesa, depth measurements were also collected with GPS snow-depth probes (a.k.a. MagnaProbes) at ~1-m intervals. During summer 2017, all reference stake positions were surveyed with <10 cm accuracy to improve overall snow depth location accuracy. During the campaign, 193 transects were measured over three weeks at Grand Mesa and 40 were collected over two weeks in Senator Beck Basin, representing more than 27,000 depth values. Each day of the campaign depth measurements were written in waterproof field books and photographed by National Snow and Ice Data Center (NSIDC) participants. The data were later transcribed and prepared for extensive quality assessment and control. Common issues such as protocol errors (e.g., survey in reverse direction), notebook image issues (e.g., halo in the center of digitized picture), and data-entry errors (sloppy writing and transcription errors) were identified and fixed on a point-by-point basis. In addition, we strove to produce a georeferenced product of fine quality, so we calculated and interpolated coordinates for every depth measurement based on surveyed stakes and the number of measurements made per transect. The product has been submitted to NSIDC in csv format. To educate data users, we present the study design and processing steps that have improved the quality and usability of this product. Also, we will address measurement and design uncertainties, which are different in open vs. forest areas.

Brucker, L.↗

Passive microwave observations of the Wedell Sea during austral winter and early spring

The results of multispectral passive microwave observations (6.7 to 90-GHz) are presented from the cruises of the FS Polarstern in the Weddell Sea from July to December 1986. This paper includes primarily the analysis of radiometric observations taken at ice station sites. Averaged emissivity spectra for first-year (FY) ice were relatively constant throughout the experiment and were not statistically different from FY ice signatures in the Arctic. Detailed ice characterization was carried out at each site to compare the microwave signatures of the ice with the physical properties. Absorption optical depths of FY ice were found to be sufficiently high that only the structure in the upper portions of the ice contributed significantly to interstation emissivity variations. The emissivities at 90-GHz, e(90), had the greatest variance. Both e(90) at vertical polarization and GR(sub e)(90, 18.7)(defined as (e(sub V)(90)-e(sub V)(18.7))/e(sub V)(90 + e(sub V)(18.7)) depended on the scattering optical depth which is a function of the snow grain diameter and layer thickness. The variance showed a latitude dependence and is probably due to an increase in the strength of snow metamorphism nearer the northern edge of the ice pack. The contribution of variations of near-surface brine volume to the emissivity was not significant over the range of values encountered at the station sites. Emissivity spectra are presented for a range of thin ice types. Unsupervised principal component analysis produced three significant eigenvectors and showed a separation among four different surface types: open water, thin ice, FY ice, and FY ice with a thick snow cover. A comparison with SMMR satellite data showed that average ice concentrations derived from the ship's ice watch log were consistent with the satellite concentrations. The surface based emissivities for FY ice were also compared with emissivities calculated from scanning multichannel microwave radiometer (SMMR) satellite radiances. Best agreement was found at 6.7 and 10-GHz, while at 18 and 37-GHz, SMMR emissivities were slightly lower than surface based results. For the three lower frequencies agreement was found within a confidence limit of 95% and for 37-GHz within about 90%.

Grenfell, T. C.↗

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980-47 2015 using the National Climate Assessment-Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah Version 3.3 land surface model at 0.125 degreesx0.125 degrees forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth and irrigation products. Mean annual trends are reported using the nonparametric Mann-Kendall test at p<0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends, and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm/yr in the Upper Great Plains and Northeast to -1 to -9 in the West and South; net radiation flux trends range from +0.05 to +0.20 W/m(exp2)/yr in the East to -0.05 to -0.20 in the West; U.S.-wide temperature trends average about +0.03K/yr. Trends in soil moisture, snow cover, latent and sensible heat fluxes and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from West to East. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days/yr while trend detection increased 11%; Trends in latent heat flux were hardly affected, RMSE decreasing only 0.17 to 0.16 W/m(exp 2)/yr, while trend detection increased 2%; NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gaged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Snow depth↗

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980–2015 using the National Climate Assessment Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah version 3.3 land surface model at 0.125° 3 0.125° forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth, and irrigation products. Mean annual trends are reported using the nonparametric Mann–Kendall test at p < 0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm/yr in the upper Great Plains and Northeast to -1 to -9 mm/yr in the West and South, net radiation flux trends range from +0.05 to +0.20 W/sq. m yr in the East to -0.05 to -0.20 W/sq. m yr in the West, and U.S.-wide temperature trends average about +0.03K/yr. Trends in soil moisture, snow cover, latent and sensible heat fluxes, and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from west to east. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days/yr while trend detection increased 11%. Trends in latent heat flux were hardly affected, with RMSE decreasing only from 0.17 to 0.16 W/sq. m yr, while trend detection increased 2%. NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gauged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Michael F. Jasinski↗

Maximum entropy modeling to identify physical drivers of shallow snowpack heterogeneity using unpiloted aerial system (UAS) lidar

Understanding the spatial variability of the snowpack is valuable for hydrologists and ecologists seeking to predict hydrological processes in a cold region. Snow distribution is a function of interactions among static variables, such as terrain, vegetation, and soil properties, and dynamic meteorological variables, such as solar radiation, wind speed and direction, and soil moisture. However, identifying the dominant physical drivers responsible for spatial patterns of the snowpack, particularly for ephemeral, shallow snowpacks, has been challenging due to the lack of the high-resolution snowpack and physical variables with high vertical accuracy as well as inherent limitations in traditional approaches. This study uses an Unpiloted Aerial System (UAS) lidar-based snow depth and static variables (1-m spatial resolution) to analyze field-scale spatial structures of snow depth and apply the Maximum Entropy (MaxEnt) model to identify primary controls over open terrain and forests at the University of New Hampshire Thompson Farm Research Observatory, New Hampshire, United States. We found that, among nine topographic and soil variables, plant functional type and terrain roughness contribute up to 80% and 76% of relative importance in the MaxEnt framework to predict locations of deeper or shallower snowpacks, respectively, across a mixed temperate forested and field landscape. Soil variables, such as organic matter and saturated hydraulic conductivity, were also important controls (up to 70% and 81%) on snow depth spatial variations for both open and forested landscapes suggesting spatial variations in soil variables under snow can control thermal transfer among soil, snowpack, and surface-atmosphere. This work contributes to improving land surface and snow models by informing parameterization of the sub-grid scale snow depths, down-scaling remotely sensed snow products, and understanding field scale snow states.

Eunsang Cho↗