Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean 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.

At least 973 records · Page 54

Four-Dimensional Oceanic and Atmosperic Data Assimilation with Tropical Rainfall Measuring Mission Data

An oceanic data assimilation system which allows to utilize the forthcoming Tropical Rainfall Measuring Mission (TRMM) data has been developed and applied to the Pacific Ocean to produce the velocity field. The assimilated data will be indispensable to examine the effects of rainfall and its variability on the structure and circulation of the tropical oceans and to assess the impact of global warming due to the increase of carbon dioxide on the ocean circulation system and the marine pollution caused by oil spill and ocean damping of radionuclide. The data will also provide the verification for the oceanic and ocean-atmosphere coupled General Circulation Models (GCM's). The system consists of oceanic GCM, analysis scheme and data. In the system the flow field has been determined to be physically consistent with the observed density field and the sea surface winds derived from the Special Sensor Microwave Imagery (SSM/I) data which drive the ocean current. The time integration has been performed for five years until the flow field near the surface attained the steady state starting from the rest ocean with observed temperature and salinity fields, and the SSM/I surface wind velocity. The resultant flow field showed high producibility of the system. Especially the flow near the ocean surface agreed well with available observed data. The system, for the first time, succeeded to produce the eastward subtropical current which has been discovered in the joint investigation on Kuroshio current (CSK) in the 1960s. To verify the quality of the flow field a trajectory analysis has been carried out and compared with the Algos buoy data. BRIEF DESCRIPTION OF THE DATA ASSIMILATION SYSTEM ## Oceanic GCM and analysis scheme--The basic equations are much the same as used for the GCM's, except for the Newtonian damping terms introduced into the prediction equations for the potential temperature and salinity to maintain these fields as observed. The C grid of 2'lat. by 2'long. in horizontal and the 11 vertical levels are applied to the entire Pacific Ocean. At the east and west ocean boundaries the periodic boundary conditions are applied creating fictitious ocean there. The SMAC Method is used to increase the accuracy of mass conservation. * Data--The JODC temperature and salinity data obtained from 1906 to 1988 are used in the system between Long.100'E. and 60'W. The surface wind data are derived from the SSM/I data by Dr-R. Atlas of NASA/GSFC. The data set contains every 6 hours data from July 1987 to June 1989 on the grid of 2'lat. by 2.5'long. The averaged for the whole period and then interpolated into the 2'lat. by 2'long. grid data are used to force the system. The sea bottom topography data was based on the General Bathymetric Chart of the Ocean (GEBCO) supplied by the Canadian Hydrographic Service under contract with the International Hydrographic Organization and International Oceanographic Commission of UNESCO.

Takano, Kenji↗

Online Simulations of Global Aerosol Distributions in the NASA GEOS-4 Model and Comparisons to Satellite and Ground-Based Aerosol Optical Depth

We have implemented a module for tropospheric aerosols (GO CART) online in the NASA Goddard Earth Observing System version 4 model and simulated global aerosol distributions for the period 2000-2006. The new online system offers several advantages over the previous offline version, providing a platform for aerosol data assimilation, aerosol-chemistry-climate interaction studies, and short-range chemical weather forecasting and climate prediction. We introduce as well a methodology for sampling model output consistently with satellite aerosol optical thickness (AOT) retrievals to facilitate model-satellite comparison. Our results are similar to the offline GOCART model and to the models participating in the AeroCom intercomparison. The simulated AOT has similar seasonal and regional variability and magnitude to Aerosol Robotic Network (AERONET), Moderate Resolution Imaging Spectroradiometer, and Multiangle Imaging Spectroradiometer observations. The model AOT and Angstrom parameter are consistently low relative to AERONET in biomass-burning-dominated regions, where emissions appear to be underestimated, consistent with the results of the offline GOCART model. In contrast, the model AOT is biased high in sulfate-dominated regions of North America and Europe. Our model-satellite comparison methodology shows that diurnal variability in aerosol loading is unimportant compared to sampling the model where the satellite has cloud-free observations, particularly in sulfate-dominated regions. Simulated sea salt burden and optical thickness are high by a factor of 2-3 relative to other models, and agreement between model and satellite over-ocean AOT is improved by reducing the model sea salt burden by a factor of 2. The best agreement in both AOT magnitude and variability occurs immediately downwind of the Saharan dust plume.

Colarco, Peter↗

Mapping Aerosol Lidar Ratios Over Ocean Using Constrained Retrievals and A Global Aerosol Model

The current NASA Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) algorithms assign one lidar ratio (i.e., extinction-to-backscatter ratio; LR) value globally for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms through the development of regional and seasonal LR climatologies. In this work, aerosol LRs are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type (e.g., marine). The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. In this presentation, we show that the twelve-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine aerosols correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). For example, smaller SSVFs (< 65%) and larger LRs (> 55 sr), are found near land masses. This is indicative of the influence of over-land aerosols (e.g., pollution and biomass burning smoke). In the remote oceans (i.e., regions likely less impacted by non-sea salt aerosols), larger SSVFs (> 95%) and smaller LRs (< 25 sr) are found. The developed relationship between the GOCART SSVFs and MODIS AOD constrained LRs (polynomial fit intersect values of ~56 sr for SSVF of 0% and ~21 sr for SSVF of 100%) is used to produce model-assisted climatological LR maps on seasonal scales. Additionally, we show maps of inferred LRs from constrained retrievals using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) product and comparisons with those from the MODIS analyses. The technique demonstrated in this study not only benefits CALIPSO in the near-term, but similar methods can be applied to future spaceborne elastic backscatter lidars with collocated passive sensors, such as those associated with the upcoming NASA Atmosphere Observing System (AOS).

