Full Water Column Research Articles

Ensemble initialization of the oceanic component of a coupled model through bred vectors at seasonal-to-interannual timescales

Abstract. We evaluate the ensemble spread at seasonal-to-interannual timescales for two perturbation techniques implemented in the ocean component of a coupled model: (1) lagged initial conditions as commonly used for decadal predictions; (2) bred vectors as commonly used for weather and seasonal forecasting. We show that relative to an uninitialized reference simulation the implementation for bred vectors can improve the ensemble spread compared to lagged initialization at timescales from one month up to three years. As bred vectors have so far mostly been used at short timescales, we initially focus on the implementation of the bred vectors in the ocean component. We introduce a depth-dependent vertical rescaling norm, accounting for the vertical dependence of the variability, and extending the commonly used upper-ocean rescaling norm to the full water column. We further show that it is sufficient for the (sub-surface) ocean to breed temperature and salinity (i.e., scalar quantities), and rely on the governing physics to carry the temperature and salinity perturbations to the flow field. Using these bred vectors with a rescaling interval of 12 months, we initialize hindcast simulations and compare them to hindcast simulations initialized with lagged initial conditions. We quantify the ensemble spread by analyzing Talagrand diagrams and spread–error ratios. For both temperature and salinity, the lagged initialized ensemble is particularly under-dispersive for the first few months of predictable lead time. The ensemble initialized with bred vectors improves the spread for temperature and salinity for the 0–700 m and 1000–3500 m means, compared to the lagged ensemble at lead times of several months to one year. As the lead time increases to years, the differences between the two ensemble initialization techniques become more difficult to discern. While the results need to be confirmed in an initialized framework, the present analysis represents a first step towards improved ensemble generation at the transition from seasonal to interannual timescales, in particular at lead times up to one year.

Biogeochemical implications from dissolved rare earth element and Nd isotope distributions in the Gulf of Alaska

Dissolved rare earth element (REE) concentrations and Nd isotope compositions were measured for surface waters and full water column profiles of the Gulf of Alaska (GoA), and compared to water mass properties and circulation in order to better understand the mechanisms controlling the input and transport of REEs in the ocean. The REEs display a typical open-ocean range of concentrations (i.e., La: 12–66pM; Lu: 0.2–2.5pM) and depth distributions (i.e., surface ocean depletion and enrichment with water depth). Nd isotope signatures are highly radiogenic, as expected for the North Pacific margin (ranging from −3.8 to +0.2εNd). The most radiogenic values were found in the coastal waters but also in the cores of eddies, indicating efficient export of REEs from the margins and across the mixed layer. This is the first time that distinct Nd isotope distributions in near surface waters can be directly assigned to offshore eddy transport. A distinct mid-depth (∼2200m) Nd isotope signal was found that most likely reflects advection of a water mass that formed through past down-welling in the Northern Pacific. Subsurface Nd isotope compositions appear to behave conservatively and can be explained through a REE distribution model proposed here. This model is based on multivariate analysis of the REEs and invokes two distinct “pools” of dissolved REEs: a “passive pool” complexed by carbonate ions, and a “bio-reactive pool” that is microbially manipulated. The latter “pool” is only significant in the upper water column and most likely reflects the indirect effects of microbial cycling of iron. Our model of the open ocean REE distribution contributes to explaining the conservative nature of Nd isotopes and provides a mechanism linking surface ocean and pore water REE dynamics.

Analysis of the vertical structure of deep ocean noise using measurements from the SPICEX and PhilSea experiments

There are open research questions about the vertical structure of low-frequency ambient noise in deep water. For example Gaul et al.'s [IEEE JOE (2007)] analysis of the Church Opal data set showed that noise decreases substantially (on the order of 20 dB) below the critical depth, whereas other researchers have reported more modest reductions [Morris, J. Acoust. Soc. Am. (1978)]. Two deep water experiments provided a unique opportunity to measure ambient noise using large vertical arrays. In 2004–2005, SPICEX used two arrays to sample a North Pacific environment. One array was centered on the sound channel axis, and the other array had hydrophones above and below the critical depth. In 2010–2011, the PhilSea experiment deployed a single array with 150 hydrophones spanning the full water column. Both experiments made repeated short measurements (each 2–3 min long) of the field at the arrays. This talk compares the ambient noise observed during SPICEX and PhilSea with results reported in the literature. Since these data sets contain receptions over the period of a year, we focus on the seasonal dependence of the noise field. In addition to investigating noise level as a function of depth, we consider wind dependence and vertical directionality.

