Transformation of Lambert Conic Conformal Coordinates from a Global Datum to a Local Datum

A detailed gravimetric geoid around Japan has been computed on the basis of 30’ × 30’ block mean free‐air gravity anomalies and GSFC GEM‐8 geopotential coefficient set. The 30’ × 30’ block means were read from various gravity maps around Japan, and the block means have been compiled into the JHDGF‐1 gravity file. Since the gravity file is restricted around Japan (see Figure 1), additional gravity data are needed to perform the Stokes’ integration in the cap with radius ψ₀ = 20°. The 1° × 1° block gravity means have been used outside the JHDGF‐1 region. The remarkable features of the gravimetric geoid occur over the trench areas. The geoidal dents over the trenches amount to ‐20~ ‐25 m in comparison with the geoidal heights in the land areas of Japan. The mean error of the 30’ × 30’ detailed gravimetric geoid obtained is estimated to be around 1.4 m, and the relative undulation of the geoid between the distance of a few hundred kilometers may be more accurate.     

Mean sea level determination from satellite altimetry

The primary experiment on the Geodynamics Experimental Ocean Satellite‐3 (GEOS‐3) is the radar altimeter. This experiment's major objective is to demonstrate the utility of measuring the geometry of the ocean surface, i.e., the geoid. Results obtained from this experiment so far indicate that the planned objectives of measuring the topography of the ocean surface with an absolute accuracy of ±5 m can be met and perhaps exceeded. The GEOS‐3 satellite altimeter measurements have an instrument precision in the range of ±25 cm to ±50 cm when the altimeter is operating in the “short pulse”; mode. After one year's operations of the altimeter, data from over 5,000 altimeter passes have been collected. With the mathematical models developed and the altimeter data presently available, mapping of local areas of ocean topography has been realized to the planned accuracy levels and better. This paper presents the basic data processing methods employed and some interesting results achieved with the early data. Plots of mean sea surface heights as inferred by the altimeter measurements are compared with a detailed 1^o × 1° gravimetric geoid.     

Comparison of geosat and gravimetric geoid profiles in the North Sea

Sea surface height profiles derived from 2‐year, repeat track, Geosat altimeter data have been compared with a regional gravimetric geoid in the western North Sea, computed using a geopotential model and terrestrial gravity data. The comparison encompasses 18 Geosat profiles covering a 750 × 850 km area of the North Sea. After a second‐order polynomial was used to model the long‐wavelength differences which cannot be clearly separated over an area of this size, results show agreement to better than ±3 cm for wavelengths between approximately 20 and 750 km. In regions where terrestrial gravity data were not available to improve the geoid, similar comparisons with the OSU91A geopotential model alone show differences of up to ±6 cm. This illustrates the importance of incorporating local gravity data in regional geoid computations, and partly validates the regional gravimetric geoid solution and Geosat sea surface profiles in the western North Sea. It is concluded that, in marine areas where the sea surface topography is known to be small in magnitude, Geosat sea surface profiles can act as an independent control on gravimetric geoids in the medium‐wavelength range.     

Residual depth and residual geoid anomalies of the South Atlantic ocean

Abstract

We have obtained the residual depth and residual geoid anomalies of the South Atlantic Ocean and interpreted them in terms of upper mantle processes. Starting from the 5’ X 5’ SYNBAPS data, we computed the 1° X 1° mean water depth. We digitized a recent sediment thickness map of the South Atlantic and corrected for the isostatic sediment loading effects. A plot of the corrected basement depth against crustal ages shows a good match to the depth‐age curve of the plate model. We then subtracted the predicted plate model depths from the corrected basement depths and obtained the 1° X 1° residual depth. The residual depth anomalies have positive values over the topographic highs and relative negative values over the ocean basins. A prominent asymmetry is observed between the residual depth over the Argentine Basin and that over the Cape Basin.