Travis Toth↗

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗

Geosat Data Assimilation with Application to the Eastern North Atlantic

An attempt is made to determine the three-dimensional ocean circulation from satellite altimeter measurements by assimilating Geosat sea surface height data into an eddy-resolving QuasiGeostrophic (QG) model of the eastern North Atlantic Ocean. Results are tested against independent information from hydrographic field observations and moored current meter data collected during the Geosat ERM. The comparison supports the concept of inferring aspects of the three-dimensional flow field from sea surface height observations by combining altimetric measurements with the dynamics of ocean circulation models. A Holland-type QG model with open boundaries is set up on a 2000 km X 2000 km domain of the eastern North Atlantic between 25 deg. and 45 deg. N, 32 deg. and 8 deg. W. By using a simple nudging technique, about two years of Geosat altimeter data are assimilated into the model every five days as space-time objective analyses on the model grid. The error information resulting from the analysis is used during the assimilation procedure to account for data uncertainties. Results show an intense eddy field, which in the surface layer interacts with a meandering Azores Front. Compared to Geosat, the model leads to smoothed fields that follow the observations. Model simulations are significantly correlated with hydrographic data from March 1988 and June 1989, both close to the surface and in the subsurface. Good agreement is also found between the model velocity fields and moored current meter data in the top two model layers. The agreement is visually weak in the bottom layer, although a coherence analysis reveals an agreement between the model simulation and current meter data over the full water column at periods exceeding 80 days.

Stammer, Detlef↗

Interdecadal Variations in the Alaska Gyre

