Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A new method for the identification of non-Gaussian line profiles in elliptical galaxies

A new parameterization for the line profiles of elliptical galaxies, the Gauss-Hermite series, is proposed. This approach expands the line profile as a sum of orthogonal functions which minimizes the correlations between the errors in the parameters of the fit. This method also make use of the fact that Gaussians provide good low-order fits to observed line profiles. The method yields measurements of the line strength, mean radial velocity, and the velocity dispersion as well as two extra parameters, h3 and h4, that measure asymmetric and symmetric deviations of the line profiles from a Gaussian, respectively. The new method was used to derive profiles for three elliptical galaxies which all have asymmetric line profiles on the major axis with symmetric deviations from a Gaussian. Results confirm that elliptical galaxies have complex structures due to their complex formation history.

Van Der Marel, Roeland P.↗
NASA NTRS1993

Low-frequency synoptic-eddy activity in the Pacific storm track

The GCM presently used to examine the nature of LF variations in synoptic-eddy activity over the North Pacific, when compared with observations, yields realistic time-mean and LF synoptic-eddy forcing of the 200-mb zonal wind. In the time mean, this forcing exhibits an east-west dipole structure that is consistent with the general view of the life cycle of baroclinic waves; these show marked upward and eastward propagation in the western and central Pacific, as well as meridional propagation to the east. Composites of the hemispheric distribution of synoptic-eddy forcing in the GCM, based on the extremes of the dominant Pacific empirical orthogonal function, show a seesaw behavior with enhanced eddy forcing in the North Pacific basin which is associated with suppressed forcing in the North Atlantic basin (and vice versa).

Higgins, R. W.↗
NASA NTRS1993

Variability in the Gulf of Alaska from Geosat altimetry data

Satellite altimetry was used to examine annual and interannual variability in the Gulf of Alaska region. Crossover data from the Geosat Geodetic Mission (GM) and collinear data from the Exact Repeat Mission (ERM) were processed separately to form sea surface height anomalies at grid points. The time series from the GM and ERM were then combined to produce a 3.75-year data set. The time series from the ERM data set agreed fairly well with hydrographic dynamic heights at several locations, with an average correlation of 0.70 between the two data sets. The combined and ERM altimetric data sets were analyzed using empirical orthogonal functions (EOFs). These revealed variability that occurs primarily on annual and interannual time scales. A comparison with EOF analysis of the atmospheric pressure field during the same time periods showed that the annual variation in pressure seemed to be reflected in both the combined and ERM altimetric data sets. The amplitude time series of the first mode in the combined data set was very similar to the North Pacific pressure index during the 1985-1989 time frame. The maximum correlation was at a lag of 250 days. Finally, an interannual mode was present in all three data sets which was closely linked to the baroclinic variations measured by the hydrographic data.

Bhaskaran, Shyam↗
NASA NTRS1994

A comparison between general circulation model simulations using two sea surface temperature datasets for January 1979

