Search NASA⌕ Search

SEARCH · Search NASA

Results for “soil heterogeneity”

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.

139 records · Page 8

Mineralogy of the Martian Surface: Crustal Composition to Surface Processes

The main results have been published in the refereed literature, and thus this report serves mainly to summarize the main findings and indicate where the detailed papers may be found. Reflectance spectroscopy has been an important tool for determining the mineralogic makeup of the near surface materials on Mars. Analysis of the spectral properties of the surface have demonstrated that these attributes are heterogeneous from the coarse spatial but high spectral resolution spectra obtained with telescopes to the high spatial but coarse spectral resolution Viking data (e.g. Arvidson et al., 1989; McEwen et al., 1989). Low albedo materials show strong evidence for the presence of igneous rock forming minerals while bright materials are generally interpreted as representing heavily altered crustal material. How these materials are physically and genetically related has important implications for understanding martian surface properties and processes, weathering histories and paths, and crustal composition. The goal of this research is to characterize the physical and chemical properties of low albedo materials on Mars and the relationship to intermediate and high albedo materials. Fundamental science questions to be pursued include: (1) the observed distributions of soil, rock, and dust a function of physical processes or weathering and (2) different stages of chemical and physical alteration fresh rock identified. These objectives will be addressed through detailed analyses and modelling of the ISM data from the Phobos-2 mission with corroborating evidence of surface composition and properties provided by data from the Viking mission.

Mustard, John F.↗

Integrating the Interconnections Between Groundwater and Land Surface Processes Through the Coupled NASA Land Information System and ParFlow Environment

Understanding the interactions between the atmosphere, the land, and the subsurface is fundamental to hydrology and is critical for a better assessment of the impacts of climate change and human management on hydrological systems. However, many land surface models simplify the subsurface hydrology and thereby these interactions. In this study, we couple the land surface model Noah-MP included in the NASA Land Information System (LIS) with the integrated hydrologic model ParFlow (ParFlow-LIS) using the Earth System Modeling Framework (ESMF) and the National United Operational Prediction Capability (NUOPC). This coupling improves the simulation of water and energy cycle processes by adding the three-dimensional variably saturated and heterogeneous flow in the subsurface using sophisticated and nonlinear physics-based equations as well as the advances in satellite remote sensing-based data assimilation of the land surface, thereby benefiting the integrated hydrologic modeling and data assimilation community. We use the High Plains aquifer, located in the central United States, as a testbed to evaluate the coupled ParFlow-LIS system. The new ParFlow-LIS system accounts for the effects of topographically driven flows on the land surface, producing more fine-scale patterns of land surface states and fluxes than standalone LIS. In addition, ParFlow-LIS enables the consideration of the effect of subsurface water storage on evapotranspiration. This is particularly important in areas and times with dry soils, such as during drought conditions or in the presence of a cone of depression due to pumping.

land surface modeling↗

Informing Improvements in Freeze/Thaw State Classification Using Subpixel Temperature

Freeze/thaw (FT) processes at the earth’s surface can have a considerable effect on global carbon, energy, and hydrologic cycles. Therefore, an accurate representation of FT is valuable to adequately monitor and model these processes. In this study, we assess the relationship between satellite-based FT products and modeled surface and soil temperatures over North America. In addition, hourly land surface temperature (LST) from the Geostationary Operational Environmental Satellite (GOES) system is also compared to FT classifications. Utilizing the higher spatial resolution temperatures (5 km), we assess subgrid-scale variability and its relationship to coarser microwave FT classifications (>25 km). We also examine product agreement and subpixel characteristics across the land cover, climate, and topography. FT classifications are shown to vary widely depending on these variables, leading to an ambiguous definition of frozen and thawed states. Our results suggest that current products can characterize FT transitions with consistent subfreezing surface characteristics in far northern regions (>50 °N). However, uncertainty associated with FT classifications is shown to increase considerably as latitude decreases. Our results also suggest that fractional FT products, utilizing data inputs, such as LST, would provide a considerable improvement in mountainous regions with high intergrid cell heterogeneity, in regions characterized by ephemeral FT events (i.e., regions <40 °N), as well as during freeze and thaw onset periods. This study also provides insight to improving the representation of surface FT state by providing a clearer definition of the subpixel scale temperature characteristics that govern existing frozen classifications.