We have obtained the 1° X 1° residual geoid of the South Atlantic by subtracting both the long wavelength features and the geoid variations due to the plate cooling from the 1° X 1° Seasat altimeter derived geoid. The long wavelength features are represented by the degree and order 10 geoid of GEM1OC, and the geoid variations due to the plate cooling effects are predicted by the plate model geoid‐age relationship. The residual geoid anomalies also show an asymmetry although weaker than the case of the residual depth, between the Argentine and Cape basins.

By taking the 5° X 5° averages, we removed possible plate bending effects on the depth and geoid anomalies. We then compared those two data sets with respect to the reported hot spot tracks in the South Atlantic. The residual depth generally shows positive values over the hot spot tracks, whereas the residual geoid does not show any sign of them. A prominent asymmetric feature of depth and geoid anomalies is observed over the Argentine and Cape basins. This asymmetry is probably caused by hotter and less dense materials beneath the Cape Basin. Hot spots or other mantle upwellings may be the heat sources.     

Combined satellite altimetry and shipborne gravimetry data processing

The recovery of quantities related to the gravity field (i.e., geoid heights and gravity anomalies) is carried out in a test area of the central Mediterranean Sea using 5' × 5' marine gravity data and satellite altimeter data from the Geodetic Mission (GM) of ERS‐J. The optimal combination of the two heterogeneous data sources is performed using (1) the space‐domain least‐squares collocation (LSC) method, and (2) the frequency‐domain input‐output system theory (IOST). The results derived by these methods agree at the level of 2 cm in terms of standard deviation in the case of the geoid height prediction. The gravity anomaly prediction results by the same methods vary between 2.18 and 2.54 mGal in terms of standard deviation. In all cases, the spectral techniques have a much higher computational efficiency than the collocation procedure. In order to investigate the importance of satellite altimetry for gravity field modeling, a pure gravimetric geoid solution, carried out in a previous study for our lest area by the fast collocation approach (FCOL), is used in comparison with the combined geoid models. The combined solutions give more accurate results, at the level of about 15 cm in terms of standard deviation, than the gravimetric geoid solution, when the geoid heights derived by each method are compared with TOPEX altimeter sea surface heights (SSHs). Moreover, nonisotropic power spectral density functions (PSDs) can be easily used by IOST, while LSC requires isotropic covariance functions. The results show that higher prediction accuracies are always obtained when using a priori nonisotropic information instead of isotropic information.     

Simulation of the spatiotemporal variability of the World Ocean sea surface hight by the INM climate models

The results of simulations of the World Ocean sea surface hight (SSH) in by various versions of the Climate Model of the Institute of Numerical Mathematics, Russian Academy of Sciences, are compared with the CNES-CLS09 fields of the mean dynamic topography (deviation of the ocean level from the geoid). Three models with different ocean blocks are considered which slightly differ in numerical schemes and have various horizontal spatial resolution, i.e., the INMCM4 model, which participated in the Climate Model Intercomparison Project (CMIP Phase 5, resolution of 1° × 1/2°); the INMCM5 model, which participates in the next project, CMIP6 (resolution of 1/2° × 1/4°); and the advanced INMCM-ER eddy-resolving model (resolution of 1/6° × 1/8°). It is shown that an increase in the spatial resolution improves the reproduction of ocean currents (with Agulhas and Kuroshio currents as examples) and their variability. A probable cause of relatively high errors in the reproduction of the SSH of Southern and Indian oceans is discussed.     

High resolution gravimetric geoid heights and gravimetric vertical deflections of Europe including marine areas

Gravimetric geoid heights and gravimetric vertical deflections have been detemined for Europe including the Mediterranean Sea, North Sea, Norwegian Sea, Baltic Sea and parts of the North Atlantic Ocean in a 12′×20′ grid. The computation has been carried out by least squares spectral combination using closed integral formulas, combining 104 000 mean free air gravity anomalies in 6′×10′ blocks, 12 000 mean free air gravity anomalies in 1⁰×1⁰ blocks and the sherical harmonic model GEM9. The precision of the computed geoid heights has been estimated to ±1 m, the precision of the computed vertical deflections has been estimated to ±2″. Comparisons of the gravimetric geoid heights and vertical deflections with a number of other solutions have been carried out, confirming the precision estimation.     