Simulations with the UCLA atmospheric general circulation model (AGCM) using two different global sea surface temperature (SST) datasets for January 1979 are compared. One of these datasets is based on Comprehensive Ocean-Atmosphere Data Set (COADS) (SSTs) at locations where there are ship reports, and climatology elsewhere; the other is derived from measurements by instruments onboard NOAA satellites. In the former dataset (COADS SST), data are concentrated along shipping routes in the Northern Hemisphere; in the latter dataset High Resolution Infrared Sounder (HIRS SST), data cover the global domain. Ensembles of five 30-day mean fields are obtained from integrations performed in the perpetual-January mode. The results are presented as anomalies, that is, departures of each ensemble mean from that produced in a control simulation with climatological SSTs. Large differences are found between the anomalies obtained using COADS and HIRS SSTs, even in the Northern Hemisphere where the datasets are most similar to each other. The internal variability of the circulation in the control simulation and the simulated atmospheric response to anomalous forcings appear to be linked in that the pattern of geopotential height anomalies obtained using COADS SSTs resembles the first empirical orthogonal function (EOF 1) in the control simulation. The corresponding pattern obtained using HIRS SSTs is substantially different and somewhat resembles EOF 2 in the sector from central North America to central Asia. To gain insight into the reasons for these results, three additional simulations are carried out with SST anomalies confined to regions where COADS SSTs are substantially warmer than HIRS SSTs. The regions correspond to warm pools in the northwest and northeast Pacific, and the northwest Atlantic. These warm pools tend to produce positive geopotential height anomalies in the northeastern part of the corresponding oceans. Both warm pools in the Pacific produce large-scale circulation anomalies with a pattern that resembles that obtained using COADS SSTs as well as EOF 1 of the control simulation; the warm pool in the Atlantic does not. These results suggest that the differences obtained with COADS SSTs and HIRS SSTs are mostly due to the differences in the datasets over the northern Pacific. There was a blocking episode near Greenland in late January 1979. Both simulations with warm SST anomalies over the northwest and northeast Pacific show a tendency toward increased incidence of North Atlantic blocking; the simulation with warm SST anomalies over the northwest Atlantic shows a tendency toward decreased incidence. These results suggest that features in both SST datasets that do not have a counterpart in the other dataset contribute signficantly to the differences between the simulated and observed fields. The results of this study imply that uncertainties in current SST distributions for the world oceans can be as important as the SST anomalies themselves in terms of their impact on the atmospheric circulation. Caution should be exercised, therefore, when linking anomalous circulation and SST patterns, especially in long-range prediction.

Ose, Tomoaki↗
NASA NTRS1994

Ozone and stratospheric height waves for opposite phases of the QBO

The stratospheric quasi-biennial oscillation (QBO) provides an important source of interannual variations in the Northern Hemisphere. O'sullivan and Salby (1990) related extra-tropical eddy transport with the phase of the tropical QBO. When the tropical wind is easterly, the zero wind line is shifted into the winter hemisphere. Enhanced wave activity in middle latitudes acts to weaken the polar vortex. When the tropical wind is in the westerly phase the situation reverses. Heights at 30 mb and ozone configurations are contrasted in this paper for these two QBO phases. When the winter vortex deforms due to the amplification of planetary waves 1 and 2, extends westward and equatorward, the complementary band of low vorticity air spirals in toward the pole from lower latitudes. Sometimes, these planetary waves break (Juckes and McIntyre, 1987) and an irreversible mixing of air takes place between high and mid-latitudes. Global ozone patterns, as obtained form satellite observations, appear to be affected by planetary wave breaking (Leovy et al. 1985). This mixing results on regions with uniform ozone and potential vorticity. In the Southern Hemisphere (SH), Newman and Randel (1988) using Total Ozone Mapping Spectrometer (TOMS) data and the NMC analyses have found strong spatial correlation between the October mean temperature in the lower stratosphere and total ozone for the 1979 through 1986 years. Recently Nogues-Paegle et al.(1992) analyzed SH ozone and height data from 1986 to 1989. They found that leading empirical orthogonal functions (EOFs) for both ozone and 50 mb heights exhibit zonal wave 1 and 2 and that the correlations between ozone and 50 mb principal components (PCs) are high. The results were found to be consistent with a linear planetary wave advecting a passive tracer. In this paper, the dominant patterns of variability for 30 mb NMC heights and TOMS total ozone are obtained for the winter to summer transition (January to May) in the Northern Hemisphere (NH) for the years 1987-1990.

Mo, Kingtse C.↗
NASA NTRS1993

Nearly steady flows in GONG prototype data