Earth observing systems↗

Comparison of Local Scale Measured and Modeled Brightness Temperatures and Snow Parameters from the CLPX 2003 by Means of a Dense Medium Radiative Transfer Theory Model

Microwave remote sensing offers distinct advantages for observing the cryosphere. Solar illumination is not required, and spatial and temporal coverage are excellent from polar-orbiting satellites. Passive microwave measurements are sensitive to the two most useful physical quantities for many hydrological applications: physical temperature and water content/state. Sensitivity to the latter is a direct result of the microwave sensitivity to the dielectric properties of natural media, including snow, ice, soil (frozen or thawed), and vegetation. These considerations are factors motivating the development of future cryospheric satellite remote sensing missions, continuing and improving on a 26-year microwave measurement legacy. Perhaps the biggest issues regarding the use of such satellite measurements involve how to relate parameter values at spatial scales as small as a hectare to observations with sensor footprints that may be up to 25 x 25 km. The NASA Cold-land Processes Field Experiment (CLPX) generated a dataset designed to enhance understanding of such scaling issues. CLPX observations were made in February (dry snow) and March (wet snow), 2003 in Colorado, USA, at scales ranging from plot scale to 25 x 25 km satellite footprints. Of interest here are passive microwave observations from ground-based, airborne, and satellite sensors, as well as meteorological and snowpack measurements that will enable studies of the effects of spatial heterogeneity of surface conditions on the observations. Prior to performing such scaling studies, an evaluation of snowpack forward modelling at the plot scale (least heterogeneous scale) is in order. This is the focus of this paper. Many forward models of snow signatures (brightness temperatures) have been developed over the years. It is now recognized that a dense medium radiative transfer (DMRT) treatment represents a high degree of physical fidelity for snow modeling, yet dense medium models are particularly sensitive to snowpack structural parameters such as grain size, density, and depth---parameters that may vary substantially within a snowpack. Microwave radiometric data and snow pit measurements collected at the Local-Scale Observation Site (LSOS) during the third Intensive Observation Period (IOP3) of the CLPX have been used to test the capabilities of a DMRT model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP). The radiometric measurements were made by the University of Tokyo s Ground Based Microwave Radiometer-7 (GBMR-7) system. We evaluate the degree to which a DMRT-based model can accurately reproduce the GBMR-7 brightness temperatures at different frequencies and incidence angles.

Tedescol, Marco↗

Observed Scaling Behavior of Snow and Soil Moisture Signatures from the CLPX and NAFE Campaigns

Seasonal snow and soil moisture (SM) are important elements of the Earth's water cycle, acting as storage elements and by modulating water and energy exchanges between the land surface and the atmosphere and space. Accurate estimates of snow water equivalent (SWE) and SM also play an important role in weather and climate models over a range of scales. There are many advantages for the retrieval of snow parameters using passive microwave sensors. However, the current scale mismatch between many high-resolution applications and available satellite snow sensor footprints (approx.25km) presents challenges. There is also a scale mismatch for planned L-band soil moisture satellite measurements with their large (50-100 km) passive microwave footprints. Combined radiometer-radar schemes are being explored as possible aids, but this is still an area of ongoing research. Retrieval accuracies are limited in part by complex, and potentially nonlinear relationships between SWE or SM and microwave emission signatures. We will use airborne and satellite observations from two field campaigns to examine the scaling behavior of brightness temperatures and SWE and SM across scales from 100 meters to 10's of kilometers. The results will provide examples of the appropriate spatial sampling scales for new sensors, and guidance for downscaling schemes. The two campaigns are NASA's Cold Land Processes Experiment- 1 (CLPX- 1 : Colorado, 2002 & 2003), and the Australian National Airborne Field Experiment (NAFE, Nov, 2005). The analyses will also illustrate the roles and spatial scales of the underlying phenomena (e.g., vegetation, topography) that control subpixel heterogeneity.

Kim, Edward J.↗

Modeling and Prediction of Wildfire Hazard in Southern California, Integration of Models with Imaging Spectrometry