Climatic dynamic topography variations in the Alaska gyre during the period 1968-1990 are described with an objective analysis of more than 12000 STD and XBT stations, and COADS wind stress data Interannual the dynamic height and SST variations were correlated and were consistent with recently described large-scale climatic shifts in the North Pacific. The gyre was centered more to the east, circulation appeared stronger, and SST was lower during the early to mid-1970s than during the 1980s. The Aleutian low (NP and PNA indices) intensified during the interim, but the response did not appear as a gyre spinup. Instead, the associated wind stress anomalies forced a slowly varying dynamic height anomaly across the eastern and northern part of the gyre through Ekman convergence, which had the effect of displacing the gyre's low somewhat to the WSW in the 1980s. The wind curl spectrum was white, and the slow oceanic response was modeled as stochastic-forced climate variability with a simple first-order Markov autoregression process. Forcing was assumed to be Ekman pumping of the pycnocline, and the damping coefficient was estimated from the data to be approx. 1 yr. A hindcast with observed winds gave estimated dynamic height patterns similar to those observed, with a canonical correlation of 0.79 at 99% confidence. This response was weak in the western half of the gyre, where slow baroclinic variability may have been influenced by long Rossby wave propagation. A simple autoregression simulation using artificial white noise forcing shows the evolution of decadal variations similar in nature to those observed. This result, along with the low frequency correlation between dynamic height and SST, suggests that the upper-ocean climatic variability in this region is primarily wind forced.

Lagerloef, Gary S.E.↗

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

Tracing the Path of Carbon Export in the Ocean though DNA Sequencing of Individual Sinking Particles

Surface phytoplankton communities were linked with the carbon they export into the deep ocean by comparing 18 S rRNA gene sequence communities from surface seawater and individually isolated sinking particles. Particles were collected in sediment traps deployed at locations in the North Pacific subtropical gyre and the California Current. DNA was isolated from individual particles, bulk-collected trap particles, and the surface seawater. The relative sequence abundance of exported phytoplankton taxa in the surface water varied across functional groups and ecosystems. Of the sequences detected in sinking particles, about half were present in large (>300 μm), individually isolated particles and primarily belonged to taxa with small cell sizes (<50 μm). Exported phytoplankton taxa detected only in bulk trap samples, and thus presumably packaged in the smaller sinking size fraction, contained taxa that typically have large cell sizes (>500 m). The effect of particle degradation on the detectable 18 S rRNA gene community differed across taxa, and differences in community composition among individual particles from the same location largely reflected differences in relative degradation state. Using these data and particle imaging, we present an approach that incorporates genetic diversity into mechanistic models of the ocean's biological carbon pump, which will lead to better quantification of the ocean’s carbon cycle.

Carbon Expert↗

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

During our participation in the NASA Crustal Dynamics Project under NASA contract NAS-27339 and grant NAG5-814 for the period 1982-1991, we published or submitted for publication 30 research papers and 52 abstracts of presentations at scientific meetings. In addition, five M.I.T. Ph.D. students (Eric Bergman, Steven Bratt, Dan Davis, Jeanne Sauber, Anne Sheehan) were supported wholly or in part by this project during their thesis research. Highlights of our research progress during this period include the following: application of geodetic data to determine rates of strain in the Mojave block and in central California and to clarify the relation of such strain to the San Andreas fault and Pacific-North American plate motions; application of geodetic data to infer post seismic deformation associated with large earthquakes in the Imperial Valley, Hebgen Lake, Argentina, and Chile; determination of the state of stress in oceanic lithosphere from a systematic study of the centroid depths and source mechanisms of oceanic intraplate earthquakes; development of models for the state of stress in young oceanic regions arising from the differential cooling of the lithosphere; determination of the depth extent and rupture characteristics of oceanic transform earthquakes; improved determination of earthquake slip vectors in the Gulf of California, an important data set for the estimation of Pacific-North American plate motions; development of models for the state of stress and mechanics of fold-and-thrust belts and accretionary wedges; development of procedures to invert geoid height, residual bathymetry, and differential body wave travel time residuals for lateral variations in the characteristic temperature and bulk composition of the oceanic upper mantle; and initial GPS measurements of crustal deformation associated with the Imperial-Cerro Prieto fault system in southern California and northern Mexico. Full descriptions of the research conducted on these topics may be found in the Semi-Annual status Reports submitted regularly to NASA over the course of this project and in the publications listed.