5’ Detailed gravimetric geoid in the northwestern Atlantic ocean

A 5’ detailed gravimetric geoid has been computed for the northwest Atlantic Ocean as ground truth for the GEOS‐3 satellite altimeter experiment. Comparisons of this geoid with satellite derived geoceiver station heights show an r.m.s. difference of 1.2 m. Initial comparisons with GEOS‐III altimeter derived geoid profiles have indicated a relative agreement of generally better than 2 m.     

Geoid from geopotential model in the Taiwan area

The geoid undulation on GRS80 in the Taiwan area at half‐degree grid points has been calculated using the reduced 30’ × 30’ block mean gravity anomalies and the OSU91A geopotential coefficient set up to degree and order 360. The OSU91A results have been used to compare with WGS84, CEM10C, and OSU86F geoid undulations determined in 18 first‐order triangulation stations of the Taiwan Geodetic Datum 1980 (TGD80). Comparisons have also been made between these free‐air anomalies determined from OSU91A, and terrestrial gravity anomalies. It has been found that the average difference between the OSU91A model‐derived, and 243 actual point free‐air anomalies is 16.8 ± 48.0 mgal. It has also been found that more reliable and dense terrestrial gravity data are needed, both for terrestrial observations and for the OSU91A model, to achieve the very high‐precision geoid on GRS80 in the area of study.     

Bathymetric and magnetic survey of California Seamount (17°40′N × 124°W)

A bathymetric and magnetic survey of the California Seamount region (17°40′N × 124°00′W) shows that existing charts are in error. California Seamount is a peak extending to within 454 m (248 fathoms) of the surface. Its true location is 17°41′N × 124°01′W, 25 km southwest of the charted position. Near the old charted position there is an elongated feature which extends to within 1818 m (994 fathoms) of the surface. Both features are located on the Clarion Fracture Zone.     

CALCULATION AND ANALYSIS OF THE HEIGHT ANOMALY AND THE DEFLECTION OF VERTICAL IN AN OCEAN AREA

By the analysis of a practical calculation, this paper describes, for the first time in China, the gravimetric method on the calculation of the height anomaly and the deflection of vertical in the ocean by Stokes' and Vening Meinesz's formula. There are 84 calculation points distributed uniformlic in a calculated area of 2°×2° in the Mid-Pacific. In the course of the calculation, the gravimetric data measured by us, the 1° ×1° mean gravity data "published in other countries and the 25-ordeic gravitational coefficients of GEM8 were used. The results (Fig. 2b) show that the calculated area is an uplift of the geoid, with a mean height anomaly of 42 m, the maximum being 45 m and the minimum 39 m. In the whole calculated area, the variation of the deflection of vertical is rather small, with the maximum 7″·1 and the minimum -0″·2. The major causes of the calculation errors are pointed out and the calculation results are compaired with the data from the satellite altimeter.     

On the Impact of Mispointing Error and Hamming Filtering on Altimeter Waveform Retracking and Skewness Retrieval