Large urban wildfires throughout southern California have caused billions of dollars of damage and significant loss of life over the last few decades. Rapid urban growth along the wildland interface, high fuel loads and a potential increase in the frequency of large fires due to climatic change suggest that the problem will worsen in the future. Improved fire spread prediction and reduced uncertainty in assessing fire hazard would be significant, both economically and socially. Current problems in the modeling of fire spread include the role of plant community differences, spatial heterogeneity in fuels and spatio-temporal changes in fuels. In this research, we evaluated the potential of Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and Airborne Synthetic Aperture Radar (AIRSAR) data for providing improved maps of wildfire fuel properties. Analysis concentrated in two areas of Southern California, the Santa Monica Mountains and Santa Barbara Front Range. Wildfire fuel information can be divided into four basic categories: fuel type, fuel load (live green and woody biomass), fuel moisture and fuel condition (live vs senesced fuels). To map fuel type, AVIRIS data were used to map vegetation species using Multiple Endmember Spectral Mixture Analysis (MESMA) and Binary Decision Trees. Green live biomass and canopy moisture were mapped using AVIRIS through analysis of the 980 nm liquid water absorption feature and compared to alternate measures of moisture and field measurements. Woody biomass was mapped using L and P band cross polarimetric data acquired in 1998 and 1999. Fuel condition was mapped using spectral mixture analysis to map green vegetation (green leaves), nonphotosynthetic vegetation (NPV; stems, wood and litter), shade and soil. Summaries describing the potential of hyperspectral and SAR data for fuel mapping are provided by Roberts et al. and Dennison et al. To utilize remotely sensed data to assess fire hazard, fuel-type maps were translated into standard fuel models accessible to the FARSITE fire spread simulator. The FARSITE model and BEHAVE are considered industry standards for fire behavior analysis. Anderson level fuels map, generated using a binary decision tree classifier are available for multiple dates in the Santa Monica Mountains and at least one date for Santa Barbara. Fuel maps that will fill in the areas between Santa Barbara and the Santa Monica Mountains study sites are in progress, as part of a NASA Regional Earth Science Application Center, the Southern California Wildfire Hazard Center. Species-level maps, were supplied to fire managing agencies (Los Angeles County Fire, California Department of Forestry). Research results were published extensively in the refereed and non-refereed literature. Educational outreach included funding of several graduate students, undergraduate intern training and an article featured in the California Alliance for Minorities Program (CAMP) Quarterly Journal.

Roberts, Dar A.↗

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong↗

Validation of Land-Surface Mosaic Heterogeneity in the GEOS DAS

The Mosaic Land-surface Model (LSM) has been included into the current GEOS Data Assimilation System (DAS). The LSM uses a more advanced representation of physical processes than previous versions of the GEOS DAS, including the representation of sub-grid heterogeneity of the land-surface through the Mosaic approach. As a first approximation, Mosaic assumes that all similar surface types within a grid-cell can be lumped together as a single'tile'. Within one GCM grid-cell, there might be 1 - 5 different tiles or surface types. All tiles are subjected to the grid-scale forcing (radiation, air temperature and specific humidity, and precipitation), and the sub-grid variability is a function of the tile characteristics. In this paper, we validate the LSM sub-grid scale variability (tiles) using a variety of surface observing stations from the Southern Great Plains (SGP) site of the Atmospheric Radiation Measurement (ARM) Program. One of the primary goals of SGP ARM is to study the variability of atmospheric radiation within a G,CM grid-cell. Enough surface data has been collected by ARM to extend this goal to sub-grid variability of the land-surface energy and water budgets. The time period of this study is the Summer of 1998 (June I - September 1). The ARM site data consists of surface meteorology, energy flux (eddy correlation and bowen ratio), soil water observations spread over an area similar to the size of a G-CM grid-cell. Various ARM stations are described as wheat and alfalfa crops, pasture and range land. The LSM tiles considered at the grid-space (2 x 2.5) nearest the ARM site include, grassland, deciduous forests, bare soil and dwarf trees. Surface energy and water balances for each tile type are compared with observations. Furthermore, we will discuss the land-surface sub-grid variability of both the ARM observations and the DAS.

Bosilovich, Michael G.↗

Real Time Land-Surface Hydrologic Modeling Over Continental US