Solomon, Sean C.↗

The Southern Ocean Sea Surface Temperature Response to Ozone Depletion: A Multi-Model Comparison

The effect of the Antarctic ozone hole extends downwards from the stratosphere, with clear signatures in surface weather patterns including a positive trend in the Southern Annular Mode (SAM). Several recent studies have used coupled climate models to investigate the impact of these changes on Southern Ocean sea surface temperature (SST), notably motivated by the observed cooling from the late 1970s. Here we examine the robustness of these model results through comparison of both previously published and new simulations. We focus on the calculation of ‘climate response functions' (CRFs), transient responses to an instantaneous step-change in ozone concentrations. The CRF for most models consists of a rapid cooling of SST, followed by a slower warming trend. However, inter-model comparison reveals large uncertainties, such that even the sign of the impact of ozone depletion on historical SST, when reconstructed from the CRF, remains unconstrained. Comparison of these CRFs with SST responses to a hypothetical step-change in the SAM, inferred through lagged linear regression, shows broadly similar results. Causes of uncertainty are explored by examining relationships between model climatologies and their CRFs. The inter-model spread in CRFs can be reproduced by varying a single subgrid-scale mixing parameter within a single model. Antarctic sea-ice CRFs are also calculated: these do not generally exhibit the two-time-scale behavior of SST, suggesting a complex relationship between the two. Finally, by constraining model climatology-response relationships with observational values, we conclude that ozone depletion in unlikely to have been the primary driver of the observed SST cooling trend.

Climate response functions (CRF)↗

Global-Scale Hydrology: Simple Characterization of Complex Simulation

Atmospheric general circulation models (AGCMS) are unique and valuable tools for the analysis of large-scale hydrology. AGCM simulations of climate provide tremendous amounts of hydrological data with a spatial and temporal coverage unmatched by observation systems. To the extent that the AGCM behaves realistically, these data can shed light on the nature of the real world's hydrological cycle. In the first part of the seminar, I will describe the hydrological cycle in a typical AGCM, with some emphasis on the validation of simulated precipitation against observations. The second part of the seminar will focus on a key goal in large-scale hydrology studies, namely the identification of simple, overarching controls on hydrological behavior hidden amidst the tremendous amounts of data produced by the highly complex AGCM parameterizations. In particular, I will show that a simple 50-year-old climatological relation (and a recent extension we made to it) successfully predicts, to first order, both the annual mean and the interannual variability of simulated evaporation and runoff fluxes. The seminar will conclude with an example of a practical application of global hydrology studies. The accurate prediction of weather statistics several months in advance would have tremendous societal benefits, and conventional wisdom today points at the use of coupled ocean-atmosphere-land models for such seasonal-to-interannual prediction. Understanding the hydrological cycle in AGCMs is critical to establishing the potential for such prediction. Our own studies show, among other things, that soil moisture retention can lead to significant precipitation predictability in many midlatitude and tropical regions.

Koster, Randal D.↗

The Impact of the Assimilation of Aquarius Sea Surface Salinity Data in the GEOS Ocean Data Assimilation System

Ocean salinity and temperature differences drive thermohaline circulations. These properties also play a key role in the ocean-atmosphere coupling. With the availability of L-band space-borne observations, it becomes possible to provide global scale sea surface salinity (SSS) distribution. This study analyzes globally the along-track (Level 2) Aquarius SSS retrievals obtained using both passive and active L-band observations. Aquarius alongtrack retrieved SSS are assimilated into the ocean data assimilation component of Version 5 of the Goddard Earth Observing System (GEOS-5) assimilation and forecast model. We present a methodology to correct the large biases and errors apparent in Version 2.0 of the Aquarius SSS retrieval algorithm and map the observed Aquarius SSS retrieval into the ocean models bulk salinity in the topmost layer. The impact of the assimilation of the corrected SSS on the salinity analysis is evaluated by comparisons with insitu salinity observations from Argo. The results show a significant reduction of the global biases and RMS of observations-minus-forecast differences at in-situ locations. The most striking results are found in the tropics and southern latitudes. Our results highlight the complementary role and problems that arise during the assimilation of salinity information from in-situ (Argo) and space-borne surface (SSS) observations

