The ECCO routine global ocean data assimilation system
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The Proudman Oceanographic Laboratory (POL), the Institute of Oceanographic Sciences Deacon Laboratory (IOSDL), and the British National Space Centre (BNSC) plan to use radar altimeter data from the TOPEX/POSEIDON mission in five main areas of research. These are: (1) investigation of variations in Southern Ocean sea levels and circulation; (2) related numerical modeling and data assimilation; (3) global and regional ocean tides; (4) North Sea and northeast Atlantic studies; and (5) winds and waves. Each one of these topics is described briefly below while more information can be extracted from the United Kingdom response to the TOPEX/POSEIDON Announcement of Opportunity.
The nature of global, large-scale sea level variability is studied using a numerical ocean circulation model.
Numerical experiments using simulated altimeter data were conducted in order to examine the assimilation of altimeter-derived sea surface heights into numerical ocean circulation models. A reduced-gravity, primitive equation circulation model of the Gulf of Mexico was utilized; the Gulf of Mexico was chosen because of its amenability to modeling and the ability of low vertical-mode models to reproduce the observed dynamical features of the Gulf circulation. The simulated data were obtained by flying an imaginary altimeter over the model ocean and sampling the model sea surface just as real altimeter would observe the true ocean. The data were used to initialize the numerical model and the subsequent forecast was compared to the true numerical solution. Results indicate that for a stationary, circular eddy, approximately three to four tracks (either ascending or descending) across the eddy are sufficient to ensure adequate spatial resolution.
Wavelet analysis of NASA scatterometer (NSCAT) backscatter and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) radiance data can be used to obtain daily sea ice drift information for the Arctic region. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars. Comparisons with ice motion derived from ocean buoys give good quantitative agreement. Both comparison results from NSCAT and SSM/I are compatible, and the results from NSCAT can definitely complement that from SSM/I when there are cloud or surface effects. Then three sea-ice drift daily results from NSCAT, SSM/I, and buoy data can be merged as a composite map by some data fusion techniques. The ice flow streamlines are highly correlated with surface air pressure contours. Examples of derived ice-drift maps in December 1996 illustrate large-scale circulation reversals over a period of four days. A method for deriving divergence and shear at the large-scale has been developed and comparison between buoys and satellite results shows a good agreement. These calibrated/validated results indicate that NSCAT, SSM/I merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model. For demonstration purpose, the ice velocities derived from satellite data are compared with the ice velocities derived from a coupled ice-ocean interaction model. The comparison reveals that the general circulation patterns of the two are quite similar but the ice velocity differences between the two are quite significant. In order to quantify the wind effects on ice motion, empirical orthogonal functions (EOF) are used in the principal component analysis for both ice motion and pressure field. Some preliminary results of sea-ice motion from QuikScat will also be presented.
The objectives of this study are the following: (1) to assess the quality of the sea level observed by altimetry on a global scale, especially in tropical regions where the atmospheric effects are the most critical, and to prepare the altimeter data for their assimilation into tropical Oceanic General Circulation Models; and (2) to validate and calibrate the altimeter-derived sea state parameters, to assimilate these parameters into numerical models, and to estimate the altimeter's sea state and electromagnetic biases.
Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.
A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.
Ocean circulation is a critical factor in determining the Earth's climate. Satellite altimetry has been proven a powerful technique for measuring the height of the sea surface for the study of global ocean circulation dynamics. A major objective of my research is to investigate the utility of altimeter data for ocean circulation studies. The 6 years' data record of TOPEX/POSEIDON have been analyzed to study the spatial and temporal characteristics of large-scale ocean variability. A major result obtained in 1998 is the discovery of large-scale oscillations in sea level with a period of 25 days in the Argentine Basin of the South Atlantic Ocean (see diagram). They exhibit a dipole pattern with counterclockwise rotational propagation around the Zapiola Rise (centered at 45S and 317E), a small seamount in the abyssal plain of the basin. The peak-to-trough amplitude is about 10 cm over a distance of 500-1000 km. The amplitude of these oscillations has large seasonal-to-interannual variations. The period and rotational characteristics of these oscillations are remarkably similar to the observations made by two current meters deployed near the ocean bottom in the region. What TOPEX/POSEIDON has detected apparently are manifestations of the movement of the entire water column (barotropic motion). The resultant transport variation is estimated to be about 50 x 10(exp 6) cubic M/S, which is about 50% of the total water transport in the region. Preliminary calculations suggest that these oscillations are topographically trapped waves. A numerical model of the South Atlantic is used to investigate the nature of and causes for these waves. A very important property of sea surface height is that it is directly related to the surface geostrophic velocity, which is related to deep ocean circulation through the density field. Therefore altimetry observations are not only useful for determining the surface circulation but also for revealing information about the deep ocean. Another objective of my research is to investigate the extent of information one can infer about the deep ocean from altimetric observations. The approach is data assimilation by ocean models. The objective is to develop a modeling/data assimilation system that will produce estimates of the 3-dimensional state of the ocean for the entire duration of the TOPEX/POSEIDON Mission.
The use of model-based global data sets of atmospheric circulation for studying fundamental dynamical and physical processes is discussed, focusing on limitations of the available model-based data sets. Data from the Global Weather Experiment in 1979 were analyzed by two authorized level IIIb data centers in 1980 and in 1981. The analyses led to difference in data-sparse regions such as the tropics. Study areas which can be addressed by an internally-consistent long-term multivariate data set for the atmospheric circulation are considered, including mean climate, forcing for the ocean models, global hydrological cycle, atmospheric energetics, intraseasonal variability, land surface processes, and structure and variability of vertical velocity, divergence, and diabatic heating. It is concluded that the most comprehensive technique for integrating space and in situ observations to produce this type of data set would be a four-dimensional data assimilation system with a realistic physical model of the type employed in operational numerical weather prediction.
We will use TOPEX/POSEIDON altimetry in combination with ongoing and planned efforts, including data from the European Remote Sensing satellite (ERS-1), in situ measurements designed specifically to complement satellite altimetry, and assimilation of these several data types into eddy-resolving numerical models in order to understand the dynamics and energetics of midlatitude western boundary currents (specifically, the Gulf Stream and the Kuroshio Extension). A better understanding of the recirculation of midlatitude gyres can best be undertaken in the format of such regional, eddy-resolving experiments. Such regional programs will enhance and be enhanced by the larger scale circulation studies of the World Ocean Circulation Experiment (WOCE) as well as by global-scale eddy-resolving models that we will develop prior to the TOPEX/POSEIDON mission. This effort includes participation on the TOPEDX/POSEIDON Science Team.
There are four distinct but related tasks in this investigation. First, we want to combine estimates of the surface dynamic topography derived independently from altimetry and from hydrography in an optimal way (weighted least squares with a priori estimates) to produce an optimal estimate (with error bounds) of the large-scale circulation. Second, we want to study the large-scale (as opposed to the mesoscale) variability using crossover differences as well as exact repeat tracks. The third task involves the study of mesoscale variability in the Kuroshio Extension, the California Current, and eventually the Pacific Ocean and the whole globe, wherein surface transport estimates are provided for both the Kuroshio Extension and the Gulf Stream. The final task is a study of how altimeter data and other data can be assimilated into numerical models.
Reanalyses are regularly gridded, retrospective depictions of the physical earth system, which are produced through the correction of a short-term forecast to available observations. In the Arctic, reanalyses are particularly well suited to marshal the sparse observing network to provide a plausible, multivariate representation of conditions. Atmospheric reanalyses such as MERRA-2 (NASA Modern-Era Retrospective analysis for Research and Applications, version 2) and ocean reanalyses such as SODA3 (Univ. Maryland Simple Ocean Data Assimilation version 3) are widely used in Arctic research for diagnostic studies of circulation, model evaluation, and as boundary conditions for a variety of process models. Here, we provide examples that illustrate the utility of reanalyses for providing information on the spatial and temporal scales of recent, rapid changes in the Arctic. Recent trends in Arctic surface temperatures, surface melt over Greenland and Arctic glaciers, and evolving freshwater conditions in the Arctic Ocean are examples where reanalyses can provide information that cannot easily be obtained via other means. These examples provide information on the scale, magnitude, and the uncertainty of recent Arctic change and provide a context for future scenarios. We further quantify uncertainties in key reanalyses variables and approaches for addressing these issues.
Models are foundational for estimating states of the earth's climate system, both as tools to extrapolate information in time and space, and as observation 'operators' used to relate what is analyzed and predicted to what is observed. Expanding the simulation approach further, observing system simulation experiments (OSSEs) are designed to mimic the complete process of analyzing the climate state by replacing real observations with entirely simulated ones determined from a model-based depiction of nature. OSSEs provide a framework to 'fly' simulated satellite instruments through a synthetic atmosphere and investigate the trade-spaces of measurements for various satellite configurations and sampling strategies, and assess their measurement impact on modeling and forecasting capabilities. Such a tool is a crucial but as yet unfulfilled need for future mission selection and design. The components of a state-of-the-art OSSE system are being assembled at the Global Modeling and Assimilation Office (GMAO, Code 610.1) at NASA/GSFC, leveraging on the GMAO's existing modeling and data assimilation infrastructure for numerical weather prediction (NWP). The OSSE framework is based on the GMAO's Goddard Earth Observing System atmospheric general circulation model, version 5 (GEOS-5) and the Gridpoint Statistical Interpolation (GSI) observational analysis scheme, combined with the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) model developed by the Atmospheric Chemistry and Dynamics Branch (Code 613.3). This system is an evolving, key component of Goddard's planned development of an Integrated Earth System Analysis (IESA) capability, which will bring together into a single, fully interactive system Goddard's modeling and assimilation efforts in atmosphere, ocean and chemistry and aerosols to provide a comprehensive analysis and prediction system for weather and climate In addition to providing a state-of-the-art capability for assimilating current observation types, GEOS-5, and the future IESA, provide the capability to identify the need for, and assess the potential impact of, future observing systems under consideration for improving weather and climate prediction.
Approximately 240 Tg of mineral dust aerosol are transported annually from Saharan Africa to the Atlantic Ocean. Dust affects the Earth radiation budget, and plays direct (through scattering and absorption of radiation) and indirect (through modification of cloud properties and environment) roles in climate. Deposition of dust to the surface provides an important nutrient source to terrestrial and oceanic ecosystems. Dust is additionally a contributor to adverse air quality. Among the tools toward understanding the lifecycle and impacts of mineral dust aerosols are numerical models. Important constraints on these models come from quantitative satellite observations, like those from the space-based Moderate Resolution Imaging Spectroradiometer (MODIS). In particular, Kauhan et al. [2005] used MODIS aerosol observations to infer transport and deposition fluxes of Saharan dust over the Atlantic, Caribbean, and Amazonian basins. Those observations are used here to constrain the transport of dust and its interannual variability simulated in the NASA GEOS-5 general circulation model and data assimilation system. Significant uncertainty exists in the MODIS-derived fluxes, however, due to uncertainty in the wind fields provided by meteorological analyses in this region. That same uncertainty in the wind fields is manifest in our GEOS-5 simulations of dust distributions. Here we use MODIS observations to investigate the seasonality and location of the Saharan dust plume and explore through sensitivity analysis of our model the meteorological controls on the dust distribution, including dust direct radiative effects and sub-gridscale source and sink processes.
A Global Ocean Observing System (GOOS) should be established now with international coordination (1) to address issues of global change, (2) to implement operational ENSO forecasts, (3) to provide the data required to apply global ocean circulation models, and (4) to extract the greatest value from the one billion dollar investment over the next ten years in ocean remote sensing by the world's space agencies. The objectives of GOOS will focus on climatic and oceanic predictions, on assessing coastal pollution, and in determining the sustainability of living marine resources and ecosystems. GOOS will be a complete system including satellite observations, in situ observations, numerical modeling of ocean processes, and data exchange and management. A series of practical and economic benefits will be derived from the information generated by GOOS. In addition to the marine science community, these benefits will be realized by the energy industries of the world, and by the world's fisheries. The basic oceanic variables that are required to meet the oceanic and predictability objectives of GOOS include wind velocity over the ocean, sea surface temperature and salinity, oceanic profiles of temperature and salinity, surface current, sea level, the extent and thickness of sea ice, the partial pressure of CO2 in surface waters, and the chlorophyll concentration of surface waters. Ocean circulation models and coupled ocean-atmosphere models can be used to evaluate observing system design, to assimilate diverse data sets from in situ and remotely sensed observations, and ultimately to predict future states of the system. The volume of ocean data will increase enormously over the next decade as new satellite systems are launched and as complementary in situ measuring systems are deployed. These data must be transmitted, quality controlled, exchanged, analyzed, and archived with the best state-of-the-art computational methods.