Doppler velocity images obtained with the GONG prototype instrument were analyzed to measure the nearly steady photospheric flows. The data consists of 88 images each of velocity, intensity, and modulation obtained at 20:00 UT on 88 days from July 1992 to February 1994. Each velocity image was temporally filtered to remove the p-mode oscillations, masked to exclude active regions, and then analyzed using spherical harmonics and orthogonal functions as described by Hathaway (1992). The spectral coefficients show very consistent results for the entire time interval with some evidence of year-to-year variations. The rotation profile agrees well with previous results and exhibits a north-south asymmetry that reverses sign during the 20 month interval. The residual rotation velocities exhibit structures with amplitudes of approximately 5 m/s that may be related to torsional oscillations. The meridional circulation is directed from the equator toward the poles with a peak velocity in the photosphere of approximately 50 m/s. The higher order components are very weak but indicate a divergent flow from the mid-latitudes (opposite that found for the June 1989 data). The convective limb shift is well fit by a 3rd order polynomial. The convection spectrum has a prominent peak at spherical harmonic degrees of l approximately 150 with very little signal in the low degree modes. Analysis of this signal shows that there is no evidence for giant cell convection at the level of approximately 10 m/s for all modes up to l = 32.

Hathaway, David H.↗
NASA NTRS1994

Principal modes of atmospheric circulation anomalies associated with global angular momentum fluctuations

This paper provides a description of the variability of global atmospheric angular momentum (GAM) and its relationship with principal modes of three-dimensional atmospheric circulation anomalies. The data used are 5-day mean global wind fields from the European Centre for Medium-Range Weather Forecasts initialized dataset for 1980-1989. Significant seasonal variation of GAM is observed with maxima in April and November and a minimum during late July. The amplitude of the annual cycle is largest in the upper troposphere and decreases toward the surface. Although the lower tropospheric contribution to the total angular momentum is relatively small, its annual cycle is out of phase with those of the upper atmosphere and GAM. Also identified is a distinct semiannual component, with double peaks appearing in April and November. This signal is most noticeable in the upper troposphere above the 300-mb level. The principal modes of zonal-mean angular momentum and meridional circulation anomalies and their coupled modes are obtained by using empirical orthogonal function analysis and singular value decomposition. It is shown that the leading modes of the angular momentum and meridional circulation are coupled with each other and are responsible for much of the variability in GAM. The coupled modes represent fluctuations of upper-level subtropical zonal flow, which are linked to the modulation of Hadley circulation intensity in both hemispheres. It is found that GAM is highly correlated with the first eigenvector of upper-level streamfunction anomalies, which consists of a superrotational flow in the tropics and subtropics, except over the central Pacific where a 'blocked' flow with two subtropical anticyclonic circulation cells straddling the equator is found. Much of the blocked flow is due to the establishment of dipole anomalies in the velocity potential with centers over the central Pacific and the Maritime Continent on the interannual time scale. On the intraseasonal time scale, GAM fluctuation is dominated by superrotational flow in the tropics, with the blocked flow present to a much lesser extent. The associated velocity potential anomaly has a weak dipole structure with centers over the Indian Ocean and the eastern Pacific. The implications of the above results on the total angular momentum balance of the earth-atmosphere system are also discussed.

Kang, In-Sik↗
NASA NTRS1994

Inference of subsurface thermohaline structure from fields measurable by satellite

Satellites now provide global measurements of the ocean's surface height and temperature. Ocean climatologies for the northwest Pacific and northwest Atlantic Oceans that relate sea surface height, sea surface temperature, day of the year, latitude, and longitude to temperature and salinity profiles were produced using least-squares regression. These analyses use over 33,000 profiles of historical temperature and salinity data and are considerably streamlined and compacted by expressing each profile in terms of empirical orthogonal functions. Evaluations and error analyses of the climatologies, including a comparison to the navy's Generalized Digital Environmental Model, were performed and differences between the regions are discussed. Two sample vertical sections are shown to be closely reproduced with the climatologies. Climatologies based on surface height and temperature are found to offer considerable improvement over climatologies based only on position.

Carnes, Michael R.↗
NASA NTRS1994

Southern hemisphere low level wind circulation statistics from the Seasat scatterometer

