ALBERT

Description: The oceans have a major impact on global geophysical processes of the Earth. Non-tidal changes in oceanic currents and ocean-bottom pressure have been shown to be a major source of polar motion excitation and also measurably change the length of the day. The changing mass distribution of the oceans causes the Earth's gravitational field to change and causes the center-of-mass of the oceans to change which in turn causes the center-of-mass of the solid Earth to change. The changing mass distribution of the oceans also changes the load on the oceanic crust, thereby affecting both the vertical and horizontal position of observing stations located near the oceans. Recognizing the important role that non-tidal oceanic processes play in Earth rotation dynamics and terrestrial reference frame definition, the International Earth Rotation Service has recently created a Special Bureau for the Oceans in order to facilitate research into these and other solid Earth geophysical processes affected by the oceans.

Description: TOPEX/Poseidon has been collecting altimeter data continuously since October 1992. Altimeter data have been used to produce maps of sea surface height, geostrophic velocity, significant wave height, and wind speed. This information is of proven use to mariners as well as to the scientific community. Uses of the data include commercial and recreational vessel routing, ocean acoustics, input to geographic information systems developed for the fishing industry, identification of marine mammal habitats, fisheries management, and monitoring ocean debris. As with sea surface temperature data from the Advanced Very High Resolution Radiometer (AVHRR) in the late 1980s and early 1990s, altimeter data from TOPEX/Poseidon and ERS-1 and -2 are in the process of being introduced to the marine world for operational maritime use. It is anticipated that over the next few years companies that specialize in producing custom products for shipping agencies, fisheries and yacht race competitors will be incorporating altimeter data into their products. The data are also being incorporated into weather and climate forecasts by operational agencies both in the US and Europe. This paper will discuss these products, their uses, operational demonstrations and means of accessing the data.

Description: Oceanic internal tides are internal waves with tidal periodicities. They are ubiquitous throughout the ocean, although generally more pronounced near large bathymetric features such as mid-ocean ridges and continental slopes. The internal vertical displacements associated with these waves can be extraordinarily large. Near some shelf breaks where the surface tides are strong, internal displacements (e.g., of an isothermal surface) can exceed 200 meters. Displacements of 10 meters in the open ocean are not uncommon. The associated current velocities are usually comparable to or larger than the currents of the surface tide. On continental shelves internal tides can occasionally generate packets of internal solitons, which are detectable in remote sensing imagery. Other common nonlinear features are generation of higher harmonics (e.g., 6-hr waves) and wave breaking. Internal tides are known to be an important energy source for mixing of shelf waters. Recent research suggests that they may also be a significant energy source for deep-ocean mixing.

Description: We study the tracer subgrid term in isopycnal coordinates, S(sub I). We employ two ingredients: the experimental data on vertical spectra of ocean turbulence measured by Gargett et al.(1981) and the stochastic approach recently developed by Dukowicz and Smith (1997). Our result confirms that S(sub I) is made of two parts: an advection and a diffusion term. However, the tracer bolus velocity u** consists of two terms u** = u(sub 1) + u(sub 2) while in the GM model there is only a term related to u(sub 1) which is shown to be: u(sub 1) = k(bar-q)(sup -1)(delta)(sub rho) where bar-q is the thickness weighted average potential vorticity, a result in agreement with the recent suggestions by Treguier et al. (1997), Lee et al. (1997) and Greatbatch (1998). The second component u(sub 2) IS new. We compute it in the geostrophic approximation using the Gargett et al. data (1981) on ocean vertical turbulence. We find that u(sub 2) much greater than u(sub 1) and that u(sub 2) is orthogonal to u(sub 1).

Description: Sea surface height (SSH) from altimeter observations from 1992 on and from modeling results is investigated to determine the modes of variability and the linkages to the state of oceanic circulation in the North Atlantic. First the altimeter and model simulated SSH are analyzed using the empirical orthogonal function (EOF) analysis. They are found to share a similar leading mode where the center of action is along the Gulf Stream and North Atlantic Current with opposite sign anomalies in the subpolar gyre and in the slope waters along the Eastern Seaboard. The time series of the leading EOF mode from the altimeter data shows that between winters of 1995 and 1996, SSH over the Gulf Stream decreased by about 12cm which change is reproduced by the model simulation. Based on the relationship from the model simulations between the time series of the SSH EOF1 and meridional heat transport, it is suggested that associated with this SSH change in 1995-96, the overturning has slowed down from its heights in the early 90's. Furthermore, it is shown that decadal variability in the leading SSH mode originates from the thermal forcing component. This adds confidence to the qualitative relationship between the state of overturning/meridional heat transport and SSH in the limited area described by the EOF1. SSH variability in the eastern side of the North Atlantic basin, outside the western boundary current region, is determined by local and remote (Rossby waves) wind stress curl forcing.

Description: The global stratospheric ozone-layer depletion results In an increase in biologically harmful ultraviolet (UV) radiation reaching the surface and penetrating to ecologically significant depths in natural waters. Such an increase can be estimated on a global scale by combining satellite estimates of UV irradiance at the ocean surface from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument with the SeaWIFS satellite ocean-color measurements in the visible spectral region. In this paper we propose a model of seawater optical properties in the UV spectral region based on the Case I water model in the visible range. The inputs to the model are standard monthly SeaWiFS products: chlorophyll concentration and the diffuse attenuation coefficient at 490nm. Penetration of solar UV radiation to different depths in open ocean waters is calculated using the RT (radiative transfer) quasi-single scattering approximation (QSSA). The accuracy of the QSSA approximation in the water is tested using more accurate codes. The sensitivity study of the underwater UV irradiance to atmospheric and oceanic optical properties have shown that the main environmental parameters controlling the absolute levels of the UVB (280-320nm) and DNA-weighted irradiance underwater are: solar-zenith angle, cloud transmittance, water optical properties, and total ozone. Weekly maps of underwater UV irradiance and DNA-weighted exposure are calculated using monthly-mean SeaWiFS chlorophyll and diffuse attenuation coefficient products, daily SeaWiFS cloud fraction data, and the TOMS-derived surface UV irradiance daily maps. The final products include global maps of weekly-average UVB irradiance and DNA-weighted daily exposures at 3m and 10m, and depths where the UVB irradiance and DNA-weighted dose rate at local noon are equal to 10% of their surface values.

Description: Seeking to improve upon the visualization of the Antarctic Circumpolar Wave (ACW) , we compare a 16-year sequence of 6-month winter averages of Antarctic sea ice extents and concentrations with those of adjacent sea surface temperatures (SSTs). Here we follow SSTs around the globe along the maximum sea ice edge rather than in a zonal band equatorward of it. The results are similar to the earlier ones, but the ACWs do not propagate with equal amplitude or speed. Additionally in a sequence of 4 polar stereographic plots of these SSTs and sea ice concentrations, we find a remarkable correlation between SST minima and sea ice concentration maxima, even to the extent of matching contours across the ice-sea boundary, in the sector between 900E and the Palmer Peninsula. Based on these observations, we suggest that the memory of the ACW in the sea ice is carried from one Austral winter to the next by the neighboring SSTS, since the sea ice is nearly absent in the Austral summer.

Description: Detailed knowledge of ocean circulation and its transport properties is prerequisite to an understanding of the earth's climate and of important biological and chemical cycles. Results from two recent experiments, THETIS-2 in the Western Mediterranean and ATOC in the North Pacific, illustrate the use of ocean acoustic tomography for studies of the large scale circulation. The attraction of acoustic tomography is its ability to sample and average the large-scale oceanic thermal structure, synoptically, along several sections, and at regular intervals. In both studies, the acoustic data are compared to, and then combined with, general circulation models, meteorological analyses, satellite altimetry, and direct measurements from ships. Both studies provide complete regional descriptions of the time-evolving, three-dimensional, large scale circulation, albeit with large uncertainties. The studies raise serious issues about existing ocean observing capability and provide guidelines for future efforts.

Description: The distribution and variation of oxygen isotopes in seawater are calculated using the Goddard Institute for Space Studies global ocean model. Simple ecological models are used to estimate the planktonic foraminiferal abundance as a function of depth, column temperature, season, light intensity, and density stratification. These models are combined to forward model isotopic signals recorded in calcareous ocean sediment. The sensitivity of the results to the changes in foraminiferal ecology, secondary calcification, and dissolution are also examined. Simulated present-day isotopic values for ecology relevant for multiple species compare well with core-top data. Hindcasts of sea surface temperature and salinity are made from time series of the modeled carbonate isotope values as the model climate changes. Paleoclimatic inferences from these carbonate isotope records are strongly affected by erroneous assumptions concerning the covariations of temperature, salinity, and delta (sup 18)O(sub w). Habitat-imposed biases are less important, although errors due to temperature-dependent abundances can be significant.

Description: Under a data purchase agreement with Goddard Space Flight Center, Orbital Sciences Corporation has been able to contract building of the Sea-Viewing Wide-Field-of-View Sensor (SeaWIFS). Orbital Sciences was then able commercialized the data that the satellite produces. These data are used to create daily fish finding maps, allowing fishing fleets to focus on locations where many commercially important surface feeding fish, like tuna and swordfish, congregate. In agriculture and forestry, SeaWIFS' images offer an alternative to direct on-site inspection or expensive serial photography.

Description: Analytical Spectral Devices was able to commercialize the "engine" of SeaSpec. The Dual CCD developed for SeaSpec is used in our FieldSpecO Dual UV/VNIR/CCD. This product is very popular for coastal research. It allows for a small suitcase sized spectrometer to be placed in a boat and the attached fiber optic cable can be lowered 30 feet below the surface allowing researchers to perform the same studies as SeaSpec just not at the deep-water depths. Most of the people who inquire about SeaSpec end up purchasing a FieldSpec Dual CCD.SeaSpec was designed to withstand harsh marine environments and depths up to 200 m. While it meets these requirements it is the only instrument of its kind and we consider it a prototype. The first version was released in April 1997. It was tested, we made modification and returned the product in March 1998.

Description: The effect of random wave fields on passive tracer spatial variations is studied.

Description: This study examines the role of the Indonesian Throughflow (ITF) variations on the El Nino/La Nina events simulated by a coupled ocean-atmosphere model of the tropical Pacific.

Description: Satellite data have revealed overall decreases in the Arctic sea ice cover since the late 1970s, although with substantial interannual variability. The ice reductions are likely tied to an overall warming in the Arctic region over the same time period, although both the warming and the ice reductions could be connected to large-scale oscillations within the system. Should the ice reductions continue, consequences to the Arctic ecosystems and climate could be considerable.

Description: The measurement of global sea surface height made by the TOPEX/POSEIDON satellite has provided the first synoptic view of large-scale oceanic variability at the intraseasonal scales from weeks to months.

Description: Satellite altimetry and hydrographic observations are used to characterize the mesoscale eddy field in the Sargasso Sea near Bermuda and to address the role of physical processes on the supply of new nutrients to the euphotic zone. The observed sea level anomaly (SLA) field is dominated by the occurrence of westward propagating features with SLA signatures as large as 25 cm, Eulerian temporal scales of roughly a month, lifetimes of several months, spatial scales of approximately 200 km, and a propagation of approximately 5 cm/s . Hydrographic estimates of dynamic height anomaly (referenced to 4000 dbar) are well correlated with satellite SLA (r(sup 2) = 0.65), and at least 85% of the observed dynamic height variability is associated with the first baroclinic mode of motion. This allows us to apply the satellite observations to remotely estimate isopycnal displacements and the flux of nutrients into the euphotic zone due to eddy pumping. Eddy pumping is the process by which mesoscale eddies induce isopycnal displacements that lift nutrient- replete waters into the euphotic zone, driving new primary production. A kinematic approach to the estimation of the eddy pumping results in a flux of 0.24+/-0.1 mol N/sq m/yr (including a scale estimate for the small contribution due to 18 deg water eddies). This flux is more than an order of magnitude larger than the diapycnal diffusive flux as well as scale estimates for the vertical transport due to isopycnal mixing along sloping isopycnal surfaces. Eddy pumping and wintertime convection are the two dominant mechanisms transporting new nutrients into the euphotic zone, and the sum of all physical new nutrient supply fluxes effectively balances previous geochemical estimates of annual new production for this site. However, if biological transports (e.g., nitrogen fixation, etc.) are significant, the new nitrogen supply budget will be in excess of geochemical new production estimates. This suggests that the various physical and biological transport fluxes, as well as geochemical inferences of new production, still need to be reconciled and many outstanding questions remain.

