Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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 1,243 records · Page 69

State Estimation for Distribution Networks with Asynchronous Sensors Using Stochastic Descent: Preprint

This paper investigates the problem of state estimation for distribution networks with asynchronous sensors comprising of a mix of smart meters and phasor measurement units (PMUs) with multiple sampling and reporting rates. We consider two independent scenarios of state estimation and tracking, with either voltages or currents as states. With these two sets, we investigate estimation under (a) full data, assuming all measurements are available and (b) limited data, where an online algorithmic approach is adopted to estimate the possibly time-varying states by processing measurements as and when available. The proposed algorithm, inspired by the classical Stochastic Gradient Descent (SGD) approach updates the states based on the previous estimate and the newly available measurements. Finally, we demonstrate the estimation and tracking efficacy through numerical simulations on the IEEE-37 test network, while also highlighting how estimation with currents as states leads to faster convergence.

asynchronous sensors↗

Enhanced orbital gyrocompassing by the optical flow sensed by an Earth-pointing camera

A new method for improving the orbital gyrocompassing process involving the attitude angle estimation of an earth-pointing satellite in low-Earth orbit uses an electro-optical sensor for direct measurement of the satellite azimuth angle. Simulations have shown that this additional measurement drastically reduces the estimator convergence time, especially when the sun sensor is rendered ineffective, e.g., by high solar elevations. The azimuth-sensing method is based on estimation of the image shift between successive picture frames of an on-board, Earth-pointing, charge-coupled device (CCD) full-matrix camera. The shift-estimation algorithm is based on minimizing a cost function which expresses mean-squared differences in brightness patterns of selected areas of the two frames. An extensive evaluation program with a computer-controlled 2-axis light table and actual satellite images has demonstrated high robustness for a wide range of variation of parameters including image texture content; camera focal length; sampling rate; and number of pixels processed. It was shown to be possible to estimate the azimuth angle within 0.1-0.2 degrees, for a suitably chosen parameter set.

Topaz, Leora↗

Analysis of Physical Properties of Dust Suspended in the Mars Atmosphere

Methods for iteratively determining the infrared optical constants for dust suspended in the Mars atmosphere are described. High quality spectra for wavenumbers from 200 to 2000 1/cm were obtained over a wide range of view angles by the Mariner 9 spacecraft, when it observed a global Martian dust storm in 1971-2. In this research, theoretical spectra of the emergent intensity from Martian dust clouds are generated using a 2-stream source-function radiative transfer code. The code computes the radiation field in a plane-parallel, vertically homogeneous, multiply scattering atmosphere. Calculated intensity spectra are compared with the actual spacecraft data to iteratively retrieve the optical properties and opacity of the dust, as well as the surface temperature of Mars at the time and location of each measurement. Many different particle size distributions a-re investigated to determine the best fit to the data. The particles are assumed spherical and the temperature profile was obtained from the CO2 band shape. Given a reasonable initial guess for the indices of refraction, the searches converge in a well-behaved fashion, producing a fit with error of less than 1.2 K (rms) to the observed brightness spectra. The particle size distribution corresponding to the best fit was a lognormal distribution with a mean particle radius, r(sub m) 0.66 pm, and variance, omega(sup 2) = 0.412 (r(sub eff) = 1.85 microns, v(sub eff) =.51), in close agreement with the size distribution found to be the best fit in the visible wavelengths in recent studies. The optical properties and the associated single scattering properties are shown to be a significant improvement over those used in existing models by demonstrating the effects of the new properties both on heating rates of the Mars atmosphere and in example spectral retrieval of surface characteristics from emission spectra.

Snook, Kelly↗

Stochastic Mixing Model with Power Law Decay of Variance