An adequate conceptual definition of the geoid is essential for the unambiguous combination of satellite tracking data, satellite al‐timetry, and surface gravity measurements to obtain sea surface topography. The factors influencing the selection of a particular level surface of the earth's gravity field include the purpose(s) for which the geoid is to be used at the 5‐cm level, and the types of data to be used in achieving these objectives. The principal reasons for high precision determinations of the shape of the geoid are: the determination of sea surface topography for applications in oceanography; and the unification of leveling datums with a resolution equivalent to that of first order geodetic leveling. A conceptual definition of the geoid acceptable to oceanographers would be: The geoid for a selected epoch of measurement is that level surface of the earth's gravity field in relation to which the average non‐tidal (or quasi‐stationary) sea surface topography is zero as sampled globally in ocean regions. In the geodetic context, it would be convenient, though not essential, to modify this definition in such a way that the global sea surface topography had zero mean as sampled for evaluations of the geodetic boundary value problem. In either case, a basis exists for unifying all leveling datums serving areas in excess of 10⁶ km², using either gravity anomaly data for the regions or precise determinations of position at first order bench marks. Unfavorable signal‐to‐noise ratios can pose problems when dealing with datums serving smaller areas. Elevation and gravity data banks must be correctly referenced to leveling datums prior to use in sea surface topography determinations. A recent attempt to upgrade the Australian gravity anomaly data bank indicates that all current data banks of this type are inadequate for the task. It is unlikely that time variations in the radial position of the geoid as conceptually defined above, will exceed ±5 cm per century, provided the rate of earth expansion was less than 1 part in 10¹⁰ yr^‐l and there is no dramatic change in the present rate of secular change in Mean Sea Level.     

Hydrography and circulation of the bay of Bengal during withdrawal phase of the southwest Monsoon

Hydrographic data were collected from 3 to 10 September 1996 along two transects; one at 18° N and the other at 90° E. The data were used to examine the thermohaline, circulation and chemical properties of the Bay of Bengal during the withdrawal phase of the southwest monsoon. The surface salinity exhibited wide spatial variability with values as low as 25.78 at 18° N / 87° E and as high as 34.79 at 8° N / 90° E. Two high salinity cells (S > 35.2) were noticed around 100 m depth along the 90° E transect. The wide scatter in T-S values between 100 and 200 m depth was attributed to the presence of the Arabian Sea High Salinity (ASHS) water mass. Though the warm and low salinity conditions at the sea surface were conducive to a rise in the sea surface topography at 18° N / 87° E, the dynamic height showed a reduction of 0.2 dyn.m. This fall was attributed to thermocline upwelling at this location. The geostrophic currents showed alternating flows across both the transects. Relatively stronger and mutually opposite currents were noticed around 25 m depth across the 18° N transect with velocity slightly in excess of 30 cm s⁻¹. Similar high velocity (> 40 cm s⁻¹) pockets were also noticed to extend up to 30 m depths in the southern region of the 90° E transect. However, the currents below 250 m were weak and in general < 5 cm s⁻¹. The net geostrophic volume transports were found to be of the order of 1.5 × 10⁶ m³ s⁻¹ towards the north and of 6 × 10⁶ m³ s⁻¹ towards west across the 18° N and 90° E transects respectively. The surface circulation patterns were also investigated using the trajectories of drifting buoys deployed in the eastern Indian Ocean around the same observation period. Poleward movement of the drifting buoy with the arrival of the Indian Monsoon Current (IMC) at about 12° N along the eastern rim of the Bay of Bengal has been noticed to occur around the beginning of October. The presence of an eddy off the southeast coast of India and the IMC along the southern periphery of the Bay of Bengal were also evident in the drifting buoy data.     

Tidal estimation in the pacific with application to SEASAT altimetry

Abstract

The techniques for computing the eigenfunctions of the velocity potential (Proudman functions) set out in Sanchez et al. (1986) in relation to the Atlantic‐Indian Ocean are here applied to the Pacific Ocean, using a 6° × 6° grid of 510 points (455 points for the associated stream functions). Normal modes are computed from the first 150 Proudman functions and have natural periods from 43.9 hours downwards. Tidal syntheses are derived from these modes by direct application of the (frictionless) dynamic equations and by least‐squares fitting of Proudman functions to the dynamically interpolated tide‐gauge data of Schwiderski (1983). The modes contributing most energy to the principal harmonic tidal constituents are different in the two computations; their natural periods are typically in the range of 9–16 hours for semidiurnal, 14–43 hours for diurnal tides. The RMS of fit for Proudman functions is in all cases better than the corresponding value for the same number of spherical harmonics.

