Engineering PapersSearch

Engineering topics

Au, A. Y.

Publications and source records attributed to Au, A. Y..

Global Biomass Variation and its Geodynamic Effects, 1982-1998

Redistribution of mass near Earth's surface alters its rotation, gravity field, and geocenter location. Advanced techniques for measuring these geodetic variations now exist, but the ability to attribute the observed modes to individual Earth system processes has been hampered by a shortage of reliable global data on such processes, especially hydrospheric processes. To address one aspect of this deficiency, 17 yrs of monthly, global maps of vegetation biomass were produced by applying field-based relationships to satellite-derived vegetation type and leaf area index. The seasonal variability of biomass was estimated to be as large as 5 kg m(exp -2). Of this amount, approximately 4 kg m(exp -2) is due to vegetation water storage variations. The time series of maps was used to compute geodetic anomalies, which were then compared with existing geodetic observations as well as the estimated measurement sensitivity of the Gravity Recovery and Climate Experiment (GRACE). For gravity, the seasonal amplitude of biomass variations may be just within GRACE'S limits of detectability, but it is still an order of magnitude smaller than current observation uncertainty using the satellite-laser-ranging technique. The contribution of total biomass variations to seasonal polar motion amplitude is detectable in today's measurement, but it is obscured by contributions from various other sources, some of which are two orders of magnitude larger. The influence on the length of day is below current limits of detectability. Although the nonseasonal geodynamic signals show clear interannual variability, they are too small to be detected.

Rodell, M.

Seasonal and Non-seasonal Sea Level Variations by Exchange of Water with Land Hydrology

The global ocean exchanges a large amount of water, seasonally or non-seasonally, with land hydrology. Apart from the long-term melting of ice sheets and glaciers, the water is exchanged directly as land runoff R, and indirectly via atmosphere in the form of precipitation minus evapo-transpiration P-E. On land, the hydrological budget balance is soil moisture S = P-E-R. The runoff R has been difficult to monitor; but now by combining the following two data sets one can obtain a global estimate, subject to the spatial and temporal resolutions afforded by the data: (1) The space gravity mission GRACE yields monthly S estimate on a spatial scale larger than approx. 1000 km over the last 2.5 years; (2) The atmospheric circulation model output, such as from NCEP, provides proxy estimates for P-E at monthly and approx. 200 km resolutions. We will discuss these estimates and the effects on the global ocean water budget and hence sea level.

Chao, Benjamin F.

Time-Variable Gravity Signal due to Extratropic Pacific Water Mass Redistribution

Using the satellite-laser-ranging (SLR) data, Cox and Chao [2002] reported the detection of a large post-1998 anomaly (in the form of a positive jump) in the time series of Earth s lowest-degree gravity harmonic 52, or the dynamic oblateness. Among several groups now examining the mass redistribution in the global geophysical fluids in search of the cause(s), we report here a temporally coinciding anomalies found in the extratropic north + south Pacific basins. Clearly seen in the leading EOFPC mode for extratropic Pacific, these anomalies occurred in sea-surface height, sea-surface temperature, and temperature- and salinity-depth profiles. We based our analysis on two different data sources: TOPEX/Poseidon altimetry, and the ECCO ocean general circulation model output assimilating T/P data. The magnitude of these changes, when converted to equivalent J2 change, appears to be a few times too small to explain the observed J2 directly. These findings, and the fact that the anomalies occurred following the strong 1997-98 El Nino, suggest strong geophysical connection of the interannual-to-decadal variation of 52 with the Pacific Decadal Oscillation (PDO) and the ultimate global-change processes that cause PDO. More work is underway, and additional independent data sources are examined, paying close attention to the fact that the J2 anomaly has been reversing back to normal since 2001. These include: (1) cryospheric contributions (melting of glaciers and ice sheets); (2) land hydrological contributions; (3) polar sea influences ( e g , via deep flow); (4) fluid flow in Earth's core; (5) time-variable gravity signals from SLR in higher harmonic degree/order, including J3,J4, (2,1), and (2,2) coefficients, considering their lower signal-to-noise ratios; (6) Earth rotation data in terms of length-of-day and polar motion.

Chao, B. F.

Time-Variable Gravity Signal Due to Extratropic Pacific Water Mass Redistribution