Description: A mesoscale resolution biogeochemical survey was carried out in the vicinity of the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-series Study (BATS) site during the summer of 1996. Real-time nowcasting and forecasting of the flow field facilitated adaptive sampling of several eddy features in the area. Variations in upper ocean nutrient and pigment distributions were largely controlled by vertical isopycnal displacements associated with the mesoscale field. Shoaling density surfaces tended to introduce cold, nutrient-rich water into the euphotic zone, while deepening isopycnals displaced nutrient-depleted water downward. Chlorophyll concentration was generally enhanced in the former case and reduced in the latter. Eddy-induced upwelling at the base of the euphotic zone was affected by features of two different types captured in this survey-, (1) a typical mid-ocean cyclone in which doming of the main thermocline raised the near-surface stratification upward and (2) a mode water eddy composed of a thick lens of 18C water, which pushed up the seasonal thermocline and depressed the main thermocline. Model hindcasts using all available data provide a four-dimensional context in which to interpret temporal trends at the BATS site and two other locations during the 2 weeks subsequent to the survey. Observed changes in near-surface structure at the BATS site included shoaling isopycnals, increased nutrient availability at the base of the euphotic zone, and enhanced chlorophyll concentration within the cuphotic zone. These trends are explicable in terms of a newly formed cyclone that impinged upon the site during this time period. These observations reveal that eddy upwelling has a demonstrable impact on the way in which the nitrate-density relationship changes with depth from the aphotic zone into the euphotic zone. A similar transition is present in the BATS record, suggesting that eddy-driven upwelling events are present in the time series of upper ocean biogeochemical properties. The variability in main thermocline temperature and nitrate in this synoptic spatial survey spans the range observed in these quantities in the 10-year time series available at BATS to date (1988-1998).

Description: A mesoscale eddy advected across the Bermuda Testbed Mooring site over a 30-day period centered on July 14, 1995. Temperature and current measurements along with biogeochemical measurements were used to characterize the biological response of the upper ocean associated with the introduction of nitrate into the euphotic layer due to the doming of isotherms associated with the eddy. Complementary shipboard data showed an anomalous water mass, which extended from a depth of approximately 50 to 1000 m, manifesting as a cold surface expression and warm anomaly at depth. Although mesoscale eddies are frequently observed in the Sargasso Sea, the present observations are particularly unique because of the high-temporal-resolution measurements of the new instrumentation deployed on the mooring. Analyzers that measure nitrate plus nitrite were placed at depths of 80 and 200 m and bio-optical sensors were located at depths of 20, 35, 45, 71, and 86 m. Peak nitrate values of nearly 3.0 micro-M at 80 m and chlorophyll a values of 1.4 mg/cubic m at 71 m were observed, a well as a 25- to 30-meter shoaling of the 1% light level depth. A Doppler shift from the inertial period (22.8 hours) to 25.2 hours was observed in several time series records due to the movement of the eddy across the mooring. Inertial pumping brought cold, nutrient-rich waters farther into the euphotic zone than would occur solely by isothermal lifting. Silicic acid was depleted to undetectable levels owing to the growth of diatoms within the eddy. The chlorophyll a values associated with the eddy appear to be the largest recorded during the eight years of the ongoing US JGOFS Bermuda Atlantic Time Series Study program.

Description: Ocean surface remote sensing techniques often rely on scattering or emission linked to shorter- scale gravity-capillary ocean wavelets. However, it is increasingly apparent that slightly longer wavelengths of O(10 to 500 cm) are vital components in the robust sea surface description needed to link varied global remote sensing data sets. This paper describes a sensor suite developed to examine sea surface slope variations in the field using an aircraft flying at very low altitude (below 30 m) and will also provide preliminary measurements detailing changes in slope characteristics versus sea state and friction velocity. Two-dimensional surface slope is measured using simultaneous range measurements from three compact short-range laser altimeters mounted in an equilateral triangle arrangement with spacing of about 1 m. In addition, all three lasers provide independent wave elevation profiles after GPS-aided correction for aircraft altitude. Laser range precision is 1 cm rms while vertical motion correction is 15 cm rms. The measurements are made along-track at approximately 1 m intervals setting the spatial scale of the measurement to cover waves of intermediate to long scale. Products available for this array then include surface elevation, two-dimensional slope distribution, and the cross- and along-track 1-D slope distributions. To complement the laser, a down-looking mm-wave radar scatterometer is centered within the laser array to measure radar backscatter simultaneously with the laser slope. The radar's footprint is nominally 1 m in diameter. Near-vertical radar backscatter is inversely proportional to the small-scale surface slope variance and to the tilt of the underlying (laser-measured) surface facet. Together the laser and radar data provide information on wave roughness from the longest scales down to about 1 cm. These measurements are complemented by aircraft turbulence probe data that provides robust surface flux information.

Description: The importance of Fe biogeochemistry has stimulated interest in Fe isotope fractionation. Recent studies using thermal ionization mass spectrometry (TIMS) and a "double spike" demonstrate the existence of biogenic Fe isotope effects. Here, we assess the utility of multiple-collector inductively-coupled plasma mass spectrometry(MC-ICP-MS) with a desolvating sample introduction system for Fe isotope studies, and present data on Fe biominerals produced by a thermophilic bacterium. Additional information is contained in the original extended abstract.

Description: A mesoscale resolution biogeochemical survey was carried out in the vicinity of the US Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study (BATS) site during the summer of 1996. Real-time nowcasting and forecasting of the flow field facilitated adaptive sampling of several eddy features in the area. Variations in upper ocean nutrient and pigment distributions were largely controlled by vertical isopycnal displacements associated with the mesoscale field. Shoaling density surfaces tended to introduce cold, nutrient-rich water into the euphotic zone, while deepening isopycnals displaced nutrient-depleted water downward. Chlorophyll concentration was generally enhanced in the former case and reduced in the latter. Eddy-induced upwelling at the base of the euphotic zone was affected by features of two different types captured in this survey: (1) a typical mid-ocean cyclone in which doming of the main thermocline raised the near-surface stratification upward; and (2) a mode water eddy composed of a thick lens of 18 C water, which pushed up the seasonal thermocline and depressed the main thermocline. Model hindcasts using all available data provide a four-dimensional context in which to interpret temporal trends at the BATS site and two other locations during the two weeks subsequent to the survey. Observed changes in near-surface structure at the BATS site included shoaling iscpycnals, increased nutrient availability at the base of the euphotic zone, and enhanced chlorophyll concentration within the euphotic zone. These trends are explicable in terms of a newly formed cyclone that impinged upon the site during this time period. These observations reveal that eddy upwelling has a demonstrable impact on the way in which the nitrate-density relationship changes with depth from the aphotic zone into the euphotic zone. A similar transition is present in the BATS record, suggesting that eddy-driven upwelling events are present in the time series of upper ocean biogeochemical properties. The variability in main thermocline temperature and nitrate in this synoptic spatial survey spans the range observed in these quantities in the 10-year time series available at BATS to date (1988-1998).

Description: Satellite altimetry and hydrographic observations are used to characterize the mesoscale eddy field in the Sargasso Sea near Bermuda and to address the role of physical processes on the supply of new nutrients to the euphotic zone. The observed sea level anomaly (SLA) field is dominated by the occurrence of westward propagating features with SLA signatures as large as 25 cm, Eulerian temporal scales of roughly a month, lifetimes of several months, spatial scales of approximately 200 km, and a propagation of approximately 5 cm/s. Hydrographic estimates of dynamic height anomaly (referenced to 4000 dbar) are well correlated with satellite SLA (r(exp 2) = 0.65), and at least 85% of the observed dynamic height variability is associated with the first baroclinic mode of motion. This allows us to apply the satellite observations to remotely sensed estimate isopycnal displacements and the flux of nutrients into the euphotic zone due to eddy pumping. Eddy pumping is the process by which mesoscale eddies induce isopycnal displacements that lift nutrient-replete waters into the euphotic zone, driving new primary production. A kinematic approach to the estimation of the eddy pumping results in a flux of 0.24 +/- 0.1 mol N/sq m (including a scale estimate for the small contribution due to 18 deg water eddies). This flux is more than an order of magnitude larger than the diapycnal diffusive flux as well as scale estimates for the vertical transport due to isopycnal mixing along sloping isopycnal surfaces. Eddy pumping and wintertime convection are the two dominant mechanisms transporting new nutrients into the euphotic zone, and the sum of all physical new nutrient supply fluxes effectively balances previous geochemical estimates of annual new production for this site. However, if biological transports (e.g., nitrogen fixation, etc.) are significant, the new nitrogen supply budget will be in excess of geochemical new production estimates. This suggests that the various physical and biological transport fluxes, as well as geochemical inferences of new production, still need to be reconciled and many outstanding questions remain.

Description: The question of whether the Earth ever passed through a magma ocean stage is of considerable interest. Geochemical evidence strongly suggests that the Moon had a magma ocean and the evidence is mounting that the same was true for Mars. Analyses of martian (SNC) meteorites have yielded insights into the differentiation history of Mars, and consequently, it is interesting to compare that planet to the Earth. Three primary features of Mars contrast strongly to those of the Earth: (i) the extremely ancient ages of the martian core, mantle, and crust (about 4.55 b.y.); (ii) the highly depleted nature of the martian mantle; and (iii) the extreme ranges of Nd isotopic compositions that arise within the crust and depleted mantle. The easiest way to explain the ages and diverse isotopic compositions of martian basalts is to postulate that Mars had an early magma ocean. Cumulates of this magma ocean were later remelted to form the SNC meteorite suite and some of these melts assimilated crustal materials enriched in incompatible elements. The REE pattern of the crust assimilated by these SNC magmas was LREE enriched. If this pattern is typical of the crust as a whole, the martian crust is probably similar in composition to melts generated by small degrees of partial melting (about 5%) of a primitive source. Higher degrees of partial melting would cause the crustal LREE pattern to be essentially flat. In the context of a magma ocean model, where large degrees of partial melting presumably prevailed, the crust would have to be dominated by late-stage, LREE-enriched residual liquids. Regardless of the exact physical setting, Nd and W isotopic evidence indicates that martian geochemical reservoirs must have formed early and that they have not been efficiently remixed since. The important point is that in both the Moon and Mars we see evidence of a magma ocean phase and that we recognize it as such. Several lines of theoretical inference point to an early Earth that was also hot and, perhaps, mostly molten. The Giant Impact hypothesis for the origin of the Moon offers a tremendous input of thermal energy and the same could be true for core formation. And current solar system models favor the formation of a limited number of large (about 1000 km) planetesimals that, upon accreting to Earth, would cause great heating, being lesser versions of the Giant Impact. Several lines of geochemical evidence do not favor this hot early Earth scenario. (i) Terrestrial man-tle xenoliths are sometimes nearly chondritic in their major element compositions, suggesting that these rocks have never been much molten. Large degrees of partial melting probably promote differentiation rather than homogenization. (ii) Unlike the case of Mars, the continental crust probably did not form as a highly fractionated residual liquid from a magma ocean (about 99% crystallization), but, rather, formed in multiple steps. [The simplest model for the formation of continental crust is complicated: (a) about 10% melting of a primitive mantle, making basalt; (b) hydrothermal alteration of that basalt, converting it to greenstone; and (c) 10% partial melting of that greenstone, producing tonalite.] This model is reinforced by the recent observation from old (about 4.1 b.y.) zircons that the early crust formed from an undepleted mantle having a chondritic Lu/Hf ratio. (iii) If the mantle were once differentiated by a magma ocean, the mantle xenolith suite requires that it subsequently be homogenized. The Os isotopic compositions of fertile spinel lherzolites place constraints on the timing of that homogenization. The Os isotopic composition of spinel lherzolites approaches that of chondrites and correlates with elements such as Lu and Al. As Lu and Al concentrations approach those of the primitive mantle, Os isotopic compositions approach chondritic. The Re and Os in these xenoliths were probably added as a late veneer. Thus, the mantle that received the late veneer must have been nearly chondritic in terms of its major elements (excluding Fe). If the mantle that the veneer was mixed into was not al-ready homogenized, then Os isotopes should not correlate with incompatible elements such as Al. Consequently, either early differentiation of the mantle did not occur or the homogenization of this differentiation must have occurred before the late veneer was added. The timing of the late veneer is itself uncertain but presumably postdated core formation at about 4.45 b.y. and did not postdate the 3.8-3.9 b.y. late bombardment of the Moon. This timing based on siderophile elements is consistent with the Hf isotopic evidence cited above. If the Earth, Moon and Mars had magma oceans, the Earth subsequently rehomogenized whereas the Moon and Mars did not. The simplest solution to this observation is that homogenization of igneous differentiates was never necessary on Earth, either because the hypothetical magma ocean never occurred or because this event did not produce mantle differentiation.