Analyses of remotely sensed low-level wind vector data over the Southern Ocean are performed. Five-day averages and monthly means are created and the month-to-month variability during the winter (July-September) of 1978 is investigated. The remotely sensed winds are compared to the Australian Bureau of Meteorology (ABM) and the National Meteorological Center (NMC) surface analyses. In southern latitudes the remotely sensed winds are stronger than what the weather services' analyses suggest, indicating under-estimation by ABM and NMC in these regions. The evolution of the low-level jet and the major stormtracks during the season are studied and different flow regimes are identified. The large-scale variability of the meridional flow is studied with the aid of empirical orthogonal function (EOF) analysis. The dominance of quasi-stationary wave numbers 3,4, and 5 in the winter flows is evident in both the EOF analysis and the mean flow. The signature of an exceptionally strong blocking situation is evident in July and the special conditions leading to it are discussed. A very large intraseasonal variability with different flow regimes at different months is documented.

Levy, Gad↗
NASA NTRS1995

Decadal variability of the Tropical Atlantic Ocean Surface Temperature in shipboard measurements and in a Global Ocean-Atmosphere model

Sea surface temperature (SST) variability was investigated in a 200-yr integration of a global model of the coupled oceanic and atmospheric general circulations developed at the Geophysical Fluid Dynamics Laboratory (GFDL). The second 100 yr of SST in the coupled model's tropical Atlantic region were analyzed with a variety of techniques. Analyses of SST time series, averaged over approximately the same subregions as the Global Ocean Surface Temperature Atlas (GOSTA) time series, showed that the GFDL SST anomalies also undergo pronounced quasi-oscillatory decadal and multidecadal variability but at somewhat shorter timescales than the GOSTA SST anomalies. Further analyses of the horizontal structures of the decadal timescale variability in the GFDL coupled model showed the existence of two types of variability in general agreement with results of the GOSTA SST time series analyses. One type, characterized by timescales between 8 and 11 yr, has high spatial coherence within each hemisphere but not between the two hemispheres of the tropical Atlantic. A second type, characterized by timescales between 12 and 20 yr, has high spatial coherence between the two hemispheres. The second type of variability is considerably weaker than the first. As in the GOSTA time series, the multidecadal variability in the GFDL SST time series has approximately opposite phases between the tropical North and South Atlantic Oceans. Empirical orthogonal function analyses of the tropical Atlantic SST anomalies revealed a north-south bipolar pattern as the dominant pattern of decadal variability. It is suggested that the bipolar pattern can be interpreted as decadal variability of the interhemispheric gradient of SST anomalies. The decadal and multidecadal timescale variability of the tropical Atlantic SST, both in the actual and in the GFDL model, stands out significantly above the background 'red noise' and is coherent within each of the time series, suggesting that specific sets of processes may be responsible for the choice of the decadal and multidecadal timescales. Finally, it must be emphasized that the GFDL coupled ocean-atmosphere model generates the decadal and multidecadal timescale variability without any externally applied force, solar or lunar, at those timescales.

Mehta, Vikram M.↗
NASA NTRS1994

Variability in pigment concentration in warm-core rings as determined by coastal zone color scanner satellite imagery from the Mid-Atlantic Bight

A time series of coastal zone color scanner (CZCS) derived chlorophyll (CZCS-chl) and sea surface temperature (SST) satellite imagery was developed for the Mid-Atlantic Bight (MAB). Warm-core rings (WCR) were identified by both the warmer SST signal as well as the low pigment concentrations of their cores. The variation in pigment concentrations and SST observed in satellite imagery over the geographic range and life span of four WCRs is investigated. The hypotheses are that pigment concentration increase during the lifetime of the WCR is a response to processes such as convective overturn, upwelling, edge enhancement due to increased vertical mixing, active convergence, or lateral exchange. Empirical orthogonal function analysis (EOF) is used to investigate the relationship between SST and pigment patterns observed in the presence of a WCR. The first two EOF modes explain more than 80% of the variability observed in all four WCRs and in both (SST and pigment) data sets. The results of this study show that, at the synoptic scales of staellite data, the variability observed in the WCRs is greater at the periphery of the rings. These results show that advective entrainment, rather than processes at ring center (e.g., shoaling of the pycnocline/nutricline in response to frictional decay) or at the periphery due to other processes such as vertical mixing, is the mechanism responsible for the observed variability.

