Search NASA⌕ Search

SEARCH · Search NASA

Results for “sea state estimation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Intercomparison of Satellite Derived Gravity Time Series with Inferred Gravity Time Series from TOPEX/POSEIDON Sea Surface Heights and Climatological Model Output

The upcoming GRACE mission promises to open a window on details of the global mass budget that will have remarkable clarity, but it will not directly answer the question of what the state of the Earth's mass budget is over the critical last quarter of the 20th century. To address that problem we must draw upon existing technologies such as SLR, DORIS, and GPS, and climate modeling runs in order to improve our understanding. Analysis of long-period geopotential changes based on SLR and DORIS tracking has shown that addition of post 1996 satellite tracking data has a significant impact on the recovered zonal rates and long-period tides. Interannual effects such as those causing the post 1996 anomalies must be better characterized before refined estimates of the decadal period changes in the geopotential can be derived from the historical database of satellite tracking. A possible cause of this anomaly is variations in ocean mass distribution, perhaps associated with the recent large El Nino/La Nina. In this study, a low-degree spherical harmonic gravity time series derived from satellite tracking is compared with a TOPEX/POSEIDON-derived sea surface height time series. Corrections for atmospheric mass effects, continental hydrology, snowfall accumulation, and ocean steric model predictions will be considered.

Cox, C.↗
NASA NTRS2018

Comparison of MERRA-2 and ECCO-V4 Ocean Surface Heat Fluxes: Consequences of Different Forcing Feedbacks on Ocean Circulation and Implications for Climate Data Assimilation

Forcing ocean models with reanalysis data is a common practice in ocean modeling. As part of this practice, prescribed atmospheric state variables and interactive ocean SST (Sea Surface Temperature) are used to calculate fluxes between the ocean and the atmosphere. When forcing an ocean model with reanalysis fields, errors in the reanalysis data, errors in the ocean model and errors in the forcing formulation will generate a different solution compared to other ocean reanalysis solutions (which also have their own errors). As a first step towards a consistent coupled ocean-atmosphere reanalysis, we compare surface heat fluxes from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to heat fluxes from a state-of-the-art oceanic reanalysis, the Estimating the Circulation and Climate of the Ocean Version 4, Release 2 (ECCO-v4). Then, we investigate the errors associated with the MITgcm (Massachusetts Institute of Technology general circulation model) ocean model in its ECCO-v4 ocean reanalysis configuration (1992-2011) when it is forced with MERRA- 2 atmospheric reanalysis fields instead of with the ECCO-v4 adjoint optimized ERA-interim state variables. This is done by forcing ECCO-v4 ocean with and without feedbacks from MERRA-2 related to turbulent fluxes of heat and moisture and the outgoing long wave radiation. In addition, we introduce an intermediate forcing method that includes only the feedback from the interactive outgoing long wave radiation. The resulting ocean circulation is compared with ECCO-v4 reanalysis and in-situ observations. We show that, without feedbacks, imbalances in the energy and the hydrological cycles of MERRA-2 (which are directly related to the fact it was created without interactive ocean) result in considerable SST drifts and a large reduction in sea level. The bulk formulae and interactive outgoing long wave radiation, although providing air-sea feedbacks and reducing model-data misfit, strongly relax the ocean to observed SST and may result in unwanted features such as large change in the water budget. These features have implications in a desired forcing recipe to be used. The results strongly and unambiguously argue for next generation data assimilation climate studies to involve fully coupled systems.

MERRA-2↗
NASA NTRS1991

The effect of the degree of wave development on the sea state bias in radar altimetry measurement

The extent by which the wave development affects the sea state bias (SSB) in the Geosat altimeter height measurements is investigated. It is shown that the altimetric SSB is not a linear function of the wave height but that it also depends on other factors of wave development, particularly on the wave age (defined as the ratio of the phase speed of the dominant ocean waves to the ocean wind speed). The wave age can be estimated crudely, as the so-called 'pseudo-wave age', on the basis of the significant wave height H(1/3) and the wind speed measured by the altimeter. The analysis of 2.7 years' worth of Geosat data showed a general trend for the dependence of the SSB on the pseudo-wave age, which agrees well with the theoretical prediction that, for a given H(1/3), the SSB decreases as the pseudo-wave age increases.

Fu, Lee-Lueng↗
NASA NTRS2009

A Decade of Satellite Ocean Color Observations