Description: Mechanisms controlling the variation in sea surface temperature warm event in the equatorial Pacific were investigated through ocean model simulations. In addition, the mechanisms of the climatological SST cycle were investigated. The dominant mechanisms governing the seasonal cycle of SST vary significantly across the basin. In the western Pacific the annual cycle of SST is primarily in response to external heat flux. In the central basin the magnitude of zonal advection is comparable to that of the external heat flux. In the eastern basin the role of zonal advection is reduced and the vertical mixing is more important. In the easternmost equatorial Pacific the vertical entrainment contribution is as large as that of vertical diffusion. The model estimate of the vertical mixing contribution to the mixed layer heat budget compared well with estimates obtained by analysis of observations using the same diagnostic vertical mixing scheme. During 1994- 1995 the largest positive SST anomaly was observed in the mid-basin and was related to reduced latent heat flux due to weak surface winds. In the western basin the initial warming was related to enhanced external heating and reduced cooling effects of both vertical mixing and horizontal advection associated with weaker than usual wind stress. In the eastern Pacific where winds were not significantly anomalous throughout 1994-1995, only a moderate warm surface anomaly was detected. This is in contrast to strong El Nino events where the SST anomaly is largest in the eastern basin and, as shown by previous studies, the anomaly is due to zonal advection rather than anomalous surface heat flux. The end of the warm event was marked by cooling in July 1995 everywhere across the equatorial Pacific.

Description: A mesoscale eddy advected across the Bermuda Testbed Mooring site over a 30-day period centered on July 14, 1995. Temperature and current measurements along with biogeochemical measurements were used to characterize the biological response of the upper ocean associated with the introduction of nitrate into the euphoric layer due to the doming of isotherms associated with the eddy. Complementary shipboard data showed an anomalous water mass, which extended from a depth of approximately 50 to 1000 m, manifesting as a cold surface expression and warm anomaly at depth. Although mesoscale eddies are frequently observed in the Sargasso Sea, the present observations are particularly unique because of the high-temporal-resolution measurements of the new instrumentation deployed on the mooring. Analyzers that measure nitrate plus nitrite were placed at depths of 80 and 200 m and bio-optical sensors were located at depths of 20, 35, 45, 71, and 86 m. Peak nitrate values of nearly 3.0 microns at 80 m and chlorophyll alpha values of 1.4 mg/cu m at 71 m were observed, as well as a 25- to 30-meter shoaling of the 1% light level depth. A Doppler shift from the inertial period (22.8 hours) to 25.2 hours was observed in several time series records due to the movement of the eddy across the mooring. Inertial pumping brought cold, nutrient-rich waters farther into the euphotic zone than would occur solely by isothermal lifting. Silicic acid was depleted to undetectable levels owing to the growth of diatoms within the eddy. The chlorophyll alpha values associated with the eddy appear to be the largest recorded during the 8 years of the ongoing U.S. JGOFS Bermuda Atlantic Time Series Study (BATS) program.

Description: The question of whether the Earth ever passed through a magma ocean stop is of considerable interest. Geochemical evidence strongly suggests that the Moon had a magma ocean and the evidence is mounting that the same was true for Mars. Analyses of mar (SNC) meteorites have yielded insights into the differentiation history of Mars, and consequently, it is interesting to compare that planet to the Earth. Three primary features of An contrast strongly to those of the Earth: (1) the extremely ancient ages of the martian core, mantle, and crust (approx. 4.55 b.y.); (2) the highly depleted nature of the martian mantle; and (3) the extreme ranges of Nd isotopic compositions that arise within the crust and depleted mantle.

Description: We have performed a series of experiments using an ocean model to study the sensitivity of tropical Pacific Ocean to variations in precipitation induced freshwater fluxes. Variations in these fluxes arise from natural causes on all time scales. In addition, estimates of these fluxes are uncertain because of differences among measurement techniques. The model used is a quasi-isopycnal model, covering the Pacific from 40 S to 40 N. The surface forcing is constructed from observed wind stress, evaporation, precipitation, and surface temperature (SST) fields. The heat flux is produced with an iterative technique so as to maintain the model close to the observed climatology, but with only a weak damping to that climatology. Climatological estimates of evaporation are combined with various estimates of precipitation to determine the net surface freshwater flux. Results indicate that increased freshwater input decreases salinity as expected, but increases temperatures in the upper ocean. Using the freshwater flux estimated from the Microwave Sounding Unit leads to a warming of up to 0.6 C in the western Pacific over a case with zero net freshwater flux. SST is sensitive to the discrepancies among different precipitation observations, with root-mean-square differences in SST on the order of 0.2-0.3 C. The change in SST is more pronounced in the eastern Pacific, with differences of over 1 C found among the various precipitation products. Interannual variation in precipitation during El Nino events leads to increased warming. During the winter of 1982-83, freshwater flux accounts for about 0.4 C (approximately 10-15% of the maximum warming) of the surface warming in the central-eastern Pacific. Thus, the error of SST caused by the discrepancies in precipitation products is more than half of the SST anomaly produced by the interannual variability of observed precipitation. Further experiments, in which freshwater flux anomalies are imposed in the western, central, and eastern Pacific, show that the influence of net freshwater flux is also spatially dependent. The imposition of freshwater flux in the far western Pacific leads to a trapping of salinity anomaly to the surface layers near the equator. An identical flux imposed in the central Pacific produces deeper and off-equatorial salinity anomalies. The contrast between these two simulations is consistent with other simulations of the western Pacific barrier layer information.

Description: Surface height, sea surface temperature and surface wind observations from satellites have given a detailed time sequence of the initiation and evolution of the 1997/98 El Nino. The data have beet complementary to the subsurface TAO moored data in their spatial resolution and extent. The impact of satellite observations on seasonal prediction in the tropical Pacific using a coupled ocean-atmosphere general circulation model will be presented.

Description: Traditionally, data from random ocean waves are presented in spectral functions. The spectra are the result of Fourier analysis. Fourier spectral analysis has dominated data analysis for, at least, the last hundred years. It has been the standard method for is examining the global amplitude-frequency distributions. Although Fourier transform valid under extremely general conditions, there are some crucial restrictions for the Fourier spectral analysis. The system must be linear, and the data must be stationary- otherwise, the resulting spectrum will make little physical sense. The stationarity requirement is also a common required criterion for most of other available data analysis methods. Nevertheless, few, if any, natural phenomena are linear and stationary. To compound these complications is the imperfection of our probes or numerical schemes the interactions of the imperfect probes even with a perfect linear system can make the final data nonlinear. Furthermore, all the available data are usually of finite duration. Under these conditions, Fourier analysis is of limited use, For lack of alternatives, however, Fourier analysis is still used to process such data. The loose application of Fourier analysis and the insouciant adoption of the stationary and linear assumptions may lead to misleading conclusions. Ocean waves are know to be nonlinear, and the wind system generating the wave field are seldom stationary- As a result, the traditional examination of the spectral form hardly made physical sense. A new method for analyzing nonlinear and nonstationary data has been developed. The key part is the Empirical Mode Decomposition (EMD) method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF) that serve as the basis of the representation of the data, This decomposition method is adaptive, and, therefore, highly efficient. The IMFs admit well-behaved Hilbert transforms, and yield instantaneous energy and frequency as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Among the main conceptual innovations is the introduction of the instantaneous frequencies for complicated data sets, which eliminate the need of spurious harmonics to represent nonlinear and nonstationary signals. The spectral forms of the ocean waves are very different. This new method offers the first physical meaningful view of the wave spectrum. Data from laboratory and field will be presented to illustrate the differences.