The land surface component of the hydrological cycle is fundamental to the overall functioning of the atmospheric and climate processes. Spatially and temporally variable rainfall and available energy, combined with land surface heterogeneity cause complex variations in all processes related to surface hydrology. The characterization of the spatial and temporal variability of water and energy cycles are critical to improve our understanding of land surface-atmosphere interaction and the impact of land surface processes on climate extremes. Because the accurate knowledge of these processes and their variability is important for climate predictions, most Numerical Weather Prediction (NWP) centers have incorporated land surface schemes in their models. However, errors in the NWP forcing accumulate in the surface and energy stores, leading to incorrect surface water and energy partitioning and related processes. This has motivated the NWP to impose ad hoc corrections to the land surface states to prevent this drift. A proposed methodology is to develop Land Data Assimilation schemes (LDAS), which are uncoupled models forced with observations, and not affected by NWP forcing biases. The proposed research is being implemented as a real time operation using an existing Surface Vegetation Atmosphere Transfer Scheme (SVATS) model at a 40 km degree resolution across the United States to evaluate these critical science questions. The model will be forced with real time output from numerical prediction models, satellite data, and radar precipitation measurements. Model parameters will be derived from the existing GIS vegetation and soil coverages. The model results will be aggregated to various scales to assess water and energy balances and these will be validated with various in-situ observations.

Houser, Paul R.↗

Surface Coatings on Lunar Volcanic Glasses

We are undertaking a detailed study of surface deposits on lunar volcanic glass beads. These tiny deposits formed by vapor condensation during cooling of the gases that drove the fire fountain eruptions responsible for the formation of the beads. Volcanic glass beads are present in most lunar soil samples in the returned lunar collection. The mare-composition beads formed as a result of fire-fountaining approx.3.4-3.7 Ga ago, within the age range of large-scale mare volcanism. Some samples from the Apollo 15 and Apollo 17 landing sites are enriched in volcanic spherules. Three major types of volcanic glass bead have been identified: Apollo 15 green glass, Apollo 17 orange glass, and Apollo 17 "black" glass. The Apollo 15 green glass has a primitive composition with low Ti. The high-Ti compositions of the orange and black glasses are essentially identical to each other but the black glasses are opaque because of quench crystallization. A poorly understood feature common to the Apollo 15 and 17 volcanic glasses is the presence of small deposits of unusual materials on their exterior surfaces. For example, early studies indicated that the Apollo 17 orange glasses had surface enrichments of In, Cd, Zn, Ga, Ge, Au, and Na, and possible Pb- and Zn-sulfides, but it was not possible to characterize the surface features in detail. Technological advances now permit us to examine such features in detail. Preliminary FE-TEM/X-ray studies of ultramicrotome sections of Apollo 15 green glass indicate that the surface deposits are heterogeneous and layered, with an inner layer consisting of Fe with minor S and an outer layer of Fe and no S, and scattered Zn enrichments. Layering in surface deposits has not been identified previously; it will be key to defining the history of lunar fire fountaining.

Wentworth, Susan J.↗

Microchemical and Structural Evidence for Space Weathering in Soils from Asteroid Itokawa

The chemistry, microstructure and optical properties of grains on the surfaces of airless bodies are continu-ously modified due to their interactions predominantly with solar energetic ions and micrometeorite impacts. Collectively known as space weathering, this phenomenon results in a discrepancy between remotely sensed spectra from asteroids and those ac-quired directly from meteorites. The return of pristine samples from the asteroid Itokawa provides insight into surface processes on airless bodies and will help in correlating remote sensing data with laboratory analysis of meteorites. Samples and Methods: We examined Itokawa samples RA-QD02-0042-01 and RA-QD-02-0042-02, ultramicrotomed sec-tions of a singular grain prepared by the Hayabusa sample cura-tion team. We analyzed these slices using a 200 keV JEOL 2010F transmission electron microscope (TEM) at Arizona State Uni-versity and a 200 keV JEOL 2500SE TEM at NASA JSC. Both field emission TEMs are equipped with energy-dispersive X-ray spectrometers (EDS) and scanning TEM (STEM) detectors. Results and Discussion: TEM observations reveal that the sectioned grain predominantly consists of a single crystal of low-Ca orthopyroxene, with subsidiary smaller regions of olivine, Fe-Ni sulfide, and Fe-Ni metal. EDS-spectrum imaging and high-resolution TEM (HRTEM) show local, nanocrystalline regions of the outermost 2 to 5 nm of the pyroxene are composed of an Fe-Mg-S-rich and Si- and O-depleted layer that is underlain by a 2- to 5-nm thick amorphous zone enriched in Si. These layers occur in multiple microtome slices and have uniform thicknesses. We also observe localized 'islands' of material on the surface of the pyroxene which HRTEM imaging indicates are amorphous and EDS measurements show are compositionally heterogeneous. A 10- to 60-nm thick partially amorphous zone occurs below the compositionally distinct rim. While this this zone is associated with the compositionally heterogeneous outer layer, it also occurs as a local stand-alone feature on the exterior rim of the grain. Ar-eas of the pyroxene grain rim also exhibit a vesicular texture. The TEM data indicate a complex history of space weather-ing for samples RA-QD02-0042-01 and -02. The outermost layer of nanocrystalline material with varied composition is consistent with previously suggested [3-4] chemical and structural pro-cessing by solar wind ions, with a possible additional role for im-pact vapor deposition [3-4]. The amorphous and compositionally distinct islands on the surface of this grain, similar to lunar glasses, suggest formation through vapor deposition via micrometeor-ite impact events. In comparison, the amorphization and vesicula-tion textures are likely a product of radiation damage from the solar wind. The depth and degree of amorphization, in conjunction with model calculations, will help provide an upper limit on exposure time for these particles.