Before fitting the Proudman functions to the altimetry from the 3‐day repeat cycle of Seasat, the data are processed by novel methods. The geoid component is eliminated by taking collinear differences at a fixed time‐lag of two repeat cycles. Orbit errors are reduced by extracting the 1 rev^?1 component at every ascending node; this component varies slowly and nonlinearly in time. The spatial fitting process includes M₂ and O₁ frequencies, both of which emerge with significant and realistic tidal mapping, but residual noise in the data limits the number of Proudman functions to about 50–60 before showing signs of “over‐fitting.”; Fitting the same data by spherical harmonics gives marginally lower predicted variance for the same number of parameters.     

Observations on the hydrology of bay of Islands,New Zealand

Observations were made on several hydrological features of Bay of Islands during 1970 to 1971. The topography of Bay of Islands allows for a gradual change from estuarine to oceanic conditions within the harbour.

Monthly means of sea surface temperatures ranged between 15°c and 23°c, and some sea stratification was observed during summer. Salinities in the main basin were about 35.5‰; water transparency ranged from 2 m to 6 m by Secchi disc in the estuaries, to more than 15m in outer basin areas. Dissolved oxygen content was high, usually exceeding 100% saturation in surface waters. Water circulation within the bay appears to consist of an anti‐clockwise movement of at least the surface water, induced by a north‐west moving current, possibly derived from the East Auckland Current.

The observations suggest that, except; for the estuarine areas, Bay of Islands is hydrologically a fairly homogeneous, well‐mixed body of water.     

Ocean temperature climate off North‐East New Zealand

The ocean temperature field off the north‐east coast of New Zealand is studied to quantify the annual cycle and reveal the intra‐ and inter‐annual variability. The data used are repeat expendable bathythermograph (XBT) sections between Auckland and either Suva or Honolulu which have been collected quarterly since 1986. These sections give temperature measurements between the surface and 800 m and Auckland and 30°S from 1986 to August 1999. The mean and annual cycle are compared with those from the NOAA World Ocean Atlas (WOA98). The results are similar; however WOA98 lacks the horizontal resolution to fully discern the East Auckland Current and North Cape Eddy, while the XBT analysis lacks the temporal resolution to discern higher frequency intra‐annual signals. The temperature variability in the mixed layer is dominated by the annual cycle, which accounts for 80–90% of the variance. The amplitude of the annual cycle diminishes rapidly with depth, from 2.8°C at the surface, to c. 0.1°C at 180 m. The phase of the annual cycle is retarded with depth, with peak temperatures occurring in February at the surface and in June/July at 180 m. Removing the annual cycle from the time series reveals the more subtle inter‐ and intra‐annual variability. This variability is of the order of 1°C in the upper 50 m, decreasing to 0.3°C at 400–500 m. The surface layer was cold between 1991 and 1994 (c. 0.7°C cooler than average), and 0.7°C warmer than average in 1999. The deeper ocean shows a different signal, being up to 0.3°C cooler in 1990–92, 0.3°C warmer in 1998, and c. 0.2°C warmer than average in 1999. The inter‐annual mixed layer variability is highly correlated with the Southern Oscillation Index and also with inter‐annual terrestrial air temperature and wind measurements from northern New Zealand. In contrast, at higher intra‐annual frequencies, the mixed layer variability is not correlated with air and wind measurements. At these higher frequencies, the air temperature is better correlated with the sea surface temperature (SST) than with the bulk mixed layer temperature.     

南极绕极流的经向输运

南极绕极流(ACC)是南大洋中最显著的流动，流量超过130×106m3/s (Nowlin et al.，1986)。传统认为，由于以东向运动为主的ACC的存在极大地阻碍了南大洋中上层的南北向物质和能量的交换，绕极流区的经向输运是非常小的。但是近些年的研究发现，穿过ACC的通量并不是可以忽略不计的，它对维持南极和亚南极区的动力和热力平衡起着重要作用，在全球气候系统中也有着深刻的影响(Doos et al.，1994)。     