Description: The goal of this research is to monitor the health and vigor of coral reef ecosystems, and their sensitivity to natural and anthropogenic climate changes. To achieve these lofty goals, this research is investigating the feasibility of using spaceborne high-resolution spectrometers (on the US Landsat, French Systeme Probatoire pour l'Observation de la Terre [SPOT] and/or the Indian Resources Satellite [IRS 1C & 1D] spacecraft) to first map the aerial extent of coral reef systems, and second separate the amount of particular corals. If this is successful, we could potentially provide a quantum leap in our understanding of coral reef systems, as well as provide much needed baseline data to measure future changes in global coral reef ecosystems. In collaboration with Tomas Tomascik, Yann Morel, and other colleagues, a series of experiments were planned to coordinate in situ coral observations, high-resolution spaceborne imagery (from Landsat, SPOT, and, possibly, IRS IC spacecraft), and NASA Space Shuttle photographs and digital images. Our eventual goal is to develop "coral health algorithms" that can be used to assess time series of imagery collected from satellite sensors (Landsat since 1972, SPOT since 1986) in concert with in situ observations. The bad news from last year was that from 1997 to mid- 1998, the extreme cloudiness over southeast Asia due to prolonged smoke from El Nino-related fires and the economic chaos in this region frustrated both our space and reef-based data collection activities. When this volatile situation stabilizes, we will restart these activities. The good news was that in collaboration with Al Strong at the National Oceanic and Atmospheric Administration (NOAA) we had an exciting year operationally using the NOAA's Advanced Very High Resolution Radiometer sensor derived sea surface temperature products to warn of coral "bleaching" at many locations throughout the tropics. Data from NOAA's satellites showed that during the El Nino of 1997 and the first half of 1998, more ocean area in the tropics experienced exceptionally high sea surface temperatures, or "hot spots," than have been observed in any fall year since the El Nino of 1982. From January to July, the coral bleaching events were concentrated in the Southern Hemisphere (during its warm season). Since July 1998, the reports of extensive coral bleaching have again spread into regions of the Northern Hemisphere following abnormally high sea surface temperatures, especially around the Philippines and throughout the Caribbean Basin, Bahamas, Bermuda and Florida Keys. These El Nino induced events clearly demonstrated that corals are the "canaries of the marine ecosystem," highly sensitive to short-term natural climate events (El Nino), and should be monitored as measures of longer-term environmental and climate change. Additional information is contained in the original.

Description: An important term in the mass balance of the Arctic Ocean sea ice is the ice export. We estimated the winter sea ice export through the Fram Strait using ice motion from satellite passive microwave data and ice thickness data from moored upward looking sonars. The average winter area flux over the 18-year record (1978-1996) is 670,000 square km, approximately 7% of the area of the Arctic Ocean. The winter area flux ranges from a minimum of 450,000 sq. km in 1984 to a maximum of 906,000 sq km in 1995. The daily, monthly and interannual variabilities of the ice area flux are high. There is an upward trend in the ice area flux over the 18-year record. The average winter volume flux over the winters of October 1990 through May 1995 is 1745 cubic km ranging from a low of 1375 cubic km in 1990 to a high of 2791 cubic km in 1994. The sea-level pressure gradient across the Fram Strait explains more than 80% of the variance in the ice flux over the 18-year record. We use the coefficients from the regression of the time-series of area flux versus pressure gradient across the Fram Strait and ice thickness data to estimate the summer area and volume flux. The average 12-month area flux and volume flux are 919,000 sq km and 2366 cubic km. We find a significant correlation (R =0.86) between the area flux and positive phases of the North Atlantic Oscillation (NAO) index over the months of December through March. Correlation between our six years of volume flux estimates and the NAO index gives R =0.56. During the high NAO years, a more intense Icelandic low increases the gradient in the sea-level pressure by almost 1 mbar across the Fram Strait thus increasing the atmospheric forcing on ice transport. Correlation is reduced during the negative NAO years because of decreased dominance of this large-scale atmospheric pattern on the sea-level pressure gradient across the Fram Strait. Additional information is contained in the original.

Description: As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Description: A data assimilation is carried out to detect the variability of the Japan Sea circulation in the range from a few days to several years and from eddy to basin scale. The model applied in this study is the same 1/6 degree GFDL MOM1 as Kim and Yoon (1999) but is driven by ECMWF daily wind stress, heat and fresh water fluxes. The satellite altimeter data of TOPEX/POSEIDON, ERS-1 (phase C and G) and -2 are assimilated by an approximate Kalman filter (Fukumori and M.Rizzoli, 1995). The approximation is made by seeking an asymptotic steady error covariance matrix (Fukumori et al., 1993) and by introducing a coarser grid model for the innovation (data-model misfit). The coarse grid model is defined on 1/2 degree horizontal resolution and consists of the barotropic stream function, first baroclinic displacement and velocity amplitudes. The assimilated estimates explain about 6cm sea level variability of the data (approximately 12cm in the southern part), which is much larger than the previous reduced-gravity model and TOPEX/POSEIDON altimeter data assimilation (Hirose et al., 1999). The impacts of the T/P and ERS data on the filtered estimates are comparable. The result also shows high correlation to subsurface water temperatures measured by CTD. Many of the mesoscale eddies/disturbances travel east-northeastward with the advection speed of 1-3cm/s though most of them generated in the western region can not pass over the Oki Spur. The quasi-biennial variability found by Hirose and Ostrovskii (1999) did not show clear propagation pattern. The shallow Oki Spur may work as a "western boundary" to this signal. This is more plausible estimation than by the R-G model which has no bottom topography.

Description: The variation in the sun output between the maximum and minimum activity over a 11-yr solar cycle is only about 0.1% of the average. Some people believe that even such a small variability in the sun's output is enough to cause our climate change. However, some others have thrown doubts upon such a relationship. Some of the arguments related to this controversial issue can be found in an earlier report by National Research Council (1982). There is a new wave of study about possible solar-climate relationship in the past decade, and the controversial discussion continues. One group still believes in the importance of the sun's role in climate change, while another still discount it with different degrees.

Description: Maintenance of the large-scale thermohaline circulation has long presented a problem to oceanographers. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon (T/P) satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Description: A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.

Description: Satellite altimetry, particularly Topex/Poseidon, has markedly improved our knowledge of global tides, thereby allowing significant progress on some longstanding problems in geodynamics. This paper reviews some of that progress. Emphasis is given to global-scale problems, particularly those falling within the mandate of the new IERS Special Bureau for Tides: angular momentum, gravitational field, geocenter motion. For this discussion I use primarily the new ocean tide solutions GOT99.2, CSR4.0, and TPXO.4 (for which G. Egbert has computed inverse-theoretic error estimates), and I concentrate on new results in angular momentum and gravity and their solid-earth implications. One example is a new estimate of the effective tidal Q at the M_2 frequency, based on combining these ocean models with tidal estimates from satellite laser ranging. Three especially intractable problems are also addressed: (1) determining long-period tides in the Arctic [large unknown effect on the inertia tensor, particularly for Mf]; (2) determining the global psi_l tide [large unknown effect on interpretations of gravimetry for the near-diurnal free wobble]; and (3) determining radiational tides [large unknown temporal variations at important frequencies]. Problems (2) and (3) are related.

Description: For the Southern Ocean Waves Experiment (SOWEX), conducted in June 1992 out of Hobart, Tasmania, the 36 GHz (8.3 mm) NASA Scanning Radar Altimeter (SRA) was shipped to Australia and installed on a CSIRO Fokker F-27 research aircraft instrumented to make comprehensive surface layer measurements of air-sea interaction fluxes. The sea surface mean square slope (mss), which is predominantly caused by the short waves, was determined from the backscattered power falloff with incidence angle measured by the SRA in the plane normal to the aircraft heading. On each flight, data were acquired at 240 m altitude while the aircraft was in a 7 deg roll attitude, interrogating off-nadir incidence angles from -15 deg through nadir to +29 deg. The aircraft turned azimuthally through 810 deg in this attitude, mapping the azimuthal dependence of the backscattered power falloff with incidence angle. Two sets of turning data were acquired on each day, before and after the aircraft measured wind stress at low altitude (12 m to 65 m). Wave topography and backscattered power for mss were also acquired during those level flight segments whenever the aircraft altitude was above the SRA minimum range of 35 m. A unique feature of this experiment was the use of a nadir-directed low-gain horn antenna (35 deg beamwidth) to acquire azimuthally integrated backscattered power data versus incidence angle before and after the turn data.

Description: The GEOS-1 general circulation model has been used to compute atmospheric torques on the oceans and solid Earth for the period 1980-1995. The time series for the various torque components have been analyzed by means of Fourier transform techniques. It was determined that the wind stress torque over land is more powerful than the wind stress torque over water by 55\%, 42\%, and 80\t for the x, y, and z components respectively. This is mainly the result of power in the high frequency range. The pressure torques due to polar flattening, equatorial ellipticity, marine geoid, and continental orography were computed. The orographic or "mountain torque" components are more powerful than their wind stress counterparts (land plus ocean) by 231\% (x), 191\% (y), and 77\% (z). The marine pressure torques due to geoidal undulations are much smaller than the orographic ones, as expected. They are only 3\% (x), 4\% (y), and 5\% (z) of the corresponding mountain torques. The geoidal pressure torques are approximately equal in magnitude to those produced by the equatorial ellipticity of the Earth. The pressure torque due to polar flattening makes the largest contributions to the atmospheric'torque budget. It has no zonal component, only equatorial ones. Most of the power of the latter, between 68\% and 69 %, is found in modes with periods under 15 days. The single most powerful mode has a period of 361 days. The gravitational torque ranks second in power only to the polar flattening pressure torque. Unlike the former, it does produce a zonal component, albeit much smaller (1\ ) than the equatorial ones. The gravitational and pressure torques have opposite signs, therefore, the gravitational torque nullifies 42\% of the total pressure torque. Zonally, however, the gravitational torque amounts to only 6\% of the total pressure torque. The power budget for the total atmospheric torque yields 7595 and 7120 Hadleys for the equatorial components and 966 Hadleys for the zonal. The x-component exhibits a large mean value (1811 H), mainly the result of polar flattening pressure torque acting on the ocean surfaces. Atmospheric torque modes with periods of 408, 440, and 476 days appear in the spectrum of the equatorial components.

Description: Changes in the ocean angular momentum (OAM) components about the equatorial axes, either due to fluctuations in currents or bottom pressure (mass redistribution), can induce movements of the Earth's pole of rotation, commonly referred to as polar motion or wobble. Output from a 1 deg resolution ocean model is used to calculate the effective equatorial OAM functions chi(sub 1, sup O) and chi(sub 2, sup O), corresponding to polar motion excitation about the equatorial axis pointing to the Greenwich and 90 deg E meridians, respectively. Time series of chi(sup O) are combined with similar atmospheric series chi(sup A), computed from the National Centers for Environmental Prediction/National Center for Atmospheric Research reanalyses, to interpret the observed low-frequency polar motion excitation for the period 1985-1996. Results indicate that the oceans are a very important excitation source for the Chandler (approximately 433 days), annual, and semiannual wobbles, providing for much better amplitude and phase agreement with the observed excitation at these periods, in comparison with what is obtained when only the atmosphere is considered. Both oceanic mass and motion terms are found to be important but with mass signals having somewhat larger amplitudes. The role of regional variability in ocean currents and bottom pressure in contributing to chi(sup O) signals is quantified. Midlatitude regions (approximately 30 deg - 70 deg) figure prominently as places of strong local oceanic excitation signals. The North Pacific basin is found to be generally important for chi(sup O) excitation, while the Southern Ocean is important for both chi(sub 1, sup O) and chi(sub 2, sup O). The largest positive covariances of local with global chi(sup O) signals occur in the Kuroshio region near the western boundary of the North Pacific for chi(sub 1, sup O) and southwest of Australia for chi(sub 2, sup O).

Description: Bottom sediment resuspension frequency, duration and extent (% of bottom sediments affected) were characterized for the fifteen month period from September 1995 to January 1997 for the Barataria Basin, LA. An empirical model of sediment resuspension as a function of wind speed, direction, fetch and water depth was derived from wave theory. Water column turbidity was examined by processing remotely sensed radiance information from visible and near-IR AVHRR imagery. Based on model predictions, wind induced resuspension occurred during all seasons of this study. Seasonal characteristics for resuspension reveal that late fall, winter and early spring are the periods of most frequent and intense resuspension. Model predictions of the critical wind speed required to induce resuspension indicate that winds of 4 m/s (averaged over all wind directions resuspend approximately 50% of bottom sediments in the water bodies examined. Winds of this magnitude (4 m/s) occurred for 80% of the time during the late fall, winter and early spring and for approximately 30% of the time during the summer. More than 50% of the bottom sedimets are resuspended throughout the year, indicating the importance of resuspension as a process affecting sediment and biogeochemical fluxes in the Barataria Basin.

Description: Changes in ocean angular momentum about the polar axis (M) are related to fluctuations in zonal currents (relative component M(sub r)) and latitudinal shifts in mass (planetary component M(sub Omega)). Output from a 1 deg ocean model is used to calculate global M(sub r), M(sub Omega), and M time series at 5-day intervals for the period January 1985-April 1996. The annual cycle in M(sub r), M(sub Omega), and M is larger than the semiannual cycle, and M(sub Omega) amplitudes are nearly twice those of M(sub r). Year-to-year modulation of the seasonal cycle is present, but interannual variability is weak. The spectrum of M is red (background slope between omega(sup (-1) and omega(sup -2)) at subseasonal periods, implying a white or blue spectrum for the external torque on the ocean. Comparisons with previous studies indicate the importance of direct atmospheric forcing in inducing subseasonal M signals, relative to instabilities and other internal sources of rapid oceanic signals. Regional angular momentum estimates show that seasonal variability tends to be larger at low latitudes but there are many local maxima due to the spatial structure of zonal current and mass variability. At seasonal timescales, latitudes approximately 20 S - 10 N contribute substantial variability to M(sub Omega), while signals in M(sub r) can be traced to Antarctic Circumpolar Current transports and associated circulation. Variability in M is found to be small when compared with similar time series for the atmosphere and the solid Earth, but ocean signals are significantly coherent with atmosphere-solid Earth residuals, implying a measurable oceanic impact on length-of-day variations.

Description: The solar radiation-based calibration (SRBC) of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) was performed on 1 November 1993. Measurements were made outdoors in the courtyard of the instrument manufacturer. SeaWiFS viewed the solar irradiance reflected from the sensor's diffuser in the same manner as viewed on orbit. The calibration included measurements using a solar radiometer designed to determine the transmittances of principal atmospheric constituents. The primary uncertainties in the outdoor measurements are the transmission of the atmosphere and the reflectance of the diffuser. Their combined uncertainty is about 5 or 6%. The SRBC also requires knowledge of the extraterrestrial solar spectrum. Four solar models are used. When averaged over the responses of the SeaWiFS bands, the irradiance models agree at the 3.6% level, with the greatest difference for SeaWiFS band 8. The calibration coefficients from the SRBC are lower than those from the laboratory calibration of the instrument in 1997. For a representative solar model, the ratios of the SRBC coefficients to laboratory values average 0.962 with a standard deviation of 0.012. The greatest relative difference is 0.946 for band 8. These values are within the estimated uncertainties of the calibration measurements. For the transfer-to-orbit experiment, the measurements in the manufacturer's courtyard are used to predict the digital counts from the instrument on its first day on orbit (August 1, 1997). This experiment requires an estimate of the relative change in the diffuser response for the period between the launch of the instrument and its first solar measurements on orbit (September 9, 1997). In relative terms, the counts from the instrument on its first day on orbit averaged 1.3% higher than predicted, with a standard deviation of 1.2% and a greatest difference of 2.4% or band 7. The estimated uncertainty for the transfer-to-orbit experiment is about 3 or 4%.

Description: This report documents the scientific activities during the first Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Bio-Optical Algorithm Round-Robin (SeaBOARR-98) experiment, which took place from 5-17 July 1998, at the Acqua Alta Oceanographic Tower (AAOT) in the northern Adriatic Sea off the coast of Italy. The ultimate objective of the SeaBOARR activity is to evaluate the effect of different measurement protocols on bio-optical algorithms using data from a variety of field campaigns. The SeaBOARR-98 field campaign was concerned with collecting a high quality data set of simultaneous in-water and above-water radiometric measurements. The deployment goals documented in this report were to: a) use four different surface glint correction methods to compute water-leaving radiances, L W (lambda), from above-water data; b) use two different in-water profiling systems and three different methods to compute L W (lambda) from in-water data (one making measurements at a fixed distance from the tower, 7.5 m, and the other at variable distances up to 29 m away); c) use instruments with a common calibration history to minimize intercalibration uncertainties; d) monitor the calibration drift of the instruments in the field with a second generation SeaWiFS Quality Monitor (SQM-II), to separate differences in methods from changes in instrument performance; and e) compare the L W (lambda) values estimated from the above-water and in-water measurements. In addition to describing the instruments deployed and the data collected, a preliminary analysis of the data is presented, and the kind of follow-on work that is needed to completely assess the estimation of L W (lambda) from above-water and in-water measurements is discussed.