Thompson, M. S.↗

Patterns of Canopy and Surface Layer Consumption in a Boreal Forest Fire from Repeat Airborne Lidar

Fire in the boreal region is the dominant agent of forest disturbance with direct impacts on ecosystem structure, carbon cycling, and global climate. Global and biome-scale impacts are mediated by burn severity, measured as loss of forest canopy and consumption of the soil organic layer. To date, knowledge of the spatial variability in burn severity has been limited by sparse field sampling and moderate resolution satellite data. Here, we used pre- and post-fire airborne lidar data to directly estimate changes in canopy vertical structure and surface elevation for a 2005 boreal forest fire on Alaskas Kenai Peninsula. We found that both canopy and surface losses were strongly linked to pre-fire species composition and exhibited important fine-scale spatial variability at sub-30m resolution. The fractional reduction in canopy volume ranged from 0.61 in lowland black spruce stands to 0.27 in mixed white spruce and broad leaf forest. Residual structure largely reflects standing dead trees, highlighting the influence of pre-fire forest structure on delayed carbon losses from above ground biomass, post-fire albedo, and variability in understory light environments. Median loss of surface elevation was highest in lowland black spruce stands (0.18 m) but much lower in mixed stands (0.02 m), consistent with differences in pre-fire organic layer accumulation. Spatially continuous depth-of-burn estimates from repeat lidar measurements provide novel information to constrain carbon emissions from the surface organic layer and may inform related research on post-fire successional trajectories. Spectral measures of burn severity from Landsat were correlated with canopy (r = 0.76) and surface (r = -0.71) removal in black spruce stands but captured less of the spatial variability in fire effects for mixed stands (canopy r = 0.56, surface r = -0.26), underscoring the difficulty in capturing fire effects in heterogeneous boreal forest landscapes using proxy measures of burn severity from Landsat.

Alonzo, Michael↗

Contributions of GRACE to Climate Monitoring