Cox and Chao [2002] reported the detection of a large anomaly in the form of a positive "jump" in the time series of Earth's lowest-degree gravity harmonic J2, or the dynamic oblateness, during 1998. This prompted us to examine the mass redistribution in the global oceans. We report here a seesaw of the sea-surface height (SSH) in the extratropic north + south Pacific basins -- the leading (nonseasonal) EOF/PC mode in SSH derived from the 10-year TOPEX/Poseidon altimetry data in the extratropic Pacific region. The mode underwent a step-like jump with time evolution that match remarkably well with the observed J2 anomaly. However, the magnitude is several times too small to explain the observed J2, even if assuming the SSH jump was all mass-induced (as opposed to any steric effect which causes no time-variable gravity signal). If one accepts the notion that this extratropic Pacific seesaw is part of the geophysical process that produced the observed 1998 J2 anomaly, then this finding suggests strong geophysical connection of the interannual-to-decadal variation of J2 with the Pacific Decadal Oscillation (PDO), as the time series of the above EOF/PC mode is actually a formally defined PDO Index series.

Chao, Benjamin F.

Gravitational Signal of Mass Redistribution Due to Interannual Meteorological Oscillations in Atmosphere and Ocean

Interannual meteorological oscillations (ENSO, QBO, NAO, etc.) have demonstrable influences on Earth's rotation. Here we study their effects on global gravitational field, whose temporal variations are being studied using SLR (satellite laser ranging) data and in anticipation of the new space mission GRACE. The meteorological oscillation modes are identified using the EOF (empirical orthogonal function)/PC (principal component) decomposition of surface fields (in which we take care of issues associated with the area-weighting and non-zero mean). We examine two fields, one for the global surface pressure field for the atmosphere obtained from the NCEP reanalysis (for the past 40 years), one for the surface topography field for the ocean from the Topex/Poseidon (T/P) data (for the past 8 years). We use monthly maps, and remove the mean-monthly ("climatology") values from each grid point, hence focusing only on non-seasonal signals. The T/P data were first subject to a steric correction where the steric contribution to the ocean surface topography was removed according to output of the numerical POCM model. The respective atmospheric and oceanic contributions to the gravitational variation, in terms of harmonic Stokes coefficients, are then combined mode-by-mode. Since the T/P data already contain the oceanic response to overlying atmospheric pressure, no regards to the inverted-barometer behavior for the ocean need be considered. Results for the lowest-degree Stokes coefficients can then be compared with space geodetic observations including the Earth's rotation and the SLR data mentioned above, to identify the importance of each meteorological oscillations in gravitational variation signals.

Chao, B. F.

Interannual Atmospheric Oscillations and Their Gravitational Effects

Past studies have examined the effects of the interannual atmospheric oscillations (IAO), or often called teleconnection patterns, on the Earth's rotation, such as El Nino/Southern Oscillation (ENSO), the quasi-biennial oscillation, and to a lesser extent the North Atlantic Oscillation (NAO). The present study focuses on IAO's effects on the gravitational field and geocenter; the purpose is to be able to identify and isolate the contributions of each prominent IAO in relation to the total contribution of the atmosphere, in terms of their magnitudes, their geographical patterns, and their interannual time history. We use the 40-year NCEP reanalysis of the monthly, global atmospheric surface pressure field as our basic data set. The method we apply to isolate the IAOs is the empirical orthogonal function (EOF) decomposition which is widely used in meteorological investigations. We do the EOF analysis first on IAOs' seasonal signals (by "collapsing" the 40-year series into 12 mean-months for each grid point) and obtain estimates for their respective contributions. Then we remove these seasonal signals from the data to focus on the (broad-band) interannual EOFs. We examine ENSO, NAO, North Pacific Oscillation, and other less prominent IAOs that can be identified from our data set, and compute their respective contributions to the variation of global gravitational field and geocenter motion. Finally, we compare the results with the available observational data, and discuss the implications w.r.t. the upcoming space missions such as GRACE.

Chao, Benjamin F.

Non-Tidal Non-Seasonal Oceanic Mass Redistribution Estimated from the TOPEX/Poseidon Observation

Topex/Poseidon altimetry data are used to estimate the non-tidal mass redistribution as a function of space-time. The goal is to study the contribution of ocean circulations in the geodynamic effects including Earth's rotational and gravitational variations. We examine the non-seasonal anomalies at monthly sampling rate over the T/P span of eight years, concentrating especially on interannual variabilities. Since the sea-surface height data obtained from altimetry is the combined effect of steric change (primarily thermal effect) and the mass flux, and because the former has zero contribution to the geodynamic effects, one needs to do a so-called steric correction by removing from the altimetry data the steric contributions. We achieve it using multiyear monthly sea-surface temperature maps, together with monthly "climatology" for mixed-layer depth maps for mean-months of the year. We analyze both sets of the altimetry data and the steric correction data using the empirical orthogonal function/principal component analysis (in which we take care of issues associated with the area-weighting and non-zero mean), and examine the most important modes, either globally or regionally. In particular, the ENSO in the tropical Pacific/Indian Oceans exhibits the most prominent pattern. The net mass transport after the steric correction can then be compared with: (1) ocean general circulation model outputs for the same period of time (such as from POCM-4B); (2) non-atmospheric Earth rotation variations obtained from space geodesy data and atmospheric angular momentum data; (3) non-atmospheric low-degree gravitational variations from satellite-laser-ranging observations and global atmospheric data.