Salinty↗

El Nino and the Global Ocean Observing System

Until a decade ago, an often-quoted expression in oceanography is that very few observations are recorded throughout the ocean. Now, the sentiment is no longer valid in the uppermost 10% of the tropical Pacific Ocean nor at the surface of the global ocean. One of the remarkable legacies of the 1985-1994 Tropical Oceans Global Atmosphere (TOGA) Program is an in situ marine meteorological and upper oceanographic measurement array throughout the equatorial Pacific to monitor the development and maintenance of El Nino episodes. The TOGA Observing System, which initially consisted of moored- and drifting-buoy arrays, a network of commercial ships, and coastal and island stations, now includes a constellation of satellites and data-assimilating models to simulate subsurface oceanographic conditions. The El Nino and La Nina tropical Pacific Ocean observing system represents the initial phase of an integrated global ocean observing system. Remarkable improvements have been made in ocean model simulation of subsurface currents, but some problems persist. For example, the simulation of the South Equatorial Current (SEC) remains an important challenge in the 2S-2N Pacific equatorial wave guide. During El Nino the SEC at the equator is reduced and sometimes the direction is reversed, becoming eastward. Both conditions allow warm water stored in the western Pacific to invade the eastern region, creating an El Nino episode. Assimilation of data is a tenet of faith to correct simulation errors caused by deficiencies in surface fluxes (especially wind stress) and parameterizations of subgrid-scale physical processes. In the first of two numerical experiments, the Pacific SEC was simulated with and without assimilation of subsurface temperature data. Along the equator, a very weak SEC occurred throughout the eastern Pacific, independent of assimilation of data. However, as displayed in the diagram, in the western Pacific there was no satisfactory agreement between the two simulations. To help determine reliability of the simulated SEC in the western Pacific, current measurements recorded during the 9-19 October 1994 voyage of the French research vessel L'Atalante are also shown in the diagram. With data assimilation, the simulated SEC was in much better agreement with L'Atalante observations. The simulated SEC with data assimilation was far from perfect, in part because of the sparsity of subsurface temperature observations. In the next experiment, TOPEX/POSEIDON sea surface height data in combination with subsurface temperatures will be assimilated to assess further improvement of the simulation of the SEC.

Halpern, David↗

Mapping Aerosol Lidar Ratios Over Ocean using MODIS AOD Constrained Retrievals and GOCART Model Simulations

After 17 years, the NASA Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission ceased science operations in August 2023. For the final CALIPSO data products release (Version 5), the CALIPSO project seeks to improve the accuracy of its aerosol extinction by advancing knowledge of aerosol lidar ratios (i.e., extinction-to-backscatter ratios; LRs) for various aerosol types. The current algorithm assigns one LR value globally for each of the seven tropospheric aerosol types. The CALIPSO team aims to improve the retrieval algorithm through the development of regional and seasonal LR climatologies for the same aerosol types. In this study, aerosol LRs are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data over oceans during daytime. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type. The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. This presentation will reveal findings that the 12-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine and dusty marine aerosols correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). For example, smaller SSVFs (< 65%) and larger LRs (> 55 sr), are found near land masses (Fig. 1). This indicates the influence of advected anthropogenic aerosols (e.g., pollution and biomass burning smoke). In the remote oceans (i.e., regions likely less impacted by non-sea salt aerosols), the SSVFs are larger (> 95%) and the LRs are smaller (< 25 sr) (Fig. 1). A polynomial fit of the MODIS AOD constrained LRs to the corresponding GOCART SSVFs (intersect values of ~58 sr for SSVF of 0% and ~21 sr for SSVF of 100%) is further used to produce model-assisted climatological LR maps on seasonal scales. Additionally, we will show results of a LR validation analysis for which we compare the revised CALIPSO AODs obtained by applying the seasonal/regional constrained LRs against CALIPSO Version 4.51 Ocean Derived Column Optical Depth (ODCOD). While the majority of the presentation will focus on LRs for CALIOP-classified marine and dusty marine aerosols, an overview of LR results will show preliminary results for other aerosol types over ocean, such as dust and elevated smoke. The technique demon-strated in this study highlights the benefits not only to the final planned CALIPSO data release in 2025, but similar methods can be applied to future spaceborne elastic backscatter lidars with collocated passive sensors (e.g., such as those associated with NASA’s proposed Atmosphere Observing System).