The distribution of neodymium isotopes and concentrations in the Eastern Equatorial Pacific: Water mass advection versus particle exchange

The radiogenic isotope composition of the Rare Earth Element (REE) neodymium (Nd) is a powerful water mass proxy for present and past ocean circulation. The processes controlling the Nd budget of the global ocean are not quantitatively understood and in particular source and sink mechanisms are still under debate.In this study we present the first full water column data set of dissolved Nd isotope compositions and Nd concentrations for the Eastern Equatorial Pacific (EEP), where one of the globally largest Oxygen Minimum Zones (OMZ) is located. This region is of particular interest for understanding the biogeochemical cycling of REEs because anoxic conditions may lead to release of REEs from the shelf, whereas high particle densities and fluxes potentially remove the REEs from the water column. Data were obtained between 1°40′N and 16°S along a nearshore and an offshore transect. Near surface zonal current bands, such as the Equatorial Undercurrent (EUC) and the Subsurface Countercurrent (SSCC), which are supplying oxygen-rich water to the OMZ are characterized by radiogenic Nd isotope signatures (εNd=−2). Surface waters in the northernmost part of the study area are even more radiogenic (εNd=+3), most likely due to release of Nd from volcanogenic material. Deep and bottom waters at the southernmost offshore station (14°S) are clearly controlled by advection of water masses with less radiogenic signatures (εNd=−7) originating from the Southern Ocean. Towards the equator, however, the deep waters show a clear trend towards more radiogenic values of up to εNd=−2. The northernmost station located in the Panama basin shows highly radiogenic Nd isotope signatures in the entire water column, which indicates that particle scavenging, downward transport and release processes play an important role. This is supported by relatively low Nd concentrations in deep waters (3000–6000m) in the EEP (20pmol/kg) compared to locations in the Northern and Central Pacific (40–60pmol/kg), which suggests enhanced removal of Nd in the EEP.

Emptying of Large-Scale Pipeline by Pressurized Air

AbstractEmptying of an initially water-filled horizontal PVC pipeline driven by different upstream compressed air pressures and with different outflow restriction conditions, with motion of an air-water front through the pressurized pipeline, is investigated experimentally. Simple numerical modeling is used to interpret the results, especially the observed additional shortening of the moving full water column due to formation of a stratified water-air “tail.” Measured discharges, water-level changes, and pressure variations along the pipeline during emptying are compared using control volume (CV) model results. The CV model solutions for a nonstratified case are shown to be delayed as compared with the actual measured changes of flow rate, pressure, and water level. But by considering water-column mass loss due to the water-air tail and residual motion, the calibrated CV model yields solutions that are qualitatively in good agreement with the experimental results. A key interpretation is that the long air...

Autonomous Depth Adjustment for Underwater Sensor Networks: Design and Applications

To fully understand the ocean environment requires sensing the full water column. Utilizing a depth adjustment system on an underwater sensor network provides this while also improving global sensing and communications. This paper presents a depth adjustment system for waters up to 50 m deep that connects to the aquanode sensor network nodes. We performed experiments characterizing and demonstrating the functionality of the depth adjustment system. We discuss the application of this device in improving acoustic communication and also verify the functionality of a decentralized depth adjustment algorithm that optimizes the placement of the nodes for collecting sensing data.

Concept Tests for a New Wire Flying Vehicle Designed to Achieve High Horizontal Resolution Profiling in Deep Water

Abstract Efficiently profiling the water column to achieve both high vertical and horizontal resolution from a moving vessel in deep water is difficult. Current solutions, such as CTD tow-yos, moving vessel profilers, and undulating tow bodies, are limited by ship speed or water depth. As a consequence, it is difficult to obtain oceanographic sections with sufficient resolution to identify many relevant scales over the deeper sections of the water column. This paper presents a new concept for a profiling vehicle that slides up and down a towed wire in a controlled manner using the lift created by wing foils. The wings provide a novel low-power method of propulsion along the cable by using the free stream velocity of the wire moving through the water in similar fashion to a sailboat sailing up wind. Scale model tests show a wide range of achievable profiling glide slopes for tow cable angles between vertical and 45°, and effective isolation of cable strum vibration from the towed vehicle body. The concept is not depth limited and will offer two-dimensional resolution that meets or exceeds current undulating tow bodies over the full water column. Additionally, this system could be used simultaneously with many other deep towed instrument packages to produce complementary datasets.