The NASA/German Gravity Recovery and Climate Experiment (GRACE) was launched in March 2002. Rather than looking downward, GRACE continuously monitors the locations of and precise distance between twin satellites which orbit in tandem about 200 km apart. Variations in mass near Earth's surface cause heterogeneities in its gravity field, which in turn affect the orbits of satellites. Thus scientists can use GRACE data to map Earth's gravity field with enough accuracy to discern month to month changes caused by ocean circulation and redistribution of water stored on and in the land. Other gravitational influences, such as atmospheric circulation, post-glacial rebound, and solid earth movements are either independently determined and removed or are negligible on a monthly to sub-decadal timescale. Despite its coarse spatial (>150,000 sq km at mid-latitudes) and temporal (approx monthly) resolutions, GRACE has enabled significant advancements in the oceanic, hydrologic, and cryospheric science, and has great potential for climate monitoring, because it is the only global observing system able to measure ocean bottom pressures, total terrestrial water storage, and ice mass changes. The best known GRACE results are estimates of Greenland and Antarctic ice sheet loss rates. Previously, scientists had estimated ice mass losses using ground and satellite based altimetry and surface mass balance estimates based on snowfall accumulation and glacier discharge. While such measurements are still very useful for their spatial detail, they are imperfectly correlated with large-scale ice mass changes, due to snow and ice compaction and incomplete spatial coverage. GRACE enables scientists to generate monthly time series of Greenland and Antarctic ice mass, which have confirmed the shrinking of the polar ice sheets, one of the most obvious and indisputable manifestations of climate change. Further, GRACE has located and quantified hot spots of ice loss in southeastern Greenland and western Antarctica. For 2002 to present, the rate of ice mass loss has been 200 to 300 GT/yr in Greenland and 70 to 210 GT/yr in Antarctica, and some scientists are suggesting that the rates are accelerating. Similarly, GRACE has been used to monitor mass changes in alpine glaciers. Tamisiea et al. first characterized glacier melt along the southern coast of Alaska, more recently estimated to be occurring at a rate of 84 GT/yr. Chen et al. estimated that Patagonian glaciers are melting at a rate of 28 GT/yr, and estimated that the high mountains of central Asia lose ice at a rate of 47 GT/yr. Tapley et al. and Wahr et al. presented the first GRACE based estimates of changes in column-integrated terrestrial water storage (TWS; the sum of ground-water, soil moisture, surface waters, snow, ice, and water stored in vegetation) at continental scales. Since then, dozens of studies have shown that GRACE based estimates of regional to continental scale TWS variations agree with independent information, and some innovative uses of GRACE data have been developed. Rodell et al. (2004) and Swenson and Wahr (2006) demonstrated that by combining GRACE derived terrestrial water storage changes with observations of precipitation and runoff in a river basin scale water budget, it was possible to produce new estimates of evapotranspiration and atmospheric moisture convergence, essential climate variables that are difficult to estimate accurately. Similarly, GRACE has been used to constrain estimates of global river discharge and the contribution of changes in TWS to sea level rise. Crowley et al. observed a negative correlation between interannual TWS anomalies in the Amazon and the Congo River basin. Yeh et al. and Rodell et al. estimated regionally averaged groundwater storage variations based on GRACE and auxiliary observations. Rodell et al. and Tiwari et al. applied that method to quantify massive groundwater depletion in northern India caused by over reliance on aquifers for irration, and Famiglietti et al. found a similar situation in California's Central Valley. Zaitchik et al. and Lo et al. described approaches to use GRACE to constrain hydrological models, enabling integration of GRACE data with other observations and achieving much higher spatial and temporal resolutions than GRACE alone. Such approaches are now supporting applications including drought and water resources monitoring. Oceanography has likewise benefitted from the independent nature of GRACE observations. One application is measurement of the mass component of sea level rise, which complements radar altimetry and in situ measurements. GRACE also measures ocean bottom pressures (OBP), which help to refine understanding and modeling of ocean circulation and the ocean's fresh water budget, among other things. For example, Hayakawa et al. showed that GRACE observes OBP patterns absent from the background models of oceanic variability. Morison et al. used GRACE to describe important decadal scale shifts in circulation and an ongoing trend of freshening of the western Arctic, important indicators of climate variability. The research of Song and Zlotnicki and Chambers and Willis on GRACE-derived ocean bottom pressures in the sub-polar gyre led to the discovery of an ENSO teleconnection and a long-term change in OBP in the North Pacific sub-polar gyre that was not predicted by an ocean model. Further, Chambers and Willis were able to identify an internal redistribution of mass between Atlantic and Pacific Oceans lasting at least six years, which was not predicted by ocean models and was the first direct evidence of sustained mass transport from one ocean basin to another on periods longer than a year. Boening et al. observed a record increase in OBP over part of the southeastern Pacific in late 2009 and early 2010, primarily caused by wind stress curl associated with a strong and persistent anticyclone and likely related to the concurrent Central Pacific El Nino. GRACE has far surpassed its 5-year design lifetime, but it will likely succumb to the aging of batteries and instrument systems sometime in the next few years. NASA has begun initial development of a follow-on to GRACE with very similar design, which could launch as soon as 2016 and would provide continuity in the data record while improving resolution slightly. Higher resolution time variable gravity missions are also on the drawing board .

Rodell, Matthew↗