Description: A number of ocean models of different complexity have been used to study changes in the oceanic mass field and angular momentum and their relation to the variable Earth rotation and gravity field. Time scales examined range from seasonal to a few days. Results point to the importance of oceanic signals in driving polar motion, in particular the Chandler and annual wobbles. Results also show that oceanic signals have a measurable impact on length-of-day variations. Various circulation features and associated mass signals, including the North Pacific subtropical gyre, the equatorial currents, and the Antarctic Circumpolar Current play a significant role in oceanic angular momentum variability.

Description: The purpose of this technical memorandum is to provide current documentation of the Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies (SIMBIOS) Project activities, NASA Research Announcement (NRA) research status, satellite data processing, data product validation, and field calibration. This documentation is necessary to ensure that critical information is related to the scientific community and NASA management. This critical information includes the technical difficulties and challenges of combining ocean color data from an array of independent satellite systems to form consistent and accurate global bio-optical time series products. This technical report is not meant as a substitute for scientific literature. Instead, it will provide a ready and responsive vehicle for the multitude of technical reports issued by an operational project.

Description: We develop and test a 1-point closure turbulence model with the following features: 1) we include the salinity field and derive the expression for the vertical turbulent diffusivities of momentum K(sub m) , heat K(sub h) and salt K(sub s) as a function of two stability parameters: the Richardson number R(sub i) (stratification vs. shear) and the Turner number R(sub rho) (salinity gradient vs. temperature gradient). 2) to describe turbulent mixing below the mixed layer (ML), all previous models have adopted three adjustable "background diffusivities" for momentum, heat and salt. We propose a model that avoids such adjustable diffusivities. We assume that below the ML, the three diffusivities have the same functional dependence on R( sub i) and R(sub rho) as derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking.measured by Gargett et al. The procedure frees the model from adjustable background diffusivities and indeed we employ the same model throughout the entire vertical extent of the ocean. 3) in the local model, the turbulent diffusivities K(sub m,h,s) are given as analytical functions of R(sub i) and R(sub rho). 5) the model is used in an O-GCM and several results are presented to exhibit the effect of double diffusion processes. 6) the code is available upon request.

Description: Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable underestimate of the extent of turbulent mixing and thus to an incorrect mixed layer depth. The predicted temperature and salinity profiles (vs. depth) are presented and compared with those of the Kolmogorov-Petruvsky-Piskunuv (KPP) model and Levitus data.

Description: Goddard Ocean Tide model GOT99.2 is a new solution for the amplitudes and phases of the global oceanic tides, based on over six years of sea-surface height measurements by the TOPEX/POSEIDON satellite altimeter. Comparison with deep-ocean tide-gauge measurements show that this new tidal solution is an improvement over previous global models, with accuracies for the main semidiurnal lunar constituent M2 now below 1.5 cm (deep water only). The new solution benefits from use of prior hydrodynamic models, several in shallow and inland seas as well as the global finite-element model FES94.1. This report describes some of the data processing details involved in handling the altimetry, and it provides a comprehensive set of global cotidal charts of the resulting solutions. Various derived tidal charts are also provided, including tidal loading deformation charts, tidal gravimetric charts, and tidal current velocity (or transport) charts. Finally, low-degree spherical harmonic coefficients are computed by numerical quadrature and are tabulated for the major short-period tides; these are useful for a variety of geodetic and geophysical purposes, especially in combination with similar estimates from satellite laser ranging.

Description: The NOAA/NASA Pathfinder program was created by the Earth Observing System (EOS) Program Office to determine how existing satellite-based data sets can be processed and used to study global change. The data sets are designed to be long time-series data processed with stable calibration and community consensus algorithms to better assist the research community. The Ocean Altimeter Pathfinder Project involves the reprocessing of all altimeter observations with a consistent set of improved algorithms, based on the results from TOPEX/POSEIDON (T/P), into easy-to-use data sets for the oceanographic community for climate research. Details are currently presented in two technical reports: Report# 1: Data Processing Handbook Report #2: Data Set Validation This report describes the validation of the data sets against a global network of high quality tide gauge measurements and provides an estimate of the error budget. The first report describes the processing schemes used to produce the geodetic consistent data set comprised of SEASAT, GEOSAT, ERS-1, TOPEX/ POSEIDON, and ERS-2 satellite observations.

Description: Seasonal and long-term sediment transport characteristics were examined using surficial sediment 7-Be inventories and the down core distribution of excess 210-Pb. Data were collected in the Barataria Basin, LA over the fifteen month period from September 1995 to January 1997. Seasonal sediment transport rates based on 7-Be inventories ranged from -1.6E3 to 1.42E4 g/m2/yr, whereas decadal sediment burial rates based on excess 210-Pb ranged from 3.83E2 to 2.00E3 g/m2/yr, respectively. Seasonal transport characteristics vary with location in the basin and appear to be largely controlled by seasonal weather patterns and the associated winds. It appears that, at less sheltered locations, long term rates of sediment burial are controlled by frontal passages and the associated strong northerly and southerly component winds; whereas at fetch limited locations burial rates are likely controlled by stronger weather events such as tropical storms and hurricanes.

Description: This report documents the fifth Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Intercalibration Round-Robin Experiment (SIRREX-5), which was held at the National Institute of Standards and Technology (NIST) on 23-30 July 1996. The agenda for SIRREX-5 was established based on recommendations made during SIRREX-4. For the first time in a SIRREX activity, instrument intercomparisons were performed at field sites, which were near NIST. The goals of SIRREX-5 were to continue the emphasis on training and the implementation of standard measurement practices, investigate the calibration methods and measurement chains in use by the oceanographic community, provide opportunities for discussion, and intercompare selected instruments. As at SIRREX-4, the day was divided between morning lectures and afternoon laboratory exercises. A set of core laboratory sessions were performed: 1) in-water radiant flux measurements; 2) in-air radiant flux measurements; 3) spectral radiance responsivity measurements using the plaque method; 4) device calibration or stability monitoring with portable field sources; and 5) various ancillary exercises designed to illustrate radiometric concepts. Before, during, and after SIRREX-5, NIST calibrated the SIRREX-5 participating radiometers for radiance and irradiance responsivity. The Facility for Automated Spectroradiometric Calibrations (FASCAL) was scheduled for spectral irradiance calibrations for standard lamps during SIRREX-5. Three lamps from the SeaWiFS community were submitted and two were calibrated.

Description: In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Description: The apparent color of the ocean is determined by the interactions of incident light with substances or particles present in the water. The most significant constituents are free-floating photosynthetic organisms (phytoplankton) and inorganic particulates. Phytoplankton contain chlorophyll, which absorbs light at blue and red wavelengths and transmits in the green. Particulate matter can reflect and absorb light, which reduces the clarity (light transmission) of the water. Substances dissolved in water can also affect its color. Observations of ocean color from space, utilizing sensors specially designed to detect the small amount of light radiating from the sea surface, provide a global picture of the patterns of biological productivity in the world's oceans. For that reason, ocean color remote sensing data is a vital resource for biological oceanography. Unlike the limited area of the ocean that can be investigated from a research ship, data from a satellite sensor covers a large region and provides a comprehensive view of the marine environment.

Description: A coupled general ocean circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. The model is driven by climatological meteorological conditions, cloud cover, and sea surface temperature. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability, and the interactions among three functional phytoplankton groups (diatoms, chorophytes, and picoplankton) and three nutrient groups (nitrate, ammonium, and silicate). Phytoplankton groups are initialized as homogeneous fields horizontally and vertically, and allowed to distribute themselves according to the prevailing conditions. Basin-scale model chlorophyll results are in very good agreement with CZCS pigments in virtually every global region. Seasonal variability observed in the CZCS is also well represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are also in good conformance, although occasional departures are apparent. Agreement of nitrate distributions with in situ data is even better, including seasonal dynamics, except for the equatorial Atlantic. The good agreement of the model with satellite and in situ data sources indicates that the model dynamics realistically simulate phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization, and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.

Description: The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) was originally calibrated by the instrument's manufacturer, Santa Barbara Research Center (SBRC), in November 1993. In preparation for an August 1997 launch, the SeaWiFS Project and the National Institute of Standards and Technology (NIST) undertook a second calibration of SeaWiFS in January and April 1997 at the facility of the spacecraft integrator, Orbital Sciences Corporation (OSC). This calibration occurred in two phases, the first after the final thermal vacuum test, and the second after the final vibration test of the spacecraft. For the calibration, SeaWiFS observed an integrating sphere from the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) at four radiance levels. The spectral radiance of the sphere at these radiance levels was also measured by the SeaWiFS Transfer Radiometer (SXR). In addition, during the calibration, SeaWiFS and the SXR observed the sphere at 16 radiance levels to determine the linearity of the SeaWiFS response. As part of the calibration analysis, the GSFC sphere was also characterized using a GSFC spectroradiometer. The 1997 calibration agrees with the initial 1993 calibration to within +/- 4%. The new calibration coefficients, computed before and after the vibration test, agree to within 0.5%. The response of the SeaWiFS channels in each band is linear to better than 1%. In order to compare to previous and current methods, the SeaWiFS radiometric responses are presented in two ways: using the nominal center wave-lengths for the eight bands; and using band-averaged spectral radiances. The band-averaged values are used in the flight calibration table. An uncertainty analysis for the calibration coefficients is also presented.