Tidally-induced sediment transport patterns in the upper Bay of Fundy: A numerical study

The Minas Basin, the eastern end of the Bay of Fundy, is well known for its high tide ranges and strong tidal currents, which can be exploited to extract electricity power. The properties of the tidally-induced sediment transport in the Minas Basin, where significant changes in tidal processes may occur due to a recently proposed tidal power project, have been studied with a three-dimensional hydrodynamic model, an empirical bed load sediment transport model and surface sediment concentrations derived from the remotely-sensed images. The hydrodynamic model was evaluated against independent observational data, which include tidal elevation, tidal current (in the full water column and bottom layer), residual current profile and tidal asymmetry indicators. The evaluation shows that the model is in good agreement with the observations. The sediment transport includes two components, bed load and suspended particulate load. The bed load is calculated using the modelled bottom shear stress and the observed grain size data. The estimated features of bed load transport roughly agree with the observed patterns of the erosion and deposition in the Minas Basin and Cobequid Bay. The transport of the suspended load is estimated using the modelled velocity fields and the surface sediment concentration derived from remote-sensing images. The comparisons between the modelled results and the limited observations illustrate that the observed directions of suspended sediment transport are basically reproduced by the model. The modelled net suspended sediment input into the Minas Basin through Minas Passage is 2.4×10 6 m 3 yr −1, which is comparable to the observed value of 1.6×10 6 m 3 yr −1. The variations of the bed load and the suspended load in space and time are also presented. The total net transport, defined as the mean value of the sum of bed and suspended load transports during the tidal cycle, shows strong spatial variability. The magnitude of the transport flux ranges from 0.1 to 0.2 kg m −1 s −1 in Minas Channel and Minas Passage, 0.1 kg m −1 s −1 in Cobequid Bay, to 0.01 kg m −1 s −1 in the central Minas Basin and Southern Bight. In Minas Channel, the sediment transport follows the structure of the tidal residual circulation, which features a large anticlockwise gyre. The sediment in Minas Passage moves eastward and deposits into the central Minas Basin. However, the sediment from the eastern part of the Basin moves westward and deposits in the central Minas Basin as well. In the Cobequid Bay, sediment moves eastward and deposits in the upper bay.

Resistance, resilience and recovery: aquatic bacterial dynamics after water column disturbance

For lake microbes, water column mixing acts as a disturbance because it homogenizes thermal and chemical gradients known to define the distributions of microbial taxa. Our first objective was to isolate hypothesized drivers of lake bacterial response to water column mixing. To accomplish this, we designed an enclosure experiment with three treatments to independently test key biogeochemical changes induced by mixing: oxygen addition to the hypolimnion, nutrient addition to the epilimnion, and full water column mixing. We used molecular fingerprinting to observe bacterial community dynamics in the treatment and control enclosures, and in ambient lake water. We found that oxygen and nutrient amendments simulated the physical-chemical water column environment following mixing and resulted in similar bacterial communities to the mixing treatment, affirming that these were important drivers of community change. These results demonstrate that specific environmental changes can replicate broad disturbance effects on microbial communities. Our second objective was to characterize bacterial community stability by quantifying community resistance, recovery and resilience to an episodic disturbance. The communities in the nutrient and oxygen amendments changed quickly (had low resistance), but generally matched the control composition by the 10th day after treatment, exhibiting resilience. These results imply that aquatic bacterial assemblages are generally stable in the face of disturbance.

Prokaryotic dynamics and heterotrophic metabolism in a deep convection site of Eastern Mediterranean Sea (the Southern Adriatic Pit)