Chao, Benjamin F.

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.

New Estimates of Hydrological and Oceanic Excitations of Variations of Earth's Rotation, Geocenter and Gravitational Field

Hydrological mass transport in the geophysical fluids of the atmosphere-hydrosphere-solid Earth surface system can excite Earth's rotational variations in both length-of-day and polar motion. These effects can be computed in terms of the hydrological angular momentum by proper integration of global meteorological data. We do so using the 40-year NCEP data and the 18-year NASA GEOS-1 data, where the precipitation and evapotranspiration budgets are computed via the water mass balance of the atmosphere based on Oki et al.'s (1995) algorithm. This hydrological mass redistribution will also cause geocenter motion and changes in Earth's gravitational field, which are similarly computed using the same data sets. Corresponding geodynamic effects due to the oceanic mass transports (i.e. oceanic angular momentum and ocean-induced geocenter/gravity changes) have also been computed in a similar manner. We here compare two independent sets of the result from: (1) non-steric ocean surface topography observations based on Topex/Poseidon, and (2) the model output of the mass field by the Parallel Ocean Climate Model. Finally, the hydrological and the oceanic time series are combined in an effort to better explain the observed non-atmospheric effects. The latter are obtained by subtracting the atmospheric angular momentum from Earth rotation observations, and the atmosphere- induced geocenter/gravity effects from corresponding geodetic observations, both using the above-mentioned atmospheric data sets.

Chao, Benjamin F.

Temporal variations of the earth's gravitational field from satellite laser ranging to LAGEOS

Monthly values of the J2 and J3 earth gravitational coefficients were estimated using LAGEOS satellite laser ranging data collected between 1980 and 1989. Monthly variations in gravitational coefficients caused by atmospheric mass redistribution were calculated using measurements of variations in surface atmospheric pressure. Results for correlation studies of the two time series are presented. The LAGEOS and atmospheric J2 time series agree well and it appears that variations in J2 can be attributed to the redistribution of atmospheric mass. Atmospheric and LAGEOS estimates for J3 show poorer agreement, J3 estimates appear to be very sensitive to unmodeled forces acting on the satellite. Results indicate that the LAGEOS data can be used to detect small variations in the gravitational field.

Nerem, R. S.

Analysis of altimetry over inland seas

Satellite-based altimetric data taken by GEOS-3 and SEASAT over the Black Sea and Caspian Sea are analyzed and a least squares collocation technique is used to predict the geoid undulation on a .25-degree by .25-degree grid and to transform these geoid undulations to free air gravity anomalies. This project entailed processing satellite altimeter data over inland seas for recovery of area mean gravity information. Gravity information in this area of the world is not readily available, so the possibility of obtaining it from the processing of altimeter observations is attractive. The principal objective was to complete and extend analyses done in a previous study, verify those results, and document the results and techniques. A secondary objective was to improve the algorithms and results, if possible. The approach used involved editing geoid height data to remove overland data; evaluating geoid height differences at crossover points; removing orbit errors from geoid heights using crossover differences; gridding geoid height data at .25-degree by .25-degree intervals; and estimating the gravity anomalies from gridded geoid heights using the collocation technique.

Au, A. Y.

User's guide: Programs for processing altimeter data over inland seas

The programs described were developed to process GEODYN-formatted satellite altimeter data, and to apply the processed results to predict geoid undulations and gravity anomalies of inland sea areas. These programs are written in standard FORTRAN 77 and are designed to run on the NSESCC IBM 3081(MVS) computer. Because of the experimental nature of these programs they are tailored to the geographical area analyzed. The attached program listings are customized for processing the altimeter data over the Black Sea. Users interested in the Caspian Sea data are expected to modify each program, although the required modifications are generally minor. Program control parameters are defined in the programs via PARAMETER statements and/or DATA statements. Other auxiliary parameters, such as labels, are hard-wired into the programs. Large data files are read in or written out through different input or output units. The program listings of these programs are accompanied by sample IBM job control language (JCL) images. Familiarity with IBM JCL and the TEMPLATE graphic package is assumed.

Au, A. Y.