Here we present a simple stochastic mixing model based on the law of large numbers (LLN). The reason why the LLN is involved in our formulation of the mixing problem is that the random conserved scalar c = c(t,x(t)) appears to behave as a sample mean. It converges to the mean value mu, while the variance sigma(sup 2)(sub c) (t) decays approximately as t(exp -1). Since the variance of the scalar decays faster than a sample mean (typically is greater than unity), we will introduce some non-linear modifications into the corresponding pdf-equation. The main idea is to develop a robust model which is independent from restrictive assumptions about the shape of the pdf. The remainder of this paper is organized as follows. In Section 2 we derive the integral equation from a stochastic difference equation describing the evolution of the pdf of a passive scalar in time. The stochastic difference equation introduces an exchange rate gamma(sub n) which we model in a first step as a deterministic function. In a second step, we generalize gamma(sub n) as a stochastic variable taking fluctuations in the inhomogeneous environment into account. In Section 3 we solve the non-linear integral equation numerically and analyze the influence of the different parameters on the decay rate. The paper finishes with a conclusion.

Fedotov, S.↗

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Machine learning informed rational design of high entropy double perovskite oxide universal air/steam electrodes for solid oxide electrochemical cells

Due to their high efficiency and versatility, solid oxide electrochemical cells (SOCs) are poised to play a significant role in future energy conversion and storage applications. In recent years, SOCs have bifurcated into two distinct categories: traditional oxygen-ion conducting SOCs that typically operate from ∼650—850 °C and the more recent proton-conducting ceramic (PCC) SOCs that typically operate from ∼400—650 °C. Current performance and lifetime of both oxygen-ion conducting SOCs and PCCs is primarily limited by the air/steam electrode, which facilitates the oxygen reduction reaction (ORR) during fuel cell operation and must also facilitate the oxygen evolution reaction (OER) during electrolysis operation. Here, we present a newly designed high-entropy double perovskite oxide suitable as a universal ORR/OER electrode for both oxygen-ion conducting SOCs and PCCs. Machine learning methods are applied to identify chemical descriptors for highly catalytic high-entropy double perovskite oxides (AA’B 2 O 6 ) across a large compositional space. Based on the machine-learning guidance, we ultimately converge on Ba 0.9 Cs 0.1 (Ca 0.2 Gd 0.2 La 0.2 Pr 0.2 Sr 0.2 )Co 1.5 Fe 0.5 O 6 (CsBaHEO) as a universal air/steam electrode. Structure stabilization is accomplished by an equimolar five-cation high-entropy composition on the A’-site, while cesium substitution on the A-site enhances the electrical conductivity and leads to a higher oxygen vacancy concentration. This material exhibits versatility and high performance in reversible oxygen-ion SOCs, reversible PCCs, and also large-scale tubular PCCs. For example, the CsBaHEO-based PCC reaches 1018 mW∙cm −2 at 600°C, while a large-scale tubular PCC using CsBaHEO for electrolysis achieves a hydrogen production rate of 21.314 ML∙min −1 at 600 °C.

Cell↗

What Supergranule Flow Models Tell Us About the Sun's Surface Shear Layer and Magnetic Flux Transport

Models of the photospheric flows due to supergranulation are generated using an evolving spectrum of vector spherical harmonics up to spherical harmonic wavenumber l~1500. Doppler velocity data generated from these models are compared to direct Doppler observations from SOHO/MDI and SDO/HMI. The models are adjusted to match the observed spatial power spectrum as well as the wavenumber dependence of the cell lifetimes, differential rotation velocities, meridional flow velocities, and relative strength of radial vs. horizontal flows. The equatorial rotation rate as a function of wavelength matches the rotation rate as a function of depth as determined by global helioseismology. This leads to the conclusions that the cellular structures are anchored at depths equal to their widths, that the surface shear layer extends to at least 70 degrees latitude, and that the poleward meridional flow decreases in amplitude and reverses direction at the base of the surface shear layer (approx.35 Mm below the surface). Using the modeled flows to passively transport magnetic flux indicates that the observed differential rotation and meridional flow of the magnetic elements are directly related to the differential rotation and meridional flow of the convective pattern itself. The magnetic elements are transported by the evolving boundaries of the supergranule pattern (where the convective flows converge) and are unaffected by the weaker flows associated with the differential rotation or meridional flow of the photospheric plasma.

Hathaway, David↗

Web Monitoring of EOS Front-End Ground Operations, Science Downlinks and Level 0 Processing