After the successful Coastal Zone Color Scanner (CZCS, 1978-1986), demonstration that quantitative estimations of geophysical variables such as chlorophyll a and diffuse attenuation coefficient could be derived from top of the atmosphere radiances, a number of international missions with ocean color capabilities were launched beginning in the late 1990s. Most notable were those with global data acquisition capabilities, i.e., the Ocean Color and Temperature Sensor (OCTS 1996-1997), the Sea-viewing Wide Field-of-view Sensor (SeaWiFS, United States, 1997-present), two Moderate Resolution Imaging Spectroradiometers, (MODIS, United States, Terra/2000-present and Aqua/2002-present), the Global Imager (GLI, Japan, 2002-2003), and the Medium Resolution Imaging Spectrometer (MERIS, European Space Agency, 2002-present). These missions have provided data of exceptional quality and continuity, allowing for scientific inquiries into a wide variety of marine research topics not possible with the CZCS. This review focuses on the scientific advances made over the past decade using these data sets.

McClain, Charles R.↗
NASA NTRS2014

Background Error Covariance Estimation Using Information from a Single Model Trajectory with Application to Ocean Data Assimilation

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory.SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Error Covariance↗
NASA NTRS2014

Background Error Covariance Estimation using Information from a Single Model Trajectory with Application to Ocean Data Assimilation into the GEOS-5 Coupled Model

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory. SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Data Assimilation↗
NASA NTRS2013

Comparison Between NPP-VIIRS Aerosol Data Products and the MODIS AQUA Deep Blue Collection 6 Dataset Over Land

Aerosols are small particles suspended in the atmosphere and have a variety of natural and man-made sources. Knowledge of aerosol optical depth (AOD), which is a measure of the amount of aerosol in the atmosphere, and its change over time, is important for multiple reasons. These include climate change, air quality (pollution) monitoring, monitoring hazards such as dust storms and volcanic ash, monitoring smoke from biomass burning, determining potential energy yields from solar plants, determining visibility at sea, estimating fertilization of oceans and rainforests by transported mineral dust, understanding changes in weather brought upon by the interaction of aerosols and clouds, and more. The Suomi-NPP satellite was launched late in 2011. The Visible Infrared Imaging Radiometer Suite (VIIRS) aboard Suomi-NPP is being used, among other things, to determine AOD. This study compares the VIIRS dataset to ground-based measurements of AOD, along with a state-of-the-art satellite AOD dataset (the new version of the Moderate Resolution Imaging Spectrometer Deep Blue algorithm) to assess its reliability. The Suomi-NPP satellite was launched late in 2011, carrying several instruments designed to continue the biogeophysical data records of current and previous satellite sensors. The Visible Infrared Imaging Radiometer Suite (VIIRS) aboard Suomi-NPP is being used, among other things, to determine aerosol optical depth (AOD), and related activities since launch have been focused towards validating and understanding this new dataset through comparisons with other satellite and ground-based products. The operational VIIRS AOD product is compared over land with AOD derived from Moderate Resolution Imaging Spectrometer (MODIS) observations using the Deep Blue (DB) algorithm from the forthcoming Collection 6 of MODIS data

Sayer, Andrew M.↗
NASA NTRS1984

Remote sensing of air-sea interactions

A number of preliminary concepts for the measurement or inference of fluxes across the air-sea interface through remote sensing are proposed. All the methods are achievable from aircraft with state-of-the-art technology. Only one is now ready for space implementation. The focus is on cold outbreaks. Sensible (latent) heat flux is inferred from the difference between initial surface air temperature (vapor mixing ratio) and the downwind SST (and corresponding saturation mixing ratio). The downwind growth rate of the PBL as measured by lidar also provides estimates of surface heating and the cross-inversion entrainment velocity. The lidar also provides a measure of the depth of the inversion and its penetration by surface-forced convection; this permits estimates of the surface heat flux. Lidar and radiometric measurements of cloud top height and temperature provide means of deducing the temperature sounding downstream so that heating is computed with the aid of a known sounding upstream.

Atlas, D.↗
NASA NTRS1999

Detecting 1mm/Year Signals in Altimetric Global Sea Level: Effect of Atmospheric Water Vapor and Precipitation