Travis D Toth↗

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↗

Computer modeling of large asteroid impacts into continental and oceanic sites: Atmospheric, cratering, and ejecta dynamics

Numerous impact cratering events have occurred on the Earth during the last several billion years that have seriously affected our planet and its atmosphere. The largest cratering events, which were caused by asteroids and comets with kinetic energies equivalent to tens of millions of megatons of TNT, have distributed substantial quantities of terrestrial and extraterrestrial material over much or all of the Earth. In order to study a large-scale impact event in detail, computer simulations were completed that model the passage of a 10 km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics associated with impact of the asteroid into two different targets, i.e., an oceanic site and a continental site. The calcuations were designed to broadly represent giant impact events that have occurred on the Earth since its formation and specifically represent an impact cratering event proposed to have occurred at the end of Cretaceous time. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock that expanded radially outward. Behind the shock front was a region of highly shock compressed and intensely heated air. Behind the asteroid, rapid expansion of this shocked air created a large region of very low density that also expanded away from the impact area. Calculations of the cratering events in both the continental and oceanic targets were carried to 120 s. Despite geologic differences, impacts in both targets developed comparable dynamic flow fields, and by approx. 29 s similar-sized transient craters approx. 39 km deep and approx. 62 km across had formed. For all practical purposes, the atmosphere was nearly completely removed from the impact area for tens of seconds, i.e., air pressures were less than fractions of a bar out to ranges of over 50 km. Consequently, much of the asteroid and target materials were ejected upward into a near vacuum. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to add substantial solid and vaporized material to the atmosphere, but these conditions were not studied.

Roddy, D. J.↗

Evaluation of ERA-interim and MERRA Cloudiness in the Southern Oceans

The Southern Ocean cloud cover modeled by the Interim ECMWF Re-Analysis (ERA-Interim) and Modern- Era Retrospective Analysis for Research and Applications (MERRA) reanalyses are compared against Moderate Resolution Imaging Spectroradiometer (MODIS) and Multiangle Imaging Spectroradiometer (MISR) observations. ERA-Interim monthly mean cloud amounts match the observations within 5%, while MERRA significantly underestimates the cloud amount. For a compositing analysis of clouds in warm season extratropical cyclones, both reanalyses show a low bias in cloud cover. They display a larger bias to the west of the cyclones in the region of subsidence behind the cold fronts. This low bias is larger for MERRA than for ERA-Interim. Both MODIS and MISR retrievals indicate that the clouds in this sector are at a low altitude, often composed of liquid, and of a broken nature. The combined CloudSat-Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) cloud profiles confirm these passive observations, but they also reveal that low-level clouds in other parts of the cyclones are also not properly represented in the reanalyses. The two reanalyses are in fairly good agreement for the dynamic and thermodynamic characteristics of the cyclones, suggesting that the cloud, convection, or boundary layer schemes are the problem instead. An examination of the lower-tropospheric stability distribution in the cyclones from both reanalyses suggests that the parameterization of shallow cumulus clouds may contribute in a large part to the problem. However, the differences in the cloud schemes and in particular in the precipitation processes, which may also contribute, cannot be excluded.

bias↗