Garcia-Moliner, Graciela↗
NASA NTRS1993

Observations on the long-period variability of the Gulf Stream downstream of Cape Hatteras

To examine the long-period variability of the Gulf Stream, sea level residuals relative to a 2-year mean sea level in the Gulf Stream downstream of Cape Hatteras (between 75 deg W and 60 deg W longitude) are used. Residuals, as derived from Geosat altimetry between November 1986 and December 1988, were gridded in space and time at a temporal resolution of 10 days and spatial resolution of 1/4 deg. Complex empirical orthogonal function (CEOF) analysis was applied to the data set to extract the spatially correlated signal with the original data subsampled to 1/2 deg. In addition to determining the space-time scales and propagation characterisitics of the different modes, wavenumber-frequency spectral techniques were used to separate the variability into propagating and stationary components. The CEOF technique applied to the data set indicated that the first four CEOF modes accounted for 60% of the variability and were found to be above the noise leve 99% of the time. CEOF 1 was associated with westward propagation at 5 km/d at a wavelength of 2000 km and eastward propagation at 1-2 km/d centered at a 500-km wavelength. This first CEOF is in good agreement with thin-jet equivalent barotropic models which predict westward propagation for wavelengths greater than 1130 km. A deflection of the wavelike pattern at 65 deg W also indicates a possible topographic effect. A simple scaling of the effect of topography indicates that for length scales longer than the internal Rossby radius of deformation, the topographic term is at least of the same order of magnitude as the beta effect. The second CEOF was more broadbanded in wavenumber space, with eastward propagation occurring in a wavenumber-frequency band between 300 and 1400 km and 0.5 and 2.0 cycles/yr. The third CEOF is similar in structure to the first, but with less energy. CEOF 4 was clearly identifiable with higher frequencies than the first three with westward propagation at 4 km/d. The spatial location of this mode along with the westward propagation indicates possible influences from eddy-stream interactions. Thus topography, Rossby wave dynamics and eddy-stream interactions all appear to have a significant role in determining the space-time scales and propagation properties of the long-period response of sea level in the Gulf Stream.

Vazquez, Jorge↗
NASA NTRS1995

Temporal variation of meandering intensity and domain-wide lateral oscillations of the Gulf Stream

The path of the Gulf Stream exhibits two modes of variability: wavelike spatial meanders associated with instability processes and large-sale lateral shifts of the path presumably due to atmospheric forcing. The objectives of this study are to examine the temporal variation of the intensity of spatial meandering in the stream, to characterize large-scale lateral oscillations in the stream's path, and to study the correlation betwen these two dynamically distinct modes of variability. The data used for this analysis are path displacemets ofthe Gulf Stream between 75 deg and 60 deg W obtained from AVHRR-derived (Advanced Very High Resolution Radiometer) infrared images for the period April 1982 through December 1989. Meandering intensity, measured by the spatial root-mean-sqaure displacement of the stream path, displays a 9-month dominant periodicity which is persistent through the study period. The 9-month fluctuation in meandering intensity may be related to the interaction of Rosseby waves with the stream. Interannual variation of meandering intensity is also found to be significant, with meandering being mich more intense during 1985 than it was in 1987. Annual variation, however,is weak and not well-defined.The spatially averaged position of the stream, which reflects nonmeandering large-scale lateral oscillations of the stream path, is dominated by an annual cycle. On average, the mean position is farthest north in November and farthest south in April. The first empirical orthogonal function mode of the space-time path displacements represents lateral oscillatins that are in-phase over the space-time domain. Interannual oscillations are also observed and are found to be weaker than the annual oscillation. The eigenvalue of the first mode indicates that about 21.5% of the total space-time variability of the stream path can be attibuted to domain-wide lateral oscillation. The correlation between meandering intensity and domain-wide lateral oscillations is very weak.

