The Challenges of Simulating SWE Beneath Forest Canopies are Reduced by Data Assimilation of Snow Depth
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
This data package contains data and information about Global Navigation Satellite System (GNSS)-based Vegetation Optical Depth (VOD), tree sway motion, and eddy-covariance evapotranspiration (ET) data collected at the Niwot Ridge Subalpine Forest AmeriFlux site (US-NR1). The raw GNSS data were collected between May 2022 and August 2023. Other processed datasets such as tree sway motion and ET data are also included. The goal was to study the water content within a subalpine forest and, more specifically, examine the canopy evaporation process. This data archive includes all data that were used within the following Biogeosciences discussion paper that further summarizes the research objectives and conclusions:Burns, S.P., V. Humphrey, E.D. Gutmann, M.S. Raleigh, D.R. Bowling, and P.D. Blanken, 2025: Using GNSS-based vegetation optical depth, tree sway motion, and eddy-covariance to examine evaporation of canopy-intercepted rainfall in a subalpine forest. EGUsphere [preprint],https://doi.org/10.5194/egusphere-2025-1755This data archive also supplements the 30-min Lawrence Berkeley National Laboratory (LBNL) AmeriFlux dataset for US-NR1 (i.e., https://doi.org/10.17190/AMF/1246088) and updates what was in the 2020 ESS-DIVE US-NR1 archive (https://doi.org/10.15485/1671825) to include data from the years 2020-2025. More specifically, the following updates are provided: (i) five-minute statistics (means, variances, covariances) of all data measured by the US-NR1 data system between Sep 2020 and Jun 2025 in netCDF format, (ii) the electronic logbook of US-NR1 site visits, (iii) a web calendar (in HTML format) documenting activity at the site (a replica of https://urquell.colorado.edu/calendar/), (iv) photos taken at the site between years 2020 and present day (Aug 2025), and (v) several auxiliary datasets, primary related to trees near the site, soil properties, soil moisture and soil temperature, and subcanopy radiation data. The data package is setup so that the web calendar, photos, and electronic logbook can be easily accessed on a local computer using a web browser. The provided data files are in either BINEX or SBF format (for the raw GNSS data), netCDF, CSV, ASCII, or MATLAB format. To obtain a better understanding about the archive, please start by reading the following PDF which is included within the data archive:README_ESS_DIVE_USNR1_2025_readme_first.pdf.
Stable water isotopes (d18O, d2H and d-excess) are important tracers in hydrologic research to understand water partitioning between vegetation, groundwater, and runoff but are rarely applied to large watersheds with persistent snowpack and complex topopgraphy. Data were collected for the Lawrence Berkeley National Laboratory Watershed Function Science Focus Area supported by the U.S Department of Energy in the East River, CO Hydrologic Unit Code (140200010) with limited data also collected in adjacent watersheds Ohio Creek and Taylor River. Data are provided in csv and includes isotopic information for precipitation (years 2014-2016), snowmelt (years 2016-2017) and snowpits (years 2016-2020). Snowpit data contain depth resolved information at 10 cm intervals for density, snow water equivalent (SWE) and stable water isotopes. Bulk isotopic data for 86 snowpits contain depth, SWE, density and SWE-weighted isotope values. Sampling locations and elevations are provided within the data files, while kmz files are provided to view sampling locations using Google Earth software.
A concept for a solar radiometer network to provide supporting data during the Satellite Preciptiation and Cloud Experiment (SPACE) was developed. Each of the 9 prime and 10 supplementary SPACE ground sites will be equipped with an upward pointing global solar pyranometer. About half of the sites will also be equipped with upward pointing diffuse (shade ring) solar pyranometers, and a downward pointing global albedo pyranometer. These radiometers will be used to monitor the spatial and temporal variability of solar insolation and haze optical depth. The insolation data will ultimately be input to numerical models of the pre-storm and near-storm boundary layer. The optical depth data will be compared with simultaneous measurements from airborne and satellite-based passive visible radiometers and airborne lidars.
The study of sea ice using airborne remote sensing platforms provides unique capabilities to measure a wide variety of sea ice properties. These measurements are useful for a variety of topics including model evaluation and improvement, assessment of satellite retrievals, and incorporation into climate data records for analysis of interannual variability and long-term trends in sea ice properties. In this paper we describe methods for the retrieval of sea ice thickness, freeboard, and snow depth using data from a multisensor suite of instruments on NASA's Operation IceBridge airborne campaign. We assess the consistency of the results through comparison with independent data sets that demonstrate that the IceBridge products are capable of providing a reliable record of snow depth and sea ice thickness. We explore the impact of inter-campaign instrument changes and associated algorithm adaptations as well as the applicability of the adapted algorithms to the ongoing IceBridge mission. The uncertainties associated with the retrieval methods are determined and placed in the context of their impact on the retrieved sea ice thickness. Lastly, we present results for the 2009 and 2010 IceBridge campaigns, which are currently available in product form via the National Snow and Ice Data Center
Asthma hospital/emergency room (patient) data are used as the foundation for creating a health outcome indicator of human response to environmental air quality. Daily U.S. Environmental Protection Agency (EPA) Air Quality System (AQS) fine particulates (PM2.5) ground data and the U.S. National Aeronautical Space Administration (NASA) MODIS aerosol optical depth (AOD) data were acquired and processed for years of 2007 and 2008. Figure 1 shows the PM2.5 annual mean composite of all the 2007 B-spline daily surfaces. Initial models for predicting the number of weekly asthma cases within a Florida county has focused on environmental variables. Weekly maximums of PM2.5, relative humidity, and the proportions of the county with smoke and fire were the environmental variables included in the model. Cosine and sine functions of time were used to account for seasonality in asthma cases. Counties were considered to be random effects, thereby adjusting for differences in socio ]demographics and other factors. The 2007 predictions for Miami ]Dade county when using B ]splines PM2.5 are displayed in Figures 2.
The 11 th Society of Petroleum Engineers Comparative Solution Project (shortened SPE11 herein) benchmarked simulation tools for geological carbon dioxide (CO 2 ) storage. A total of 45 groups from leading research institutions and industry across the globe signed up to participate, with 18 ultimately contributing valid results that were included in the comparative study reported here. This paper summarizes the SPE11 results. A comprehensive introduction and qualitative discussion of the submitted data are provided, together with an overview of online resources for accessing the full depth of data. A global metric for analyzing the relative distance between submissions is proposed and used to conduct a quantitative analysis of the submissions. This analysis attempts to statistically resolve the key aspects influencing the variability between submissions. The study shows that the major qualitative variation between the submitted results is related to thermal effects, dissolution-driven convective mixing, and resolution of facies discontinuities. Moreover, a strong dependence on grid resolution is observed across all three versions of the SPE11. However, our quantitative analysis suggests that the observed variations are predominantly influenced by factors not documented in the technical responses provided by the participants. We therefore identify that unreported variations due to human choices within the process of setting up, conducting, and reporting on the simulations underlying each SPE11 submission are at least as impactful as the computational choices reported.
The establishment of a limited spatial data base for the U.S. eastern seaboard vicinity is discussed along with the demonstration of computer assisted analysis techniques for investigating spatial patterns and relationships among ship survey data and remotely sensed sea surface temperature. Ship survey variables included concentrations of two zooplankton, two icthyoplankton, and two fish species, in addition to physical data of depth to bottom and surface and bottom water temperatures. Continuous spatial distributions of these data were created by both weighted nearest neighbor and iterative smoothing interpolation techniques. Maps of surface water temperature were created by digitizing GOES satellite images. All mapped data were spatially registered by conversion of latitude and longitude coordinates to rotated Lambert conic conformal rectilinear coordinates and stored in grid format of approximately one hundred square kilometers per cell. The analysis of these data include the generation of statistical summaries and maps describing the joint occurrence among variables.
Climate forcing by anthropogenic aerosols is a significant but highly uncertain factor in global climate change. Only satellites can offer the global coverage essential to reducing this uncertainty; however, satellite measurements must be coupled with correlative, in situ measurements both to constrain the aerosol optical properties required in satellite retrieval algorithms and to provide chemical identification of aerosol sources. This grant funded the third year of a three-year project which seeks to develop methodologies for combining spaceborne lidar with in-situ aerosol data sets to improve estimates of direct aerosol climate forcing. Progress under this one-year grant consisted in analysis and publication of field studies using a new in-situ capability for measuring aerosol 180 deg backscatter and the extinction-to-backscatter ratio. This new measurement capacity allows definitive lidar/in-situ comparisons and improves our ability to interpret lidar data in terms of climatically relevant quantities such as the extinction coefficient and optical depth. Analyzed data consisted of measurements made along the coast of Washington State, in Central Illinois, over the Indian Ocean, and in the Central Pacific. Thus, this research, combined with previous measurements by others, is rapidly building toward a global data set of extinction-to-backscatter ratio for key aerosol types. Such information will be critical to interpreting lidar data from the upcoming PICASSO-CENA, or P-C, satellite mission. Another aspect of this project is to investigate innovative ways to couple the lidar-satellite signal with target in-situ measurements toward a direct determination of aerosol forcing. This aspect is progressing in collaboration with NASA Langley's P-C lidar simulator.
An investigation was designed for the Skylab satellite to determine the feasibility of remote sensing of the soil moisture content of the surface from various microwave sensors. Skylab data for the experiment were collected during passes 5, 10, 16, and 38 across the two test sites selected in eastern Kansas and western Texas. Pass 38 covered both test sites giving five data sets for the analysis. As Skylab data were being taken the moisture content of the soil was sampled by ground crews for each 2.5 centimeter depth from the surface to 15 centimeters at interval of about six kilometers along two different routes along the test sites. This resulted in a total of 2250 soil moisture samples corresponding to different locations and six different depths. Skylab data were collected by passive microwave radiometers at wavelengths of 2.1 and 21 centimeters by the S193 and S194 microwave sensors. An active microwave system also collected scatterometer data at a wavelength of 2.1 centimeters. The analysis of microwave data has revealed that the longer wavelength L-Band passive radiometer gives the best correlation with soil moisture content of the upper 2.5 centimeter depth of soil.
These data accompany the published paper Lewis et al., 202X and are from a brackish coastal wetland in situ soil warming experiment (GENX) equipped with automated flux chambers. Methane (CH4) and carbon dioxide (CO2) fluxes were measured in 12 automated chambers using custom-built automated chambers connected to an LI-7810 CH4/CO2 analyzer. The chambers are 1.5 m tall and contain the dominant vegetation species of the site (Schoenoplectus americanus, Spartina patens, and Distichlis spicata). The chambers are also distributed across a soil warming gradient, ranging from ambient to 6°C above ambient, that was started in February 2022. This dataset contains the following files: (1) CH4 and CO2 fluxes from each chamber for March to November 2025, statistics for each flux, and environmental data (water depth, salinity, air temperature) at the time of the flux measurement; (2) 15-minute soil temperature data for each chamber; (3) Aboveground vegetation biomass (total and by species) and stem counts and dimensions for S. americanus; (4) Elevation for each chamber. All data processing code is available on Github.
This work investigates multi-agent path planning in strong, dynamic currents using thousands of highly under-actuated vehicles. We address the specific task of path planning for a global network of ocean-observing floats. These submersibles are typified by the Argo global network consisting of over 3000 sensor platforms. They can control their buoyancy to float at depth for data collection or rise to the surface for satellite communications. Currently, floats drift at a constant depth regardless of the local currents. However, accurate current forecasts have become available which present the possibility of intentionally controlling floats' motion by dynamically commanding them to linger at different depths. This project explores the use of these current predictions to direct float networks to some desired final formation or position. It presents multiple algorithms for such path optimization and demonstrates their advantage over the standard approach of constant-depth drifting.
The stratospheric aerosol measurement 2, stratospheric aerosol and gas experiment (SAGE) 1, and SAGE 2 series of solar occultation satellite instruments were designed for the study of stratospheric aerosols and gases and have been extensively validated in the stratosphere. They are also capable, under cloud-free conditions, of measuring the extinction due to aerosols in the troposphere. Such tropospheric extinction measurements have yet to be validated by appropriate lidar and in situ techniques. In this paper published atmospheric aerosol optical depth measurements, made from high-altitude observatories during volcanically quiet periods, have been compared with optical depths calculated from local SAGE 1 and SAGE 2 extinction profiles. Surface measurements from three such observatories have been used, one located in Hawaii and two within the continental United States. Data have been intercompared on a seasonal basis at wave-lenths between 0.5 and 1.0 micron and found to agree within the range of measurement errors and expected atmospheric variation. The mean rms difference between the optical depths for corresponding satellite and surface measured data sets is 29%, and the mean ratio of the optical depths is 1.09.
ERTS-1 MSS data taken on October 10, 1972 of the Little Bahama Bank are being used to demonstrate the use of ERTS-1 data for mapping of shallow water features for the purpose of upgrading world navigation charts. Marked reflectance differences occur for the shallow water areas in Bands 4, 5, and 6. Digital processing of two adjacent data tapes within the ERTS frame covering an area of about 40 by 40 miles has been completed. Correlation of depth measurements to 5 meters has been successful. A mathematical model for depth measurements using ratio of voltages in Band 4 and 5 has been successfully developed and is being tested for accuracy. Additional studies for areas near Puerto Rico and in Northern Lake Michigan will be undertaken. Satellite data will also provide geographical evidence for verifying existence or nonexistence of doubtful shoal waters now appearing on world charts and considered to be hazardous to shipping.
The Mars Science Laboratory (MSL) was protected during its Mars atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA). PICA is a lightweight carbon fiber/polymeric resin material that offers excellent performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of atmospheric entry and material response. MEDLI recorded, among others, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). These measurements have been widely used in the literature as a validation benchmark for state-of-the-art ablation codes. The objective of this work is to perform an inverse estimate of the MSL heatshield material properties and aerothermal environment during Mars entry from the MISP flight data.
Introduction: NASA’s next mission to Mars, the Mars 2020, will use the same heatshield of the Mars Science Laboratory (MSL) for thermal protection during entry, descent and landing. The heatshield is a tiled system made of Phenolic Impregnated Carbon Ablators (PICA) blocks [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers excellent performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of atmospheric entry and material response. MEDLI recorded, among others, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to compare the thermal response of the MSL heatshield to the MISP flight data. In preparation to Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against MEDLI flight data.
The Mars Science Laboratory (MSL) was protected during its Mars atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA). PICA is a lightweight carbon fiber/polymeric resin material that offers excellent performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of atmospheric entry and material response. MEDLI recorded, among others, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to showcase the capability of the Design, Analysis, and Optimization of Thermal Protection Materials (DAOTPM) software. DAO-TPM is a Python based framework that works as a link between mission design, aerothermal and radiative environment computation, Thermal Protection Systems (TPS) microstructure analysis, material response and optimization tools. The toolbox has a Graphical User Interface (GUI) that allows the user to build as well as run the various software and utilities used to design, analyze and optimize a heatshield during atmospheric entry.
Abstract Wind-wave development is governed by the fetch- or duration-limited growth principle that is expressed as a pair of similarity functions relating the dimensionless elevation variance (wave energy) and spectral peak frequency to fetch or duration. Combining the pair of similarity functions, the fetch or duration variable can be removed to form a dimensionless function of elevation variance and spectral peak frequency, which is interpreted as the wave energy evolution with wave age. The relationship is initially developed for quasi-neural stability and quasi-steady wind forcing conditions. Further analyses show that the same fetch, duration, and wave-age similarity functions are applicable to unsteady wind forcing conditions, including rapidly accelerating and decelerating mountain gap wind episodes and tropical cyclone (TC) wind fields. Here it is shown that with the dimensionless frequency converted to dimensionless wavenumber using the surface wave dispersion relationship, the same similarity function is applicable in all water depths. Field data collected in shallow to deep waters and mild to TC wind conditions and synthetic data generated by spectrum model computations are assembled to illustrate the applicability. For the simulation work, the finite-depth wind-wave spectrum model and its shoaling function are formulated for variable spectral slopes. Given wind speed, wave age, and water depth, the measured and spectrum-computed significant wave heights and the associated growth parameters are in good agreement in forcing conditions from mild to TC winds and in all depths from deep ocean to shallow lake. Significance Statement This paper presents a growth function and spectrum model to describe wind-wave development in all water depths. Their applicability covers a wide range of wind forcing conditions including steady, accelerating, decelerating, and tropical cyclone events. Support for the unified spectrum model and growth function is presented with field observations and numerical computations.