Several research efforts exist to use Topography Experiment (TOPEX)/ Projet d'Observatorie de Surveillance et d'Etudes Integrees de la Dynamique des Oceans (Poseidon) (T/P) to detect changes in global sea level possibly associated with climate change. This requires much better than 1 mm/yr accuracy, something that none of the instruments in T/P [or the European Remote Sensing (ERS-2) satellite, or the U.S. Navy's Geosat Follow-On (GFO) satellite] were designed for. This work focuses on the ability of the T/P microwave radiometer (TMR) to retrieve the path delay due to atmospheric water vapor along the altimeter's path with accuracy in the time changes below 1 mm/yr on global average. In collaboration with Stephen Keihm of JPL and Christopher Ruf of Pennsylvania State University, we compared TMR path delay (PD) estimates with atmospheric precipitable water (PW) from the Special Sensor Microwave Imager (SSMI) aboard the Defense Meteorological Satellite Program (DMSP) series of satellites for 1992-1998 to selected radiosondes, and we also looked at the brightness temperatures measured by TMR in the lowest 1% of the histogram. The conclusion is that TMR had a slow instrumental drift, associated with the 18-GHz channel, which causes an approximate underestimation of water vapor at a rate equivalent to 1.2 mm/yr in path delay between 1992 and 1996; this effect stopped and no drift is detected in 1997. The same study concluded that there is no detectable scale error (one which is proportional to measured vapor) in TMR. In related work, carried out with graduate student Damien Cailliau, we investigated the relative abilities of TMR, SSMI and the UP dual-frequency radar altimeter to detect rain, relative to a climatology of shipborne observations. Rain is a crucial but poorly measured variable in studies of the climate system, and a dedicated mission, Tropical Rainfall Measuring Mission (TRMM), was recently launched to measure it. However, the climatologies built over the past 10 years used passive radiometers, such as SSMI, or infrared observation from National Oceanic and Atmospheric Administration (NOAA) geostationary satellites. We concluded that the dual-frequency altimeter does an excellent job at retrieving rain on scales much smaller than the passive radiometer, and a better job at retrieving total precipitation (rain as well as snow) at high latitudes. This work improves the sea level observations by allowing better flagging of observations contaminated by rain, and more importantly, provides an alternative way to measure rainfall over the oceans.

Zlotnicki, Victor↗
NASA NTRS1986

Role of seasonality in the evolution of climate during the last 100 million years

A simple climate model has been used to calculate the effect of past changes in the land-sea distribution on the seasonal cycle of temperatures during the last 100 million years. Modeled summer temperature decreased over Greenland by more than 10 C and over Antarctica by 5 to 8 C. For the last 80 million years, this thermal response is comparable in magnitude to estimated atmospheric carbon dioxide effects. Analysis of paleontological data provides some support for the proposed hypothesis that large changes due to seasonality may have sometimes resulted in an ice-free state due to high summer temperatures rather than year-round warmth. Such 'cool' nonglacials may have prevailed for as much as one-third of the last 100 million years.

Crowley, T. J.↗
NASA NTRS2014

Assessing the Impact of Laurentide Ice-sheet Topography on Glacial Climate

Simulations of past climates require altered boundary conditions to account for known shifts in the Earth system. For the Last Glacial Maximum (LGM) and subsequent deglaciation, the existence of large Northern Hemisphere ice sheets caused profound changes in surface topography and albedo. While ice-sheet extent is fairly well known, numerous conflicting reconstructions of ice-sheet topography suggest that precision in this boundary condition is lacking. Here we use a high-resolution and oxygen-isotopeenabled fully coupled global circulation model (GCM) (GISS ModelE2-R), along with two different reconstructions of the Laurentide Ice Sheet (LIS) that provide maximum and minimum estimates of LIS elevation, to assess the range of climate variability in response to uncertainty in this boundary condition.We present this comparison at two equilibrium time slices: the LGM, when differences in ice-sheet topography are maximized, and 14 ka, when differences in maximum ice-sheet height are smaller but still exist. Overall, we find significant differences in the climate response to LIS topography, with the larger LIS resulting in enhanced Atlantic Meridional Overturning Circulation and warmer surface air temperatures, particularly over northeastern Asia and the North Pacific. These up- and downstream effects are associated with differences in the development of planetary waves in the upper atmosphere, with the larger LIS resulting in a weaker trough over northeastern Asia that leads to the warmer temperatures and decreased albedo from snow and sea-ice cover. Differences between the 14 ka simulations are similar in spatial extent but smaller in magnitude, suggesting that climate is responding primarily to the larger difference in maximum LIS elevation in the LGM simulations. These results suggest that such uncertainty in ice-sheet boundary conditions alone may significantly impact the results of paleoclimate simulations and their ability to successfully simulate past climates, with implications for estimating climate sensitivity to greenhouse gas forcing utilizing past climate states.

glaciology↗
NASA NTRS1999

Ocean-Atmosphere Interaction in Climate Changes

The diagram, which attests the El Nino teleconnection observed by the NASA Scatterometer (NSCAT) in 1997, is an example of the results of our research in air-sea interaction - the core component of our three-part contribution to the Climate Variability Program. We have established an interplay among scientific research, which turns spacebased data into knowledge, a push in instrument technology, which improves observations of climate variability, and an information system, which produces and disseminates new data to support our scientific research. Timothy Liu led the proposal for advanced technology, in response to the NASA Post-2002 Request for Information. The sensor was identified as a possible mission for continuous ocean surface wind measurement at higher spatial resolution, and with the unique capability to measure ocean surface salinity. He is participating in the Instrument Incubator Program to improve the antenna technology, and is initiating a study to integrate the concept on Japanese missions. He and his collaborators have set up a system to produce and disseminate high level (gridded) ocean surface wind/stress data from NSCAT and European missions. The data system is being expanded to produce real-time gridded ocean surface winds from Quikscat, and precipitation and evaporation from the Tropical Rain Measuring Mission. It will form the basis for a spacebased data analysis system which will include momentum, heat and water fluxes. The study on 1997 El Nino teleconnection illustrates our interdisciplinary and multisensor approach to study climate variability. The diagram shows that the collapse of trade wind and the westerly wind anomalies in the central equatorial Pacific led to the equatorial ocean warming. The equatorial wind anomalies are connected to the anomalous cyclonic wind pattern in the northeast Pacific. The anomalous warming along the west coast of the United States is the result of the movement of the pre-existing warm sea surface temperature anomalies with the cyclonic wind anomalies toward the coast. The results led to a new study which identifies decadal ocean variations in the Northeast Pacific. Three studies of oceanic responses to wind forcing caused by the seasonal change of monsoons, the passage of a typhoon, and the 1997 El Nino, were successfully conducted. Besides wind forcing, we continue to examine new techniques for estimating thermal and hydrologic fluxes, through the inverse ocean mixed-layer model, through divergence of atmospheric water transport, and by direct retrieval from radiances observed by microwave radiometers. Greenhouse warming has been linked to water vapor measured by two spaceborne sensors in two studies. In the first study, strong baroclinicity and deep convection were found to transport water vapor to the upper atmosphere and increase greenhouse trapping over the storm tracks of the North Pacific and Atlantic. In another study, the annual cycle of greenhouse warming were related to sea surface temperature (SST) and integrated water vapor, and the latitudinal dependence of the magnitudes and phases of the annual cycles were compared.

Liu, W. Timothy↗
NASA NTRSDate not supplied

Blue Carbon Stock Estimation of US Coastal Marshes from New York to Georgia How Much Do We Have to Lose?

Valued for their myriad multifunctional roles in preserving habitat for biodiversity, cycling nutrients, protecting shorelines, improving water quality, and sequestering blue carbon, coastal and riverine marshes continue to disappear not only due to land use change but also sea level rise. Blue carbon’s role in the carbon cycle has been underappreciated until the last few decades, and the carbon stocks for U.S. Atlantic coastal margins have previously been estimated for a depth of 0.5-1 meter. We select several North American East Coast databases that have been utilized to study sea level rise and instead use them to estimate marsh depth from New York to Georgia; in local marshes that lack depth data, we present new depth data and perform several full-depth transects. This alongside marsh area data from the National Wetlands Inventory (NWI) enables us to estimate marsh volume. We then utilize an average carbon density of 27 kg/m 3 by Holmquist et al. 2018 to calculate blue carbon stocks. We present the inventories of peat depth by state, marsh area by state, and total marsh volume by state, followed by the blue carbon stocks. The peat depth inventory of eastern U.S. marshes reveals that riverine marshes are significantly deeper than coastal marshes. While both are important, we highlight the jeopardy of rapidly eroding urban marshes due to the heavy metals these marshes have sequestered for centuries as they contaminate the surrounding estuary. Uncertainties and caveats include carbon density values, loss of marsh since the NWI assessment, and depth averages. The high cost of marsh disappearance underscores the urgent need for their preservation.

Coastal marshes↗
NASA NTRS1991

Effective gain measurements in chromium-doped forsterite

Effective gain cross section in tetravalent chromium-doped forsterite laser crystal was measured over the 1180-1330 nm spectral range. The experiment was performed using two collinear laser beams in a pump-and-probe arrangement. The peak-gain cross section from this measurement is estimated to be 1.9 x 10 to the -19th sq cm at 1215 nm, which is comparable to the value of about 2 x 10 to the -19th sq cm predicted by fluorescence linewidth and lifetime measurements. These results indicate that excited-state absorption is not a major loss mechanism in tetravalent chromium-doped forsterite.

Petricevic, V.↗
NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗
NASA NTRS2013

The GLAS Algorithm Theoretical Basis Document for Precision Orbit Determination (POD)

The Geoscience Laser Altimeter System (GLAS) was the sole instrument for NASA's Ice, Cloud and land Elevation Satellite (ICESat) laser altimetry mission. The primary purpose of the ICESat mission was to make ice sheet elevation measurements of the polar regions. Additional goals were to measure the global distribution of clouds and aerosols and to map sea ice, land topography and vegetation. ICESat was the benchmark Earth Observing System (EOS) mission to be used to determine the mass balance of the ice sheets, as well as for providing cloud property information, especially for stratospheric clouds common over polar areas. The GLAS instrument operated from 2003 to 2009 and provided multi-year elevation data needed to determine changes in sea ice freeboard, land topography and vegetation around the globe, in addition to elevation changes of the Greenland and Antarctic ice sheets. This document describes the Precision Orbit Determination (POD) algorithm for the ICESat mission. The problem of determining an accurate ephemeris for an orbiting satellite involves estimating the position and velocity of the satellite from a sequence of observations. The ICESatGLAS elevation measurements must be very accurately geolocated, combining precise orbit information with precision pointing information. The ICESat mission POD requirement states that the position of the instrument should be determined with an accuracy of 5 and 20 cm (1-s) in radial and horizontal components, respectively, to meet the science requirements for determining elevation change.

Rim, Hyung Jin↗
NASA NTRS2003

Optimal Observations for Variational Data Assimilation

An important aspect of Ocean state estimation is the design of an observing system that allows the efficient study of climate aspects in the ocean. A solution of the design problem is presented here in terms of optimal observations that emerge as nondimensionalized singular vectors of the modified data resolution matrix. The actual computation is feasible only for scalar quantities in the limit of large observational errors. In the framework of a lo resolution North Atlantic primitive equation model it is demonstrated that such optimal observations when applied to determining the strength of the volume and heat transport across the Greenland-Scotland ridge, perform significantly better than traditional section data. On seasonal to inter-annual time-scales optimal observations are located primarily along the continental shelf and information about heat-transport, wind stress and stratification is being communicated via boundary waves and advective processes. On time-scales of about a month, sea surface height observations appear to be more efficient in reconstructing the cross-ridge heat transport than hydrographic observations. Optimal observations also provide a tool for understanding how the ocean state is effected by anomalies of integral quantities such as meridional heat transport.

Koehl, Armin↗
NASA NTRS2003

Collaborative Research: Atlantic Ocean Tropical/Subtropical Processes from Seasonal to Decadal Time Scales: Model/Data, Model/Model Comparison and Model/Data Synthesis Through Assimilation

The effort of this first year of research has been focused on the assimilation of TOPEX/Poseidon altimetric data into a primitive equation model of the Atlantic tropical/subtropical circulation. A reduced-rank, stationary Kalman filter has been constructed to assimilate the altimetric sea surface height anomaly (SHA) into the model. The goal is to assess how the inter-hemispheric transports between the Atlantic subtropics and tropics are affected by the assimilation and how the subsurface thermocline structure , and its variability ,is dynamically constrained by the SHA. The model is a reduced-gravity primitive equation GCM of the upper Atlantic Ocean between 30 S and 30 N. The assimilation scheme is an approximation to the extended Kalman filter in which the error covariances of the state estimates are calculated only in a reduced- dimension subspace. The subspace is defined by the leading empirical orthogonal functions calculated from an unconstrained model calculation. Both an identical twin experiment using simulated SHA observations and assimilation of the real TOPEX data were performed. Results from the twin experiments demonstrate the ability of the method to constrain the ocean circulation and the subsurface temperature structure. The impact on the subsurface temperature structure of TOPEX assimilation was assessed using data from expandable bathythermographs. This showed a substantial improvement in the estimated temperature variability only within 13 degrees in latitude around the equator. The impact of TOPEX SHA assimilation on zonally integrated meridional transport across different latitudes was also estimated. Again within 13 degrees from the equator both the mean amplitude and interannual variability of the surface and subsurface transports were significantly enhanced, while the transports were insensitive to the assimilation in the subtropics.

Malannotte-Rizzoli, Paola↗