Lee, Tong↗
NASA NTRS1995

Dynamics of the Antarctic Circumpolar Current. Evidence for Topographic Effects from Altimeter Data and Numerical Model Output

Geosat altimeter data and numerical model output are used to examine the circulation and dynamics of the Antarctic Circumpolar Current (ACC). The mean sea surface height across the ACC has been reconstructed from height variability measured by the altimeter, without assuming prior knowledge of the geoid. The results indicate locations for the Subantarctic and Polar Fronts which are consistent with in situ observations and indicate that the fronts are substantially steered by bathymetry. Detailed examination of spatial and temporal variability indicates a spatial decorrelation scale of 85 km and a temporal e-folding scale of 34 days. Empirical Orthogonal Function analysis suggests that the scales of motion are relatively short, occuring on 1000 km length-scales rather than basin or global scales. The momentum balance of the ACC has been investigated using output from the high resolution primitive equation model in combination with altimeter data. In the Semtner-Chervin quarter-degree general circulation model topographic form stress is the dominant process balancing the surface wind forcing. In stream coordinates, the dominant effect transporting momentum across the ACC is bibarmonic friction. Potential vorticity is considered on Montgomery streamlines in the model output and along surface streamlines in model and altimeter data. (AN)

Gille, Sarah T.↗
NASA NTRS1993

Interannual Sea Level Variations in the Tropical Indian Ocean from Geosat and Shallow Water Simulations

Sea level variations of the Indian Ocean north of 20 deg S are analyzed from Geosat satellite altimeter data over April 1985-September 1989. These variations are compared and interpreted with numerical simulations derived from a reduced gravity model forced by FSU observed winds over the same period. After decomposition into complex empirical orthogonal functions, the low-frequency anomalies are described by the first two modes for observations as well as for simulations. The sums of the two modes contain 34% and 40% of the observed and simulated variances, respectively. Averaged over the basin, the observed and simulated sea level changes are correlated by 0.92 over 1985-1988. The strongest change happens during the El Ninio 1986-1987: between winter 1986 and summer 1987 the basin-averaged sea level rises by approx. 1 cm. These low-frequency variations can partly be explained by changes in the Sverdrup circulation. The southern tropical Indian Ocean between 1O deg and 20 deg S is the domain where those changes are strongest: the averaged sea level rises by approx. 4.5 cm between winter 1986 and winter 1987. There, the signal propagates southwestward across the basin at a speed similar to free Rossby waves. Sensitivity of observed anomalies is examined over 1987-1988, with different orbit ephemeris, tropospheric corrections, and error reduction processes. The uncertainty of the basin-averaged sea level estimates is mostly due to the way the orbit error is reduced and reaches approx. 1 cm. Nonetheless, spatial correlation is good between the various observations and better than between observations and simulations. Sensitivity of simulated anomalies to the wind uncertainty, examined with Former Soviet Union (FSU) and European Center for Medium-Range Weather Forecasting (ECMWF) forcings over 1985-1988, shows that the variance of the simulations driven by ECMWF is 52% smaller, as FSU winds are stronger than ECMWF. Results show that the wind strength also affects the dynamic response of the ocean: anomalies propagate westward across the basin more than twice as fast with FSU than with ECMWF. It is found that the discrepancy is larger between ECMWF and FSU simulations than between observations and FSU simulations.

Perigaud, Claire↗
NASA NTRS1999

The Prospect for Remote Sensing of Cirrus Clouds with a Submillimeter-Wave Spectrometer