We report on investigations of prokaryotic abundance, biomass, extracellular enzymatic activity, prokaryotic heterotrophic production and respiration in the full water column (∼1200m) of a deep convection site (the Southern Adriatic Pit), carried out on six cruises in 2006–2008. Prokaryotic abundance (PA) varied vertically and temporally and ranged from 1.2 to 20.4×105cellml−1. Cell volumes, generally increased with depth; the lowest mean cell volume was observed in a period with no active convective process (Feb-07) and the highest in a period of stratification (Jun-08) following the convection process occurred in Feb-08. Prokaryotic biomass decreased with the depth and was related with both seasonal cycles of organic matter and hydrological processes. The picophytoplankton ranged in the upper layer (UL) from 0.089 to 10.71×104cellml−1. Cells were also recorded till 500m depth in Feb-08 and this finding could be linked to water convection occurred in the Southern Adriatic Pit in that month. In UL the variations of enzymatic activities as well as leucine-aminopeptidase/ß-glucosidase ratio showed a seasonal trend probably linked to the productive processes of the photic layer. An inverse relation between alkaline phosphatase activity (APA) and phosphate concentrations was found (APA=0.0003PO4−1.7714, R2=0.333, P<0.05). Generally cell-specific enzymatic activities increased with depth as did cell-specific carbon dioxide production rates, while cell-specific prokaryotic heterotrophic production had an opposite trend. High values of prokaryotic growth efficiency registered in the deep layers in Nov-06 reflected a supply of preformed C transported within the deep water masses. Overall, in 2007 when no convective phenomenon was observed, the variability of prokaryotic metabolism was governed by the seasonal cycle of the organic matter, while in Nov-06 and Jun-08 the dynamics of deep water ventilation influenced the trend along the water column of many microbial parameters. The yearly trophic balance of the study site appeared to move towards autotrophy only in UL, whilst in the whole water column, the prokaryotic carbon demand exceeded POC availability rained down from euphotic zone. This mismatch was balanced by the DOC entrapped in the “younger waters” of new formation that alters the normal flux of the biological pump and fuels the deep marine biota in this area of deep water convection.

Modal depth function estimation using time-frequency analysis

Acoustic propagation in shallow water is characterized by a set of depth-dependent modes, the modal depth functions, which propagate in range according to their horizontal wavenumbers. For inversion purposes, modal depth function estimation in shallow water is an issue when the environment is not known. Classical methods that provide blind mode estimation rely on the singular value decomposition of the received field at different frequencies over a vertical array of transducers. These methods require that the vertical array spans the full water column. This is obviously a strong limitation for the application of such methods in an operational context. To overcome these shortcomings, this study proposes to replace the spatial diversity constraint by a frequency diversity condition, and thus considers the case of a field emanating from an impulsive source. Indeed, because of the discrete nature of the wavenumber spectrum and due to their dispersive behavior, the modes are separated in the time-frequency domain. This phenomenon enables the design of a modal filtering scheme for signals received on a single receiver. In the case of a vertical receiver array, the modal contributions can be isolated for each receiver even when using a partial water column spanning array. This method thus eliminates the receiving constraints of classical methods of modal depth function estimation, although it imposes the use of an impulsive source. The developed algorithm is benchmarked on numerical simulations and validated on laboratory experimental data recorded in an ultrasonic waveguide. Practical applications are also discussed.

Current and Density Observations of Packets of Nonlinear Internal Waves on the Outer New Jersey Shelf

AbstractClosely spaced observations of nonlinear internal waves (NLIWs) were made on the outer continental shelf off New Jersey in June 2009. Nearly full water column measurements of current velocity were made with four acoustic Doppler current profilers (ADCPs) that were moored about 5 km apart on the bottom along a line approximately normal to the bathymetry between water depths of 67 and 92 m. Density profiles were obtained from two vertical strings of temperature and conductivity sensors that were deployed near each of the interior ADCP moorings. In addition, a towed ScanFish provided profiles and fixed-level records of temperature and salinity through several NLIW packets near the moorings. Several case studies were selected to describe the propagation of the NLIWs. One to three solitary waves of depression were observed in five selected packets. There were also occurrences of multiple-phase dispersive wave packets. The average propagation speed corrected for advection of the observed waves was 0.51 ...

Current and Density Observations of Packets of Nonlinear Internal Waves on the Outer New Jersey Shelf

Abstract Closely spaced observations of nonlinear internal waves (NLIWs) were made on the outer continental shelf off New Jersey in June 2009. Nearly full water column measurements of current velocity were made with four acoustic Doppler current profilers (ADCPs) that were moored about 5 km apart on the bottom along a line approximately normal to the bathymetry between water depths of 67 and 92 m. Density profiles were obtained from two vertical strings of temperature and conductivity sensors that were deployed near each of the interior ADCP moorings. In addition, a towed ScanFish provided profiles and fixed-level records of temperature and salinity through several NLIW packets near the moorings. Several case studies were selected to describe the propagation of the NLIWs. One to three solitary waves of depression were observed in five selected packets. There were also occurrences of multiple-phase dispersive wave packets. The average propagation speed corrected for advection of the observed waves was 0.51 ± 0.09 m s−1. The waves were directed primarily shoreward (~northwestward) along the mooring line with average wavelengths and periods of about 300 m and 10 min, respectively. Wave amplitudes and energies decreased with decreasing water depth. The observed wave parameters can be locally described by a two-layer Korteweg–de Vries (KdV) model, except for the decreasing amplitudes, which may be due to shear-induced dissipation and/or bottom drag. The various complementary observations utilized in this study present a unique description of NLIWs.

