On the Transfer Function between Surface Fields and the Geostrophic Stream Function in the Mediterranean Sea
|
|
- Abigayle Blair
- 6 years ago
- Views:
Transcription
1 1406 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 On the Transfer Function between Surface Fields and the Geostrophic Stream Function in the Mediterranean Sea JORDI ISERN-FONTANET, MAHESH SHINDE, AND CRISTINA GONZÁLEZ-HARO Institut Catala de Ciencies del Clima, Barcelona, Spain (Manuscript received 28 August 2013, in final form 19 February 2014) ABSTRACT The real transfer function and the phase shift between sea surface height (SSH) and sea surface buoyancy (SSB) were estimated from the output of a realistic eddy-resolving model of the Mediterranean Sea circulation. The analysis of their temporal evolution unveiled the existence of a clear seasonal cycle closely related to that of the mixed layer depth. The phase shifts between SSH and SSB attain their minimum for deep mixed layers, which is different from zero. Besides, the spectral slope of the transfer function at scales shorter than 100 km fluctuates between k 21 and k 22. For deep mixed layers, it is close to k 21, as predicted by the surface quasigeostrophic (SQG) solution. At longer wavelengths, it is approximately constant under the different environmental conditions in all of the subbasins analyzed with the exception of the Gulf of Lions. The capability to observe sea surface temperature (SST) from satellites motivated the extension of this analysis to SST and SSH. Results showed a similar qualitative behavior but with larger phase shifts. In spite of the presence of a phase shift, even for deep mixed layers, results revealed that it is still possible to reconstruct surface dynamics from SST using a transfer function, provided that the mixed layer is deep enough. For the present study, a threshold value of 70 m was enough to identify the appropriate environmental conditions. In addition, the results revealed that a precise estimation of the transfer function significantly improves the reconstruction of the flow in comparison with the application of the classical SQG solution. 1. Introduction The Mediterranean Sea is a hot spot for climate change and plays a major role in the climate variability over Europe. This has led to efforts being made to better characterize its surface currents and to develop numerical models of its circulation. However, the dynamics in most of the Mediterranean are characterized by a Rossby radius of deformation of the order of km, which requires spatial resolutions higher than the available observations of surface velocities (;100 km). Besides, current observations of sea surface temperature (SST) do have the appropriate spatial resolutions. However, the lack of a direct connection between SST and surface currents has proven to be an obstacle difficult to overcome for their quantitative estimation from SST satellite observations. The diagnosis of surface velocities from SST images is based on the exploitation of the properties of upper-ocean Corresponding author address: Jordi Isern-Fontanet, Institut Catala de Ciencies del Clima, Doctor Trueta 203, Barcelona, Spain. jisern@ic3.cat dynamics. One of such properties is that small-scale patterns tend to be advected by larger-scale velocities in a turbulent flow. Then, mesoscale currents (;20 km) can be inferred, in principle, by tracking submesoscale patterns (1 10 km) in consecutive cloud-free SST images (e.g., Emery et al. 1986; Bowen et al. 2002). Another possibility consists of the exploitation of heat conservation and the inversion of the temperature equation to retrieve ocean currents using a sequence of SST images (e.g., Kelly 1989; Ostrovskii and Piterbarg 1995; Chen et al. 2008). A third approach, which has emerged during the last years, is based on the property that SST and surface currents have a similar complex phase in the upper ocean under the appropriate conditions (e.g., Isern-Fontanet et al. 2006a; LaCasce and Mahadevan 2006; Isern-Fontanet et al. 2008; Isern- Fontanet and Hascoët 2014). Indeed, Lapeyre and Klein (2006) put forward that potential vorticity (PV) in the upper ocean is correlated to the sea surface buoyancy (SSB), which locks the phase between SSB and the streamfunction. Contrary to the previous approaches, this property allows diagnosing velocities from a single SST field. DOI: /JPO-D Ó 2014 American Meteorological Society
2 MAY 2014 I S E R N - F O N T A N E T E T A L Indeed, Lapeyre and Klein (2006) and LaCasce and Mahadevan (2006) proposed that the mesoscale and submesoscale dynamics of the upper-ocean layers could be modeled using an effective version of the surface quasigeostrophic (SQG) equations. This implies that the surface streamfunction can be retrieved directly from SSB by convolving it with a kernel proportional to k 21, where k is the modulus of the two-dimensional wavevector. This kernel is usually known as the transfer function in signal processing theory. Other solutions can be found depending on the geometrical conditions (finite depth vs semi-infinite ocean; Tulloch and Smith 2006), the stratification (constant vs exponential; LaCasce 2012), or the relative contribution of the interior PV to respect the surface boundary condition (Lapeyre and Klein 2006; Klein et al. 2010). All of these solutions are also characterized by the phase locking between the streamfunction and SSB. It is worth mentioning that although the relevant dynamical variable in the PV inversion problem is the SSB, the field that can be remotely observed is the SST. Therefore, some additional assumptions regarding the alignment between SST and buoyancy gradients (Klein et al. 1998; Ferrari and Paparella 2003) have to be made if the objective is to retrieve ocean dynamics from SST. Because these solutions imply that SSB (or SST) shares the phase with the streamfunction, it is necessary to verify this condition before applying these transfer functions and identify under which conditions the assumption of zero phase shift holds. The validity of the SQG approach implies an energy spectrum with slopes of k 25/3 (Blumen 1978) and a spectrum of k 211/3 for sea surface height (SSH). This has motivated the study of its validity through the analysis of altimetric observations of SSH. Le Traon et al. (2008) showed that spectral slopes in the mesoscale band in high energy areas are significantly different from a k 25 law, predicted by quasigeostrophic (QG) turbulence, but they very closely follow the k 211/3 slope, which indicates that the surface QG is a much better dynamical framework than the QG turbulence theory to describe the ocean surface dynamics. However, results in low energy areas revealed shallower slopes (Xu and Fu 2011). Sasaki and Klein (2012) further investigated such areas using numerical simulations and showed that even in low energy areas SSH had spectra of the type k 24, close to SQG predictions, and suggested that instrumental noise could hide such slopes. Furthermore, they showed that the scales for which the k 24 slope could be seen depend on the energy level. In particular, for low energy areas with sea level variance on the order of 25 cm 2, the k 24 slope can be observed at scales shorter than 150 km. These studies have been complemented by the analysis of high-frequency radar just off the U.S. West Coast, which found spectral slopes equivalent to a SSH k 24 slope for scales smaller than 100 km (Kim et al. 2011). Nevertheless, this does not necessarily imply that velocities can be diagnosed from SST. On one side, Armi and Flament (1985) pointed out that different dynamical regimes could have similar spectral slopes and stressed the importance of the complex phases. On the other side, Klein and Hua (1990) showed that the emergence and evolution of the SST mesoscale variability depends on the mixed layer (ML) dynamics, which can introduce a phase shift between SST and currents. The Mediterranean Sea is a semienclosed basin with an excess of evaporation mainly compensated by the entrance of the relatively fresher waters from the Atlantic Ocean through the Strait of Gibraltar. Satellite images of SST have been widely used to investigate the circulation patterns in the Mediterranean Sea and track coherent vortices (e.g., Puillat et al. 2002; Hamad et al. 2005). However, the difficulty to quantitatively estimate surface velocities from SST patterns has limited its use to qualitative analysis. Consequently, quantitative studies have been based on altimetric (e.g., Larnicol et al. 2002; Pujol and Larnicol 2005; Isern-Fontanet et al. 2006c) and in situ measurements (e.g., Millot 1999; Malanotte-Rizzoli et al. 1999; Poulain et al. 2012, and references therein). The picture that emerges from both satellite and in situ observations reveals that incoming fresh waters circulate anticlockwise along the coast in the western Mediterranean Sea and the Levantine Basin (Millot and Taupier-Letage 2005; Poulain et al. 2012). The circulation in the Ionian Sea is still under debate, and different circulation schemata have been proposed (e.g., Millot and Taupier-Letage 2005; Pinardi et al. 2005; Poulain et al. 2012). The major discrepancies concern the existence of a jet that crosses the Ionian Sea, which is supported by the observations provided by surface drifters (Millot and Gerin 2010; Poulain et al. 2012), although some researchers consider it an artifact (Millot and Gerin 2010). In addition to these patterns, the instability of the inflow and local wind action often generate coherent vortices in several parts of the basin that tend to propagate following rather well-defined patterns in some areas of the Mediterranean Sea (Isern- Fontanet et al. 2006b). Although the use of SST has proven to be in good agreement with in situ measurements (e.g., Taupier-Letage et al. 2003; Isern-Fontanet et al. 2004) and very useful to get insight about ocean circulation in the area, there are some potential difficulties. By the end of summer, surface resident waters are warmer than incoming Atlantic water, while in winter the Atlantic water is warmer than the resident Mediterranean waters. This has two main implications
3 1408 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 for the diagnosis of surface velocities from SST: some features may not be visible during certain periods of the year, and the alignment between salinity and temperature gradients may vary (Taupier-Letage 2008). The main objective of this study is to characterize the transfer function between SSB and SSH and better determine the conditions under which SSB shares the complex phase with SSH in the Mediterranean Sea. Because SSB cannot be determined from direct observations, the study has been extended to include the analysis of SST to clearly identify the capability to reconstruct surface velocities from satellite observations of SST. The paper is organized as follows. Section 2 defines the theoretical framework used in this study. Section 3 describes the numerical simulations used for this study and the characteristics of the areas analyzed. Section 4 describes our results and section 5 discusses them. 2. Theoretical framework The principle of invertibility of PV (Hoskins et al. 1985) allows the diagnosis of the 3D dynamics of a balanced flow from the knowledge of PV in the ocean interior and density on the vertical boundaries. If we assume that the flow is in QG equilibrium, the problem consists of inverting the equation = 2 c 1 z! f0 2 c N 2 5 Q, (1) z where Q(x, z) is the PV anomaly, c(x, z) is the streamfunction of the flow, f 0 is the local Coriolis frequency, N(z) is the Brunt V ais al a frequency, x 5 (x, y), and $ [ ( x, y ). The appropriate boundary conditions can be derived from the hydrostatic equation. In particular, for the ocean surface they can be written as c f 0 z 5 b s, (2) s where b s (x) is the SSB and the subscript s stands for surface fields. This problem can be split into two solutions as proposed by Bretherton (1966) and Lapeyre and Klein (2006): an interior solution c i (x, z) obtained assuming zero surface buoyancy and nonzero interior PV, that is, b s 5 0 and Q 6¼ 0; and a surface solution c s (x, z) obtained assuming nonzero surface buoyancy and zero interior PV, that is, b s 6¼ 0 and Q 5 0. Then, the total solution to the inversion problem is the sum of both contributions, that is, c(x, z) 5 c s (x, z) 1 c i (x, z). (3) a. The transfer function formalism Lapeyre and Klein (2006) showed that the large-scale forcing in density and PV can lead to the property that the interior PV mesoscale anomalies are correlated to the surface buoyancy anomalies in the upper ocean. In that case, the PV anomaly can be separated as Q(x, z) j(z)b s (x), (4) with j(z) being a function that specifies the amplitude of PV anomaly. Then, the total solution to the PV inversion problem [Eqs. (1) and (2)] has the form ^c(k, z) 5 ~ F b (k, z) ^b s (k), (5) where the hat stands for the Fourier transform, k 5 (k x, k y ) is the wavevector, and F ~ b (k, z) is a transfer function that depends on the modulus of the wavevector k 5 kkk and the depth z. This transfer function can also be separated into interior and surface contributions, ~F bi (k, z) and F ~ bs (k, z), respectively. A simple analytical solution for F ~ bs (k, z) and F ~ bi (k, z) can be obtained for a finite depth H with constant stratification N 0 [ f 0 n 0. Then, the surface transfer function is given by ~F bs (k, z) 5 1 cosh[n 0 (H 1 z)k], (6) n 0 f 0 k sinh(n 0 kh) (Tulloch and Smith 2006) from which the classical SQG solution ~F bs (k, z) 5 1 n 0 f 0 k exp(n 0kz) (7) is recovered in the limit of infinite depth (Held et al. 1995). Notice that the surface behavior of Eq. (6) at short scales (large k) is similar to the classical SQG solution ~ F bs (k,0); k 21, while at large scales (small k)itis steeper, ~ Fbs (k,0); k 22. Besides, the interior solution for constant stratification and finite depth is given by j(z) ~F bi (k, z) 52! f 0 k (8) n 2 0 H2 (e.g., Klein et al. 2010), which corresponds to the baroclinic mode. At long scales, this solution tends to a constant, while at short scales it behaves as ~ F bi (k,0); k 22. The relative dominance of each solution can be separated
4 MAY 2014 I S E R N - F O N T A N E T E T A L by a critical wavelength that depends on the large-scale properties of the flow (Lapeyre 2009; Klein et al. 2010). b. The transfer function at the ocean surface Equation (5) implies that the streamfunction and surface buoyancy are in phase. However, previous studies have shown that there may exist a phase shift between the streamfunction and surface buoyancy (Isern-Fontanet et al. 2008; Hausmann and Czaja 2012). To take this into account, we introduce a phase shift 1 Du b (k) intoeq.(5) and focus at z 5 0. Then, ^c(k) 5 F b (k)e idu b (k) ^bs (k), (9) where i is the complex unit, c(x) is now the surface streamfunction, and F b (k) is the total transfer function at the ocean surface, that is, F b (k) [ ~ F bi (k,0)1 ~ F bs (k, 0). (10) The transfer function F b (k) is a real function that controls the amplitude of the resulting streamfunction and only depends on the modulus of the wavevector. Besides, the geostrophic streamfunction at the ocean surface is proportional to the SSH h(x), that is, c(x) 5 g f 0 h(x). (11) This implies that the transfer functions for surface buoyancy can be calculated as F b (k) g f 0 hj^hji k hj ^b s ji k, (12) where angle brackets with a subscript k indicate that the average is taken over those wavevectors with the same modulus. On the other side, the phase shift can be obtained as " Du b (k) 5 cos 21 <(^h ^b # s *) j^hjj ^b. (13) s j Equivalently, a transfer function F T (k) and phase shift Du T (k) can be defined for the SST anomaly T s (x). Then, Eq. (9) becomes ^c(k) 5 F T (k)e idu T (k) ^Ts (k). (14) 1 In our notation, a phase shift with no imaginary part has been absorbed into the definition of the transfer function, that is, F b (k) can take negative values. Notice that, if salinity would not contribute to surface buoyancy variability then, F T (k) } F b (k) and Du T (k) 5Du b (k). 3. Numerical simulations a. The Mediterranean Forecasting System In this study, we used the nowcasts provided by the Mediterranean Forecasting System (MFS) running at the Istituto Nazionale di Geofisica e Vulcanologia (INGV) ( The model is an implementation of the Nucleus for European Modelling of the Ocean Ocean Parallelise (NEMO OPA), version 3.2 (Madec et al. 1998; Madec 2008) that covers the entire Mediterranean Sea and part of the Atlantic Ocean with a spatial resolution of 1 / /168 and 71 unevenly spaced vertical levels (Tonani et al. 2008; Oddo et al. 2009). In the Atlantic, the model is nested within the monthly-mean climatological fields computed from the daily output of the 1 / /48 global model (Drevillon et al. 2008). The model is forced by momentum, water, and heat fluxes interactively computed by bulk formulae using the 6-h, 0.58 horizontalresolution operational analyses from the European Centre for Medium-Range Weather Forecasts (ECMWF) and the model-predicted surface temperatures (Tonani et al. 2008). In addition, evaporation is derived from the latent heat flux, precipitation is taken from monthly-mean Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) data (Xie and Arkin 1997), and river runoff is composed of monthly-mean climatological data for the seven major rivers in the area (Fekete et al. 1999; Raicich 1996). The flow through the Dardanelles Strait is parameterized as a river with monthly climatological net inflow rates taken from Kourafalou and Barbopoulos (2003). The system assimilates sea level anomaly (SLA), sea surface temperature, in situ temperature profiles from Voluntary Observing Ship (VOS) expendable bathythermographs (XBTs), in situ temperature and salinity profiles from Argo floats, and in situ temperature and salinity profiles from CTD using the three-dimensional variational data assimilation (3DVAR) scheme (Dobricic and Pinardi 2008). The background error correlation matrix is estimated from the temporal variability of parameters in a historical model simulation. Background error correlation matrices vary seasonally and in 13 regions of the Mediterranean that have different physical characteristics (Dobricic et al. 2007). The mean dynamic topography used for the assimilation of SLA has been computed by Dobricic (2005). Satellite SST data are used for the correction of surface heat fluxes.
5 1410 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 FIG. 1. (from top to bottom) RMS of the MLD, RMS of wind stress, and SLV. Squares indicate the location and extension of the four areas analyzed: 1) Algerian Basin, 2) Ionian Sea, 3) Levantine Basin, and 4) Gulf of Lions. b. Regions under study This study has been focused on four regions of the Mediterranean Sea representative of its circulation: the Algerian Basin, the Ionian Sea, the Levantine Basin, and the Gulf of Lions (Fig. 1). The dynamics and environmental conditions present during the period analyzed were characterized by the sea level variance (SLV), the root-mean-square (RMS) of wind stress, and the RMS of the mixed layer depth (MLD; see Table 1). Here, the MLD was estimated using the classical difference criterion with a reference depth of 10 m, as in de Boyer Montegut et al. (2004) and D Ortenzio et al. (2005). Notice that the SLV values shown in Table 1 correspond to low energy areas. The Algerian Basin is a highly energetic area characterized by the presence of eddies with sizes ranging from 30 to 150 km generated by the instability of the
6 MAY 2014 I S E R N - F O N T A N E T E T A L TABLE 1. Spatially averaged (indicated by the angle brackets) SLV hs 2 hi, RMS of wind stress, and RMS of MLD for the four boxes seen in Fig. 1. No. Area hs 2 h i (cm2 ) ht 2 i 1/2 (N m 22 ) hh 2 i 1/2 (m) 1 Algerian Basin Ionian Sea Levantine Basin Gulf of Lions fresh waters of Atlantic origin that flow along the Algerian coast. These eddies tend to follow well-defined patterns characterized by a recirculation loop approximately enclosed within the box used in our analysis (Isern-Fontanet et al. 2006b). The presence of such eddies produced the largest spatially averaged SLV (Table 1). On the contrary, the Ionian Sea had the lowest SLV of all four boxes, mainly concentrated in the southern part of the box. Its dynamics are characterized by a lower concentration of intense eddies (Isern- Fontanet et al. 2006b) and the presence of a jet crossing the box from east to west (e.g., Poulain et al. 2012). The Levantine Basin had dynamical characteristics that resemble those of the Algerian Basin. It is characterized by the presence of energetic eddies generated in part by the instability of the waters that circulate counterclockwise along the coast. However, the length of the time series here analyzed was not long enough to produce more homogeneous patterns of SLV. Indeed, the box in this area was dominated by the presence of a highly energetic eddy in its southwestern corner, which is evident in Fig. 1. Finally, the Gulf of Lions, which is an area of deep-water formation, was characterized by the strongest wind forcing, very deep MLD, and quite active dynamics, mainly in its southwestern part due to the contribution of Algerian eddies. In addition to these eddies, this area is characterized by the presence of smaller vortices (;50 km) generated by the instability of the Liguro Provençal Current that follows the continental shelf (e.g., Rubio et al. 2005). Salinity gradients play a major role in the dynamics of the Mediterranean Sea. This made it necessary to explore the contribution of salinity and temperature to surface buoyancy. In particular, the alignment and compensation in density between salinity and temperature gradients were analyzed using the complex ratio given by Ferrari and Paparella (2003): R 5 a( x T s 1 i y T s ) b( x S s 1 i y S s ), (15) where a, 0 is the thermal expansion coefficient, b. 0 is the haline contraction coefficient, and S s is the sea surface salinity (SSS). The phase f of the complex density ratio R quantifies the degree of alignment between SST and SSS gradients, while its magnitude quantifies the relative strength of the gradients. For f 5 0radorf 56p, there is a thermohaline alignment. Figure 2 shows the probability density functions (PDF) of the FIG. 2. PDF of the alignment between temperature and salinity gradients. Dotted line corresponds to the Algerian Basin, dashed line to the Ionian Sea, dashed dotted line to the Levantine Basin, and solid line to the Gulf of Lions.
7 1412 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 alignment between salinity and temperature gradients observed in the simulations. The resulting PDFs were bimodal with peaks at f 5 0 and f 5 p, indicating a strong tendency toward the alignment. However, the relative importance of each peak depended on the area analyzed. In the western Mediterranean Sea (Algerian Basin and the Gulf of Lions), it was more likely to find alignments around f 5 0, while in the eastern Mediterranean Sea (Ionian Sea and Levantine Basin), it was more likely to find alignments around f 5 p. The probability of finding a misalignment between salinity and density gradients larger than p/4 was estimated to be approximately 26%. If only MLs deeper than 70 m were considered, such probability reduced to values between 21% and 23%. Alignment properties in all boxes were similar, although in the Ionian Basin the probability of finding a misalignment was found to be slightly larger. The probability of having compensated fronts, that is, jrj 5 1andf 5 0 rad, was estimated to be less than 0.1%. 4. Results a. Phase shift between SSB and SSH The temporal evolution of the phase shift between SSB and SSH for different k was investigated for all boxes using a spectral correlation coefficient r(k) based on Eq. (13): r(k) [ hcosdu b i k. (16) Results revealed a clear seasonal cycle with the largest shifts between 30 and 200 km during summer (May September), which was the period of the year with the shallowest ML and lower winds (see the example in Fig. 3). Because wind stress was highly variable, its similarity with the MLD was more evident once the daily winds were smoothed using a moving mean with an amplitude of 45 days. For those days that had smaller phase shifts, a higher degree of homogeneity across k in the western Mediterranean Sea (Gulf of Lions and Algerian Basin) than in the eastern Mediterranean (not shown) could be observed. Indeed, for the Ionian Sea and the Levantine Basin, it was evident that the phase shifts tended to be smaller for wavelengths larger than 100 km. To complement the spectral analysis, a streamfunction c b (x) was built combining the phases of SSB with the spectrum of SSH for the spectral band between 20 and 300 km, that is, ^c b [ g j^hj f 0 j ^b ^b s. (17) s j By construction, such streamfunction only differed from Eq. (11) by the relative position of flow patterns. Both streamfunctions were compared, once properly filtered with the same cutoff wavelengths, using the lineal correlation between vorticities as a measure (Fig. 4). Results also revealed a clear seasonal cycle characterized by higher correlations during winter. However, correlations never reached 1, suggesting the presence of a phase shift between SSB and SSH. This is in agreement with results shown in Fig. 3, which showed that r(k), 1. Results also unveiled differences between basins and from year to year. Indeed, correlations were larger or equal to 0.9 in the Algerian and Levantine Basins between December and March, while in the Ionian Basin they were around 0.8 during winter In the Gulf of Lions, the period of correlations higher or equal to 0.9 was shorter. The interannual variability was evident, for example, in the Ionian Basin, where higher correlations were obtained from January to April 2009 than during the same period in Interestingly, a difference between these two periods was the depth of the ML, which was thicker in The difficulty in reaching correlations higher than 0.9 still indicated the existence of a phase shift between SSB and SSH. To estimate it, daily phase shifts associated with spatial correlations larger or equal to 0.9 were averaged (Fig. 5). Results revealed that even for high spatial correlations, there was a phase shift between SSB and SSH of the order of ;p/6. Although it was quite similar in all boxes for wavelengths shorter than 100 km, there were significant differences at longer wavelengths. In particular, the Ionian Sea and Levantine Basin exhibited smaller phase shifts between 100 and 250 km than the boxes in the western Mediterranean. Unsurprisingly, higher threshold values of correlation led to lower phase shifts. b. The transfer function between SSB and SSH In addition to the phase shift, the characteristics and temporal evolution of the transfer function were investigated. The transfer function was computed from daily SSB and SSH fields using Eq. (12). Results shown in Fig. 3 revealed that it was characterized by negative slopes particularly for wavelengths shorter than 100 km, and a clear seasonal cycle characterized by steeper slopes for wavelengths shorter than 100 km during summer. A qualitative correspondence between the patterns observed in the transfer function and those observed for the temporal evolution of the phase shift could be found. In general, all boxes exhibited similar features with the exception of the Gulf of Lions, which had less marked patterns (not shown).
8 MAY 2014 I S E R N - F O N T A N E T E T A L FIG. 3. (left) Time series for the mean MLD hhi (thick dashed line), mean wind stress modulus hti (thin solid line), and smoothed mean wind stress modulus (thick solid line). (middle) Spectral correlation between SSB and SSH. (right) Modulus of the transfer function between SSB and SSH. Results correspond to the box in the Levantine Basin. The transfer function associated with the smallest phase shifts was estimated by averaging the observations for those days with vorticity correlations larger or equal to 0.9 as before (Fig. 5). From results, two different bands in most basins could be clearly identified: one characterized by an SQG-like slope, that is, F b (k) ; k 21, between 20 and 100 km (see Table 2); and another characterized by a plateau, that is, F b (k) ; C, for wavelengths larger than 100 km. Nevertheless, the situation was slightly different in the Gulf of Lions, where the F b (k) ; k 21 extended along the whole range of wavelengths analyzed. c. Dependence on environmental conditions The strong seasonal dependence of the phase shift and transfer function pushed us to investigate which parameters better characterized such evolution. We focused on the spatial averages of the MLD hhi, wind stress hti, and the stratification just below the ML n,
9 1414 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 FIG. 4. Spatial correlations between the vorticity diagnosed from SSH and the vorticity diagnosed from the streamfunctions derived from SSB (black line) and SST (gray line) using the spectral amplitudes of SSH. because these parameters play a key role in the deepening of the ML (Klein and Hua 1988). Results (Fig. 6) showed a clear dependence on the MLD, while the dependence on wind was much less clear. Notice, however, that when averaging wind using a running mean with an amplitude of 45 days, results are similar to those for the MLD. It is worth mentioning that in the Gulf of Lions deep ML did not always imply high correlations, as is evident by the set of circles with MLD larger than 70 m and correlations lower than 0.8 seen in Fig. 6. On the other side, the results were less clear using the stratification than with the MLD. The averaged phase shift and transfer function were explored using the MLD as a criterion to classify the daily results. In particular, the mean transfer function was obtained under different MLD ranges, as is shown in Fig. 7. For MLD deeper than 70 m, the results were similar to those obtained selecting correlations larger or equal to 0.9 (Fig. 5), except for the Gulf of Lions, which exhibited larger phase shifts and more of a departure from the mean transfer functions than that observed for the other boxes. The reduction of the MLD threshold value used to select which days were averaged lead to an increase of the phase shift and a change of the transfer function for wavelengths shorter than 100 km from F b (k) ; k 21 to F b (k) ; k 22. The plateau for wavelengths longer than 100 km did not change. The kinetic energy spectra estimated from SSH was compared to the spectra of surface buoyancy in order to further understand the variation of the transfer function. Notice that the kinetic energy spectra is related to SSH and SSB through E k2 j^hj f 0 2 g k2 jf b j 2 j ^b 2 s j. (18) This implies that the kinetic energy spectra and surface buoyancy spectra will have the same slope if F b (k) ; k 21. As for the transfer functions, spectra were averaged based on the MLD (Fig. 8). For MLD deeper than 70 m, surface buoyancy and kinetic energy had the same slopes at wavelengths shorter than 100 km, as expected by the
10 MAY 2014 I S E R N - F O N T A N E T E T A L FIG. 5. (left) Spectral dependence of the mean phase shift between SSB and SSH. (right) Mean transfer function modulus. The mean has been taken over those days with correlations of vorticity larger or equal than 0.9. Dotted line corresponds to the Algerian Basin, dashed line to the Ionian Sea, dashed dotted line to the Levantine Basin, and solid line to the Gulf of Lions. k 21 slope of the corresponding transfer function (Fig. 7). For shallow MLD, the main difference corresponded to the kinetic energy spectra, which was steeper than for deep MLD at wavelengths shorter than 100 km approximately. In both situations the spectra of surface buoyancy was qualitatively similar, which originated the ;k 22 slope in the transfer function. On the contrary, the kinetic energy had a shallower spectrum at larger wavelengths for the different environmental conditions. d. Flow reconstruction from SST The same analysis applied to surface buoyancy was applied to SST. Phase shifts and transfer function temporal evolutions were qualitatively similar to those observed for surface buoyancy with the main difference that, in general, phases at wavelengths shorter than 20 km tended to be orthogonal and the annual cycle stronger (not shown). On the contrary, the mean transfer function estimated for different MLD showed the same patterns as for buoyancy, that is, a F T (k) ; k 21 behavior for wavelengths shorter than 100 km and F T (k) ; C for wavelengths longer than 100 km. The observed transfer function in the Gulf of Lions also revealed a SQG-like behavior for the whole range of wavelengths analyzed. As before, a transfer function c T (x)equivalent to Eq. (17) was defined: ^c T [ g j^hj f 0 j ^T ^T s j s. (19) Resulting vorticity correlations time series (Fig. 4) were systematically lower than for SSB, although both followed a similar seasonal cycle. Correlations in the Ionian Sea were quite small (less than 0.8) in comparison to the other boxes. The capability to diagnose surface dynamics from SST was investigated, focusing on the reconstruction capabilities for the surface streamfunction, surface velocity v 5 (u, y), and vorticity and its dependence on the MLD. In particular, the surface streamfunction was derived from SST and SSH using Eq. (19), focusing on two different bands: the band with SQG-like behavior ( km) and the band with a flat slope ( km). Figure 9 shows the dependence of the averaged correlations with the MLD. Averaged correlations for all variables showed the same qualitative behavior as in Fig. 6: a rapid increase of correlation with the MLD until reaching a maximum value for MLD between 50 and 70 m (Fig. 9). In addition, the capability to reconstruct dynamical variables was better for wavelengths larger than 100 km than for shorter wavelengths, indicating a smaller phase shift at these wavelengths in agreement with the results reported in previous sections. 5. Discussion It has been shown that the transfer function between the streamfunction and the surface buoyancy does not follow the classical SQG solution for large scales, with TABLE 2. Spectral slopes between 20 and 100 km of the transfer function shown in Fig. 5 and the associated uncertainties in the four boxes analyzed. No. Area Slope 1 Algerian Basin Ionian Sea Levantine Basin Gulf of Lions
11 1416 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 FIG. 6. Scatterplot of the correlations between the vorticity diagnosed from SSB and model s vorticity for the band between 20 and 300 km and different environmental parameters: ML depth hhi, mean wind stress hti, and stratification below the ML n. the exception of the Gulf of Lions, and the existence of a phase shift between the streamfunction and surface buoyancy has been evidenced. Indeed, results have shown that the amplitude of the mean transfer function for both, surface density and SST, has an SQG-like response for wavelengths shorter than 100 km and deep ML in all the regions analyzed. On the contrary, it has almost flat amplitude for longer wavelengths in most boxes. The qualitative comparison of these results with the analytical solution obtained assuming a finite depth and a constant stratification [Eqs. (6) and (8)] suggests a dominance of the interior solution at large scales and of the surface solution at smaller scales, in agreement with theoretical studies (Lapeyre and Klein 2006; Lapeyre 2009). The above results were observed mainly during winter, when a deep ML can be found. On the contrary, observations for shallow ML reveal a different response characterized by steeper slopes in the amplitude of the transfer function at short wavelengths, that is, F b (k) ; k 22. A recent study by Callies and Ferrari (2013) showed that the interior QG solution dominates the flow at scales larger than 20 km in highly energetic areas, such as the Gulf Stream (s 2 h. 400 cm2 ). But in low energy areas, such as the eastern subtropical North Pacific (s 2 h ; 21 cm2 ), results were inconsistent with the dominance of both surface and interior solutions. The sea level variance of the boxes analyzed in this study ranges between s 2 h ; 21 cm2 and s 2 h ; 34 cm2 (Table 1), which is closer to the variance in the eastern subtropical North Pacific than in the Gulf Stream area. Consequently, we would expect that the dynamics here observed would be closer to that of the subtropical North Pacific. However, contrary to the observations done by Callies and Ferrari (2013), our results are statistically consistent with the predictions of the SQG solution with constant stratification, at least for scales shorter than 100 km. The slightly higher variability in the Mediterranean Sea could be at the origin of these differences. Besides, our results are based on the numerical simulation of the Mediterranean circulation, which may not include all the physical process that contributes to the energy. This implies the need to extend our study to real data in order to confirm the results here found and compare them with the observations of Callies and Ferrari (2013). The deviation of the streamfunction from the SQG model, particularly at larger scales, has an impact on the performance of the flow reconstruction from satellite SST. The use of the SQG model as in Isern-Fontanet et al. (2006b, 2008) would imply an overestimation of the relative importance of the large scales to respect small scales. Indeed, the comparison of the performance of the SQG reconstruction with the full knowledge of the transfer function represented dramatic improvements. In particular, the correlation for vorticity obtained applying the SQG approach and using Eq. (19) for the range of wavelengths between 40 and 200 km were compared. Results shown in Fig. 10 confirm such an improvement, which can be so important that correlations can increase from values below 0.5 to values on the
12 MAY 2014 I S E R N - F O N T A N E T E T A L FIG. 7. (left) Spectral dependence of the mean phase shift between surface buoyancy and SSH. (right) Mean transfer function modulus. (from top to bottom) Observations for MLD shallower than 20 m, observations for MLD between 20 and 70 m, and results for MLD deeper than 70 m. Dotted line corresponds to the Algerian Basin, dashed line to the Ionian Sea, dashed dotted line to the Levantine Basin, and solid line to the Gulf of Lions. order of 0.8. These results suggest the introduction of an heuristic transfer function given by ( Cg k 2g 1 if k # k c F T (k) 5 C g k g 2g 2 1 c k 2g (20) 2 if k $ k c and based on the family of two-dimensional fluid models proposed by Pierrehumbert et al. (1994), where C g is a constant, and g 1 and g 2 are the spectral slopes in the two bands separated by the wavevector k c. The four parameters of this model could be estimated from simultaneous satellite observations of SST and SSH.
13 1418 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 FIG. 8. Mean spectra of SSB (dotted line), SST (dashed line), and SSH (solid line) for different MLD in the Levantine Basin. The spectra of SSH have been multiplied by k 2 for comparison with SST and SSB. Indeed, current observations of SSH and SST have synergetic characteristics. On one side, along-track altimetric measurements of SSHs are very well suited to quantify across-track currents and have been widely used to estimate turbulence spectra in the ocean (e.g., Le Traon et al. 1990; Stammer 1997; Xu and Fu 2011), but their sampling geometry and the number of available platforms strongly constrain the reconstruction of twodimensional streamfunctions (Pascual et al. 2006). On the other side, SST measurements provided by infrared radiometers are able to locate ocean structures with small phase shifts, if the environmental conditions are appropriate. This suggests to derive the transfer function between simultaneous SST and SSH measurements to FIG. 9. Correlation between reconstructed fields from SST using the optimal transfer function and fields provided by the model obtained by binning the MLD and taking the median in each bin. Solid line corresponds to wavelengths longer than 100 km and dashed line to wavelengths between 20 and 100 km. (from top to bottom) Streamfunction c, vorticity z, zonal velocity u, and meridional velocity y.
14 MAY 2014 I S E R N - F O N T A N E T E T A L FIG. 10. Scatterplot of the correlations for vorticity diagnosed from SST using the SQG model z SQG T and using the full knowledge of the transfer function z T for the band between 40 and 200 km in the different basins. obtain its best estimation. Furthermore, along-track SSH measurements could then be used to quantify if the reconstruction is good enough to be used for scientific or operational applications by correlating acrosstrack velocities derived from SSH with those diagnosed from the combination of SSH and SST. Our results have shown that the different response in the transfer function originates from a different slope of SSH (Fig. 8). Such a result should be verified with satellite observations of SSH. However, the main difficulty is the low signal-tonoise ratio of available altimetric measurements in the Mediterranean Sea, which could hide such behavior. Another approach to exploit the complementary behavior between SST and SSH is to use the information provided by each field to reconstruct the flow at different scales. Indeed, our results reveal that surface velocities can be reconstructed from SST fields in the band between 20 and 100 km with correlations on the order of 0.8, if the ML is deep enough (Fig. 9). Furthermore, the transfer function for this band is close to the SQG predictions (Fig. 7). This suggests that the spatial resolution of SSH observations could be improved by merging the low-resolution altimetric measurements (l. 100 km) with the high-resolution (20, l, 100 km) streamfunctions derived from SST using the SQG model as a first approach or incorporating the along-track SSH spectrum. This idea is similar to the underlying approach proposed by Gaultier et al. (2013). In particular, their procedure is based on extracting high-resolution information from the SST gradients and finite size Lyapunov exponents computed from altimetry (e.g., d Ovidio et al. 2009) to correct the low-resolution altimetric maps. The underlying concept to retrieve the high-resolution information about the velocity field has some similarities with the technique proposed by Turiel et al. (2005) and Isern-Fontanet et al. (2007). Nevertheless, the method to reconstruct the velocities from this information was different. As it has been outlined in the introduction, there are alternative methods to derive surface currents from SST. Among different techniques, the maximum cross correlation (MCC) method is by far the most widely used (Emery et al. 1986; Bowen et al. 2002; Barton 2002; Notarstefano et al. 2008). The basic idea consists of computing spatial correlations over windows between consecutive images to estimate their displacement. In general, the MCC method acts poorly in regions of uniform SST and presents some difficulties to estimating the velocity along the front of the scalar field (e.g., Zavialov et al. 2002). Consequently, the resulting velocity fields tend to be sparse, which requires merging of the estimated velocities with other measurements such as altimetry to retrieve the two-dimensional velocity field (Wilkin et al. 2002). Other approaches used to estimate velocities from SST are the constrained optical flow methods to solve the heat conservation equation (e.g., Kelly 1989; Vigan et al. 2000) orvariationalfilterandinterpolation techniques (Afanasyev et al. 2002). All these techniques rely on the availability of cloud-free very highresolution (;1 km) sequences of images over short enough time periods. This latter restriction is imposed by the lack of absolute conservation of SST. These requirements impose strong restrictions on their applicability and geographically and seasonally limit the regions over which velocities can be estimated. On the contrary, our approach based on the combination of SST and SSH is able to provide dense velocity fields from an SST snapshot. The diagnosis of subsurface fields requires separating the interior and surface contributions due to their different vertical variations. Besides, the empirical determination of the transfer function would be only able to provide the total contribution implying the need to find a criterion to separate both contributions. A possible solution could be the method of Wang et al. (2013), who propose to first determine the surface contribution using surface buoyancy (or SST) through the inversion of the problem given by Eqs. (1) and (2) with Q 5 0 and b s 6¼ 0 and then find the amplitude of the barotropic and first baroclinic mode by matching to the surface streamfunction (or SSH) resulting from the subtraction of the surface component. The underlying assumption needed to reconstruct the flow from SST is that Du T (k) 5 0, that is, the surface
15 1420 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 44 streamfunction or the SSH has to be in phase with surface buoyancy or SST. However, our results have shown that this is not the case. Even under the most favorable environmental conditions there is a phase shift present, although it is small enough to allow the reconstruction of the flow with high correlations. A similar result can be deduced from the study of Isern-Fontanet et al. (2008), who showed that the spectral correlations were never equal to 1 (see their Fig. 2). Assuming that the correct transfer function can be deduced from satellite observations, the phase shift represents the major limitation to reconstruct the flow from SST measurements. Lapeyre and Klein (2006) suggested that PV and SSB are in phase [Eq. (4)], implying that the solution of the PV problem can be written as the convolution between a real transfer function and SSB. However, the failure of this hypothesis could introduce a phase shift between SSH and SSB. Following Lapeyre and Klein (2006), the linear correlation between PV below the ML and SSB was computed and then averaged for the different MLD used to analyze the results (Fig. 11). Results clearly revealed that correlations were higher for the deep ML than for the shallow ML, but even under the most favorable environmental conditions, values larger than 0.5 were constrained to the upper 100 m. Such observations could explain the phase shift between SSH and SSB observed in the numerical simulations because the former responds not only to SSB but also to the interior PV. Some studies have investigated and characterized the phase shift between SST and SSH and the mechanisms associated with it. Hausmann and Czaja (2012) investigated the phase shift between SST and SSH associated with mesoscale eddies in the global ocean. They showed that SST anomalies tend to be westward shifted with respect to the eddy core for both cyclonic and anticylonic eddies and poleward (equatorward) for anticyclonic (cyclonic) eddies, recalling how vortices propagate in the ocean (Cushman-Roisin and Beckers 2011). A qualitative exploration of our simulations also revealed a systematic shift associated with eddies, although our analysis did not allow us to determine any preferred direction. Notice that the propagation paths of Mediterranean vortices form complex patterns due to the interaction between them and the constrains imposed by the topography (Isern-Fontanet et al. 2006b). Besides, the strong dependence on the MLD of the observed phase shifts suggests that it can play a nonnegligible role. Indeed, Klein and Hua (1990) investigated the role of the ML dynamics in the emergence of SST patterns and identified two different regimes. One regime is characterized by the ML deepening that generated the spatial variability of the SST that was a linear combination of subsurface temperatures and relative FIG. 11. Mean spatial correlation between surface buoyancy anomaly and PV anomaly in the Levantine Basin for MLD deeper than 70 m (solid line) and MLD shallower than 20 m (dashed line). vorticity depending on the strength of the ML deepening. The other regime, which can be found after the wind stops, is characterized by the advection of the SST by the flow. This process generates fronts and energetic smaller scales. Our results have shown that it is possible to reconstruct surface currents from SST observations in the Mediterranean Sea. However, the need to minimize the phase shift between SST and SSH implies that the reconstruction is constrained to winter, when less cloudfree images are available. This would force the use of the geostationary satellite to minimize the impact of clouds. Besides, our results have shown that the reconstruction in the Ionian Basin was not possible during most of the analyzed period. Because the performance of the reconstruction from SSB was similar to the results in the other areas, it points to the lower capability of SST to identify SSB patterns. 6. Conclusions The conditions under which the SST can be used to diagnose surface velocities in the Mediterranean Sea have been identified. It has been shown that the reconstruction of upper-ocean dynamics from SST works better for larger scales than shorter scales and for ML deeper than approximately 70 m. In addition, results have confirmed that the transfer function has a SQG-like behavior for wavelengths shorter than 100 km. However, in most of the studied areas the transfer function differs
Three-dimensional reconstruction of oceanic mesoscale currents from surface information
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jc004692, 2008 Three-dimensional reconstruction of oceanic mesoscale currents from surface information Jordi Isern-Fontanet, 1 Guillaume Lapeyre,
More informationProject of Strategic Interest NEXTDATA. Deliverables D1.3.B and 1.3.C. Final Report on the quality of Reconstruction/Reanalysis products
Project of Strategic Interest NEXTDATA Deliverables D1.3.B and 1.3.C Final Report on the quality of Reconstruction/Reanalysis products WP Coordinator: Nadia Pinardi INGV, Bologna Deliverable authors Claudia
More informationDo altimeter wavenumber spectra agree with interior or surface. quasi-geostrophic theory?
Do altimeter wavenumber spectra agree with interior or surface quasi-geostrophic theory? P.Y. Le Traon*, P. Klein*, Bach Lien Hua* and G. Dibarboure** *Ifremer, Centre de Brest, 29280 Plouzané, France
More informationGlobal Variability of the Wavenumber Spectrum of Oceanic Mesoscale Turbulence
802 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 41 Global Variability of the Wavenumber Spectrum of Oceanic Mesoscale Turbulence YONGSHENG XU AND LEE-LUENG FU Jet Propulsion Laboratory,
More informationEddy-resolving Simulation of the World Ocean Circulation by using MOM3-based OGCM Code (OFES) Optimized for the Earth Simulator
Chapter 1 Atmospheric and Oceanic Simulation Eddy-resolving Simulation of the World Ocean Circulation by using MOM3-based OGCM Code (OFES) Optimized for the Earth Simulator Group Representative Hideharu
More information3. Midlatitude Storm Tracks and the North Atlantic Oscillation
3. Midlatitude Storm Tracks and the North Atlantic Oscillation Copyright 2006 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without permission. EFS 3/1 Review of key results
More informationNorth Atlantic circulation in three simulations of 1/12, 1/25, and 1/50
North Atlantic circulation in three simulations of 1/12, 1/2, and 1/ Xiaobiao Xu and Eric Chassignet Center for ocean-atmospheric prediction studies Florida State University Motivation Numerical models
More informationWhat Vertical Mode Does the Altimeter Reflect? On the Decomposition in Baroclinic Modes and on a Surface-Trapped Mode
NOVEMBER 2009 L A P E Y R E 2857 What Vertical Mode Does the Altimeter Reflect? On the Decomposition in Baroclinic Modes and on a Surface-Trapped Mode GUILLAUME LAPEYRE Laboratoire de Météorologie Dynamique/IPSL,
More informationVertical velocities in the upper ocean from glider and altimetry data 1
Vertical velocities in the upper ocean from glider and altimetry data 1 In this poster we show results on the combination of new glider technology data with altimetry observations to diagnose vertical
More informationCoastal Ocean Circulation Experiment off Senegal (COCES)
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Coastal Ocean Circulation Experiment off Senegal (COCES) Pierre-Marie Poulain Istituto Nazionale di Oceanografia e di Geofisica
More informationNon-linear patterns of eddy kinetic energy in the Japan/East Sea
Non-linear patterns of eddy kinetic energy in the Japan/East Sea O.O. Trusenkova, D.D. Kaplunenko, S.Yu. Ladychenko, V.B. Lobanov V.I.Il ichev Pacific Oceanological Institute, FEB RAS Vladivostok, Russia
More informationA new global surface current climatology, with application to the Hawaiian Island region. Rick Lumpkin
A new global surface current climatology, with application to the Hawaiian Island region Rick Lumpkin (Rick.Lumpkin@noaa.gov) Drogue presence reanalysis Left: time-mean zonal currents from drifters and
More informationLecture 1. Amplitude of the seasonal cycle in temperature
Lecture 6 Lecture 1 Ocean circulation Forcing and large-scale features Amplitude of the seasonal cycle in temperature 1 Atmosphere and ocean heat transport Trenberth and Caron (2001) False-colour satellite
More informationFeature resolution in OSTIA L4 analyses. Chongyuan Mao, Emma Fiedler, Simon Good, Jennie Waters, Matthew Martin
Feature resolution in OSTIA L4 analyses Chongyuan Mao, Emma Fiedler, Simon Good, Jennie Waters, Matthew Martin GHRSST XVIII, Qingdao, China, 5-9 June 2017 Talk outline Introduction NEMOVAR in OSTIA Methods
More informationComparison Figures from the New 22-Year Daily Eddy Dataset (January April 2015)
Comparison Figures from the New 22-Year Daily Eddy Dataset (January 1993 - April 2015) The figures on the following pages were constructed from the new version of the eddy dataset that is available online
More informationOcean currents from altimetry
Ocean currents from altimetry Pierre-Yves LE TRAON - CLS - Space Oceanography Division Gamble Workshop - Stavanger,, May 2003 Introduction Today: information mainly comes from in situ measurements ocean
More informationOn the relationship between the water mass pathways and eddy variability in the Western Mediterranean Sea
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2005jc003174, 2007 On the relationship between the water mass pathways and eddy variability in the Western Mediterranean Sea E. K. Demirov 1 and
More informationthat individual/local amplitudes of Ro can reach O(1).
Supplementary Figure. (a)-(b) As Figures c-d but for Rossby number Ro at the surface, defined as the relative vorticity ζ divided by the Coriolis frequency f. The equatorial band (os-on) is not shown due
More informationThe Planetary Circulation System
12 The Planetary Circulation System Learning Goals After studying this chapter, students should be able to: 1. describe and account for the global patterns of pressure, wind patterns and ocean currents
More informationOptimal Spectral Decomposition (OSD) for GTSPP Data Analysis
Optimal Spectral Decomposition (OSD) for GTSPP Data Analysis Peter C Chu (1),Charles Sun (2), & Chenwu Fan (1) (1) Naval Postgraduate School, Monterey, CA 93943 pcchu@nps.edu, http://faculty.nps.edu/pcchu/
More informationAPPENDIX B PHYSICAL BASELINE STUDY: NORTHEAST BAFFIN BAY 1
APPENDIX B PHYSICAL BASELINE STUDY: NORTHEAST BAFFIN BAY 1 1 By David B. Fissel, Mar Martínez de Saavedra Álvarez, and Randy C. Kerr, ASL Environmental Sciences Inc. (Feb. 2012) West Greenland Seismic
More informationThe general circulation: midlatitude storms
The general circulation: midlatitude storms Motivation for this class Provide understanding basic motions of the atmosphere: Ability to diagnose individual weather systems, and predict how they will change
More informationModeling and Parameterizing Mixed Layer Eddies
Modeling and Parameterizing Mixed Layer Eddies Baylor Fox-Kemper (MIT) with Raffaele Ferrari (MIT), Robert Hallberg (GFDL) Los Alamos National Laboratory Wednesday 3/8/06 Mixed Layer Eddies Part I: Baroclinic
More informationTransient and Eddy. Transient/Eddy Flux. Flux Components. Lecture 3: Weather/Disturbance. Transient: deviations from time mean Time Mean
Lecture 3: Weather/Disturbance Transients and Eddies Climate Roles Mid-Latitude Cyclones Tropical Hurricanes Mid-Ocean Eddies Transient and Eddy Transient: deviations from time mean Time Mean Eddy: deviations
More informationThe feature of atmospheric circulation in the extremely warm winter 2006/2007
The feature of atmospheric circulation in the extremely warm winter 2006/2007 Hiroshi Hasegawa 1, Yayoi Harada 1, Hiroshi Nakamigawa 1, Atsushi Goto 1 1 Climate Prediction Division, Japan Meteorological
More informationPRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Event Response
PRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK June 2014 - RMS Event Response 2014 SEASON OUTLOOK The 2013 North Atlantic hurricane season saw the fewest hurricanes in the Atlantic Basin
More informationOCEANIC SUBMESOSCALE SAMPLING WITH WIDE-SWATH ALTIMETRY. James C. McWilliams
. OCEANIC SUBMESOSCALE SAMPLING WITH WIDE-SWATH ALTIMETRY James C. McWilliams Department of Atmospheric & Oceanic Sciences Institute of Geophysics & Planetary Physics U.C.L.A. Recall the long-standing
More informationNear-Surface Dispersion and Circulation in the Marmara Sea (MARMARA)
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Near-Surface Dispersion and Circulation in the Marmara Sea (MARMARA) Pierre-Marie Poulain Istituto Nazionale di Oceanografia
More informationDiagnosis of vertical velocities in the upper ocean from high resolution sea surface height
Click Here or Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L12603, doi:10.1029/2009gl038359, 2009 Diagnosis o vertical velocities in the upper ocean rom high resolution sea surace height P. Klein,
More informationOcean Mixing and Climate Change
Ocean Mixing and Climate Change Factors inducing seawater mixing Different densities Wind stirring Internal waves breaking Tidal Bottom topography Biogenic Mixing (??) In general, any motion favoring turbulent
More informationThermohaline variability of the North Ionian and South Adriatic Sea in
Thermohaline variability of the North Ionian and South Adriatic Sea in 2013-2016 Κλικ για προσθήκη κειμένου Zoi kokkini, P. - M. Poulain, R. Gerin, E. Mauri and G. Notarstefano 6th Euro-Argo Users Meeting,
More information2013 ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Cat Response
2013 ATLANTIC HURRICANE SEASON OUTLOOK June 2013 - RMS Cat Response Season Outlook At the start of the 2013 Atlantic hurricane season, which officially runs from June 1 to November 30, seasonal forecasts
More informationToward Environmental Predictions MFSTEP. Executive summary
Research Project co-funded by the European Commission Research Directorate-General 5 th Framework Programme Energy, Environment and Sustainable Development Contract No. EVK3-CT-2002-00075 Project home
More informationPS4a: Real-time modelling platforms during SOP/EOP
PS4a: Real-time modelling platforms during SOP/EOP Mistral Tramontane Bora Etesian Major sites of dense water formation Major sites of deep water formation influence of coastal waters Chairs: G. Boni,
More informationRecent warming and changes of circulation in the North Atlantic - simulated with eddy-permitting & eddy-resolving models
Recent warming and changes of circulation in the North Atlantic - simulated with eddy-permitting & eddy-resolving models Robert Marsh, Beverly de Cuevas, Andrew Coward & Simon Josey (+ contributions by
More informationUpper Ocean Circulation
Upper Ocean Circulation C. Chen General Physical Oceanography MAR 555 School for Marine Sciences and Technology Umass-Dartmouth 1 MAR555 Lecture 4: The Upper Oceanic Circulation The Oceanic Circulation
More information2. Baroclinic Instability and Midlatitude Dynamics
2. Baroclinic Instability and Midlatitude Dynamics Midlatitude Jet Stream Climatology (Atlantic and Pacific) Copyright 26 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without
More informationGeneral Comment on Lab Reports: v. good + corresponds to a lab report that: has structure (Intro., Method, Results, Discussion, an Abstract would be
General Comment on Lab Reports: v. good + corresponds to a lab report that: has structure (Intro., Method, Results, Discussion, an Abstract would be a bonus) is well written (take your time to edit) shows
More informationCHAPTER 2 DATA AND METHODS. Errors using inadequate data are much less than those using no data at all. Charles Babbage, circa 1850
CHAPTER 2 DATA AND METHODS Errors using inadequate data are much less than those using no data at all. Charles Babbage, circa 185 2.1 Datasets 2.1.1 OLR The primary data used in this study are the outgoing
More informationNOTES AND CORRESPONDENCE. Effect of Sea Surface Temperature Wind Stress Coupling on Baroclinic Instability in the Ocean
1092 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 37 NOTES AND CORRESPONDENCE Effect of Sea Surface Temperature Wind Stress Coupling on Baroclinic Instability in the Ocean MICHAEL A.
More informationDr Marc Lucas CLS Toulouse, France.
Dr Marc Lucas CLS Toulouse, France. Oceanology International 15th March 2012 Slide 1 Ocean depiction In the past: Information mainly comes from in situ measurements (ADCP) Now: The role of satellite data
More informationSIO 210: Dynamics VI (Potential vorticity) L. Talley Fall, 2014 (Section 2: including some derivations) (this lecture was not given in 2015)
SIO 210: Dynamics VI (Potential vorticity) L. Talley Fall, 2014 (Section 2: including some derivations) (this lecture was not given in 2015) Variation of Coriolis with latitude: β Vorticity Potential vorticity
More informationApplications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations
Applications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations Miyazawa, Yasumasa (JAMSTEC) Collaboration with Princeton University AICS Data
More informationModeling of deep currents in the Japan/East Sea
Modeling of deep currents in the Japan/East Sea Olga Trusenkova V.I.Il ichev Pacific Oceanological Institute, FEB RAS Vladivostok, Russia PICES 2014 Annual Meeting, 16-26 October 2014, Korea, Yeosu Deep
More informationCapabilities of Ocean Mixed Layer Models
Capabilities of Ocean Mixed Layer Models W.G. Large National Center for Atmospheric Research Boulder Co, USA 1. Introduction The capabilities expected in today s state of the art models of the ocean s
More informationOn the relationship between the water mass pathways and mesoscale variability in the Western Mediterranean Sea
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, On the relationship between the water mass pathways and mesoscale variability in the Western Mediterranean Sea E.K. Demirov Department of Physics
More informationGlobal observations of large oceanic eddies
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L15606, doi:10.1029/2007gl030812, 2007 Global observations of large oceanic eddies Dudley B. Chelton, 1 Michael G. Schlax, 1 Roger M. Samelson, 1 and Roland A. de
More informationLecture 8. Lecture 1. Wind-driven gyres. Ekman transport and Ekman pumping in a typical ocean basin. VEk
Lecture 8 Lecture 1 Wind-driven gyres Ekman transport and Ekman pumping in a typical ocean basin. VEk wek > 0 VEk wek < 0 VEk 1 8.1 Vorticity and circulation The vorticity of a parcel is a measure of its
More informationPotential use of microwave sea surface temperatures for the estimation of ocean currents - art. no. L24608
Please note that this is an author-produced PDF of an article accepted for publication following peer review. The definitive publisher-authenticated version is available on the publisher Web site GEOPHYSICAL
More informationNote that Rossby waves are tranverse waves, that is the particles move perpendicular to the direction of propagation. f up, down (clockwise)
Ocean 423 Rossby waves 1 Rossby waves: Restoring force is the north-south gradient of background potential vorticity (f/h). That gradient can be due to either the variation in f with latitude, or to a
More informationGFDL, NCEP, & SODA Upper Ocean Assimilation Systems
GFDL, NCEP, & SODA Upper Ocean Assimilation Systems Jim Carton (UMD) With help from Gennady Chepurin, Ben Giese (TAMU), David Behringer (NCEP), Matt Harrison & Tony Rosati (GFDL) Description Goals Products
More informationEVALUATION OF THE GLOBAL OCEAN DATA ASSIMILATION SYSTEM AT NCEP: THE PACIFIC OCEAN
2.3 Eighth Symposium on Integrated Observing and Assimilation Systems for Atmosphere, Oceans, and Land Surface, AMS 84th Annual Meeting, Washington State Convention and Trade Center, Seattle, Washington,
More informationEnhancing predictability of the Loop Current variability using Gulf of Mexico Hycom
Enhancing predictability of the Loop Current variability using Gulf of Mexico Hycom Matthieu Le Hénaff (1) Villy Kourafalou (1) Ashwanth Srinivasan (1) Collaborators: O. M. Smedstad (2), P. Hogan (2),
More informationA note on the numerical representation of surface dynamics in quasigeostrophic turbulence: Application to the nonlinear Eady model
A note on the numerical representation of surface dynamics in quasigeostrophic turbulence: Application to the nonlinear Eady model Ross Tulloch and K. Shafer Smith Center for Atmosphere Ocean Science Courant
More informationWhat governs the location of the Southern Ocean deep winter mixing in CESM
NSF NCAR WYOMING SUPERCOMPUTER CENTER DOE SCIDAC FUNDED PROJECT What governs the location of the Southern Ocean deep winter mixing in CESM Justin Small Dan Whitt Alice DuVivier Matt Long Acknowledging:
More informationWind Gyres. curl[τ s τ b ]. (1) We choose the simple, linear bottom stress law derived by linear Ekman theory with constant κ v, viz.
Wind Gyres Here we derive the simplest (and oldest; Stommel, 1948) theory to explain western boundary currents like the Gulf Stream, and then discuss the relation of the theory to more realistic gyres.
More informationSUPPLEMENTARY INFORMATION
Figure S1. Summary of the climatic responses to the Gulf Stream. On the offshore flank of the SST front (black dashed curve) of the Gulf Stream (green long arrow), surface wind convergence associated with
More informationOCN/ATM/ESS 587. Ocean circulation, dynamics and thermodynamics.
OCN/ATM/ESS 587 Ocean circulation, dynamics and thermodynamics. Equation of state for seawater General T/S properties of the upper ocean Heat balance of the upper ocean Upper ocean circulation Deep circulation
More informationLarge-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling Eric D. Skyllingstad
More informationMonitoring and modeling the Eastern Mediterranean circulation and its climatic variability
UNIVESRITY OF ATHENS, OCEAN PHYSICS AND MODELING GROUP Monitoring and modeling the Eastern Mediterranean circulation and its climatic variability S. Sofianos and OPAM group EastMed Symposium, November
More informationHYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling. Naval Research Laboratory, Stennis Space Center, USA
HYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling By BIROL KARA, ALAN WALLCRAFT AND JOE METZGER Naval Research Laboratory, Stennis Space Center, USA MURAT GUNDUZ Institute
More informationAn Introduction to Coupled Models of the Atmosphere Ocean System
An Introduction to Coupled Models of the Atmosphere Ocean System Jonathon S. Wright jswright@tsinghua.edu.cn Atmosphere Ocean Coupling 1. Important to climate on a wide range of time scales Diurnal to
More informationPacific HYCOM. E. Joseph Metzger, Harley E. Hurlburt, Alan J. Wallcraft, Luis Zamudio and Patrick J. Hogan
Pacific HYCOM E. Joseph Metzger, Harley E. Hurlburt, Alan J. Wallcraft, Luis Zamudio and Patrick J. Hogan Naval Research Laboratory, Stennis Space Center, MS Center for Ocean-Atmospheric Prediction Studies,
More informationBaroclinic Rossby waves in the ocean: normal modes, phase speeds and instability
Baroclinic Rossby waves in the ocean: normal modes, phase speeds and instability J. H. LaCasce, University of Oslo J. Pedlosky, Woods Hole Oceanographic Institution P. E. Isachsen, Norwegian Meteorological
More informationQuasi-geostrophic ocean models
Quasi-geostrophic ocean models March 19, 2002 1 Introduction The starting point for theoretical and numerical study of the three dimensional large-scale circulation of the atmosphere and ocean is a vorticity
More informationRole of Horizontal Density Advection in Seasonal Deepening of the Mixed Layer in the Subtropical Southeast Pacific
ADVANCES IN ATMOSPHERIC SCIENCES, VOL. 33, APRIL 2016, 442 451 Role of Horizontal Density Advection in Seasonal Deepening of the Mixed Layer in the Subtropical Southeast Pacific Qinyu LIU and Yiqun LU
More informationIntegration of ARGO trajectories in the Mediterranean Forecasting System and impact on. the regional analysis of the Western Mediterranean circulation
Integration of ARGO trajectories in the Mediterranean Forecasting System and impact on the regional analysis of the Western Mediterranean circulation by V. Taillandier 1, S. Dobricic 2, P. Testor 3, N.
More informationA note on the numerical representation of surface dynamics in quasigeopstrophic turbulence: Application to the nonlinear Eady model
Submitted to JAS A note on the numerical representation of surface dynamics in quasigeopstrophic turbulence: Application to the nonlinear Eady model Ross Tulloch and K. Shafer Smith Center for Atmosphere
More informationWater mass transport associated with the oceanic fronts in the northwestern Pacific Ocean HIDEYUKI NAKANO (METEOROLOGICAL RESEARCH INSTITUTE)
Water mass transport associated with the oceanic fronts in the northwestern Pacific Ocean HIDEYUKI NAKANO (METEOROLOGICAL RESEARCH INSTITUTE) How is the Kuroshio-origin water distributed in the subtropical
More informationOverview of data assimilation in oceanography or how best to initialize the ocean?
Overview of data assimilation in oceanography or how best to initialize the ocean? T. Janjic Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany Outline Ocean observing system Ocean
More informationAssimilation scheme of the Mediterranean Forecasting System: operational implementation
Annales Geophysicae (2003) 21: 189 204 c European Geosciences Union 2003 Annales Geophysicae Assimilation scheme of the Mediterranean Forecasting System: operational implementation E. Demirov 1, N. Pinardi
More informationModeling the Columbia River Plume on the Oregon Shelf during Summer Upwelling. 2 Model
Modeling the Columbia River Plume on the Oregon Shelf during Summer Upwelling D. P. Fulton August 15, 2007 Abstract The effects of the Columbia River plume on circulation on the Oregon shelf are analyzed
More informationEquatorial Superrotation on Tidally Locked Exoplanets
Equatorial Superrotation on Tidally Locked Exoplanets Adam P. Showman University of Arizona Lorenzo M. Polvani Columbia University Synopsis Most 3D atmospheric circulation models of tidally locked exoplanets
More informationDynamics of the Upper Oceanic Layers in Terms of Surface Quasigeostrophy Theory
VOLUME 36 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y FEBRUARY 26 Dynamics of the Upper Oceanic Layers in Terms of Surface Quasigeostrophy Theory G. LAPEYRE AND P. KLEIN Laboratoire de Physique
More informationLecture 5: Atmospheric General Circulation and Climate
Lecture 5: Atmospheric General Circulation and Climate Geostrophic balance Zonal-mean circulation Transients and eddies Meridional energy transport Moist static energy Angular momentum balance Atmosphere
More informationEddy-induced meridional heat transport in the ocean
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L20601, doi:10.1029/2008gl035490, 2008 Eddy-induced meridional heat transport in the ocean Denis L. Volkov, 1 Tong Lee, 1 and Lee-Lueng Fu 1 Received 28 July 2008;
More informationImprovements to the Wind Driven Component of the OSCAR Surface Current Product
OVWST Meeting May 18-2, 29 p. 1 Improvements to the Wind Driven Component of the OSCAR Surface Current Product Kathleen Dohan, Gary S. E. Lagerloef, and John T. Gunn Earth & Space Research Seattle, Washington
More informationFor the operational forecaster one important precondition for the diagnosis and prediction of
Initiation of Deep Moist Convection at WV-Boundaries Vienna, Austria For the operational forecaster one important precondition for the diagnosis and prediction of convective activity is the availability
More informationCAIBEX workshop Mesoscale experiments and modelling Cape Ghir
CAIBEX workshop Mesoscale experiments and modelling Cape Ghir C. Troupin 1, P. Sangrà 2, J. Arístegui 2 1 GHER-MARE, AGO Department, University of Liège, Belgium 2 Facultad de Ciencias del Mar, Universidad
More informationO.M Smedstad 1, E.J. Metzger 2, R.A. Allard 2, R. Broome 1, D.S. Franklin 1 and A.J. Wallcraft 2. QinetiQ North America 2. Naval Research Laboratory
An eddy-resolving ocean reanalysis using the 1/12 global HYbrid Coordinate Ocean Model (HYCOM) and the Navy Coupled Ocean Data Assimilation (NCODA) scheme O.M Smedstad 1, E.J. Metzger 2, R.A. Allard 2,
More informationDecadal variability in the Kuroshio and Oyashio Extension frontal regions in an eddy-resolving OGCM
Decadal variability in the Kuroshio and Oyashio Extension frontal regions in an eddy-resolving OGCM Masami Nonaka 1, Hisashi Nakamura 1,2, Youichi Tanimoto 1,3, Takashi Kagimoto 1, and Hideharu Sasaki
More informationTHE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS
THE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS Brian J. Soden 1 and Christopher S. Velden 2 1) Geophysical Fluid Dynamics Laboratory National Oceanic and Atmospheric Administration
More informationCoastal Ocean Circulation Experiment off Senegal (COCES)
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Coastal Ocean Circulation Experiment off Senegal (COCES) Pierre-Marie Poulain Istituto Nazionale di Oceanografia e di Geofisica
More informationPUBLICATIONS. Journal of Geophysical Research: Oceans
PUBLICATIONS Journal of Geophysical Research: Oceans RESEARCH ARTICLE Key Points: Altimeter data errors Spectral analysis Dynamics and modeling Atlantic sea surface height and velocity spectra inferred
More informationDynamics and Kinematics
Geophysics Fluid Dynamics () Syllabus Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3,
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationBay of Bengal Surface and Thermocline and the Arabian Sea
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Bay of Bengal Surface and Thermocline and the Arabian Sea Arnold L. Gordon Lamont-Doherty Earth Observatory of Columbia
More informationOn the Wind Power Input and Eddy Residence Time
Hamburg Workshop 215 On the Wind Power Input and Eddy Residence Time Xiaoming Zhai Centre for Ocean and Atmospheric Sciences School of Environmental Sciences, University of East Anglia With David Marshall,
More informationGeophysics Fluid Dynamics (ESS228)
Geophysics Fluid Dynamics (ESS228) Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3, 4,
More informationFinescale Water-Mass Variability from ARGO Profiling Floats
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Finescale Water-Mass Variability from ARGO Profiling Floats Eric Kunze Applied Physics Lab, University of Washington 1013
More informationJohn Steffen and Mark A. Bourassa
John Steffen and Mark A. Bourassa Funding by NASA Climate Data Records and NASA Ocean Vector Winds Science Team Florida State University Changes in surface winds due to SST gradients are poorly modeled
More informationIsland Wakes in Shallow Water
Island Wakes in Shallow Water Changming Dong, James C. McWilliams, et al Institute of Geophysics and Planetary Physics, University of California, Los Angeles 1 ABSTRACT As a follow-up work of Dong et al
More informationPrecipitation processes in the Middle East
Precipitation processes in the Middle East J. Evans a, R. Smith a and R.Oglesby b a Dept. Geology & Geophysics, Yale University, Connecticut, USA. b Global Hydrology and Climate Center, NASA, Alabama,
More information8 Mechanisms for tropical rainfall responses to equatorial
8 Mechanisms for tropical rainfall responses to equatorial heating More reading: 1. Hamouda, M. and Kucharski, F. (2019) Ekman pumping Mechanism driving Precipitation anomalies in Response to Equatorial
More informationThe role of the midlatitude ocean in sub-seasonal prediction
The role of the midlatitude ocean in sub-seasonal prediction R. Saravanan Xiaohui Ma, Xue Liu, J. Steinweg-Woods J. Kurian, R. Montuoro, P. Chang, I. Szunyogh Yinglai Jia, Ocean University of China J.
More informationPosition variability of the Kuroshio Extension sea surface temperature front
Acta Oceanol. Sin., 2016, Vol. 35, No. 7, P. 30 35 DOI: 10.1007/s13131-016-0909-7 http://www.hyxb.org.cn E-mail: hyxbe@263.net Position variability of the Kuroshio Extension sea surface temperature front
More informationPUBLICATIONS. Journal of Geophysical Research: Oceans. Retrieving density and velocity fields of the ocean s interior.
PUBLICATIONS Journal of Geophysical Research: Oceans RESEARCH ARTICLE Key Points: The ocean s interior can be retrieved from surface data The isqg method could serve as a simplified dynamic framework The
More informationCHAPTER 2 - ATMOSPHERIC CIRCULATION & AIR/SEA INTERACTION
Chapter 2 - pg. 1 CHAPTER 2 - ATMOSPHERIC CIRCULATION & AIR/SEA INTERACTION The atmosphere is driven by the variations of solar heating with latitude. The heat is transferred to the air by direct absorption
More informationTrends in Climate Teleconnections and Effects on the Midwest
Trends in Climate Teleconnections and Effects on the Midwest Don Wuebbles Zachary Zobel Department of Atmospheric Sciences University of Illinois, Urbana November 11, 2015 Date Name of Meeting 1 Arctic
More informationApplying Basin-Scale HyCOM Hindcasts in Providing Open Boundary Conditions for Nested High-Resolution Coastal Circulation Modeling
Applying Basin-Scale HyCOM Hindcasts in Providing Open Boundary Conditions for Nested High-Resolution Coastal Circulation Modeling Ruoying He Woods Hole Oceanographic Institution December 7, 2005 Cape
More information