This paper addresses the efforts undertaken and the technology deployed to aggregate and distribute the metadata characterizing the real-time operations associated with NASA Earth Observing Systems (EOS) high-rate front-end systems and the science data collected at multiple ground stations and forwarded to the Goddard Space Flight Center for level 0 processing. Station operators, mission project management personnel, spacecraft flight operations personnel and data end-users for various EOS missions can retrieve the information at any time from any location having access to the internet. The users are distributed and the EOS systems are distributed but the centralized metadata accessed via an external web server provide an effective global and detailed view of the enterprise-wide events as they are happening. The data-driven architecture and the implementation of applied middleware technology, open source database, open source monitoring tools, and external web server converge nicely to fulfill the various needs of the enterprise. The timeliness and content of the information provided are key to making timely and correct decisions which reduce project risk and enhance overall customer satisfaction. The authors discuss security measures employed to limit access of data to authorized users only.

Cordier, Guy R.↗

The Impact of Radiative Transfer at Reduced Spectral Resolution in Large‐Eddy Simulations of Convective Clouds

Abstract Many radiative transfer schemes approximate the spectral integration over ∼10 5 to ∼10 6 wavelengths with correlated k ‐distributions methods that typically require only 10 1 –10 2 spectral integration points ( g ‐points). The exact number of g ‐points is then chosen as an optimal balance between computational costs and accuracy, normally assessed in terms of a number of radiative quantities. How this radiative accuracy propagates to simulation accuracy, however, is not straightforward. In this study, we therefore explore the sensitivity of cloud properties in large‐eddy simulations (LES) to the accuracy of radiative fluxes and heating rates. We first generate smaller sets of g ‐points from existing k ‐distributions by repeatedly combining adjacent g ‐points while maintaining the highest possible accuracy on a chosen set of radiative metrics. Next, we perform three sets of LES with varying cloud—radiation coupling pathways, and therefore different requirements for the accuracy of the radiative transfer computations, to investigate how these smaller and thus less accurate k ‐distributions affect simulation characteristics. The decrease in radiative accuracy with 3–4 times smaller k ‐distributions results in biases in cloud properties that are relative small compared to their temporal fluctuations. These results show potential for speeding up radiative transfer computations in cloud‐resolving models by reducing the resolved spectral detail. However, more statistically converged simulations and a wider set of case studies is required to fully assess the robustness of our results.

Meteorology & Atmospheric Sciences↗

Baroclinic instability of flows along sloping boundaries

In the two-layer quasi-geostrophic model with boundaries sloping perpendicular to the basic flow, the ratios of the slopes of the bottom and the top to that of the interface between the fluid layers in the basic state are important parameters in the expression of the growth rate of unstable waves. When Eady's (1949) model is extended to include sloping bottom and top boundaries, the growth rates of unstable waves depend on the ratios of the slopes of the bottom and the top to that of the isentropes of the basic state. For the Eady model with sloping bottom, an important parameter characterizing the instability is the ratio between the vertical and horizontal heat transports by the wave divided by the slope of the isentropes of the basic state. An interpretation of these ratios and their relations clarifies the stabilization of the system for large slopes, the variation of the wavelength of the most unstable wave with the bottom slope, and the destabilization of some short waves for negative bottom slopes. It is found that the most unstable wave of the system has zero vertical energy flux convergence at the sloping bottom.

Mechoso, C. R.↗

Magnetized Target Fusion Driven by Plasma Liners

Magnetized target fusion (MTF) attempts to combine the favorable attributes of magnetic confinement fusion (MCF) for energy confinement with the attributes of inertial confinement fusion (ICF) for efficient compression heating and wall-free containment of the fusing plasma. It uses a material liner to compress and contain a magnetized plasma. For practical applications, standoff drivers to deliver the imploding momentum flux to the target plasma remotely are required. Spherically converging plasma jets have been proposed as standoff drivers for this purpose. The concept involves the dynamic formation of a spherical plasma liner by the merging of plasma jets, and the use of the liner so formed to compress a spheromak or a field reversed configuration (FRC). For the successful implementation of the scheme, plasma jets of the requisite momentum flux density need to be produced. Their transport over sufficiently large distances (a few meters) needs to be assured. When they collide and merge into a liner, relative differences in velocity, density and temperature of the jets could give rise to instabilities in the development of the liner. Variation in the jet properties must be controlled to ensure that the growth rate of the instabilities are not significant over the time scale of the liner formation before engaging with the target plasma. On impact with the target plasma, some plasma interpenetration might occur between the liner and the target. The operating parameter space needs to be identified to ensure that a reasonably robust and conducting contact surface is formed between the liner and the target. A mismatch in the "impedance" between the liner and the target plasma could give rise to undesirable shock heating of the liner leading to increased entropy (thermal losses) in the liner. Any irregularities in the liner will accentuate the Rayleigh-Taylor instabilities during the compression of the target plasma by the liner.