Description: The primary objective of the Sea-viewing, Wide Field-of-view Sensor (SeaWiFS) Project is to produce water- leaving radiances with an uncertainty of 5% in clear-water regions and chlorophyll a concentrations within +/- 35% over the range of 0.05-50 mg/cu m. Any global mission, like SeaWiFS, requires validation data be submitted from a wide variety of investigators which places a significant challenge on quantifying the total uncertainty associated with the in situ measurements, because each investigator follows slightly different practices when it comes to implementing all of the steps associated with collecting field data, even those with a prescribed set of protocols. This study uses data from multiple cruises to quantify the uncertainties associated with implementing data collection procedures while utilizing differing in-water optical instruments and deployment methods. A comprehensive approach is undertaken and includes: (1) the use of a portable light source and in-water intercomparisons to monitor the stability of the field radiometers, (2) alternative methods for acquiring reference measurements, and (3) different techniques for making in-water profiles. The only system to meet the 5% radiometric objective of the SeaWiFS Project was a free-fall profiler using (relatively inexpensive) modular components, although a more sophisticated (and comparatively expensive) profiler using integral components was very close and only 1% higher. A relatively inexpensive system deployed with a winch and crane was also close, but the ship shadow contamination problem increased the total uncertainty to approximately 6.5%.

Description: The purpose of this series of technical reports is to provide current documentation of the Sensor Intercomparison and Merger for Biological and Interdisciplinary Ocean Studies (SIMBIOS) Project activities, NASA Research Announcement (NRA) research status, satellite data processing, data product validation and field calibration. This documentation is necessary to ensure that critical information is related to the scientific community and NASA management. This critical information includes the technical difficulties and challenges of combining ocean color data from an array of independent satellite systems to form consistent and accurate global bio-optical time series products. This technical report is not meant to substitute for scientific literature. Instead, it will provide a ready and responsive vehicle for the multitude of technical reports issues by an operational project.

Description: It has recently been demonstrated that the GPS signal which has reflected from the ocean surface contains useful geophysical data from which the sea surface wind speed and other parameters can be extracted. This can be used for remote sensing, similar to present day use of radar altimeters or scatterometers, but with significantly smaller instrumentation because of the utilization of the existing GPS broadcast signal for illumination. Several campaigns of aircraft experimentation have been completed demonstrating this technique and reflected GPS data has been reliably collected from 25 km altitude on a balloon. However, there has not yet been a demonstration that the reflected GPS signal can be detected from orbit with sufficient signal to noise ratio (SNR) to make useful remote sensing measurements. A technology demonstration experiment is planned for a Space Shuttle flight in the late 2000 using the Spartan 251 recoverable carrier. This experiment will be also be the first flight validation of the PiVoT GPS receiver developed in house at the Goddard Space Flight Center. The open-architecture, design of this receiver allows the software modifications to be made which control code-correlator spacing to map out the shape of the reflected signal waveform, which is the most basic data product generated by this instrumentation. A moderate gain left-hand circularly polarized antenna, constructed from an array of off-the-shelf hemispherical antennas will be used to give approximately 6dB of additional gain. Preliminary SNR predictions have been done indicating that this antenna will offer sufficient gain to record waveform measurements. A system level description of the experiment instrumentation, including the receiver, antenna and data storage and retrieval will be given. The visibility of GPS reflections over the mission duration of several hours will be studied, including the effects of the limited beamwidth of the antenna. Additionally, because of the small size of this instrumentation, several other space flight opportunities are being considered in the near term. If these missions materialize, they will be described in the paper as well.

Description: Work continues on estimating tilt modulation distortions in the wave topography measured by a scanning radar altimeter. To quantify this effect, a two-dimensional simulation has been performed in the cross-track plane only which assumed that sinusoidal waves of constant wavelength propagate in the cross-track direction (and are infinitely long-crested in the along-track direction). The initial results reported earlier for a Gaussian surface scattering model indicated that when the highest reasonable value of mss is used in the simulation (the Plant limit of 0.08), the nadir values of the ratio of the apparent to actual wave height are the same as for the omnidirectional scattering case. But as the off-nadir angle increased, the apparent wave height increased and became larger than the actual wave height by about 10 degrees off-nadir. And the shorter the wavelength, the larger the apparent wave height increase. This represented a systematic over-estimate of the wave amplitude for waves propagating in the cross-track direction. For lower values of mss the situation worsened. The simulation has been improved by incorporating actual variations of backscattered power with incidence angle measured by the SRA instead of the Gaussian model. The resulting distortion was about half that originally reported. The 3-dimensional simulation to model waves propagating at various azimuthal angles relative to the cross-track plane is still in progress. The results of this model will be verified by comparing them with Scanning Radar Altimeter (SRA) data and optimum correction procedures will be developed. An assessment of the improved 2-dimensional model indicates that the distortion will generally be small in the data taken during the Southern Ocean Waves Experiment (SOWEX) presently be analyzed. But tilt modulation effects are of great concern for SRA data collected during the 1998 hurricane season since the minimum aircraft altitude was 1.5 km and it frequently flew higher. For a given wave length and steepness, the effect increases with increasing aircraft altitude, and is larger for waves propagating in the cross-track direction, which was generally the situation during the hurricane flights.

Description: As the use of satellite data for icesheet research has accelerated, new discoveries have changed the way that we think about ice discharge, and new tools have changed the way that we measure it. We have begun to apply a combination of remote sensing and ice sheet modeling techniques to study ice flow in the recently identified northeast Greenland ice stream. This ice stream shows evidence of organized flow far into the interior and has low-slope areas of rapid flow and regions of enhanced shear that resemble the ice streams of West Antarctica. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. Our investigation is aimed at increasing our understanding of this ice stream and its origins and role in the discharge of ice from the ice sheet and the effects that discharge may have on the history and evolution of the ice in northern Greenland. The ice flow data comes from satellite radar interferometry tied to reference points. We have applied this technique to measure the velocity over the entire ice stream, which is more than 700 km in length. The velocity contours reveal distinct, well-developed shear margins, which are coincident with linear features in the Synthetic Aperture Radar (SAR) amplitude imagery. It is surprising to see such defined, enhanced flow so far inland and so close to the ice divide. The velocity data and topography data are being used to feed a modeling effort that allows us to study the dynamics and character rapid flow. Additional information is contained in the original.

Description: The primary objective of this study is to estimate the flux of heat and freshwater resulting from sea ice melt in the polar seas. The approach taken is to examine the decay of sea ice in the summer months primarily through the use of spaceborne Synthetic Aperture Radar (SAR) imagery. The improved understanding of the dynamics of the melt process can be usefully combined with ice thermodynamic and upper ocean models to form more complete models of ice melt. Models indicate that more heat is absorbed in the upper ocean when the ice cover is composed of smaller rather than larger floes and when there is more open water. Over the course of the summer, floes disintegrate by physical forcing and heating, melting into smaller and smaller sizes. By measuring the change in distribution of floes together with open water over a summer period, we can make estimates of the amount of heating by region and time. In a climatic sense, these studies are intended to improve the understanding of the Arctic heat budget which can then be eventually incorporated into improved global climate models. This work has two focus areas. The first is examining the detailed effect of storms on floe size and open water. A strong Arctic low pressure storm has been shown to loosen up the pack ice, increase the open water concentration well into the pack ice, and change the distribution of floes toward fewer and smaller floes. This suggests episodic melting and the increased importance of horizontal (lateral) melt during storms. The second focus area is related to an extensive ship-based experiment that recently took place in the Arctic called Surface Heat Budget of the Arctic (SHEBA). An icebreaker was placed purposely into the older pack ice north of Alaska in September 1997. The ship served as the base for experimenters who deployed extensive instrumentation to measure the atmosphere, ocean, and ice during a one-year period. My experiment will be to derive similar measurements (floe size, open water, temporal change) using spaceborne SAR data obtained during the summer of 1998, and compare these results with an ocean and ice model of summer melt. Additional information is contained in the original.

Description: Until a decade ago, an often-quoted expression in oceanography is that very few observations are recorded throughout the ocean. Now, the sentiment is no longer valid in the uppermost 10% of the tropical Pacific Ocean nor at the surface of the global ocean. One of the remarkable legacies of the 1985-1994 Tropical Oceans Global Atmosphere (TOGA) Program is an in situ marine meteorological and upper oceanographic measurement array throughout the equatorial Pacific to monitor the development and maintenance of El Nino episodes. The TOGA Observing System, which initially consisted of moored- and drifting-buoy arrays, a network of commercial ships, and coastal and island stations, now includes a constellation of satellites and data-assimilating models to simulate subsurface oceanographic conditions. The El Nino and La Nina tropical Pacific Ocean observing system represents the initial phase of an integrated global ocean observing system. Remarkable improvements have been made in ocean model simulation of subsurface currents, but some problems persist. For example, the simulation of the South Equatorial Current (SEC) remains an important challenge in the 2S-2N Pacific equatorial wave guide. During El Nino the SEC at the equator is reduced and sometimes the direction is reversed, becoming eastward. Both conditions allow warm water stored in the western Pacific to invade the eastern region, creating an El Nino episode. Assimilation of data is a tenet of faith to correct simulation errors caused by deficiencies in surface fluxes (especially wind stress) and parameterizations of subgrid-scale physical processes. In the first of two numerical experiments, the Pacific SEC was simulated with and without assimilation of subsurface temperature data. Along the equator, a very weak SEC occurred throughout the eastern Pacific, independent of assimilation of data. However, as displayed in the diagram, in the western Pacific there was no satisfactory agreement between the two simulations. To help determine reliability of the simulated SEC in the western Pacific, current measurements recorded during the 9-19 October 1994 voyage of the French research vessel L'Atalante are also shown in the diagram. With data assimilation, the simulated SEC was in much better agreement with L'Atalante observations. The simulated SEC with data assimilation was far from perfect, in part because of the sparsity of subsurface temperature observations. In the next experiment, TOPEX/POSEIDON sea surface height data in combination with subsurface temperatures will be assimilated to assess further improvement of the simulation of the SEC.

Description: Combining analysis of satellite data (altimeter, scatterometer, high-resolution visible and infrared images, etc.) with mathematical modeling of non-linear wave processes, we investigate various ocean wave fields (on scales from capillary to planetary), their role in ocean dynamics and turbulent transport (of heat and biogeochemical quantities), and their effects on satellite altimeter measuring accuracy. In 1998 my attention was focused on long internal gravity waves (10 to 1000 km), known also as baroclinic inertia-gravity (BIG) waves. We found these waves to be a major factor of altimeter measurements "noise," resulting in a greater uncertainty [up to 10 cm in terms of sea surface height (SSH) amplitude] in the measured SSH signal than that caused by the sea state bias variations (up to 5 cm or so). This effect still remains largely overlooked by the satellite altimeter community. Our studies of BIG waves address not only their influence on altimeter measurements but also their role in global ocean dynamics and in transport and turbulent diffusion of biogeochemical quantities. In particular, in collaboration with Prof Peter Weichman, Caltech, we developed a theory of turbulent diffusion caused by wave motions of most general nature. Applied to the problem of horizontal turbulent diffusion in the ocean, the theory yielded the effective diffusion coefficient as a function of BIG wave parameters obtainable from satellite altimeter data. This effort, begun in 1997, has been successfully completed in 1998. We also developed a theory that relates spatial fluctuations of scalar fields (such as sea surface temperature, chlorophyll concentration, drifting ice concentration, etc.) to statistical characteristics of BIG waves obtainable from altimeter measurements. A manuscript is in the final stages of preparation. In order to verify the theoretical predictions and apply them to observations, we are now analyzing Sea-viewing Wide Field of view Sensor (SeaWiFS) and Field of view Sensor (SeaWiFS) and Advanced Very High-Resolution Radiometer (AVHRR) data on sea surface temperature (SST) and chlorophyll concentration jointly with TOPEX/POSEIDON data on SSH variations.

