Using Precipitation, Soil Moisture, & Snow Cover Observations to Constrain the Land Surface Water Cycle in the NASA GEOS Modeling and Assimilation Framework
Explore the source record for details and available documents.
Rolf Reichle↗
SEARCH · Search NASA
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.
MODIS Collection 6 MAIAC Algorithm
This paper describes the latest version of the algorithm MAIAC (Multi-Angle Implementation of Atmospheric Correction) used for processing the MODIS (Moderate-resolution Imaging Spectroradiometer) Collection6 data record. Since initial publication in 2011-2012, MAIAC has changed considerably to adapt to global processing and improve cloud/snow detection, aerosol retrievals and atmospheric correction of MODIS data. The main changes include (1) transition from a 25 to 1 km scale for retrieval of the spectral regression coefficient (SRC) which helped to remove occasional blockiness at 25 km scale in the aerosol optical depth (AOD) and in the surface reflectance, (2) continuous improvements of cloud detection, (3) introduction of smoke and dust tests to discriminate absorbing fine- and coarse mode aerosols, (4) adding over-water processing, (5) general optimization of the LUT (LookUp-Table)-based radiative transfer for the global processing, and others. MAIAC provides an interdisciplinary suite of atmospheric and land products, including cloud mask (CM), column water vapor (CWV), AOD at 0.47 and 0.55 m, aerosol type (background, smoke or dust) and fine-mode fraction over water; spectral bidirectional reflectance factors (BRF), parameters of Ross-thick Lisparse (RTLS) bidirectional reflectance distribution function (BRDF) model and instantaneous albedo. For snow-covered surfaces, we provide subpixel snow fraction and snow grain size. All products come in standard HDF4 (software library) format at 1 km resolution, except for BRF, which is also provided at 500 m resolution on a sinusoidal grid adopted by the MODIS Land team. All products are provided on per-observation basis in daily files except for the BRDF/Albedo product, which is reported every 8 days. Because MAIAC uses a time series approach, BRDF/Albedo is naturally gap-filled over land where missing values are filled-in with results from the previous retrieval. While the BRDF model is reported for MODIS Land bands 1-7 and ocean band 8, BRF is reported for both land and ocean bands 1-12. This paper focuses on MAIAC cloud detection, aerosol retrievals and atmospheric correction and describes MCD19 data products and quality assurance (QA) flags.
MAIAC Algorithm↗
Advances in Land Data Assimilation at the NASA Goddard Space Flight Center
Research in land surface data assimilation has grown rapidly over the last decade. In this presentation we provide a brief overview of key research contributions by the NASA Goddard Space Flight Center (GSFC). The GSFC contributions to land assimilation primarily include the continued development and application of the Land Information System (US) and the ensemble Kalman filter (EnKF). In particular, we have developed a method to generate perturbation fields that are correlated in space, time, and across variables and that permit the flexible modeling of errors in land surface models and observations, along with an adaptive filtering approach that estimates observation and model error input parameters. A percentile-based scaling method that addresses soil moisture biases in model and observational estimates opened the path to the successful application of land data assimilation to satellite retrievals of surface soil moisture. Assimilation of AMSR-E surface soil moisture retrievals into the NASA Catchment model provided superior surface and root zone assimilation products (when validated against in situ measurements and compared to the model estimates or satellite observations alone). The multi-model capabilities of US were used to investigate the role of subsurface physics in the assimilation of surface soil moisture observations. Results indicate that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Building on this experience, GSFC leads the development of the Level 4 Surface and Root-Zone Soil Moisture (L4_SM) product for the planned NASA Soil-Moisture-Active-Passive (SMAP) mission. A key milestone was the design and execution of an Observing System Simulation Experiment that quantified the contribution of soil moisture retrievals to land data assimilation products as a function of retrieval and land model skill and yielded an estimate of the error budget for the SMAP L4_SM product. Terrestrial water storage observations from GRACE satellite system were also successfully assimilated into the NASA Catchment model and provided improved estimates of groundwater variability when compared to the model estimates alone. Moreover, satellite-based land surface temperature (LST) observations from the ISCCP archive were assimilated using a bias estimation module that was specifically designed for LST assimilation. As with soil moisture, LST assimilation provides modest yet statistically significant improvements when compared to the model or satellite observations alone. To achieve the improvement, however, the LST assimilation algorithm must be adapted to the specific formulation of LST in the land model. An improved method for the assimilation of snow cover observations was also developed. Finally, the coupling of LIS to the mesoscale Weather Research and Forecasting (WRF) model enabled investigations into how the sensitivity of land-atmosphere interactions to the specific choice of planetary boundary layer scheme and land surface model varies across surface moisture regimes, and how it can be quantified and evaluated against observations. The on-going development and integration of land assimilation modules into the Land Information System will enable the use of GSFC software with a variety of land models and make it accessible to the research community.
Reichle, Rolf↗
Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures
In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.
Haokui Xu↗
Global Assimilation of Multi-Sensor Snow Observations for Improved Characterization of Snow Processes
Snow conditions on the land surface are recognized to be key components of the global hydrological cycle as they play a critical role in the determination of local and regional climate. In many mid-latitude and high-latitude regions, the seasonal water storage and associated spring snowmelt dominate the local hydrology. The contribution to the runoff and moisture conditions from snow is vital in supporting agriculture and in determining water resources management practices. Consequently, accurate characterization of snow properties becomes important for both end-use applications and weather and climate research. Recently a joint effort between the u.S. Air Force and NASA has enabled a blended, multi-sensor snow product known as the AFWA NASA Snow Algorithm (ANSA). This global snow dataset has been generated by utilizing the Earth Observation System (EOS) Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Microwave Scanning Radiometer for EOS (AMSR-E) datasets. ANSA product includes estimates of snow cover extent, snow water equivalent (SWE) and SWE-derived snow depth fields. The MODIS-based products enable snow cover mappings under cloud-free conditions whereas the passive microwave data from AMSR-E provides measurements under cloudy conditions. These remotely-sensed snow observations are further augmented with the information from ground-based snow measurements through data fusion techniques. The resulting ANSA products are employed in the NASA Land Information System (LIS) data assimilation framework, which provides a comprehensive environment for integrating community land surface models, ground and satellite-based observations, and ensemble-based data assimilation tools. LIS incorporates the multisensor ANSA snow retrievals with the land surface model estimates to generate spatially and temporally continuous estimates of snow states, through data assimilation. A suite of experiments to assimilate ANSA snow cover, SWE and snow depth estimates with different land surface models in LIS are conducted and the resulting estimates of snow conditions are evaluated against a number of in-situ observational datasets, over several regions of the world. These evaluations are used to compare and contrast the advantages and disadvantages of these multi-sensor snow observations.
Kumar, Sujay↗
Remote sensing of ice and snow - Review and status
The interpretation of satellite data is reviewed to address the possibility of extracting geophysically useful data regarding ice and snow. The scientific and practical utility of such data is examined in terms of the required and available scales, resolutions, and data types. Satellite instrumentation for ice and snow studies includes the AVHRR, the Special Sensor Microwave Imager, and the altimeters on Geosat and the European Remote Sensing Satellite. The satellites are presently used for conducting investigations of the microwave radiometry of snow, shortwave snowcover inventories, and sensor and model validations. It is noted that these studies encompass a range of ice properties and employ approaches and scientists from a number of fields. Areas that remain to be addressed in the field include snow structure and substrate effects, fall cooling, floe flooding, and brine expulsion.
Carsey, Frank↗
Shape from spectra
We introduce a new unified atmospheric–topographic correction approach that estimates surface geometry directly from the radiance measurement. Surface topography influences the at-sensor radiance measurement, making precise topography modeling critical in applications like vegetation or snow studies in mountainous terrain. Currently, elevation maps are used to derive topographic variables such as the slope and sky-view factor. This process is error-prone since static global digital elevation models do not generally achieve the accuracy required, and even minor mismatches in spatial resolution can introduce significant artifacts in downstream processing. Here we demonstrate that it is possible to estimate topographic parameters directly from spectral data, ensuring perfect physical consistency, temporal coincidence, and spatial alignment. We present experiments estimating topographic slope in two scenes in Southern California, with data from NASA’s Next Generation Airborne Visible/Near Infrared Imaging Spectrometer (AVIRIS-NG). We compared our radiance-based estimates against high-resolution lidar datasets. Our initial validation result showed a correlation of R 2 = 0.864 (n = 160) over the homogeneous surface of Beckman Auditorium’s cone-shaped roof on the Caltech campus in Pasadena, California. We then validate the model over a larger study site near Santa Clarita, California, finding R 2 = 0.923 (n = 40,000) in a 350 x 350m area. The accuracy of our model estimates, combined with its systematic advantages over the alternative, show the potential of the approach for use in both airborne campaigns and orbital missions.
Nimrod Carmon↗
Cloud-atmosphere impacts on the central Arctic surface energy budget (Final Report)
This project investigated how clouds and atmosphere-sea ice interactions shape the central Arctic surface energy budget using observations from the MOSAiC expedition and coordinated model experiments. Analyses focused on three themes. First, a new method applied to thermistor string data revealed that snow thermal conductivity is systematically higher than assumed in most models, strongly influencing conductive heat flux and sea ice growth, especially under clear skies. Second, a full annual cycle of surface energy budget observations, combined with detailed cloud microphysical retrievals, demonstrated how cloud regimes regulate radiative, turbulent, and conductive fluxes, driving seasonal contrasts in sea ice energy balance: in winter, radiative forcing elicits strong turbulent and conductive responses, while in summer excess energy is primarily partitioned into surface melt. Third, these observations were used to evaluate and improve regional and global forecast models. While some models captured key Arctic cloud-radiation states, most exhibited persistent deficiencies in representing liquid-containing clouds, snow-on-sea-ice processes, and flux responses to radiative forcing. Together, these results provide unprecedented benchmarks for understanding coupled Arctic processes, improve parameterizations of snow, sea ice, and cloud interactions, and highlight critical pathways by which atmospheric variability drives sea ice change in a rapidly evolving Arctic system. Lastly, the project was a resounding success, resulting in 38 peer-reviewed publications and dozens of presentations, while also supporting an early career scientist.
54 ENVIRONMENTAL SCIENCES↗
Southeast Greenland Winter Precipitation Strongly Linked to the Icelandic Low Position
Greenland's largest precipitation flux occurs in its southeast (SE) region during the winter, controlled primarily by easterly winds and frequent cyclogenesis in the North Atlantic. Several studies have attempted to link SE Greenland precipitation to the North Atlantic Oscillation (NAO) but results are inconsistent. This work uses reanalysis, automatic weather station data, and regional climate model output to show that the east-west position of the Icelandic low is a better predictor of SE Greenland precipitation (average correlation of r equals -0.48 in DJF (December January February)) than climate indices such as the NAO (r equals -0.06 in DJF). In years when the Icelandic low is positioned extremely west, moisture transport increases up to approximately 40 percent (or up to 40 kilograms per meter per second) off the SE Greenland coast compared to when the low is in an extreme east position. Furthermore, in years when the Icelandic low is positioned extremely west, storm track density and intensity increase just off the SE coast of Greenland. Thus, the Icelandic low's longitudinal position dominates SE Greenland ice sheet's wintertime precipitation, a positive term in the ice sheet mass balance. Given SE Greenland's importance in the overall ice sheet mass balance, the position of the Icelandic low is therefore important for making projections of future sea level.
Arctic↗
International Symposium on Remote Sensing of Environment, 17th, University of Michigan, Ann Arbor, MI, May 9-13, 1983, Proceedings. Volumes 1, 2 & 3
The collection, processing, and analysis of remote-sensing data from ground-based, airborne, and spaceborne instruments for application to the monitoring and management of the earth and environment and resources are examined in reviews and reports, some in summary form. Subject areas covered include US policy and directions on remote sensing (RS); the future of terrestrial RS from space; RS of land, oceans, and atmosphere from a global perspective; RS in hydrological modeling; microprocessing technology; array processors; geobased information systems; artificial intelligence; the Shuttle imaging radar; and current results from Landsat-4. Among the specific topics discussed are RS application to hydrocarbon exploration, airborne gamma-radiation assessment of snow water equivalent, surface-vegetation-biomass modeling from AVHRR and Landsat data, Landsat imagery of Mediterranean pollution, fast two-dimensional filtering of thermal-scanner data, RS of severe convective storms, registration of rotated images by invariant moments, and the geometric accuracy of Landsat-4 Thematic-Mapper P-tapes.
Source record↗
Mapping Firn Saturation over Greenland using NASA’s Soil Moisture Active Passive Satellite
Mapping the spatial extent of recently identified englacial hydrological features (i.e., ice slabs and perennial firn aquifers) formed by meters-thick water-saturated firn layers over the percolation facies of the Greenland Ice Sheet using L-band microwave radiometry has recently been demonstrated. However, these initial maps are binary, and do not provide a parameter to estimate the spatial variability in the thickness and volumetric fraction of meltwater stored within the firn pore space. Here, we exploit enhanced-resolution vertical-polarization L-band brightness temperature (T^(B)_(V)) imagery (2015–2019) generated using observations collected over Greenland by NASA’s Soil Moisture Active Passive (SMAP) satellite and a simple two-layer L-band brightness temperature model. We map water-saturated firn layers via a “firn saturation” parameter, and interpret our results together with ice slab and perennial firn aquifer spatial extents, estimates of snow accumulation simulated via the Regional Atmospheric Climate Model (RACMOp2.3), and airborne radar surveys collected via NASA’s Operation IceBridge (OIB) campaigns. We find that variable firn saturation parameter values are mapped in lower snow accumulation ice slab areas in western, northern, and northeastern Greenland, where firn is colder and water-saturated firn layers seasonally refreeze as solid-ice. Higher firn saturation parameter values are mapped in higher snow accumulation perennial firn aquifer areas in southeastern, southern, and northwestern Greenland, where firn is near the melting point, and meters-thick water-saturated firn layers exist. Our results have implications for identifying expansive englacial reservoirs that store significant volumes of meltwater in locations that are vulnerable to meltwater-induced hydrofracturing and accelerated outlet glacier flow year-round.
Greenland ice sheet↗
Refinements to SSiB with an Emphasis on Snow-Physics: Evaluation and Validation Using GSWP and Valdai Data
Refinements to the snow-physics scheme of SSiB (Simplified Simple Biosphere Model) are described and evaluated. The upgrades include a partial redesign of the conceptual architecture to better simulate the diurnal temperature of the snow surface. For a deep snowpack, there are two separate prognostic temperature snow layers - the top layer responds to diurnal fluctuations in the surface forcing, while the deep layer exhibits a slowly varying response. In addition, the use of a very deep soil temperature and a treatment of snow aging with its influence on snow density is parameterized and evaluated. The upgraded snow scheme produces better timing of snow melt in GSWP-style simulations using ISLSCP Initiative I data for 1987-1988 in the Russian Wheat Belt region. To simulate more realistic runoff in regions with high orographic variability, additional improvements are made to SSiB's soil hydrology. These improvements include an orography-based surface runoff scheme as well as interaction with a water table below SSiB's three soil layers. The addition of these parameterizations further help to simulate more realistic runoff and accompanying prognostic soil moisture fields in the GSWP-style simulations. In intercomparisons of the performance of the new snow-physics SSiB with its earlier versions using an 18-year single-site dataset from Valdai Russia, the version of SSiB described in this paper again produces the earliest onset of snow melt. Soil moisture and deep soil temperatures also compare favorably with observations.
Mocko, David M.↗
The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach
Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.
Wang, Lijing [University of Connecticut, Storrs, C↗
Surface Quantitative Precipitation Estimates (SQUIRE) of Snow Water Equivalent from the Surface Atmospheric Integrated Field Laboratory
The upper Colorado River basin is the primary source of water for 40 million people. With declining snowpack in the basin, forecasting hydrological budgets in the Southwest United States is more important than ever. However, due in part, to a lack of reliable observations of precipitation in complex terrain, hydrological models struggle to assess and forecast snowpack snow water equivalent (SWE) in the upper Colorado River basin (UCRB). Therefore, the need for more reliable SWE forecasts in the UCRB motivated the U.S. Department of Energy Atmospheric Radiation Measurement Facility’s Surface Atmospheric Integrated Field Laboratory (SAIL) that occurred from June 2021 to June 2023. During SAIL, the X-band precipitation radar from Colorado State University conducted volume scans sampling the precipitation properties over the UCRB. The ARM facility developed a gridded Surface Quantitative Precipitation Estimates (SQUIRE) product from the radar observations. To do this, various daily SWE estimates from radar using the radar reflectivity factor Z e and specific differential phase K dp were compared against ground-based precipitation gauges. SWE in precipitation calculated from Wolfe and Snider’s S–Z e estimator was in best agreement with the rain gauges for the days when SWE < 12 mm. For days with SWE > 12 mm, the WSR-88D Intermountain West relationship had the best agreement with the precipitation gauges. Airborne snow depth observations show that SQUIRE captures regions of orographic enhancement in the mountains to the west and northwest of the SAIL study area, indicating that the scientific community should focus on understanding and ultimately simulating orographic atmospheric precipitation processes to improve UCRB snowpack SWE assessment and forecasting.
Hydrology↗
The effect of surface anisotropy on the accuracy of total ozone estimates from satellite observations
The total amount of ozone in a vertical column of the earth's atmosphere is being derived from satellite measurements of the intensity of ultraviolet sunlight scattered by the earth-atmosphere system. The algorithm for deriving the ozone amount utilizes the assumption that the earth's surface reflects the incident light isotropically according to Lambert's law. Natural surface reflection deviates more or less from this law. Two extreme examples of anisotropic reflection from dark ocean and from bright snow are analyzed by means of models for their effects on the derived values of ozone.
Fraser, R. S.↗
The anomalous extinction curve in the direction of Rho Ophiuchi from 950 to 1180 A
We report the interstellar extinction from 950 to 1180 A in the direction of Rho Oph. These data were obtained with a moderate-resolution sounding rocket observation of Ro Oph, and two lightly reddened B stars of identical spectral type. Unlike Voyager observations in this wavelength region, our spectral resolution is sufficient to directly measure the contribution of molecular hydrogen along the line of sight and unambiguously identify the dust contribution. The derived curves compare favorably with the extrapolation of functional fits to the Rho Oph derived from longer wavelengths. In addition, we have fitted a model of dust extinction to the derived curve and find that the best fits indicate that the line of sight exhibits a larger than average grain size and a preferential concentration of graphite relative to silicates. The sharp rise seen at the shortest wavelengths (below 1000 A) cannot be explained by standard silicate-graphite dust models and may be indicative of a nonstandard grain shape or a missing component in the extinction models.
Green, James C.↗
New Total Ozone Algorithm for Application to the Total Ozone Mapping Spectrometer
The Total Ozone Mapping Spectrometer (TOMS) series comprises four instruments providing a total of 25 years of daily global stratospheric ozone data over the sunlit portion of the Earth. A new retrieval algorithm has been developed for TOMS, designated Version 8. The algorithm is based on differential absorption across a pair of wavelength channels chosen close together to minimize the impact of wavelength dependent forward modeling errors. Version 8 enhancements include correction for the presence of tropospheric aerosols and sun glint from water surfaces, a better treatment of variability due to tropospheric ozone and temperature dependence, and an improved forward model, particularly in regions of persistent snow and ice. Among other things, the Version 8 enhancements have reduced latitudinal dependence seen previously in TOMS - Dobson comparisons, predominantly in the Southern Hemisphere's summer, when the tropospheric ozone, temperature, and snow/ice corrections are additive. The basic components of the algorithm and its impact on derived total ozone will be discussed.
Wellemeyer, C.↗
Two Successive Martian Years on the Orbit: Similarities and Differences of CO2 Seasonal Cycle from HEND/ODYSSEY Data
The three years of Mars Odyssey successful work on the martian orbit provide a lot of new information about peculiarities of long term variations of CO2 seasonal cycle. To start such analysis we have used observations of neutron albedo of Mars obtained by High Energy Neutron detector (HEND) mounted onboard Mars Odyssey spacecraft. The high latitude northern and southern regions of Mars are affected by global redistribution of atmospheric CO2 which resulted in 25% of atmospheric mass condensed on martian surface of these regions during winter period of time. The seasonal deposit is formed starting from 60N/60S latitudes and achieve its maximal thickness about 1 m at latitudes close to martian poles. Changes of CO2 deposit thickness is the reason for significant variations of neutron flux above martian poles from summer to winter seasons because CO2 frost effectively hides upper water rich surface layers from the orbit observations in neutrons and gamma-rays. This effect was used to estimate column density of CO2 deposit at different latitudes on North and South of Mars and reconstruct multidimensional model of CO2 deposit showing how snow depth varies as function of latitude, longitude and time. In this presentation we tried to make a next step in our study of martian seasonal CO2 cycle and look for similarities and differences between two successive martian years.
Litvak, M. L.↗