Experimental determination of the three-dimensional primary field of a seismic airgun array.

The Littoral Acoustic Demonstration Center (LADC) recorded acoustic and related data on three moored arrays and ship-deployed hydrophones, which together spanned the full water column to measure the 3-D acoustic field of a seismic airgun array, in September 2007. The seismic source vessel shot specified lines to give detailed solid angle and range information concerning the field of the primary arrival. The data were collected in the western Gulf of Mexico between the East Break and Alamos Canyon regions. Peak pressures, sound exposure levels, total shot energy spectra, one-third octave band analyzes, and source directivity studies are used to characterize the field. Three-dimensional maps of these quantities are generated to show dependence on emission and azimuthal angles and range. Longer lines were shot for propagation analyzes. Source signatures from NUCLEUS and GUNDALF are used in propagation modeling. [Research supported by the Joint Industry Programme through the International Association of Oil and Gas Producers.]

Interannual to decadal variability of outflow from the Labrador Sea

A decade of weak convection in the Labrador Sea associated with decreasing water mass transformation, in combination with advective and eddy fluxes into the convection area, caused significant warming of the deep waters in both the central Labrador Sea and boundary current system along the Labrador shelf break. The connection to the export of Deep Water was studied based on moored current meter stations between 1997 and 2009 at the exit of the Labrador Sea, near the shelf break at 53°N. More than 100 year‐long current meter records spanning the full water column have been analyzed with respect to high frequency variability, decaying from the surface to the bottom layer, and for the annual mean flow, showing intra‐ to interannual variability but no detectable decadal trend in the strength of the deep and near‐bottom flow out of the Labrador Sea.

Tidal and atmospheric influences on near‐surface turbulence in an estuary

Estuarine near‐surface turbulence is important for transport, mixing, and air‐water exchanges of many important constituents but has rarely been studied in detail. Here, we analyze a unique set of estuarine observations of in situ atmospheric and full water column measurements, estimated air‐sea exchanges, and acoustic measurements of several terms in the turbulent kinetic energy (TKE) budget. Observations from a 5.1 m deep site in the Hudson River estuary include dissipation at 50 cm depth (ɛ50), as well as profiles of TKE, shear production of TKE (P), and net turbulent vertical TKE transport (TD). Regressions suggest that the principal controlling factor for ɛ50 was wind (through the surface shear velocity, U*) and that the surface heat flux and tidal currents played a secondary role. For ebb spring tides, the TKE budget at 50 cm depth was closed within noise levels. Ebbs had high ɛ50 due to local shear production, which nearly balanced ɛ50. Floods had TD approaching P in the upper water column but generally weak near‐surface shear and turbulence. Examining buoyancy fluxes that impact near‐surface stratification and can indirectly control turbulence, solar heat input and tidal straining caused similar buoyancy fluxes on a sunny, calm weather day, promoting ebb tide restratification. Wind‐driven mixing was found to dominate during a fall season storm event, and strong overnight heat loss after the storm helped delay restratification afterward. These results demonstrate the utility of combining detailed air‐sea interaction and physical oceanographic measurements in future estuary studies.

Variability of the semidiurnal internal tides observed on the Timor Shelf

Semidiurnal internal tides drive strong flow variations on the shelf between Timor Island and Australia. The variability of the internal tide over seasonal to interannual time scales is studied using temperature and velocity measurements from 11 moorings, a first attempt for an observational regional synthesis of internal tide behavior. The energy density and flux carried by the semidiurnal tides were calculated from the full water column profile estimated by the Gauss‐Markov fit of the vertical modes incorporating the effects of sloping bottom topography. The estimated energy flux can be interpreted either as generated at and propagating from nearby forcing sites where strong barotropic tides flow over steep topography or as trapped by and circulating around an island topography. The relevant horizontal length scale of the topography is several tens of kilometers, which is the wavelength of the vertical first‐mode internal tides. At some of the moorings with sufficiently long records, seasonal variations of energy density and flux are found with maxima in late summer to winter, which is the same phase as the seasonal variation of the temperature stratification between 100 and 200 m depths. The relationship between the temperature stratification and the internal tide energy density is not linear. For the energy density, the amplitude of the seasonal variation as well as the standard deviation correlates well with the mean value with a correlation coefficient &gt;0.8 and a p value &lt;0.03, but the correlation between the energy density and area‐integrated barotropic tidal forcing is not significant.