Description: The US Navy's GEOSAT Follow-On (GFO) spacecraft was launched on February 10, 1998 and the primary objective of the mission was to map the oceans using a radar altimeter. The spacecraft tracking complement consisted of GPS receivers, a laser retroreflector and Doppler beacons. Since the GPS receivers have not yet returned reliable data, the only means of providing high-quality precise orbits has been though satellite laser ranging (SLR). The spacecraft has been tracked by the international satellite laser ranging network since April 22, 1998, and an average of 7.4 passes per day have been obtained from US and participating foreign stations. Since the predicted radial orbit error due to the gravity field is two to three cm, the largest contributor to the high SLR residuals (7-10 cm RMS for five day arcs) is the mismodelling of the non-conservative forces, not withstanding the development of a three-dimensional eight-panel model and an analytical attitude model for the GFO spacecraft. The SLR residuals show a clear correlation with beta-prime (solar elevation) angle, peaking in mid-August 1998 when the beta-prime angle reached -80 to -90 degrees. In this paper we discuss the tuning of the non-conservative force model, for GFO and report the subsequent addition of the GFO tracking data to the Earth gravity model solutions.

Description: The compiled, quality-controlled satellite multichannel passive-microwave record of polar sea ice now spans over 18 years, from November 1978 through December 1996, and is revealing considerable information about the Arctic sea ice cover and its variability. The information includes data on ice concentrations (percent areal coverages of ice), ice extents, ice melt, ice velocities, the seasonal cycle of the ice, the interannual variability of the ice, the frequency of ice coverage, and the length of the sea ice season. The data reveal marked regional and interannual variabilities, as well as some statistically significant trends. For the north polar ice cover as a whole, maximum ice extents varied over a range of 14,700,000 - 15,900,000 sq km, while individual regions experienced much greater percent variations, for instance, with the Greenland Sea having a range of 740,000 - 1,110,000 sq km in its yearly maximum ice coverage. In spite of the large variations from year to year and region to region, overall the Arctic ice extents showed a statistically significant, 2.80% / decade negative trend over the 18.2-year period. Ice season lengths, which vary from only a few weeks near the ice margins to the full year in the large region of perennial ice coverage, also experienced interannual variability, along with spatially coherent overall trends. Linear least squares trends show the sea ice season to have lengthened in much of the Bering Sea, Baffin Bay, the Davis Strait, and the Labrador Sea, but to have shortened over a much larger area, including the Sea of Okhotsk, the Greenland Sea, the Barents Sea, and the southeastern Arctic.

Description: Maintenance of the large-scale thermohaline circulation has long presented a problem to oceanographers. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Description: Maintenance of the large-scale thermohaline circulation has long presented an interesting problem. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Description: It was recognized prior to the launch of TOPEX/POSEIDON, that the most important source of orbit error other than the gravity field, was due to nonconservative force modelling. Accordingly, an intensive effort was undertaken to study the nonconservative forces acting on the spacecraft using detailed finite element modelling (Antreasian, 1992; Antreasian and Rosborough, 1992). However, this detailed modelling was not suitable for orbit determination, and a simplified eight plate "box-wing" model was developed that took into account the aggregate effect of the various materials and associated thermal properties of each spacecraft surface. The a priori model was later tuned post launch with actual tracking data [Nerem et al., 1994; Marshall and Luthcke, 1994; Marshall et al., 1995]. More recently, Kubitschek (1997] developed a newer box-wing model for TOPEX/POSEIDON, which included updated material properties, accounted for a solar array deflection, and modelled solar array warping due to thermal effects. We have used this updated model as a basis to retune the macromodel for TOPEX/POSEIDON, and report on preliminary results using at least 36 cycles (one year) of SLR and DORIS data in 1993.

Description: 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.

Description: The response of the Greenland ice sheet to climate change could significantly alter sea level. The ice sheet was much thicker at the last glacial maximum. To gain insight into the global change process and the future trend, it is important to evaluate the ice mass variation as a function of time and space. The Gravity Recovery and Climate Experiment (GRACE) mission to fly in 2001 for 5 years will measure gravity changes associated with the current ice variation and the solid earth's response to past variations. Our objective is to assess the separability of different change sources, accuracy and resolution in the mass variation determination by the new gravity data and possible Global Positioning System (GPS) bedrock uplift measurements. We use a reference parameter state that follows a dynamic ice model for current mass variation and a variant of the Tushingham and Peltier ICE-3G deglaciation model for historical deglaciation. The current linear trend is also assumed to have started 5 kyr ago. The Earth model is fixed as preliminary reference Earth model (PREM) with four viscoelastic layers. A discrete Bayesian inverse algorithm is developed employing an isotropic Gaussian a priori covariance function over the ice sheet and time. We use data noise predicted by the University of Texas and JPL for major GRACE error sources. A 2 mm/yr uplift uncertainty is assumed for GPS occupation time of 5 years. We then carry out covariance analysis and inverse simulation using GRACE geoid coefficients up to degree 180 in conjunction with a number of GPS uplift rates. Present-day ice mass variation and historical deglaciation are solved simultaneously over 146 grids of roughly 110 km x 110 km and with 6 time increments of 3 kyr each, along with a common starting epoch of the current trend. For present-day ice thickness change, the covariance analysis using GRACE geoid data alone results in a root mean square (RMS) posterior root variance of 2.6 cm/yr, with fairly large a priori uncertainties in the parameters and a Gaussian correlation length of 350 km. Simulated inverse can successfully recover most features in the reference present-day change. The RMS difference between them over the grids is 2.8 cm/yr. The RMS difference becomes 1.1 cm/yr when both are averaged with a half Gaussian wavelength of 150 km. With a fixed Earth model, GRACE alone can separate the geoid signals due to past and current load fairly well. Shown are the reference geoid signatures of direct and elastic effects of the current trend, the viscoelastic effect of the same trend starting from 5 kyr ago, the Post Glacial Rebound (PGR), and the predicted GRACE geoid error. The difference between the reference and inverse modeled total viscoelastic signatures is also shown. Although past and current ice mass variations are allowed the same spatial scale, their geoid signals have different spatial patterns. GPS data can contribute to the ice mass determination as well. Additional information is contained in the original.

Description: The objectives of this work are to determine the current state of mass balance of the Greenland and Antarctic Ice Sheets. Our approach combines different techniques, which include satellite synthetic-aperture radar interferometry (InSAR), radar and laser altimetry, radar ice sounding, and finite-element modeling. In Greenland, we found that 3.5 times more ice flows out of the northern part of the Greenland Ice Sheet than previously accounted for. The discrepancy between current and past estimates is explained by extensive basal melting of the glacier floating sections in the proximity of the grounding line where the glacier detaches from its bed and becomes afloat in the ocean. The inferred basal melt rates are very large, which means that the glaciers are very sensitive to changes in ocean conditions. Currently, it appears that the northern Greenland glaciers discharge more ice than is being accumulated in the deep interior, and hence are thinning. Studies of temporal changes in grounding line position using InSAR confirm the state of retreat of northern glaciers and suggest that thinning is concentrated at the lower elevations. Ongoing work along the coast of East Greenland reveals an even larger mass deficit for eastern Greenland glaciers, with thinning affecting the deep interior of the ice sheet. In Antarctica, we found that glaciers flowing into a large ice shelf system, such as the Ronne Ice Shelf in the Weddell Sea, exhibit an ice discharge in remarkable agreement with mass accumulation in the interior, and the glacier grounding line positions do not migrate with time. Glaciers flowing rapidly into the Amudsen Sea, unrestrained by a major ice shelf, are in contrast discharging more ice than required to maintain a state of mass balance and are thinning quite rapidly near the coast. The grounding line of Pine Island glacier (see diagram) retreated 5 km in 4 years, which corresponds to a glacier thinning rate of 3.5 m/yr. Mass imbalance is even more negative on Thwaites Glacier. This sector of West Antarctica probably initiated its collapse decades or centuries ago, once the embaying ice shelves in front of them started to melt because of enhanced basal melting from warmer ocean waters. Additional information is contained in the original.

Description: Various data over the period 1980-1998 have been used to improve El Nino Southern Oscillation (ENSO) simulations and forecasts performed by coupled ocean/atmosphere models. For such a study, oceanic and atmospheric satellite data are necessary. Models assume a baroclinic ocean with a mixed layer coupled to a baroclinic atmosphere to simulate the anomalous fields. The various data have been used to revisit the parameterizations of the subsurface temperature in the mixed layer, atmospheric convection, and friction in the baroclinic ocean. This new model, named Tsub.Conv, simulates 4-year oscillations like the model with the original parameterization, but with anomalous thermocline and wind anomalies that agree a lot better with reality. The atmospheric model still presents some deficiencies in the eastern Pacific and when it is replaced by a statistical atmosphere, the coupled model (Tsub.Astat) performs even better. It is found that the off-equatorial ocean and wind curl anomalies play a crucial role for all these models and that Tsub.Astat simulates an ocean heat content recharge coming from the north prior to a warm event like the ones observed during the 1985-1992 period. The ENSO forecasts over 1980-1998 are considerably improved with Tsub.Astat. One needs to add the simulation of westerly wind bursts to be able to forecast the 1997 event one year in advance. Additional information is contained in the original.

Description: Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Description: Ocean modeling plays an important role in understanding the current climatic conditions and predicting the future climate change. Modeling the ocean at eddy-permitting and/or eddy resolving resolutions (1/3 degree or higher) has a two-fold objective. One part is to represent the ocean as realistically as possible, because mesoscale eddies have an impact on the large-scale circulation. The second objective is to learn how to represent effects of mesoscale eddies without explicitly resolving them. This is particularly important for climate models which cannot be run at eddy-resolving resolutions because of the computational constraints. At JPL, a 1/6 degree latitude by 1/6 degree longitude with 37 vertical levels Atlantic Ocean model has been developed. The model is based on the Parallel Ocean Program (POP) developed at Los Alamos National Laboratory (LANL). Using the 256-processor Cray T3D, we have conducted a 40-year integration of this Atlantic eddy-resolving ocean model. A regional analysis demonstrate that many observed features associated with the Caribbean Sea eddies can be realistically simulated by this model. Analysis of this Atlantic eddy-resolving ocean model further suggests that these Caribbean Sea eddies are connected with eddies formed outside the Caribbean Sea at the confluence of the North Brazil Current (NBC) and the North Equatorial Countercurrent. The diagram of the model simulated surface current shows that the Caribbean eddies ultimately originate in the NBC retroflection region, traveling more than a year from the North Brazil coast through the Lesser Antilles into the Caribbean Sea and eventually into the Gulf of Mexico. Additional information is contained in the original.

Description: 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.

Description: We report on the lunar and solar measurements used to determine the changes in the radiometric sensitivity of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Radiometric sensitivity is defined as the output from the instrument (or from one of the instrument bands) per unit spectral radiance at the instrument's input aperture. Knowledge of the long-term repeatability of the SeaWiFS measurements is crucial to maintaining the quality of the ocean scenes derived from measurements by the instrument. For SeaWiFS bands 1 through 6 (412 nm through 670 rim), the change in radiometric sensitivity is less than 0.2% for the period from November 1997 through November 1998. For band 7 (765 nm), the change is about 1.5%, and for band 8 (865 nm) about 5%. The rates of change of bands 7 and 8, which were linear with time for the first eight months of lunar measurements, are now slowing. The scatter in the data points about the trend lines in this analysis is less than 0.3% for all eight SeaWiFS bands. These results are based on monthly measurements of the moon. Daily solar measurements using an onboard diffuser show that the radiometric sensitivities of the SeaWiFS bands have changed smoothly during the time intervals between lunar measurements. Since SeaWiFS measurements have continued past November 1998, the results presented here are considered as a snapshot of the instrument performance as of that date.