东海浮游翼足类(Pteropods)数量分布的研究

根据1997～2000年东海海域23°30'～33°00'N,118°30'～128°00'E的4个季节海洋调查资料,运用定量、定性方法,探讨了东海浮游翼足类总丰度的平面分布、季节变化及变化的动力学机制.结果表明,东海翼足类总丰度和出现频率有明显的季节变化,均为秋季最高,夏季次之,春季最低;总丰度在各个季节基本上呈东海南部高于北部、外海高于近海的分布趋势;春季的尖笔帽螺(Creseis acicula)、夏季的锥笔帽螺(Creseis virgula)、秋季的蝴蝶螺(Desmopterus papilio)和冬季的马蹄螔螺(Limacina trochiformis)是导致总丰度季节变化的最主要的种类;冬、春和夏3个季节丰度变化及4季总丰度的变化同表层或10m层水温有非常显著的线性相关关系,与底层温度及盐度的相关关系不显著.夏季翼足类高丰度区位于台湾暖流与黑潮暖流的分支处;从夏季到秋季,翼足类随着台湾暖流向北扩展,并在与长江冲淡水,闽浙沿岸水团,黄海水团等交汇处形成高丰度(大于500×10-2个/m3)和较高丰度(250×10-2～500×10-2个/m3)分布区.水温和海流是影响东海翼足类总丰度分布的主要环境因素.     

Global ocean circulation patterns based on seasat altimeter data and the GEML2 gravity field

The Seasat altimeter data has been completely adjusted by a crossing arc technique to reduce the crossover discrepancies to approximately ±30 cm in five regional adjustments. This data was then used to create sea surface heights at 1° intersections in the ocean areas with respect to the GRS80 ellipsoid. These heights excluded the direct tidal effects but included the induced permanent deformation. A geoid corresponding to these sea surface heights was computed, based on the potential coefficients of the GEML2 gravity field up to degree 6, augmented by Rapp's coefficients up to degree 180. The differences between sea surface heights and the geoid were computed to give approximate estimates of sea surface topography. These estimates are dominated by errors in both sea surface heights and geoid undulations. To optimally determine sea surface topography a spherical harmonic analysis of raw estimates was carried out and the series was further truncated at degree 6, giving estimates with minimum wavelengths on the order of 6000 km. The direction of current flow can be computed on a global basis using the spherical harmonic expansion of the sea surface topography. Ths has been done, not only for Seasat/GEML2 estimates, but also using the recent dynamic topography estimates of Levitus. The results of the two solutions are very similar and agree well with the major circulation features of the oceans.     

The Role of Multi-Mission ERS Altimetry in the Determination of the Marine Geoid in the Azores

This study concerns the determination of a regional geoid model in the North Atlantic area surrounding the Azores islands by combining multi-mission altimetry from the ERS (European Remote Sensing) satellites and surface gravity data. A high resolution mean sea surface, named AZOMSS99, has been derived using altimeter data from ERS-1 and ERS-2 35-day cycles, spanning a period of about four years, and from ERS-1 geodetic mission. Special attention has been paid to data processing of points around the islands due to land contamination on some of the geophysical corrections. A gravimetric geoid has been computed from all available surface gravity, including land and sea observations acquired during an observation campaign that took place in the Azores in October 1997 in the scope of a European and a Portuguese project. Free air gravity anomalies were derived by altimetric inversion of the mean sea surface heights. These were used to fill the large gaps in the surface gravity and combined solutions were computed using both types of data. The gravimetric and combined solutions have been compared with the mean sea surface and GPS (Global Positioning System)-levelling derived geoid undulations in five islands. It is shown that the inclusion of altimeter data improves geoid accuracy by about one order of magnitude. Combined geoid solutions have been obtained with an accuracy of better than one decimetre.