Thio, Y. C. Francis↗

Multifaceted aerosol effects on precipitation

Aerosols have been proposed to influence precipitation rates and spatial patterns from scales of individual clouds to the globe. However, large uncertainty remains regarding the underlying mechanisms and importance of multiple effects across spatial and temporal scales. Here, in this study, we review the evidence and scientific consensus behind these effects, categorized into radiative effects via modification of radiative fluxes and the energy balance, and microphysical effects via modification of cloud droplets and ice crystals. Broad consensus and strong theoretical evidence exist that aerosol radiative effects (aerosol–radiation interactions and aerosol–cloud interactions) act as drivers of precipitation changes because global mean precipitation is constrained by energetics and surface evaporation. Likewise, aerosol radiative effects cause well-documented shifts of large-scale precipitation patterns, such as the intertropical convergence zone. The extent of aerosol effects on precipitation at smaller scales is less clear. Although there is broad consensus and strong evidence that aerosol perturbations microphysically increase cloud droplet numbers and decrease droplet sizes, thereby slowing precipitation droplet formation, the overall aerosol effect on precipitation across scales remains highly uncertain. Global cloud-resolving models provide opportunities to investigate mechanisms that are currently not well represented in global climate models and to robustly connect local effects with larger scales. This will increase our confidence in predicted impacts of climate change.

54 ENVIRONMENTAL SCIENCES↗

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

The Oak Ridge National Laboratory (ORNL) Second Target Station (STS) neutron production facility is an accelerator driven pulsed neutron source that is currently being actively developed at ORNL. The neutrons are produced by proton-induced spallation reactions. A proton beam of 700 kW power is delivered to a spallation target in short, less than 1 µs long pulses, with 15 Hz repetition rate. The spallation target of ORNL STS is a rotating water-cooled tungsten target with tantalum cladding housed in a stainless-steel shroud. It is divided into 21 segments. These segments become highly activated due to spallation reactions or nuclei transmutation by the emitted neutrons. The radioactive nuclides continue to decay after ceasing operation. The decay dose rates generated from the target segments once they are removed from their operational location within the core vessel must be accurately quantified to determine the shielding configurations of remote handling tools and transport casks and to aid in planning maintenance events. To determine the shielding configurations needed for an activated target segment after ceasing operation, both the hybrid unstructured mesh (UM)/constructive solid geometry (CSG) approach that was previously utilized for STS analyses [1] and the ADVANTG code [2] were used. Even though the ADVANTG code does not include UM capability, the utilization of its advanced variance reduction technique was crucial to accelerate the extremely difficult final photon transport calculation in this analysis. This paper also describes the procedures taken to mitigate the convergence issues that occur when ADVANTG uses a source definition that does not match the source of the final Monte Carlo (MC) calculation. These convergence issues often occur because of inconsistencies between source and transport biasing parameters.

Ibrahim, Ahmad↗

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↗

Radiative heating rates near the stratospheric fountain

Radiative heating rates are computed for various sets of conditions thought to be appropriate to the stratospheric fountain region: with and without a layer of cirrus cloud between 100 and 150 mbar; with standard ozone and with decreased ozone in the lower stratosphere, again with and without the cirrus cloud; and with different temperatures in the tropopause region. The presence of the cloud decreases the radiative cooling below the cloud in the upper troposphere and increases the cooling above it in the lower stratosphere. The cloud is heated at the base and cooled at the top and thus radiatively destabilized; overall it gains energy by radiation. Decreasing ozone above the cloud also tends to cool the lower stratosphere. The net effect is a tendency for vertical convergence and horizontal divergence in the cloud region. High resolution profiles of temperature, ozone, and cloudiness within the fountain region are required in order to assess the final balance of the various processes.

Doherty, G. M.↗