Description: The Odden Ice tongue of the Greenland Sea normally forms locally In winter as frazfl-pancake ice, allowing high positive salt fluxes during freezing that leads to open ocean convection. We report observations from satellites, aircraft, ships and submarines which show that in two recent years (1987 and 1996) a late-season Odden developed composed of old ice advected by the East Greenland Current. The Impact of such Odden is different in that it is in a state of melt and serves to stabilize the surface water in the region. The history of Oddens since 1978 is reviewed to examine the frequency of both modes.

Description: The mapping of mesoscale ocean features in the coastal zone is a major potential application for satellite data. The evolution of mesoscale features such as oil slicks, fronts, eddies, and ice edge can be tracked by the wavelet analysis using satellite data from repeating paths. The wavelet transform has been applied to satellite images, such as those from Synthetic Aperture Radar (SAR), Advanced Very High-Resolution Radiometer (AVHRR), and ocean color sensor for feature extraction. In this paper, algorithms and techniques for automated detection and tracking of mesoscale features from satellite SAR imagery employing wavelet analysis have been developed. Case studies on two major coastal oil spills have been investigated using wavelet analysis for tracking along the coast of Uruguay (February 1997), and near Point Barrow, Alaska (November 1997). Comparison of SAR images with SeaWiFS (Sea-viewing Wide Field-of-view Sensor) data for coccolithophore bloom in the East Bering Sea during the fall of 1997 shows a good match on bloom boundary. This paper demonstrates that this technique is a useful and promising tool for monitoring of coastal waters.

Description: The main result of this paper is the derivation of a new expression for the tracer subgrid term in level coordinates S(l) to be employed in O-GCM. The novel feature is the proper account of the random nature of the density field which strongly affects the transformation from isopycnal to level coordinates of the variables of interest, velocity and tracer fields, their correlation functions and ultimately the subgrid terms. In deriving our result we made use of measured properties of vertical ocean turbulence. The major new results are: 1) the new subgrid expression is different from that of the heuristic GM model, 2) u++(tracer)=1/2u+(thickness), where u++ and u+ are the tracer and thickness bolus velocities. In previous models, u++ = u+, 2) the subgrid for a tracer tau is not the same as that for the density rho even when one accounts for the obvious absence of a diffusion term in the latter. The difference stems from a new treatment of the stochastic nature of the density, 3) the mesoscale diffusivity enters both locally and non-locally, as the integral over all z's from the bottom of the ocean to the level z.

Description: The daily mean heat and momentum fluxes at the surface derived from the SSM/I and Japan's GMS radiance measurements are used to study the temporal and spatial variability of the surface energy budgets and their relationship to the sea surface temperature during the COARE intensive observing period (IOP). For the three time legs observed during the IOP, the retrieved surface fluxes compare reasonably well with those from the IMET buoy, RV Moana Wave, and RV Wecoma. The characteristics of surface heat and momentum fluxes are very different between the southern and northern warm pool. In the southern warm pool, the net surface heat flux is dominated by solar radiation which is, in turn, modulated by the two Madden-Julian oscillations. The surface winds are generally weak, leading to a shallow ocean mixed layer. The solar radiation penetrating through the bottom of the mixed layer is significant, and the change in the sea surface temperature during the IOP does not follow the net surface heat flux. In the northern warm pool, the northeasterly trade wind is strong and undergoes strong seasonal variation. The variation of the net surface heat flux is dominated by evaporation. The two westerly wind bursts associated with the Madden-Julian oscillations seem to have little effect on the net surface heat flux. The ocean mixed layer is deep, and the solar radiation penetrating through the bottom of the mixed layer is small. As opposed to the southern warm pool, the trend of the sea surface temperature in the northern warm pool during the IOP is in agreement with the variation of the net heat flux at the surface.

Description: Recent papers have described 18-year trends and annual oscillations in the Arctic and Antarctic sea ice extents, areas, and enclosed open water areas based on a newly-formulated 18.2-year ice concentration time series. This time series includes data for the entire Arctic and Antarctic ice covers, as well as for previously defined subregions consisting of 5 sectors in the Antarctic and 9 regions in the Arctic. It was obtained by fine-tuning the sea ice algorithm tie points individually for each of the four sensors used to acquire the data. In this paper, we extend these analyses to an examination of the intrinsic modes of these time series, obtained by means of Empirical Mode Decomposition, with emphasis on periodicities greater than the annual cycle. Quasibiennial and quasiquadrennial oscillations observed with a different technique and reported earlier for the first 8.8 years of this time series were also observed in the present series. However, the intrinsic modes were not monochromatic; they feature frequency as well as amplitude modulation within their respective frequency bands. Modal periods of up to 18 years are observed, with important implications for the trend analyses published earlier. These results are compared with the oscillations in the Length-of-Day and North Atlantic Oscillation parameters similarly determined for the same 18.2-year period.

Description: Despite numerous previous studies, two relationships between deep convection and the sea-surface temperature (SST) of the tropics remain unclear. The first is the cause for the sudden emergence of deep convection at about 28 deg SST, and the second is its proximity to the highest observed SST of about 30 C. Our analysis provides a rational explanation for both by utilizing the Improved Meteorological (IMET) buoy data together with radar rainfall retrievals and atmospheric soundings provided by the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA-COARE). The explanation relies on the basic principles of moist convection as enunciated in the Arakawa-Schubert cumulus parameterization. Our analysis shows that an SST range of 28-29 C is necessary for "charging" the atmospheric boundary layer with sufficient moist static energy that can enable the towering convection to reach up to the 200 hPa level. In the IMET buoy data, the changes in surface energy fluxes associated with different rainfall amounts show that the deep convection not only reduces the solar flux into the ocean with a thick cloud cover, but it also generates downdrafts which bring significantly cooler and drier air into the boundary-layer thereby augmenting oceanic cooling by increased sensible and latent heat fluxes. In this way, the ocean seasaws between a net energy absorber for non-raining and a net energy supplier for deep-convective raining conditions. These processes produce a thermostat-like control of the SST. The data also shows that convection over the warm pool is modulated by dynamical influences of large-scale circulation embodying tropical easterly waves (with a 5-day period) and MJOs (with 40-day period); however, the quasi-permanent feature of the vertical profile of moist static energy, which is primarily maintained by the large-scale circulation and thermodynamical forcings, is vital for both the 28 C SST for deep convection and its upper limit at about 30 C.

Description: A general framework for calculating lunar and solar torques on the oceanic tides is developed in terms of harmonic constituents. Axial torques and their associated angular momentum and earth-rotation variations are deduced from recent satellite-altimeter and satellite-tracking tide solutions. Torques on the prograde components of the tide produce the familiar secular braking of the rotation rate. The estimated secular acceleration is approximately -1300 sec/century(sup 2) (less 4% after including atmospheric tides); the implied rate of change in the length of day is 2.28 milliseconds/century. Torques on the retrograde components of the tide produce periodic rotation variations at twice the tidal frequency. Interaction torques, e.g. solar torques on lunar tides, generate a large suite of rotation-rate variations at sums and differences of the original tidal frequencies. These are estimated for periods from 18.6 years to quarter-diurnal. At subdaily periods the angular momentum variations are 5 to 6 orders of magnitude smaller than the variations caused by ocean tidal currents.

Description: A nine-year simulation (1984-1992) of biological processes in the equatorial Pacific Warm Pool is presented. A modified version of the 4-component (phytoplankton, zooplankton, nitrate and ammonium) ecosystem model by McClain et al. (1996) is used. Modifications include use of a spectral model for computation of PAR and inclusion of fecal pellet remineralization and ammonium nitrification. The physical parameters (horizontal and vertical velocities and temperature) required by the ecosystem model were derived from an improved version of the Gent and Cane (1990) ocean general circulation model (Murtugudde and Busalacchi, 1997). Surface downwelling spectral irradiance was estimated using the clear-sky models of Frouin et al. (1989) and Gregg and Carder (1990) and cloud cover information from the International Satellite Cloud Climatology Project (ISCCP). The simulations indicate considerable variability on interannual time scales in all four ecosystem components. In particular, surface chlorophyll concentrations varied by an order of magnitude with maximum values exceeding 0.30 mg/cu m in 1988, 1989, and 1990, and pronounced minimums during 1987 and 1992. The deep chlorophyll maximum ranged between 75 and 125 meters with values occasionally exceeding 0.40 mg/cu m. With the exception of the last half of 1988, surface nitrate was always near depletion. Ammonium exhibited a subsurface maximum just below the DCM with concentrations as high as 0.5 mg-atN/cu m . Total integrated annual primary production varied between 40 and 250 gC/sq m/yr with an annual average of 140 gC/sq m/yr. Finally, the model is used to estimate the mean irradiance at the base of the mixed layer, i.e., the penetration irradiance, which was 18 Watts/sq m over the nine year period. The average mixed layer depth was 42 m.

Description: One of the primary goals of the NASA Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies (SIMBIOS) project is to develop methods for meaningful comparison and possible merging of data products from multiple ocean color missions. The Modular Optoelectronic Scanner (MOS) is a German instrument that was launched in the spring of 1996 on the Indian IRS-P3 satellite. With the successful launch of NASA's Sea-viewing Wide Field-of-view Sensor (SeaWiFS) in the summer of 1997, there are now two ocean color missions in concurrent operation and there is interest within the scientific community to compare data from these two sensors. In this paper, we describe our efforts to retrieve ocean optical properties from both SeaWiFS and MOS using consistent methods. We first briefly review the atmospheric correction, which removes more than 90% of the observed radiances in the visible, and then describe how the atmospheric correction algorithm used for the SeaWiFS data can be modified for application to other ocean color sensors. Next, since the retrieved water-leaving radiances in the visible between MOS and SeaWiFS are significantly different, we developed a vicarious intercalibration method to recalibrate the MOS spectral bands based on the optical properties of the ocean and atmosphere derived from the coincident SeaWiFS measurements. We present and discuss the MOS retrieved ocean optical properties before and after the vicarious calibration, and demonstrate the efficacy of this approach. We show that it is possible and efficient to vicariously intercalibrate sensors between one and another.

Description: The feasibility of assimilating satellite altimetry data into a global ocean general ocean general circulation model is studied. Three years of TOPEX/POSEIDON data is analyzed using a global, three-dimensional, nonlinear primitive equation model.

Description: The European Center for Medium-Range Weather Forecasts (ECMWF) and Institut Francais Pour la Recherche et 1'Exploitation de la Mer (IFREMER) European Remote-Sensing Satellite (ERS-1), named IFR2, surface wind velocity data products are compared during July 1995 over the Arabian Sea.

Description: Empirical models for the monthly and fortnightly ocean tides are estimated using sea surface height data from repeat cycles 10 to 190 of the TOPEX/POSEIDON (T/P) altimetry mission.

Description: Two issues regarding primary productivity, as it pertains to the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Program and the National Aeronautics and Space Administration (NASA) Mission to Planet Earth (MTPE) are presented in this volume. Chapter 1 describes the development of a science plan for deriving primary production for the world ocean using satellite measurements, by the Ocean Primary Productivity Working Group (OPPWG). Chapter 2 presents discussions by the same group, of algorithm classification, algorithm parameterization and data availability, algorithm testing and validation, and the benefits of a consensus primary productivity algorithm.

Description: This flow chart describes the baseline science algorithms for the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Data Processing System (SDPS). As such, it includes only processing steps used in the generation of the operational products that are archived by NASA's Goddard Space Flight Center (GSFC) Distributed Active Archive Center (DAAC). It is meant to provide the reader with a basic understanding of the scientific algorithm steps applied to SeaWiFS data. It does not include non-science steps, such as format conversions, and places the greatest emphasis on the geophysical calculations of the level-2 processing. Finally, the flow chart reflects the logic sequences and the conditional tests of the software so that it may be used to evaluate the fidelity of the implementation of the scientific algorithm. In many cases however, the chart may deviate from the details of the software implementation so as to simplify the presentation.

Description: Empirical models of the two largest constituents of the long-period ocean tides, the monthly and the fortnightly constituents, are estimated from repeat cycles 10 to 210 of the TOPEX/POSEIDON (T/P) mission.