A Global Analysis of Sverdrup Balance Using Absolute Geostrophic Velocities from Argo [Journal of Physical Oceanography]
(Journal of Physical Oceanography Via Acquire Media NewsEdge) ABSTRACT
Using observations from the Argo array of profiling floats, the large-scale circulation of the upper 2000 decibars (db) of the global ocean is computed for the period from December 2004 to November 2010. The geostrophic velocity relative to a reference level of 900 db is estimated from temperature and salinity profiles, and the absolute geostrophic velocity at the reference level is estimated from the trajectory data provided by the floats. Combining the two gives the absolute geostrophic velocity on 29 pressure surfaces spanning the upper 2000 db of the global ocean. These velocities, together with satellite observations of wind stress, are then used to evaluate Sverdrup balance, the simple canonical theory relating meridional geostrophic transport to wind forcing. Observed transports agree well with predictions based on the wind field over large areas, primarily in the tropics and subtropics. Elsewhere, especially at higher latitudes and in boundary regions, Sverdrup balance does not accurately describe meridional geostrophic transports, possibly due to the increased importance of the barotropic flow, nonlinear dynamics, and topographic influence. Thus, while it provides an effective framework for understanding the zero-order wind-driven circulation in much of the global ocean, Sverdrup balance should not be regarded as axiomatic.
(ProQuest: ... denotes formulae omitted.)
The general circulation of the ocean plays a central role in the global transport of heat, freshwater, carbon, oxygen, nutrients, and other constituents. The ocean's large-scale circulation strongly influences the distribu- tion of these important quantities and thus shapes both the global climate and patterns of marine biological production. A fundamental zero-order theory of the large- scale circulation, Sverdrup balance (Sverdrup 1947) di- rectly relates meridional geostrophic transport to wind stress. Despite its key role in our understanding of ocean dynamics, the extent to which this tenet quantitatively describes the observed global ocean circulation has re- ceived surprisingly little investigation.
One of the most likely causes for the lack of analyses of Sverdrup balance is the dearth of subsurface velocity observations. Subsurface measurements have tradi- tionally been limited in both space and time because of the resource-intensive nature of shipboard observations and the sheer magnitude of the world's oceans. In his seminal series, Reid (1989, 1994, 1997, 2003) compiled over six decades of hydrographic and tracer data to es- timate the mean geostrophic circulation of the Atlantic, Pacific, and Indian Oceans. His analysis relied on somewhat subjective choices for constraints on the total transport and reference velocities and also assumed that both the flow and mixing take place along isopycnal surfaces. Several other studies (Macdonald and Wunsch 1996; Macdonald 1998; Ganachaud and Wunsch 2000; Ganachaud 2003) have used hydrographic sections from the World Ocean Circulation Experiment (WOCE) in conjunction with inverse box models to produce quan- titative estimates of global ocean circulation. Such esti- mates, however, are sensitive to particular assumptions concerning the accuracy of the conservation equations used to constrain the flow. Although these analyses represent considerable progress in our ability to deter- mine subsurface geostrophic flow, the spatial and tem- poral coverage of the data used in these estimates remains relatively narrow.
In contrast, the Argo array of more than 3000 auton- omous profiling floats (Roemmich et al. 2004) offers an abundance of hydrographic and velocity observations with unprecedented resolution in both space and time. These instruments provide velocity estimates and mea- sure temperature and salinity in the upper 2000 decibars (db) of the ocean, relaying the data to shore-based data centers in real time. Argo floats have been deployed throughout the global ocean with a nominal horizontal spacing of 300 km. Each float collects hydrographic profiles at regular intervals, typically 10 days. In this study, absolute geostrophic velocities are determined by combining these hydrographic profiles with the direct estimates of velocity provided by the float trajectories. The resulting quantitative estimate of the large-scale circulation has significant advantages, including vastly improved spatial and temporal data coverage and no ad hoc assumption regarding the reference velocity.
The goals of this work are to 1) compute the absolute geostrophic velocity in the upper 2000 db of the global ocean during the period from 1 December 2004 to 30 November 2010 using observations from the Argo array of profiling floats and 2) use the resulting velocity esti- mate, together with satellite observations of wind stress, to examine the validity of Sverdrup balance on a global point-by-point basis. Given that this highly simplified theory of the wind-driven circulation is not expected to correspond exactly to the observed ocean everywhere, determining if and when it quantitatively predicts geo- strophic circulation is valuable. Investigating where Sverdrup balance succeeds and where it fails, and un- derstanding why, will also help further our knowledge of the dynamics that dominate different ocean regimes. Section 2 gives a brief derivation of Sverdrup balance and reviews relevant results from previous work. The data and methods used to compute the gridded velocities and evaluate Sverdrup balance are discussed in sections 3 and 4, respectively. Section 5 presents the absolute geostrophic velocity fields followed by an analysis of the validity of Sverdrup balance and the implications of these findings. Last, section 6 reviews the conclusions of this study.
Since Sverdrup (1947) first outlined the elegant for- mulation of the wind-driven circulation now known as Sverdrup balance (hereafter SB), it has become a ca- nonical part of circulation theory. The derivation of SB is based on the linear vorticity equation connecting meridional geostrophic velocity yg to the vertical gra- dient of vertical velocity w:
Following standard notation, f is the Coriolis parameter and b is its meridional gradient. The validity of (1) requires steady flow as well as small Rossby and Ekman numbers. Integrating vertically from a depth z 5 h to the base of the Ekman layer (denoted as z 5 0) gives
where Vg is the meridional component of the geo- strophic transport. Basic Ekman theory provides that the vertical velocity at the base of the Ekman layer is related to the wind stress vector t 5 (tx, ty) through
where r is the density and k is the unit vector in the vertical direction. Next we make the critical assumption that w(z 5 h) 5 0. Given this, (2) and (3) can be com- bined to produce
Expanding the right-hand side of (4) gives a useful partition in terms of the meridional component of Ek- man transport VE and what is commonly thought of as the Sverdrup transport VSv:
This statement of SB separates Vg from the wind- derived transport VSv 2 VE, while also showing that the total meridional transport V 5 Vg 1 VE equals VSv. The form of SB provided by (5) is directly tested in this study, using observations provided by the Argo array.
Relatively few studies have attempted to assess the validity of SB using observational data, and those that have are mostly limited to fairly narrow regions. Leetmaa et al. (1977) first explored a zonally integrated form of SB in the North Atlantic. Using a hydrographic transect across 258N and assuming a level of no mo- tion, they found that the zonally integrated southward transport was consistent with SB. Wunsch and Roemmich (1985), however, identified several issues with this study that called this conclusion into question. These included the extrapolation of a level of no motion through the Gulf Stream recirculation, contrary to the evidence there; inconsistency with the meridional heat transport needed to produce a substantial quantity of deep water at high latitudes; and a contradiction be- tween the sign of the bottom w obtained by integrating (1) downward from a surface Ekman pumping velocity velocity and the bottom w expected based on continuity and the slope of the seafloor. At the same latitude, Schmitz et al. (1992) subsequently concluded that the simple balance between southward interior geostrophic transport and northward Gulf Stream transport assumed in the Leetmaa study was not supported by water mass characteristics.
In the North Pacific, Qiu and Joyce (1992) determined that interannual variations in Sverdrup transport were strongly correlated with the transport of the Kuroshio and the North Equatorial Current. Hautala et al. (1994) found good agreement between the zonally integrated meridional geostrophic transport inferred from a hy- drographic section across 248N and that computed from wind data according to SB. Using altimetric data, Vivier et al. (1999) examined a time-dependent topographic SB, in which the integral in (2) is computed from the seafloor at z 5 H to the sea surface and w(z 5 H) is given by conti- nuity [w(z 5 H) 5 u ^ $H]. They found that this version of SB explained variations in sea surface height in much of the North Pacific and the subtropical South Pacific.
More recently, Wunsch (2011) carried out a global quantitative analysis of SB using the ocean state esti- mate produced by the Estimating the Circulation & Climate of the Ocean-Global Ocean Data Assimilation Experiment (ECCO-GODAE) project from a least squares fit of general circulation models to global-scale datasets available after 1992. On a point-by-point basis, the form of SB given by (2) was tested by integrating the meridional velocity, averaged over 16 years, from the depth where jwj , 1028 ms21 to the surface. This transport was then compared to f/b[w(z 52117 m)], where w(z 52117 m) was assumed to be the Ekman pumping velocity. Using that formulation, SB was found to be reasonably accurate over much of the interior subtropical and tropical oceans. In this study, we provide another global assessment of the validity of SB, based solely on data obtained as part of the ongoing Argo program together with satellite wind stress observations.
a. Oceanographic data
The oceanographic data in this study were collected by the Argo array of autonomous profiling floats and consist of upper ocean profiles of temperature T and salinity S as a function of pressure p, along with float trajectories determined by satellite positioning. The Argo array is composed of several different types of quasi-Lagrangian profiling floats. While the individual float types vary in terms of their specific components, all the floats used here were equipped with a buoyancy engine, a satellite antenna for positioning and transmission, and an instrument package containing a conductivity- temperature-depth (CTD) sensor. To ensure high qual- ity, nearly all of the data included in our study (99.9% of the total profiles) were collected by floats with SeaBird Electronics CTDs, with the remaining small quantity of data from a selected subset of floats with Falmouth Scientific, Inc. (FSI) CTD sensors.
Each float repeatedly executed a mission cycle that consisted of four phases: 1) descent to the predefined ''parking pressure'' ppark; 2) subsurface drift at ppark for a predetermined length of time; 3) ascent to the sur- face with concurrent measurements of T, S, and p; and 4) surface drift with transmission of the hydrographic data and position fixes to a satellite. The data trans- mitted from the floats were compiled at two onshore Global Data Assembly Centers (GDACs) and were made publicly available usually within 24 h. More details regarding the technical aspects and operation of the floats as well as the processing of the data are given by Roemmich et al. (2004) and the Argo Data Management site (http://www.argodatamgt.org/Documentation).
All data used in this study were obtained from the Coriolis GDAC at the French Research Institute for Exploitation of the Sea (IFREMER), France (www. coriolis.eu.org). Only data that had undergone thorough quality control testing and were designated as delayed mode were included. In addition, data from floats with known problems (i.e., those on the ''greylist'' compiled by the GDACs) were discarded. Data from marginal seas, such as the Gulf of Mexico and the Indonesian seas, were also excluded. The time period for this analysis spanned from 1 December 2004 to 30 November 2010. Prior to this time period, global Argo data coverage was sparse, and complete global delayed-mode data were not available for later years at the time of preparation.
A total of 485 704 T and S profiles from 5027 Argo floats were used in this study. The profiles had variable vertical resolution, with observations typically at 5- or 10-db increments near the surface and gradually in- creasing to 100-db increments below 1500 db. Only profiles extending to 870 db or deeper were included in this dataset. The data are for the most part evenly dis- tributed in space (Fig. 1a), although the edges of the basins are not as well covered as the interior. The tem- poral distribution of the data is skewed toward the later part of the study period (Fig. 1b) because the Argo array did not reach full coverage until early 2007.
At each profile location, T, S, and p were used to compute 1) a profile of dynamic height D as a function of p; 2) profiles of potential temperature u, S, and p as a function of potential density su; and 3) maximum p, average u, and average S of the mixed layer. To de- termine D, specific volume anomalies d(S, T, p) 5 1/r(S, T, p) 2 1/r(35, 0, p) were calculated and then integrated with respect to pressure according to
The reference pressure p0 was set to 900 db to allow easy combination with the velocity data derived from the float trajectories. Dynamic height was computed at 29 regular pressure levels spanning the upper 2000 db. Extrapola- tion of the D profile beyond the pressure range of the original T and S profiles was limited to less than 28 db at the bottom of the profile and less than 4 db at the top.
The integral in (6) was calculated by interpolating the d profile with a piecewise polynomial composed of cubic Hermite basis functions (Fritsch and Carlson 1980) and then analytically integrating the interpolant. This form of interpolating polynomial was chosen because it efficiently produces a smooth function that is monotonic on each interval and is therefore shape preserving with no superfluous ''bumps'' or ''wiggles.'' The resulting interpolated profiles were visually appealing and had minimal introduced density inversions. To further ensure the quality of the profile data, we manually inspected all interpolated profiles. Profiles that contained gaps of more than 400 db, density inversions $0.02 kg m23 db21,or unrealistic interpolations (as determined by comparison to nearby profiles) were either corrected by removing a portion of the data or discarded entirely if correction was not possible.
The deepest p of the mixed layer along with the av- erage mixed layer u and S were computed for each profile using the density algorithm of Holte and Talley (2009). Profiles for which the algorithm failed to pro- duce a mixed layer p, such as those that did not include any p , 20 db, were excluded, resulting in a dataset of 476 319 mixed layer data points. For each profile, su was computed, and then u, S, and p were interpolated onto regular su levels, again using a piecewise cubic Hermite polynomial without extrapolation. To span the full range of the Argo profiles used here, the su levels ranged from 22.4 to 27.95 kg m23 and were linearly spaced with decreasing intervals from 0.25 to 0.01 kg m23.
Trajectory data for each float were used to estimate the drift velocity at the parking pressure u(ppark) for each cycle. Only floats with a nominal ppark $ 800 db were included in this study. Drift velocities were com- puted using the time the float began its descent tdes and the time it reached the surface at the end of the sub- sequent ascent tasc, together with the corresponding diving location xdes and surfacing location xasc:
The t and x of the velocity estimate were taken as the average of tasc and tdes and the average of xasc and xdes, respectively. This analysis only used subsurface drifts for which 5 # (tasc 2 tdes) # 25 days. Shorter cycles do not adequately measure geostrophic velocity due to a stronger influence from high-frequency fluctuations, while longer cycles result in excessive smoothing in space and time.
For each cycle along the float trajectory, xdes and xasc were determined with one of two methods depending on the type of satellite positioning system used by the float, the Argos system (96.6% of the data used here) or global positioning system (GPS; 3.4%). For Argos floats, which remain at the surface for 6-12 h while obtaining multiple position fixes of varying accuracy, the procedure fol- lowed that detailed by Park et al. (2005). Both xasc and xdes are generally unknown because a variable amount of time elapses between tasc and the first position fix and between the last position fix and tdes. The float experi- ences much stronger currents at the surface than at ppark, and thus the accuracy of u(ppark) is greatly improved if xasc and xdes are estimated taking into account the sur- face drift and the elapsed time (Park et al. 2005). To accomplish this, tdes and tasc are needed. These times were either obtained directly from the GDAC when available or estimated from the float metadata and/or surface fixes using the protocols given in Park et al. (2005). When it was not possible to calculate tasc and tdes from the available data for a particular cycle (37.6% of the data), they were approximated by computing the average Dt between tasc and the first surface fix and the average Dt between tdes and the last surface fix for all other cycles of the float and then adding these to the first and last fixes, respectively. In the small number of cases where we were unable to estimate tasc and/or tdes by any of these means (3.8%), the times of the first and last surfaces fixes were used instead.
Next xasc and xdes were estimated via a least squares fit of the fixes from each surface drift to a theoretical tra- jectory x(t), derived from a linear velocity UL and an inertial oscillation of magnitude UI and phase f:
where t0 is the time of the first fix. Each surface fix was weighted by 1/^2, where ^ is the standard deviation of the fix (150, 350, and 1000 m for Argos fixes of level 3, 2, and 1, respectively). The fitted trajectory was then extrapo- lated to tdes and tasc to give xdes and xasc. Park et al. (2005) showed that this method produces velocities that are more accurate and less biased than simple linear ex- trapolation of the surface fixes. The total uncertainty of each velocity estimate computed in this manner is on the order of 0.2cms21 (Park et al. 2005).
Floats using GPS positioning obtain a single surface fix with a much higher accuracy (10 m) and remain on the surface for a significantly shorter period of time (typically ,30 min). In these cases, xasc and xdes for each surface drift were both taken as the sole position fix, and tdes and tasc were set to the corresponding time.
The mean of the p measurements recorded along the drift track gave the measured ppark for each cycle. If no in situ p values were reported, ppark was assumed to be the preprogrammed nominal ppark. Only cycles for which 780 # ppark # 2020 db were included in this da- taset. Data were manually inspected if the difference between the nominal ppark and the measured ppark ex- ceeded 100 db. When it appeared that the float had shoaled, the pressure sensor had failed, or the float had come into contact with the seafloor (as determined by comparing the float trajectory with available bathy- metric maps), the affected data were discarded. Any cycles flagged by the GDACs for likely grounding of the float were also removed from the dataset.
After estimating ppark and u(ppark), any outlying ve- locity estimates more than 10 standard deviations from the mean over the whole dataset were excluded from further analysis. To combine the velocity data with the dynamic height data, the u(ppark) estimates were ad- justed to the common reference pressure p0 (900 db). The geostrophic shear between ppark and p0 was calcu- lated by interpolating the mapped D fields (described below in section 4) to the location of u(ppark). The u(p0) was then computed by adding this shear to u(ppark). Data for which the magnitude of the velocity adjustment was greater than 20 cm s21 were eliminated from the dataset. The discarded data were all located within 58 of the equator, where the assumption of geostrophy be- comes unreliable as the Rossby number grows as 1/f.
The resulting velocity dataset, consisting of 442 211 estimates of u(p0) from 4799 floats, provides a repre- sentative sample of the geostrophic velocity at p0. The distribution in space (not shown) and time (Fig. 1b) is comparable to that of the profile data, although there are slightly fewer estimates due to the differences in the quality and requirements of each dataset.
b. Wind stress data
Wind stress data needed to compute the Sverdrup and Ekman transports in (5) were derived from observations from the National Aeronautics and Space Adminis- tration (NASA) SeaWinds scatterometer onboard the Quick Scatterometer (QuikSCAT). These data were obtained from the Centre de Recherche et d'Exploita- tion Satellitaire (CERSAT) at IFREMER in Plouzan^e (France). A gridded product was used that gave global monthly 10-m wind fields at 0.5830.58 resolution, spanning from August 1999 to October 2009. Wind stress magnitude and direction as well as wind stress curl were provided.
a. Geostrophic velocity
Using the Argo dynamic height and velocity data, monthly absolute geostrophic velocities in the Atlantic, Pacific, Indian, and Southern Oceans were computed on 29 pressure surfaces spanning the upper 2000 db for the period from December 2004 to November 2010. The procedure consisted of three basic steps. First, the D data were gridded on each p surface to produce 18318 maps of both D and velocity relative to p0 urel, which is derived from D according to
Second, the u(p0) estimates were mapped onto the same spatial and temporal grid. A geostrophic streamfunction c0 was computed simultaneously, where c0 is defined here as
Third, the simple addition of the two fields gave monthly gridded absolute geostrophic velocities at each p level, as shown by summing (7) and (8). Sim- ilarly, the streamfunction for the absolute geostrophic flow c(p) was found by combining c0 and D(p) at each level.
All mapping was carried out using a variation of the method put forth by LeTraon (1990). Essentially this technique maps a large-scale signal using the objective function fitting procedure developed by Davis (1985) and a small-scale signal using the conventional objective analysis described by Bretherton et al. (1976). The two signals are characterized by different scales in both space and time. This procedure has the advantage of requiring no a priori knowledge of the statistics of the large-scale field; however, the statistics of the small- scale signal must be specified. Rather than simply as- suming a decorrelation length scale, as is typically done, the present study uses a novel iterative approach to de- termine the statistics of the small-scale signal directly from the data. Full details of the mapping procedure including the iterative scheme are given in the appendix. For all p levels, mapping was restricted to areas where the depth of the seafloor H $ 1000 m, as determined from a box average of 1 arc-minute global relief model of Earth's surface (ETOPO1) bathymetric data (Amante and Eakins 2009). Additional grid points were excluded for p levels .1000db if p$ H.
Uncertainties in the gridded velocities arise from a number of sources including array bias, measurement errors, and sampling errors (Katsumata and Yoshinari 2010). The mapping procedure used here has the distinct advantage of producing at each grid point a quantitative estimate of the mean squared error e2. This error esti- mate combines the usual errors in the large- and small- scale signals with the error arising in the small-scale field due to uncertainty in the large-scale signal (LeTraon 1990). The total error in the absolute velocity was given by the sum of the e2 associated with the relative and reference velocities.
The gridded absolute velocities described here [Ab- solute Geostrophic Velocities from Argo (AGVA)], together with accompanying error estimates, are avail- able for use by the oceanographic community and can be obtained online (http://flux.ocean.washington.edu/agva).
Assuming hydrostatic balance, volume transport can be written as the integral of velocity between two pres- sure surfaces p1 and p2 (with p1 , p2):
where g is the acceleration due to gravity. Because this expression requires estimates of in situ r in addition to u, monthly maps of u, S,andp on su surfaces were produced using the iterative two-step mapping technique pre- sented above. These quantities were mapped on su surfaces instead of p surfaces due to the significant spurious anomalies that are introduced by averaging u and S on p surfaces (Lozier et al. 1994). Mixed layer p, u, and S were also mapped with the same method. The gridded u, S, and p on su surfaces, together with the mixed layer maps, were used to compute monthly 18318 maps of r and su as a function of p. Mean squared errors were also calculated for each mapped quantity. The gridded u and S, vertically interpolated onto the same p levels as the u fields, are also available at the website listed above.
c. Sverdrup balance
To assess the extent to which SB accounts for global meridional geostrophic transport, the two sides of (5) were computed. There is some question as to what time period is appropriate for testing SB. In other words, does SB hold over synoptic time scales or does it only apply when the wind stress and geostrophic velocities are av- eraged over a longer time scale? Sverdrup (1947) and most oceanographic textbooks imply that SB should emerge after averaging over a long interval but do not specify a particular length. Instead of averaging, some authors (e.g., Niiler and Koblinksy 1985; Vivier et al. 1999) chose to discuss a time-dependent SB. Similar to Wunsch (2011), this study will examine the longest av- eraging periods permitted with the datasets used here, namely 6 years for the Argo data and 10 years for the QuikSCAT data, with an overlap between the two of almost 5 years.
The wind-derived transport in (5) was calculated using the zonal wind stress and wind stress curl derived from QuikSCAT together with the mean mixed layer r computed from the gridded mixed layer u and S. Monthly maps of wind-derived transport were produced for the entire QuikSCAT era (from August 1999 to October 2009) and then interpolated onto the same 183 18 horizontal grid as the mapped Argo data. These were then averaged over the whole time period and over the period contemporaneous with the Argo data (from December 2004 to October 2009).
To evaluate the geostrophic transport in (5), the mapped yg and r were used to compute Vg according to (9). In practice the integration was carried out by first finding the mean over the study period (from December 2004 to November 2010) of yg(p), u(su), S(su), p(su), and the mixed layer properties. The mean was also cal- culated for the period overlapping with the QuikSCAT data. Using these mean values, hr(p)i was constructed for each grid point (where the angle brackets denote the mean). Next hyg i and hri were interpolated to common p levels with piecewise cubic Hermite polynomials. The resulting interpolants were used to calculate hyg i/ghri, which was then integrated analytically. While in princi- pal hyg i/ghri is not exactly equal to hyg/gri, the differ- ence here is considered negligible as h1/ri ' 1/hri.
Of critical importance to any analysis of SB is the choice of h, the bottom boundary of the transport in- tegral in (2). One approach is to inquire as to the expected depth of the wind-driven circulation. The ventilated thermocline theory of Luyten et al. (1983) suggests that the wind-driven transport extends to the densest out- cropping isopycnal in a given basin. Rhines and Young (1982) identify the depth of potential vorticity homog- enization as the base of the wind-driven circulation. Motivated by these considerations and the observations of Talley (1988), Hautala et al. (1994) selected the depth of the 26.5 and 27.0 isopycnals in the eastern and central subtropical North Pacific, respectively. Other authors have integrated velocity from the surface to a fixed depth (e.g., Leetmaa et al. 1977). The most straightfor- ward idea, namely, that h should be the depth where w goes to zero, arises directly from the derivation of SB. Although historical observations of w are generally unsuitable for practical application of this concept, the ocean state estimate used by Wunsch (2011) did allow h to be defined in this way.
Here, we have evaluated SB using two different methods of determining h. The first choice for h is the depth of the maximum mean isopycnal outcropping in wintertime along the line of zero Sverdrup transport, which was prompted by the ventilated thermocline theory. In the second approach, instead of selecting an a priori value for h , at each grid point the depth and corresponding pressure were found that gave the best fit to SB. If the meridional geostrophic transport agrees with the transport predicted by SB, then this depth should also be the location in the water column where w is approximately zero. Although the existence of such a depth does not by itself guarantee that SB is satisfied, this value will be used as a starting point for further analysis. In both cases, the mean maps of u, S, p, and mixed layer properties were used to determine the de- sired limits of integration.
The uncertainty associated with the computed geo- strophic transport was derived from the error estimates for yg, u, S, and p using a Monte Carlo simulation with 500 iterations. The standard error around the mean was calculated for each quantity using the uncertainties on the monthly maps. In each iteration, hyg i, hui, hSi, hpi, and the mean mixed layer properties were individually perturbed by an amount randomly drawn from a uni- form distribution on the interval (2e, e), and the trans- port was calculated as described above. The error in hyg i was decomposed into the error associated with hy(p0)i and that associated with hyrel i, and these were varied separately. At each grid point, the errors in each vari- able were assumed to be correlated vertically but independent of the errors in the other quantities. The variance of the resulting series of 500 transport values provided an uncertainty on the computed geostrophic transport.
5. Results and discussion
a. Absolute geostrophic velocities
We first present the mean absolute velocity fields at four representative pressure levels: 5, 200, 1000, and 1500 db (Figs. 2-5). In each of these plots, c(p) is con- toured at the intervals specified in the captions. Colors indicate highs in c(p) in red and lows in blue; as (8) shows, these represent anticyclonic and cyclonic flows, respectively. All estimates at latitudes #58N/S are dis- carded because the geostrophic assumption becomes problematic in this region. In addition, results are only plotted where the uncertainty in the geostrophic velocity ejuj , sjuj/8, where sjuj is the standard deviation of juj.
At the shallow, near-surface level (Fig. 2), the domi- nant features are the anticyclonic subtropical gyres found in all basins except the northern Indian, the eastward- flowing Antarctic Circumpolar Current (ACC) south of approximately 408S, and the cyclonic subpolar gyres in the North Pacific and North Atlantic. Eastward flow is observed in the northern tropical Pacific and Atlantic, associated with the North Equatorial Countercurrent in each basin. While the long-term mean flow in the north- ern Indian Ocean is weak, an analysis of the seasonal cycle (not shown) reveals a strong signal likely associated with monsoonal forcing. At this p level, the geostrophic flow shown here makes up only a part of the total velocity field, with Ekman flow constituting a significant portion, especially at lower latitudes.
At 200 db (Fig. 3), the subtropical gyres have all shifted poleward compared to their locations at 5 db. In the South Pacific, the subtropical gyre now has two distinct lobes, one to the east of New Zealand and one at about 258S. Strong velocities are still found in the western boundary current regions and the ACC, while the circulation in the tropical regions is substantially weaker compared to the near-surface level.
In the mean flow at 1000 db (Fig. 4), the North Pacific subtropical gyre has become smaller and weaker relative to the upper levels. A generally southward flow is ob- served in the North Atlantic subtropics. The subpolar gyres, while somewhat reduced in strength, continue to occupy the entire North Atlantic and North Pacific north of approximately 458N. This observation agrees well with the fact that the flow at higher latitudes is known to be less baroclinic than the flow in the sub- tropics (Talley et al. 2011). In the Southern Hemisphere, the subtropical gyres are shifted even farther poleward, with a strong recirculation east of southern Africa as- sociated with the Agulhas retroflection. Anticyclonic flows in the western South Pacific, south of Australia, and in the southern Indian Ocean appear connected, indicative of a supergyre structure (Ridgway and Dunn 2007). The ACC continues to be quite vigorous at this pressure level, but elsewhere mean flow is relatively weak.
Last, the mean velocity field at 1500 db (Fig. 5) is similar in structure to the flow at 1000 db but exhibits more small-scale features, probably influenced by the smaller signal-to-noise ratio here. The increased south- ward flow along the western edge of the North Atlantic is suggestive of a deep western boundary current. A relatively strong southward flow along the east coast of South America is also present. In the Southern Hemi- sphere the ACC dominates the flow at this level, with significant anticyclonic flows persisting east of southern Africa and New Zealand and in the Tasman Sea.
The results presented here are broadly consistent with modern ideas of the large-scale circulation of the upper ocean. The most comprehensive observation-based de- scription of the subsurface circulation up to this point was given by Reid (1989, 1994, 1997, 2003), which relied on decades of hydrographic and tracer data. Contrasting our results with these syntheses shows a high degree of similarity at the upper pressure levels in many areas of the global ocean, especially the North Pacific and North Atlantic. In the Southern Hemisphere, where the his- torical dataset was sparser, more differences between our results and Reid's are apparent. For example, sev- eral small-scale undulations in the Agulhas retroflection and in the Pacific sector of the ACC (between 1408 and 1658E and also between 1608 and 1208W) stand out in our maps, where Reid's analyses show purely zonal flows. Deeper in the water column, more substantial disagreement is found. For instance, at 1000 and 1500 db the Southern Hemisphere subtropical gyres are oriented predominantly zonally in our results, whereas Reid gives a stronger meridional component at these pressures. In general, the results from Reid show a smoother and more coherent flow at the deeper levels than is seen in the maps computed from Argo observations.
One might expect there to be dissimilarities between Reid's results and the circulation estimates in this study due to the difference in the time periods over which the data were collected. While our estimates span from December 2004 to November 2010, the data compiled by Reid for his studies were collected from 1932 to 1991. Another major difference is the amount of data used in this study (over 480 000 Argo profiles) compared to Reid (5932 ship-based stations). In addition, our results em- ploy a direct estimate of the absolute velocities at the reference level, whereas Reid relied on a subjective ex- amination of tracer data to estimate absolute flow.
b. Sverdrup balance
The wind-derived transport [VSv 2 VE in (5)], aver- aged over the entire QuikSCAT period, exhibits the expected large-scale pattern of equatorward flow in the subtropics and tropics and poleward flow in the high latitudes (Fig. 6). The line of zero wind-derived trans- port is approximately aligned with the axes of the Kuroshio and Gulf Stream Extensions in the Northern Hemisphere and with the ACC in the Southern Hemi- sphere. The largest mean wind-derived transports are lo- cated in the higher latitudes, but there are local maxima in the eastern parts of most gyres. All of the following results will use wind transport averaged from August 1999 to October 2009 and Argo velocities averaged from December 2004 to November 2010. The longer averaging periods were adopted because they offered a slight re- duction in the uncertainties and using the shorter contem- poraneous periods produced only negligible differences.
The results of our first test of SB, using a priori esti- mates of the bottom boundary of the wind-driven cir- culation, are shown in Fig. 7. The maximum wintertime mixed layer su values along the line of zero Sverdrup transport, computed from the mapped mixed layer u and S, are 26.24, 27.24, and 27.25 kg m23 for the North Pa- cific, Southern Hemisphere and northern Indian, and North Atlantic, respectively. Because the zero Sverdrup transport line in the Southern Hemisphere was located in the circumpolar ACC, the same value was used for all basins there. The normalized difference between the two sides of (5) was determined as
with Vg computed from the mean Argo-derived velocity and density fields using the depth of the given isopycnals for h and VSv 2 VE computed from the QuikSCAT wind stress fields. The difference D defined in this way is zero where SB is exactly correct and is bounded by 61. Fig- ure 7 gives the minimum D after factoring in the un- certainty on Vg. Yellow indicates perfect agreement, while those areas where the designated isopycnal was not found in the mean are shown in dark gray. Overlaid on the normalized difference are contours of c at 5 db.
Using these three isopycnals, we infer that SB pro- vides a quantitative description of meridional transport for a large fraction of the global ocean, primarily in the subtropics and tropics. However, SB fails to accurately represent meridional transports in almost all of the ACC and subpolar North Pacific and North Atlantic for this choice of h. In the North Pacific subpolar gyre, the wind- derived transport predicts more poleward flow than is found in the geostrophic transport, which is perhaps not too surprising given that su 26.24 was chosen based on a ventilated thermocline model of the subtropical gyre. Using su 27.24, the northward geostrophic trans- port in the North Atlantic Current region is larger than the wind-derived transport, but north of this area the discrepancy reverses (at least in the part of the sub- polar gyre where this isopycnal is present). For the most part, we find more poleward geostrophic trans- port than wind-derived transport in the northern part of the ACC and less in the southern part. Furthermore, meridional transports in the western boundary areas and the boundary current extensions do not adhere to SB. The poleward wind-derived transport is generally much smaller than the observed geostrophic transport in these regions. Last, the eastern boundary regions in the Pacific and Atlantic also exhibit poor agreement with SB.
Examining SB in a different way, we find the mini- mum p that produces agreement between the two sides of (5), within one standard deviation of Vg (Fig. 8). Regions where no p was found that satisfied SB, within the depth range of these velocity fields (2000 db), are shown in dark gray. While the existence of a p that sat- isfies (5) does not prove per se that SB is valid, it does indicate areas that are suitable candidates for SB. Al- though the results are somewhat noisy, the same picture emerges as in the first test: potentially good agreement with SB in the subtropics and tropics and poorer agreement elsewhere. Where agreement is possible, the wind-driven circulation appears to be confined to shal- lower pressure levels in lower latitudes than in higher latitudes, although small-scale features are prevalent, as expected for an average over 6 years (Wunsch 2010).
Even if there exists a value for h that balances the two sides of (5), is that enough to claim that SB is accurate on anything more than a purely mathematical level? Sverdrup's original derivation assumes a flat bottom ocean, entirely wind-driven, above a quiescent abyss. While in reality the deep velocity certainly is not zero due to local thermohaline and topographic forcing, we expect that for SB to be valid, the wind-driven circula- tion must be sufficiently isolated from the deep cir- culation. To ensure this, we impose as an additional condition needed to satisfy SB that yg, within 6ey, is zero for some depth range below h. Because y g is related to ?w/?z through (1) and w(z 5 h) 5 0 when SB is satisfied, this criterion guarantees that w ' 0 in this depth range and thus that the wind-driven layer is well separated from the deeper circulation. The range chosen here corresponded to 200 db, which is arbitrary, although the results changed insignificantly for any choice from about 150 to 500 db. Stipulating this additional condition sub- stantially decreases the area where SB provides an accu- rate description of the meridional geostrophic transport (Fig. 9). For the most part, the regions where SB fails under this stricter definition are in higher latitudes and western boundary regions, as well as areas in the central North Atlantic and Indian subtropical gyres.
These findings align reasonably well with the results of Wunsch (2011), who used a 16-yr average of data- assimilated model output to test the form of SB given by (2), with wE given by w(z 52117 m). In that analysis, roughly 40% of the global ocean was found to agree with SB. The agreement was largely restricted to the interior subtropical gyres and the tropics, similar to what is presented here, although differences between those re- sults and ours are apparent in the small-scale features. Considering the different datasets and methods of test- ing SB, the overall agreement suggests that these results are robust.
Given the evaluation of SB offered by Figs. 7-9, the question that arises is why does SB succeed in some regions and not others? One obvious reason for an ap- parent mismatch between the observed geostrophic and wind-derived transports in this study is that w may reach zero at a depth deeper than 2000 db, the limit of the geostrophic velocity dataset used here. Similarly, if w never actually vanishes, (2) is typically integrated from z 5 0 to the seafloor so that a kinematic constraint on the flow at the bottom, w(z 5 H) 5 u ^ $H, can be used. This is likely a contributing factor to the lack of agreement in the high latitudes and the ACC, where the flow is more barotropic. At the points in these regions where SB is found to be valid, h is typically closer to 2000 db than is generally the case in the subtropics and tropics.
One potential reason that w might not become suffi- ciently small is that there are strong vertical motions associated with flow over topography. A simple test comparing the mean depth and bottom slope between the grid points where SB succeeds, as defined by the test shown in Fig. 9, and those where no agreement is found results in statistically significant differences (p value , 0.001) between the two groups. Places where SB is valid are on average deeper and flatter than those where SB is inaccurate. While this result does not demonstrate causality, it does validate our notion that SB is most accurate for a deep, flat bottom ocean, much as Sverdrup originally envisioned.
Another possible reason for disagreement in the high latitudes is that this study uses velocities averaged over 6 years, whereas the transit time for a first baroclinic Rossby wave is much longer in the high latitudes due to the reduction in the Rossby deformation radius (Chelton et al. 1998). While it is unknown whether a considerably longer time average would improve SB in the high latitudes, given the variability of the wind stress and other factors, the time needed for linear baroclinic adjustment there would certainly preclude the 6-yr time period used here. Furthermore, many high-latitude cur- rents, such as the North Atlantic Current, are forced by both wind and buoyancy gradients. In these cases, we would expect that a separate wind-driven layer cannot reasonably be isolated from other elements of the cir- culation. Last, in the western boundary currents and their extensions, nonlinear dynamics likely play an important role. Because SB is derived from the linear vorticity balance given in (1), which neglects such ef- fects, the lack of agreement in these regions is hardly surprising.
In this study, we have exploited the fact that the Argo array of profiling floats provides an exceptional tool for the estimation of large-scale absolute geostrophic ve- locities in the upper 2000 db of the ocean with un- paralleled spatial and temporal resolution. The maps of global mean geostrophic velocities presented here cor- respond well with what has been previously reported based on traditional hydrographic data in the upper depth range, while at deeper pressures a more variable and complicated flow is observed. Limitations of the gridded velocity fields stem from the sparse coverage of the Argo array in some regions, specifically the inability to resolve boundary currents. As a purely data-based analysis, these velocity fields make an important con- tribution to our knowledge of the detailed circulation of the global upper ocean and should be useful for initial- izing and validating models and state estimations. These fields can also provide the foundation for observational analyses of meridional heat and freshwater transport and help further our understanding of the role of ocean circulation in influencing these important climatic quantities. As long as the Argo array is maintained, these fields can be continually augmented with more data and expanded to cover a longer time period.
The absolute geostrophic velocities presented here have been used to assess the extent to which Sverdrup balance, a simple but ubiquitous theory of the relation- ship between wind forcing and ocean circulation, accu- rately predicts observed geostrophic transports on a point-by-point basis. We find that within the uncer- tainties there is good agreement over large regions of the world's oceans, namely, the subtropics and tropics away from the boundaries. Thus, for a substantial part of the global ocean, Sverdrup balance provides a reason- able quantitative picture of the observed circulation. However, the depth of the wind-driven circulation, de- fined as the location where the wind-derived transport matches the meridional geostrophic transport, has con- siderable geographic variability. While the simple model described by the Sverdrup balance is a useful starting point for a theoretical treatment of the large-scale cir- culation, the results presented here imply that it should not be unreservedly accepted as truth.
Acknowledgments. The data used in this work were collected by the International Argo Program and the many national programs that have been involved. We thank each of the many agencies and individuals that have contributed to Argo over the past decade. In the United States, this work was supported by the National Oceanographic and Atmospheric Administration (NOAA) through Grant NA17RJ1232 Task 2 to the Joint Insti- tute for the Study of the Ocean and Atmosphere at the University of Washington. We are especially indebted to Dr. Steve Piotrowicz of NOAA for his able and enthu- siastic management of the U.S. Argo program since its inception. In addition, we thank the two reviewers for their insightful and helpful comments.
Amante, C., and B. W. Eakins, 2009: ETOPO1 1 arc-minute global relief model: Procedures, data sources and analysis. NOAA Tech. Memo. NESDIS NGDC-24, 19 pp.
Batchelor, G. K., 1970: The Theory of Homogeneous Turbulence. Cambridge University Press, 197 pp.
Bretherton, F. P., R. E. Davis, and C. B. Fandry, 1976: A technique for objective analysis and design of oceanographic experi- ments applied to MODE-73. Deep-Sea Res., 23, 559-582, doi:10.1016/0011-7471(76)90001-2.
Chelton, D. B., R. A. DeSzoeke, M. G. Schlax, K. E. Naggar, and N. Siwertz, 1998: Geographical variability of the first baroclinic Rossby radius of deformation. J. Phys. Ocean- ogr., 28, 433-460, doi:10.1175/1520-0485(1998)028,0433: GVOTFB.2.0.CO;2.
Davis, R. E., 1985: Objective mapping by least-squares fit- ting. J. Geophys. Res., 90 (C3), 4773-4777, doi:10.1029/ JC090iC03p04773.
-, 1998: Preliminary results from directly measuring middepth circulation in the tropical and South Pacific. J. Geophys. Res., 103 (C11), 24619-24639, doi:10.1029/98JC01913.
-, 2005: Intermediate-depth circulation of the Indian and South Pacific Oceans measured by autonomous floats. J. Phys. Oceanogr., 35, 683-707, doi:10.1175/JPO2702.1.
Fritsch, F. N., and R. E. Carlson, 1980: Monotone piecewise cubic interpolation. SIAM J. Numer. Anal., 17, 238-246, doi:10.1137/ 0717021.
Ganachaud, A., 2003: Large-scale mass transports, water mass formation, and diffusivities estimated from World Ocean Circulation Experiment (WOCE) hydrographic data. J. Geo- phys. Res., 108, 3213, doi:10.1029/2002JC001565.
-, and C. Wunsch, 2000: Improved estimates of global ocean circulation, heat transport and mixing from hydrographic data. Nature, 408, 453-457, doi:10.1038/35044048.
Gille, S. T., 2003: Float observations of the Southern Ocean. Part I: Estimating mean fields, bottom velocities, and topographic steering. J. Phys. Oceanogr., 33, 1167-1181, doi:10.1175/ 1520-0485(2003)033,1167:FOOTSO.2.0.CO;2.
Hautala, S. L., D. H. Roemmich, and W. J. Schmitz, 1994: Is the North Pacific in Sverdrup balance along 248N? J. Geophys. Res., 99 (C8), 16041-16052, doi:10.1029/94JC01084.
Holte, J., and L. Talley, 2009: A new algorithm for finding mixed layer depths with applications to Argo data and Subantarctic Mode Water formation. J. Atmos. Oceanic Technol., 26, 1920- 1939, doi:10.1175/2009JTECHO543.1.
Katsumata, K., and H. Yoshinari, 2010: Uncertainties in global mapping of Argo drift data at the parking level. J. Oceanogr., 66, 553-569, doi:10.1007/s10872-010-0046-4.
Leetmaa, A., P. Niiler, and H. Stommel, 1977: Does the Sverdrup relation account for the mid-Atlantic circulation? J. Mar. Res., 35, 1-10.
LeTraon, P. Y., 1990: A method for optimal analysis of fields with spatially-variable mean. J. Geophys. Res., 95 (C8), 13 543- 13547, doi:10.1029/JC095iC08p13543.
Lozier, M. S., M. S. McCartney, and W. B. Owens, 1994: Anoma- lous anomalies in averaged hydrographic data. J. Phys. Oce- anogr., 24, 2624-2638, doi:10.1175/1520-0485(1994)024,2624: AAIAHD.2.0.CO;2.
Luyten, J. R., J. Pedlosky, and H. Stommel, 1983: The ventilated thermocline. J. Phys. Oceanogr., 13, 292-309, doi:10.1175/ 1520-0485(1983)013,0292:TVT.2.0.CO;2.
Macdonald, A. M., 1998: The global ocean circulation: A hydro- graphic estimate and regional analysis. Prog. Oceanogr., 41, 281-382, doi:10.1016/S0079-6611(98)00020-2.
-, and C. Wunsch, 1996: An estimate of global ocean circu- lation and heat fluxes. Nature, 382, 436-439, doi:10.1038/ 382436a0.
Niiler, P. P., and C. J. Koblinksy, 1985: A local time-dependent Sverdrup balance in the eastern North Pacific Ocean. Science, 229, 754-756, doi:10.1126/science.229.4715.754.
Park, J. J., K. Kim, B. A. King, and S. C. Riser, 2005: An advanced method to estimate deep currents from profiling floats. J. Atmos. Oceanic Technol., 22, 1294-1304, doi:10.1175/ JTECH1748.1.
Qiu, B., and T. M. Joyce, 1992: Interannual variability in the midlatitude and low-latitude western North Pacific. J. Phys. Oceanogr., 22, 1062-1079, doi:10.1175/1520-0485(1992)022,1062: IVITMA.2.0.CO;2.
Reid, J. L., 1989: On the total geostrophic circulation of the South Atlantic Ocean: Flow patterns, tracers, and transports. Prog. Oceanogr., 23, 149-244, doi:10.1016/ 0079-6611(89)90001-3.
-, 1994: On the total geostrophic circulation of the North Atlantic Ocean: Flow patterns, tracers, and transports. Prog. Oceanogr., 33, 1-92, doi:10.1016/ 0079-6611(94)90014-0.
-, 1997: On the total geostrophic circulation of the Pacific Ocean: Flow patterns, tracers, and transports. Prog. Ocean- ogr., 39, 263-352, doi:10.1016/S0079-6611(97)00012-8.
-, 2003: On the total geostrophic circulation of the Indian Ocean: Flow patterns, tracers, and transports. Prog. Oceanogr., 56, 137-186, doi:10.1016/S0079-6611(02)00141-6.
Rhines, P. B., and W. R. Young, 1982: A theory of the wind- driven circulation. I. Mid-ocean gyres. J. Mar. Res., 40, 559- 596.
Ridgway, K. R., and J. R. Dunn, 2007: Observational evi- dence for a Southern Hemisphere oceanic super- gyre. Geophys. Res. Lett., 34, L13612, doi:10.1029/ 2007GL030392.
Roemmich, D., S. Riser, R. Davis, and Y. Desaubies, 2004: Au- tonomous profiling floats: Workhorse for broad-scale ocean observations. Mar. Technol. Soc. J., 38 (2), 21-29, doi:10.4031/ 002533204787522802.
Schmitz, W. J., J. D. Thompson, and J. R. Luyten, 1992: The Sverdrup circulation for the Atlantic along 248N. J. Geophys. Res., 97 (C5), 7251-7256, doi:10.1029/ 92JC00417.
Sverdrup, H., 1947: Wind-driven currents in a baroclinic ocean; with application to the equatorial currents of the eastern Pa- cific. Proc. Natl. Acad. Sci. USA, 33 (11), 318-326, doi:10.1073/ pnas.33.11.318.
Talley, L. D., 1988: Potential vorticity distribution in the North Pacific. J. Phys. Oceanogr., 18, 89-106, doi:10.1175/ 1520-0485(1988)018,0089:PVDITN.2.0.CO;2.
-, G. L. Pickard, W. J. Emery, and J. H. Swift, 2011: Descriptive Physical Oceanography: An Introduction. 6th ed. Academic Press, 560 pp.
Vivier, F., K. A. Kelly, and L. Thompson, 1999: Contributions of wind forcing, waves, and surface heating to sea surface height observations in the Pacific Ocean. J. Geophys. Res., 104 (C9), 20767-20788, doi:10.1029/1999JC900096.
Wunsch, C., 2010: Toward a midlatitude ocean frequency- wavenumber spectral density and trend determina- tion. J. Phys. Oceanogr., 40, 2264-2281, doi:10.1175/ 2010JPO4376.1.
-, 2011: The decadal mean ocean circulation and Sverdrup balance. J. Mar. Res., 69, 417-434, doi:10.1357/ 002224011798765303.
-, and D. Roemmich, 1985: Is the North Atlantic in Sverdrup balance? J. Phys. Oceanogr., 15, 1876-1880, doi:10.1175/ 1520-0485(1985)015,1876:ITNAIS.2.0.CO;2.
ALISON R. GRAY AND STEPHEN C. RISER
School of Oceanography, University of Washington, Seattle, Washington
(Manuscript received 16 October 2012, in final form 7 January 2014)
Corresponding author address: Alison R. Gray, School of Oceanography, University of Washington, Box 355351, Seattle, WA 98195.
Objective Mapping of Argo Data
The technique used in this study to map the quantities of interest [D, u(p0), u, S, p, and mixed layer properties] closely follows that developed by LeTraon (1990). This method estimates both large- and small-scale signals from scattered data. The large-scale component, for which statistics are not known and perhaps not even clearly de- finable, is computed using the objective function fitting procedure of Davis (1985, 2005). A set of functions with spatial scales greater than a given limit, chosen from a complete set of functions, is fit to the data, subject to a zero-bias constraint. The large-scale signal is then removed from the data to produce small-scale anomalies that are objectively mapped following the conventional tech- niques of Bretherton et al. (1976).
In this study, the large- and small-scale signals were defined by different spatial and temporal scales. Spa- tially, the basis functions chosen to represent the large- scale field were spherical harmonics up to order 15. Thus, the large-scale component was the part of the signal with scales *248 latitude/longitude. This order was chosen for the spherical harmonic basis functions to accurately depict the large-scale signal and avoid over- lap with the mesoscale field. Temporally, the large-scale signal was computed as a mean monthly seasonal cycle over the 6-yr study period. All data from each month, regardless of year, were used to compute a large-scale field according to the protocol outlined by LeTraon (1990) for nonsimultaneous data, whereby the covari- ance between data from different years was set to zero. The small-scale signal, in contrast, was mapped sepa- rately for each individual month. It is important to note that separation into large- and small-scale signals is strictly statistical and thus does not necessarily corre- spond to distinct physical processes.
The statistics of the small-scale signal, that is, the autocovariance of the anomalies, must be specified a priori and are used in estimating both the small-scale field and the large-scale field (for which the least squares fit to the given basis functions takes into account the correlation of the anomalies). A Gaussian form for the spatial autocovariance was used here, although other forms were tried and did not substantially affect the results. The autocovariance was assumed to depend only on separation distance s according to
for variance C0 and decorrelation length scale L.
While typical analyses chose a value for L somewhat arbitrarily, we employ a novel method that determines C0 and L in an iterative fashion directly from the data for each p or su level. An initial guess for C0 and L was used to first estimate the large-scale signal via objective function fitting and remove that signal from the data. Here, we used C0 5 1andL 5 250 km, but testing showed that the final parameters were insensitive to the choice of the initial values. The resulting anomalies were then used to calculate the observed autocovariance C using all pairs of data collected in the same month and separated by spatial lag s:
where the angle brackets denote an ensemble mean taken over all valid pairs of data, and the prime denotes anomalies from the large-scale signal. The expression (A2) is written for D, but analogous equations were used for all scalar properties being mapped. All data pairs were grouped according to separation distance into 25- km bins and averaged to give a mean autocovariance curve. Next, a least squares fit of the observed mean autocovariance curve to (A1) gave both C0 and L. The data at the smallest separation distances (#25 km) were excluded from the fit due to probable space-time alias- ing. Using the fitted autocovariance parameters, the mapping of the large-scale signal and calculation of the autocovariance were then repeated. This procedure was iterated until the parameter values converged. This novel iterative approach has been tested and validated using both analytically constructed fields and output from a high-resolution model, the results of which can be found in A. R. Gray and S. C. Riser (2014, unpublished manuscript).
To maximize the number of data pairs included in the autocovariance estimation, and hence the degrees of freedom, the autocovariance is assumed to be isotropic. While this is not strictly true, Davis (1998) found that in the South Pacific the variability is generally isotropic except within about 58 of the equator. Thus, the assump- tion of isotropy most likely introduces significant errors only in the equatorial region, where the validity of the geostrophic relationship is already questionable.
Once the covariance parameters were determined for each p or su level, the large- and small-scale signals were both mapped onto a 18318 horizontal grid at the above- mentioned temporal resolution. Observational error and subgrid-scale noise were specified to account for an additional 50% of the measured variance in the data; this value was chosen to be reasonably conservative, following Gille (2003).
To map urel and u(p0), slight adjustments were made to the basic procedure. LeTraon (1990) and Bretherton et al. (1976) outline methods for mapping multivariate data and for determining fields related to the data through linear operators. As (7) and (8) show, the rel- ative and reference components of velocity are related linearly to D and c0, respectively. Using partial de- rivatives of both the autocovariance function (A1) and the spherical harmonic basis functions allowed us to map urel from the D data and c0 from the u(p0) data. In ad- dition, because the u(p0) data were vectors, a simple application of (A2) was not possible. These data were first multiplied by f to account for the latitudinal de- pendence of (8). Then for each pair of velocities the vec- tors were decomposed into longitudinal and transverse components, which are aligned parallel and perpendicular to the axis between the two data points, respectively. Longitudinal and transverse autocovariance curves were computed following Batchelor (1970), who showed that these can be derived from the covariance function describing the velocity streamfunction given that the statistics of the field are isotropic, homogeneous, and stationary. Last, a least squares method was used to simultaneously fit the mean observed longitudinal and transverse curves to the analytical forms of those functions.
(c) 2014 American Meteorological Society
[ Back To Technology News's Homepage ]