Full-depth profiles of prokaryotes, heterotrophic nanoflagellates, and ciliates along a transect from the equatorial to the subarctic central Pacific Ocean

Studies in epipelagic waters report higher heterotrophic microbial biomass in the productive high latitudes than in the oligotrophic low latitudes; however, biogeographical data are scarce in the deep ocean. To examine the hypothesis that the observed latitudinal differences in heterotrophic microbial biomass in the epipelagic zone also occur at depth, abundance and biomass of heterotrophic prokaryotes, nanoflagellates (HNF), and ciliates were determined at depths of 5-5000 m in the central Pacific between August and September of 2005. Heterotrophic microbial biomass increased from the tropical to the subarctic region over the full water column, with latitudinal differences in prokaryotic biomass increasing from 2.3-fold in the epipelagic zone to 4.4-fold in the bathypelagic zone. However, the latitudinal difference in HNF and ciliate biomass decreased with depth. In the mesopelagic zone, the vertical attenuation rate of prokaryotic abundance, which was calculated as the linear regression slope of log-log plot of abundance versus depth, ranged from -0.55 to -1.26 and was more pronounced (steeper slope) in the lower latitudes. In contrast, the vertical attenuation rate of HNF in the mesopelagic zone (-1.06 to -1.27) did not differ with latitude. In the subarctic, the attenuation rate of HNF was 1.7 times steeper than for prokaryotes. These results suggest the accumulation of prokaryotes in the deep subarctic Pacific, possibly due to low grazing pressure. Although the vertical attenuation rate of ciliates was steepest in the bathypelagic zone, HNF abundance did not further decrease at depths below 1000 m, except for at 2000 m where HNF was lowest across the study area. Ciliate abundance ranged 0.3-0.8 cells l -1 at 4000 m, and were below the detection limit ( < 0.1 cells l -1 ) at 5000 m. To our knowledge, this study presents the first data for ciliates below 2000 m.

The three‐dimensional acoustic field of primary arrivals from a seismic airgun array.

The Source Characterization Study 2007 (SCS07) measured the three‐dimensional (3‐D) acoustic field of a seismic airgun array. The Littoral Acoustic Demonstration Center (LADC) performed the experiment, collecting acoustic and related data on three moored hydrophone arrays and one ship‐deployed hydrophone array which together spanned the full water column. Sensitive and desensitized phones were deployed at each position to extend the dynamic range. An ultra short baseline localization system was deployed with the EARS moorings to provide array shape. With postanalysis this results in time‐dependent positions for each of the acoustic sensors. Every channel is calibrated. A seismic source vessel shot a series of lines designed to give detailed angle and range information concerning the field of the primary arrival. Peak pressures, sound exposure levels, total shot energy spectra, and one‐third octave band analyses give important insights into details of the acoustic field. Images of these quantities are generated to show dependence on emission and azimuthal angles and range. 3‐D visualizations and two‐dimensional cuts through the data are shown. [Research supported by the Joint Industry Programme through the International Association of Oil and Gas Producers.]

Sound transmission through a fluctuating ocean (1979) and ocean acoustic tomography (1995): An intertwined history.

Early attempts at understanding scintillations in sound transmission through the ocean assumed that ocean fine structure is homogeneous and isotropic. It is, in fact, dominated by internal waves, which are neither homogeneous nor isotropic. Sound Transmission Through a Fluctuating Ocean (Flatté et al., Cambridge University Press, Cambridge, England 1979) combined recently developed internal-wave models with path-integral methods to predict the fluctuations of resolved acoustic multipaths. At that time, Ocean Acoustic Tomography (Munk et al., Cambridge University Press, Cambridge, England 1995), which uses ray travel times to determine large-scale ocean structure, had just been proposed. The wideband acoustic sources and receivers required to make travel-time measurements for tomography also provided the technology needed to quantify acoustic fluctuations. Conversely, internal-wave-induced scattering limits travel-time measurement precision. Tomography experiments at 25-km range provided some of the earliest tests of path-integral predictions. Subsequent tomography experiments provided data at ever-increasing ranges and decreasing frequencies. The Acoustic Thermometry of Ocean Climate project made measurements at megameter ranges with 1400-m vertical line array receivers. The next step is a modular distributed VLA capable of spanning the full water column that is under development to enable separation of acoustic modes using spatial filtering and to fully characterize deep-water acoustic time fronts.