Given the substantial radiative effects of cirrus clouds and the need to validate cirrus cloud mass in climate models, it is important to measure the global distribution of cirrus properties with satellite remote sensing. Existing cirrus remote sensing techniques, such as solar reflectance methods, measure cirrus ice water path (IWP) rather indirectly and with limited accuracy. Submillimeter/wave radiometry is an independent method of cirrus remote sensing based on ice particles scattering the upwelling radiance emitted by the lower atmosphere. A new aircraft instrument, the Far Infrared Sensor for Cirrus (FIRSC), is described. The FIRSC employs a Fourier Transform Spectrometer (FTS). which measures the upwelling radiance across the whole submillimeter region (0.1 1.0-mm wavelength). This wide spectral coverage gives high sensitivity to most cirrus particle sizes and allows accurate determination of the characteristic particle size. Radiative transfer modeling is performed to analyze the capabilities of the submillimeter FTS technique. A linear inversion analysis is done to show that cirrus IWP, particle size, and upper-tropospheric temperature and water vapor may be accurately measured, A nonlinear statistical algorithm is developed using a database of 20000 spectra simulated by randomly varying most relevant cirrus and atmospheric parameters. An empirical orthogonal function analysis reduces the 500-point spectrum (20 - 70/cm) to 15 "pseudo-channels" that are then input to a neural network to retrieve cirrus IWP and median particle diameter. A Monte Carlo accuracy study is performed with simulated spectra having realistic noise. The retrieval errors are low for IWP (rms less than a factor of 1.5) and for particle sizes (rins less than 30%) for IWP greater than 5 g/sq m and a wide range of median particle sizes. This detailed modeling indicates that there is good potential to accurately measure cirrus properties with a submillimeter FTS.

Evans, K. Franklin↗
NASA NTRS1999

Reconstruction of the 1997/1998 El Nino from TOPEX/POSEIDON and TOGA/TAO Data Using a Massively Parallel Pacific-Ocean Model and Ensemble Kalman Filter

Two massively parallel data assimilation systems in which the model forecast-error covariances are estimated from the distribution of an ensemble of model integrations are applied to the assimilation of 97-98 TOPEX/POSEIDON altimetry and TOGA/TAO temperature data into a Pacific basin version the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. in the first system, ensemble of model runs forced by an ensemble of atmospheric model simulations is used to calculate asymptotic error statistics. The data assimilation then occurs in the reduced phase space spanned by the corresponding leading empirical orthogonal functions. The second system is an ensemble Kalman filter in which new error statistics are computed during each assimilation cycle from the time-dependent ensemble distribution. The data assimilation experiments are conducted on NSIPP's 512-processor CRAY T3E. The two data assimilation systems are validated by withholding part of the data and quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The pros and cons of each system are discussed.

Keppenne, C. L.↗
NASA NTRS1999

Estimation of Non-Seasonal Oceanic Mass Redistribution and its Geodynamic Effects from the Topex/Poseidon Observation

Mass redistribution in terms of angular momentum will change the Earth's rotation under the conservation of angular momentum. It also changes the external gravitational field according to Newton's gravitational law. These small geodynamic variations can be measured very accurately by modem space geodetic techniques. We report on new estimates of the oceanic tidal angular momentum and hence the tidal variations in the Earth's rotation rate and polar motion. These are deduced from the new tidal solutions GOT99 and TPXO.3, plus M2 and Mf estimates (with error bars) from TPXO.4a. Comparisons with observed diurnal/semidiurnal tidal variations are generally similar to previous model estimates. Comparisons of models and observations of Mf polar motion are confused, with neither models nor observations agreeing among themselves or with each other. Part of this is caused by differences in Mf currents; model TPXO.4 shows current-intensive Rossby modes along western boundaries of major basins, which are indeed expected on theoretical grounds. Unfortunately, there are no reliable in situ observations of Mf currents useful for comparing models. On non-tidal oceanic mass redistribution we compute the principal components of the multi-year T/P sea surface height (SSH) field using EOF (empirical orthogonal functions) taking into consideration of the latitude-dependent area weighting. Emphasis is placed upon non-seasonal and interannual oscillations on regional scales, including ENSO and the North Atlantic Oscillation. Steric corrections to the observed SSH due to sea surface temperature changes (which has no geodynamic effects) is now under study using POCM ocean model output.

Chao, Benjamin F.↗