Search Results
You are looking at 11  20 of 23 items for :
 Author or Editor: Paola Malanotte Rizzoli x
 Article x
 Refine by Access: Content accessible to me x
Abstract
One of the most important forthcoming synoptic datasets for ocean circulation studies will he the seasurface height data provided by the TOPEX /POSEIDON satellite. The TOPEX/POSEIDON project is in the planning stage and must still decide upon the particular characteristics of the satellite track. The repeat period will be between 10 and 20 days for a variety of technical and strategic reasons. These choices win give a global coverage with spatial resolution (eastwest or northsouth separation of crossover points) in midlatitudes of roughly 2.8° of latitude and longitude for a 10day repeat orbit and 1.4° of latitude and longitude for a 20day repeat orbit. Thus, the crucial question we address in the present study is: what is the effect of changing space or time resolution or both upon the success of a numerical model in reconstructing a fourdimensional picture a the ocean circulation through the assimilation of altimetric data?
To answer this question we carry out a series of numerical experiments with a threelayer, eddyresolving quasigeostrophic model of the ocean circulation in which we systematically vary the space and time resolutions of the data available for assimilation experiments. The experiments are carded out under the “best possible” conditions for the assimilation to be successful, namely: (i) the model is “perfect"; (ii) the data have no errors; and (iii) the data are dynamically compatible with the model since they are simulated by the model itself in a control run.
We reach the following conclusions. In principle, assimilation of altimetric data with a simple relaxation (“nudging”) technique can be very successful in driving the assimilation model to the control run even in the deep layers for which no data are supplied. This is achieved with a “nearly perfect” spacetime resolution surface height dataset in which data are supplied at every model grid point and every 0.5 day in time. The residual errors after one year of continuous assimilation amount to less than 10% in all three layers. When the altimetric data are provided along tracks with a given realistic separation (but complete time information), the decrease in space resolution degrades the model estimates somewhat. With data provided at every time step but a track separation of 280 km and making use of the best choice of assimilation procedures we have found that the residual rms errors amount to about 45% after six years of continuous assimilation. While the patterns of the circulation are somewhat different from those of the control run and the flow intensifies are slightly underestimated, the correspondences between the assimilation run and the control run are considerable. When the altimetric data are provided with a realistic time sampling period (but with space resolution at every grid point), the intensity of the flow fields also are somewhat underestimated, especially in the deep layer. The assimilation procedure is again capable, however, of reproducing quite faithfully the flow patterns throughout the water column.
When the altimetric data are assimilated along the actual tracks, that is only at the track gridpoints and at the actual time of arrival, the best assimilation results achieved with TOPEX repeat periods of 10 or 20 days are about equally effective for improving the model estimates of the circulation. The residual errors after 6 years of continuous time assimilation are from 60% to 70% for both 10 and 20day repeats. Apparently, the tradeoff between space and time resolutions just about compensate for each other. The results show that under the best of conditions (small errors, good model) a single satellite makes only minor improvements in the model estimates, and it cannot reconstruct the details of the mesoscale eddy field.
It should be kept in mind that these results depend on the space and time scales of motion in the region to be studied. Moreover, the conclusions reached here depend, to an unknown extent, on the assimilation technique used. Better techniques might allow us to better differentiate between the different spacetime choices for TOPEX and to reproduce the actual oceanic circulation more faithfully.
Abstract
One of the most important forthcoming synoptic datasets for ocean circulation studies will he the seasurface height data provided by the TOPEX /POSEIDON satellite. The TOPEX/POSEIDON project is in the planning stage and must still decide upon the particular characteristics of the satellite track. The repeat period will be between 10 and 20 days for a variety of technical and strategic reasons. These choices win give a global coverage with spatial resolution (eastwest or northsouth separation of crossover points) in midlatitudes of roughly 2.8° of latitude and longitude for a 10day repeat orbit and 1.4° of latitude and longitude for a 20day repeat orbit. Thus, the crucial question we address in the present study is: what is the effect of changing space or time resolution or both upon the success of a numerical model in reconstructing a fourdimensional picture a the ocean circulation through the assimilation of altimetric data?
To answer this question we carry out a series of numerical experiments with a threelayer, eddyresolving quasigeostrophic model of the ocean circulation in which we systematically vary the space and time resolutions of the data available for assimilation experiments. The experiments are carded out under the “best possible” conditions for the assimilation to be successful, namely: (i) the model is “perfect"; (ii) the data have no errors; and (iii) the data are dynamically compatible with the model since they are simulated by the model itself in a control run.
We reach the following conclusions. In principle, assimilation of altimetric data with a simple relaxation (“nudging”) technique can be very successful in driving the assimilation model to the control run even in the deep layers for which no data are supplied. This is achieved with a “nearly perfect” spacetime resolution surface height dataset in which data are supplied at every model grid point and every 0.5 day in time. The residual errors after one year of continuous assimilation amount to less than 10% in all three layers. When the altimetric data are provided along tracks with a given realistic separation (but complete time information), the decrease in space resolution degrades the model estimates somewhat. With data provided at every time step but a track separation of 280 km and making use of the best choice of assimilation procedures we have found that the residual rms errors amount to about 45% after six years of continuous assimilation. While the patterns of the circulation are somewhat different from those of the control run and the flow intensifies are slightly underestimated, the correspondences between the assimilation run and the control run are considerable. When the altimetric data are provided with a realistic time sampling period (but with space resolution at every grid point), the intensity of the flow fields also are somewhat underestimated, especially in the deep layer. The assimilation procedure is again capable, however, of reproducing quite faithfully the flow patterns throughout the water column.
When the altimetric data are assimilated along the actual tracks, that is only at the track gridpoints and at the actual time of arrival, the best assimilation results achieved with TOPEX repeat periods of 10 or 20 days are about equally effective for improving the model estimates of the circulation. The residual errors after 6 years of continuous time assimilation are from 60% to 70% for both 10 and 20day repeats. Apparently, the tradeoff between space and time resolutions just about compensate for each other. The results show that under the best of conditions (small errors, good model) a single satellite makes only minor improvements in the model estimates, and it cannot reconstruct the details of the mesoscale eddy field.
It should be kept in mind that these results depend on the space and time scales of motion in the region to be studied. Moreover, the conclusions reached here depend, to an unknown extent, on the assimilation technique used. Better techniques might allow us to better differentiate between the different spacetime choices for TOPEX and to reproduce the actual oceanic circulation more faithfully.
Abstract
Many coastal regions in the world ocean are characterized by wellmixed conditions to shelf depth in the density field during the winter season. In these situations it is appropriate to construct a model based on the assumption that the shelf is vertically well mixed. Such a model has been constructed assuming that (i) vertical mixing of momentum is stronger than either horizontal mixing or inertial effects; (ii) the density field is also vertically well mixed, i.e., varies only in horizontal at zero order in the expansion in the vertical Peclet number, (iii) the crossshelf scale is small compared to the alongshelf scale; (iv) depth varies only in crossshelf direction. The transport streamfunction equation and the advective density equation can then be combined into a single model equation by noting that in the vertically wellmixed flow, density is conserved as it is advected along streamlines.
This model is used to study two different configurations quite common for shelf circulations in the world ocean. The first configuration considers the effect of a deep baroclinic ocean in driving the shelf circulation. Past studies show that a barotropic deepocean pressure gradient cannot drive significant shelf flow and that the continental slope effectively “insulates”the shelf from the deep ocean. Thus, the basic question we want to answer is: how does the baroclinic structure of a deep ocean flow affect its ability to penetrate the shelf and determine its circulation? The second configuration to which the model is applied is that of a coastal current driven by an alongshore buoyancy source, such as a river discharge or an alongshore jet.
In all model applications studied, the basic mechanism by which the flow is able to cross topography is bottom friction. In the first situation the flow is driven by prescribing the velocity and density at the outer edge of the shelf. Three types of shelf forcing by the deep ocean are studied: a wide inflow, a narrow inflow and forcing by a Gulf Stream ring. A specific application is made to the northern Adriatic coastal shelf forced by a dense water pool formed in wintertime in the Adriatic interior. In all cases, the main conclusion for the baroclinic deepocean inflow onto the wellmixed shelf is that, like in the barotropic case, the tendency for the flow to follow isobaths is much stronger than the degree to which bottom friction allows crossisobath motion. The deep ocean inflow forces a horizontal boundary layer against the shelf edge. The width of this boundary layer, and therefore the shelf penetration, is larger if the drag coefficient is higher, the latitude lower or the bottom slope weaker. Also, surfaceintensified deep ocean flows penetrate the shelf few strongly than bottom intensified flows.
In the two examples studied of flow entering the shelf region near shore, i.e. a coastal river outflow and a coastal jet, it is again shown that the vertical shear of the shore inflow again determines the overall flow pattern. Because of bottom friction, a light alongshore jet or a low density river, i.e., surface intensified flow, expands across the topography more slowly than a bottom intensified flow, such as a high density river or a heavy jet.
Abstract
Many coastal regions in the world ocean are characterized by wellmixed conditions to shelf depth in the density field during the winter season. In these situations it is appropriate to construct a model based on the assumption that the shelf is vertically well mixed. Such a model has been constructed assuming that (i) vertical mixing of momentum is stronger than either horizontal mixing or inertial effects; (ii) the density field is also vertically well mixed, i.e., varies only in horizontal at zero order in the expansion in the vertical Peclet number, (iii) the crossshelf scale is small compared to the alongshelf scale; (iv) depth varies only in crossshelf direction. The transport streamfunction equation and the advective density equation can then be combined into a single model equation by noting that in the vertically wellmixed flow, density is conserved as it is advected along streamlines.
This model is used to study two different configurations quite common for shelf circulations in the world ocean. The first configuration considers the effect of a deep baroclinic ocean in driving the shelf circulation. Past studies show that a barotropic deepocean pressure gradient cannot drive significant shelf flow and that the continental slope effectively “insulates”the shelf from the deep ocean. Thus, the basic question we want to answer is: how does the baroclinic structure of a deep ocean flow affect its ability to penetrate the shelf and determine its circulation? The second configuration to which the model is applied is that of a coastal current driven by an alongshore buoyancy source, such as a river discharge or an alongshore jet.
In all model applications studied, the basic mechanism by which the flow is able to cross topography is bottom friction. In the first situation the flow is driven by prescribing the velocity and density at the outer edge of the shelf. Three types of shelf forcing by the deep ocean are studied: a wide inflow, a narrow inflow and forcing by a Gulf Stream ring. A specific application is made to the northern Adriatic coastal shelf forced by a dense water pool formed in wintertime in the Adriatic interior. In all cases, the main conclusion for the baroclinic deepocean inflow onto the wellmixed shelf is that, like in the barotropic case, the tendency for the flow to follow isobaths is much stronger than the degree to which bottom friction allows crossisobath motion. The deep ocean inflow forces a horizontal boundary layer against the shelf edge. The width of this boundary layer, and therefore the shelf penetration, is larger if the drag coefficient is higher, the latitude lower or the bottom slope weaker. Also, surfaceintensified deep ocean flows penetrate the shelf few strongly than bottom intensified flows.
In the two examples studied of flow entering the shelf region near shore, i.e. a coastal river outflow and a coastal jet, it is again shown that the vertical shear of the shore inflow again determines the overall flow pattern. Because of bottom friction, a light alongshore jet or a low density river, i.e., surface intensified flow, expands across the topography more slowly than a bottom intensified flow, such as a high density river or a heavy jet.
Abstract
In this work we take a first step in the process of assimilating data into models of the ocean general circulation. The goals is not prediction but rather understanding how the data insertion process affects, and is affected by, the dynamics governing the model. The chosen model ocean is steady, weakly nonlinear and highly frictional Strong vertical friction plays the role of eddy fluxes in driving the circulation in the deep layers.
In the data insertion process we capitalize upon the two principles that (i) the available dynamical models are imperfect; (ii) oceanographic data are measured locally. Three major questions are addressed; 1) what is the influence of local data insertion in terms of improving estimates of the model general circulation? 2) how does the model dynamics affect the spreading of information from the data insertion region? 3) what can we learn about the model physics from the effects of data insertion
Density (or temperature) measurements along long hydrographic or tomographic sections or arrays are chosen as data. We vary the location of the section as well as its orientation. In our highly frictional ocean, the most effective sections are meridional, long and located at a distance from the western boundary. Model estimates are then significantly improved over the broad region extending from the data section to the western boundary itself.
Advective effects are minimal and influence the spreading of information only in the intense western boundary current. Rather, the structure of the gyre interior manifests itself through a quite important steering effect exerted by the motion in the intermediate layer upon the spread of information in the surface layer. Due to this effect the region southwest of the data section is consistently preferred for the improvement of the estimates. Simple analytical computations are carried out to rationalize the numerical results. This effect is likely to persist in more realistic, fully eddyresolving simulations in which the interfacial eddy stresses would play the role here given to vertical friction.
The dependence of spreading of information upon the internal physics and/or external forcing is used to examine what is imperfect in the model parameterizations. In a simple analytical example we scan the twodimensional parameter space defined by internal friction and wind stress amplitude. The “correct” values of the above parameters cannot be inferred by this simple scanning due to the nonuniqueness of the solution.
Abstract
In this work we take a first step in the process of assimilating data into models of the ocean general circulation. The goals is not prediction but rather understanding how the data insertion process affects, and is affected by, the dynamics governing the model. The chosen model ocean is steady, weakly nonlinear and highly frictional Strong vertical friction plays the role of eddy fluxes in driving the circulation in the deep layers.
In the data insertion process we capitalize upon the two principles that (i) the available dynamical models are imperfect; (ii) oceanographic data are measured locally. Three major questions are addressed; 1) what is the influence of local data insertion in terms of improving estimates of the model general circulation? 2) how does the model dynamics affect the spreading of information from the data insertion region? 3) what can we learn about the model physics from the effects of data insertion
Density (or temperature) measurements along long hydrographic or tomographic sections or arrays are chosen as data. We vary the location of the section as well as its orientation. In our highly frictional ocean, the most effective sections are meridional, long and located at a distance from the western boundary. Model estimates are then significantly improved over the broad region extending from the data section to the western boundary itself.
Advective effects are minimal and influence the spreading of information only in the intense western boundary current. Rather, the structure of the gyre interior manifests itself through a quite important steering effect exerted by the motion in the intermediate layer upon the spread of information in the surface layer. Due to this effect the region southwest of the data section is consistently preferred for the improvement of the estimates. Simple analytical computations are carried out to rationalize the numerical results. This effect is likely to persist in more realistic, fully eddyresolving simulations in which the interfacial eddy stresses would play the role here given to vertical friction.
The dependence of spreading of information upon the internal physics and/or external forcing is used to examine what is imperfect in the model parameterizations. In a simple analytical example we scan the twodimensional parameter space defined by internal friction and wind stress amplitude. The “correct” values of the above parameters cannot be inferred by this simple scanning due to the nonuniqueness of the solution.
Abstract
The major objective of oceanic data assimilation studies has been thus far to obtain a fourdimensional realization (space plus time) of the oceanic flow simultaneously consistent with the observations and the model dynamics. In these latest years, however, the forecasting of oceanic motions has emerged as a legitimate and important goal per se. In particular, the operational prediction of mesoscale flows and frontal systems has been the objective of recent assimilation applications in various regional systems of the World Ocean. One such effort focused on the shortterm prediction of the Gulf Stream system in the DAMEE GSR (Data Assimilation and Model Evaluation Experiments Gulf Stream Region) sponsored by the U.S. Navy. The objective of DAMEE GSR phases I and II was 1–2week forecast experiments. Phase III extended the suite of case studies by adding a 2monthlong assimilation experiment to assess the impact of longterm assimilations on model performance and forecasting skill.
In this paper the authors report the results of DAMEE GSR phase III but broaden the perspective by addressing two further issues, namely, the model sensitivity to the choice of the initial fields and the frequency of intermittent data assimilation. Two versions of the OTIS3 (Optimum Thermal Interpolation System) of the U.S. Navy Fleet Numerical Oceanography Center were available, providing slightly different distributions of temperature and salinity over the entire Gulf Stream system. They are referred to as OTIS3a, available with biweekly frequency from 4 May 1988 through 28 December 1988, in the context of a different assimilation work; and OTIS3b, provided by the DAMEE GSR phase III effort, for the 2month period 4 May–4 July 1988, with a slightly irregular frequency, weekly on the average. The main results can be summarized as follows.
The intermittent assimilation of the OTIS3b datasets with average weekly frequency profoundly improves the model forecasting skill. Without assimilation the model never beats persistence. With the assimilation, the modelpredicted Gulf Stream north wall is in excellent agreement with the verification infrared (IR) north wall, remaining always within the error bar of the IR north wall estimate, ±15km.
Two types of sensitivity experiments to the initial conditions were carried out: first, reconstruction of the initial fields with the two different OTIS3a and OTIS3b datasets but with the same initialization method; second, reconstruction of the initial fields with the same OTIS3a dataset but with two different initialization methods. The results show that the initial velocity field is much more crucial in affecting the model evolution and hence its predictive skill as it determines the stability properties of the Gulf Stream jet. Hence, it is very important to use the same dynamical initialization for velocity when starting from different distributions of temperature and salinity, as the jet profiles thus obtained will be very similar in structure and strength. This identical dynamical initialization will allow meaningful comparisons of experiments that start from slightly different density distributions.
Finally, we compare weekly assimilations of OTIS3b with biweekly and monthly assimilations of OTIS3a, initialized with the same procedure. The authors conclude that a weekly assimilation of the global OTIS3 dataset is not necessary and that a biweekly assimilation is equally effective in improving the model predictive skill.
Abstract
The major objective of oceanic data assimilation studies has been thus far to obtain a fourdimensional realization (space plus time) of the oceanic flow simultaneously consistent with the observations and the model dynamics. In these latest years, however, the forecasting of oceanic motions has emerged as a legitimate and important goal per se. In particular, the operational prediction of mesoscale flows and frontal systems has been the objective of recent assimilation applications in various regional systems of the World Ocean. One such effort focused on the shortterm prediction of the Gulf Stream system in the DAMEE GSR (Data Assimilation and Model Evaluation Experiments Gulf Stream Region) sponsored by the U.S. Navy. The objective of DAMEE GSR phases I and II was 1–2week forecast experiments. Phase III extended the suite of case studies by adding a 2monthlong assimilation experiment to assess the impact of longterm assimilations on model performance and forecasting skill.
In this paper the authors report the results of DAMEE GSR phase III but broaden the perspective by addressing two further issues, namely, the model sensitivity to the choice of the initial fields and the frequency of intermittent data assimilation. Two versions of the OTIS3 (Optimum Thermal Interpolation System) of the U.S. Navy Fleet Numerical Oceanography Center were available, providing slightly different distributions of temperature and salinity over the entire Gulf Stream system. They are referred to as OTIS3a, available with biweekly frequency from 4 May 1988 through 28 December 1988, in the context of a different assimilation work; and OTIS3b, provided by the DAMEE GSR phase III effort, for the 2month period 4 May–4 July 1988, with a slightly irregular frequency, weekly on the average. The main results can be summarized as follows.
The intermittent assimilation of the OTIS3b datasets with average weekly frequency profoundly improves the model forecasting skill. Without assimilation the model never beats persistence. With the assimilation, the modelpredicted Gulf Stream north wall is in excellent agreement with the verification infrared (IR) north wall, remaining always within the error bar of the IR north wall estimate, ±15km.
Two types of sensitivity experiments to the initial conditions were carried out: first, reconstruction of the initial fields with the two different OTIS3a and OTIS3b datasets but with the same initialization method; second, reconstruction of the initial fields with the same OTIS3a dataset but with two different initialization methods. The results show that the initial velocity field is much more crucial in affecting the model evolution and hence its predictive skill as it determines the stability properties of the Gulf Stream jet. Hence, it is very important to use the same dynamical initialization for velocity when starting from different distributions of temperature and salinity, as the jet profiles thus obtained will be very similar in structure and strength. This identical dynamical initialization will allow meaningful comparisons of experiments that start from slightly different density distributions.
Finally, we compare weekly assimilations of OTIS3b with biweekly and monthly assimilations of OTIS3a, initialized with the same procedure. The authors conclude that a weekly assimilation of the global OTIS3 dataset is not necessary and that a biweekly assimilation is equally effective in improving the model predictive skill.
Abstract
Ocean Acoustic Tomography was proposed by Munk and Wunsch as a method for making measurements of ocean variability over large areas. After the successful demonstration of the feasibility of the idea in the 1981 threedimensional Mesoscale Experiment the tomography group has proposed a new experiment to be carried out in the Eastern Pacific Ocean, on ranges as long as the subtropical gyre scale.
This paper address the question of which average quantities of importance for the ocean general circulation and ocean climate can be measured by tomography and with what accuracy. The paper focuses upon the following quantities i) measurement of the heat content vertical profile horizontally averaged over a tomographic section; ii) time variability of the average heat content, or average pycnocline displacement, at different depths; iii) measurement of the average pycnocline slope at different depths.
To answer the above question the tomographic experiment is simulated in a given model ocean, using Holland's eddyresolving general circulation quasigeostrophic model. The results of the modeling simulations can be summarized as follows.

The tomographic technique bars upon the use of inverse methods to reconstruct the interior sound speed perturbation field, or, equivalently, the heat content field. Over ranges as long as the gyre scale, the typical result of a single inversion is to provide an ocean with warm or cold biases. A simple iterative procedure allows the removal of these biases. The final estimates of the mean heat content (averaged over the tomographic section) at different depths is very good.

Through a timeevolution experiment carried out for the duration of a full year, the time evolution of the average pycnocline displacement can be monitored at various depths. Thus tomography can measure the frequency spectrum of the average pycnocline displacement in layers below the surface mixed layer in which the circulation is basically winddriven.

The initial estimate of the average heat content can be significantly improved through a better specification of the statistics of the region, like the inclusion of a spatial mean in the horizontal covariance function for the sound speed perturbation. In particular, the inclusion of an inhomogeneous covariance modeling longscale pycnocline trends allow us to estimate the average pycnocline slope at various depths. The obtained slope estimates are very good. Measurement of isopycnal slopes averaged in time could be used for βspiral calculations. Thus, simple “density” tomography would provide a tool to evaluate the absolute velocity field and not only the geostrophic velocity shear.
Abstract
Ocean Acoustic Tomography was proposed by Munk and Wunsch as a method for making measurements of ocean variability over large areas. After the successful demonstration of the feasibility of the idea in the 1981 threedimensional Mesoscale Experiment the tomography group has proposed a new experiment to be carried out in the Eastern Pacific Ocean, on ranges as long as the subtropical gyre scale.
This paper address the question of which average quantities of importance for the ocean general circulation and ocean climate can be measured by tomography and with what accuracy. The paper focuses upon the following quantities i) measurement of the heat content vertical profile horizontally averaged over a tomographic section; ii) time variability of the average heat content, or average pycnocline displacement, at different depths; iii) measurement of the average pycnocline slope at different depths.
To answer the above question the tomographic experiment is simulated in a given model ocean, using Holland's eddyresolving general circulation quasigeostrophic model. The results of the modeling simulations can be summarized as follows.

The tomographic technique bars upon the use of inverse methods to reconstruct the interior sound speed perturbation field, or, equivalently, the heat content field. Over ranges as long as the gyre scale, the typical result of a single inversion is to provide an ocean with warm or cold biases. A simple iterative procedure allows the removal of these biases. The final estimates of the mean heat content (averaged over the tomographic section) at different depths is very good.

Through a timeevolution experiment carried out for the duration of a full year, the time evolution of the average pycnocline displacement can be monitored at various depths. Thus tomography can measure the frequency spectrum of the average pycnocline displacement in layers below the surface mixed layer in which the circulation is basically winddriven.

The initial estimate of the average heat content can be significantly improved through a better specification of the statistics of the region, like the inclusion of a spatial mean in the horizontal covariance function for the sound speed perturbation. In particular, the inclusion of an inhomogeneous covariance modeling longscale pycnocline trends allow us to estimate the average pycnocline slope at various depths. The obtained slope estimates are very good. Measurement of isopycnal slopes averaged in time could be used for βspiral calculations. Thus, simple “density” tomography would provide a tool to evaluate the absolute velocity field and not only the geostrophic velocity shear.
Abstract
Increasing evidence suggests that the Gulf Stream system is the origin for much of the eddy variability observed in the northwestern Atlantic Ocean. However, the dynamical mechanisms by which eddy energy, once generated in the neighborhood of the Gulf Stream, penetrates into the midocean regions are not fully understood. Here we explore the proposition that radiation away from the Gulf Stream can be interpreted as an oscillatory response of the ocean gyre interior to the forcing associated with translation and temporal evolution of the Gulf Stream current and its meanders. Specifically, we propose as a possible mechanism responsible for pules of Rossby wave radiation the sudden growth to finite amplitude and successive amplitude pulsations of quasistationary or eastwardmoving finiteamplitude meanders. To examine this mechanism in detail, we construct a model for radiation in which the quasigeostrophic equivalent barotropic potential vorticity equation is forced to its northern boundary by assigned distributions of streamfunction and vorticity corresponding to 1) a stationary pulsating meander, 2) a slowly propagating meander, and 3) a combination of 1 and 2.
We use both analytical and numerical techniques to solve the linear boundaryforced initial value problem with and without lateral friction and with and without bottom topography. We focus not only on the steady forced response but also on the transient component of the interior Rossby radiation field. On the βplane without topography, it is shown that radiation into the far field is only possible for meander pulsation below a critical frequency. Below this frequency, the transient farfield response is large and force basin modes are necessary, together with the forced equilibrium solution, to accommodate the prescribed initial state of the system. With the addition of viscosity to damp the transient basin modes, and for frequencies higher than the critical value, the response is trapped at the northern wall and no radiation can be supported in the far field.
The inclusion of a simple parabolic relief along the full channel does not modify the qualitative behavior of the interior field. The only effect of the parabolic relief is to modify the effective value of planetary β and the crosschannel structure of the westwardpropagating Rossby modes. The above results are, however, profoundly altered with the addition of nonlinearity. In particular, new mechanisms are possible for the development of farfield radiation and for the appearance of nonlinear, coherent solutions, as will he shown in Part II of this paper.
Abstract
Increasing evidence suggests that the Gulf Stream system is the origin for much of the eddy variability observed in the northwestern Atlantic Ocean. However, the dynamical mechanisms by which eddy energy, once generated in the neighborhood of the Gulf Stream, penetrates into the midocean regions are not fully understood. Here we explore the proposition that radiation away from the Gulf Stream can be interpreted as an oscillatory response of the ocean gyre interior to the forcing associated with translation and temporal evolution of the Gulf Stream current and its meanders. Specifically, we propose as a possible mechanism responsible for pules of Rossby wave radiation the sudden growth to finite amplitude and successive amplitude pulsations of quasistationary or eastwardmoving finiteamplitude meanders. To examine this mechanism in detail, we construct a model for radiation in which the quasigeostrophic equivalent barotropic potential vorticity equation is forced to its northern boundary by assigned distributions of streamfunction and vorticity corresponding to 1) a stationary pulsating meander, 2) a slowly propagating meander, and 3) a combination of 1 and 2.
We use both analytical and numerical techniques to solve the linear boundaryforced initial value problem with and without lateral friction and with and without bottom topography. We focus not only on the steady forced response but also on the transient component of the interior Rossby radiation field. On the βplane without topography, it is shown that radiation into the far field is only possible for meander pulsation below a critical frequency. Below this frequency, the transient farfield response is large and force basin modes are necessary, together with the forced equilibrium solution, to accommodate the prescribed initial state of the system. With the addition of viscosity to damp the transient basin modes, and for frequencies higher than the critical value, the response is trapped at the northern wall and no radiation can be supported in the far field.
The inclusion of a simple parabolic relief along the full channel does not modify the qualitative behavior of the interior field. The only effect of the parabolic relief is to modify the effective value of planetary β and the crosschannel structure of the westwardpropagating Rossby modes. The above results are, however, profoundly altered with the addition of nonlinearity. In particular, new mechanisms are possible for the development of farfield radiation and for the appearance of nonlinear, coherent solutions, as will he shown in Part II of this paper.
Abstract
The improvement in the climatological behavior of a numerical model as a consequence of the assimilation of surface data is investigated. The model used for this study is a quasigeostrophic (QG) model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height that have been obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations.
Comparisons of the assimilation results with available in situ observations show a significant improvement in the degree of realism of the climatological model behavior, with respect to the model in which no data are assimilated. The remaining discrepancies in the model mean circulation seem to be mainly associated with deficiencies in the mean component of the surface data that are assimilated. On the other hand, the possibility of building into the model more realistic eddy characteristics through the assimilation of the surface eddy field proves very successful in driving components of the mean model circulation that are in relatively good agreement with the available observations. Comparisons with current meter time series during a time period partially overlapping the Geosat mission show that the model is able to “correctly” extrapolate the instantaneous surface eddy signals to depths of approximately 1500 m. The correlation coefficient between current meter and model time series varies from values close to 0.7 in the top 1500 m to values as low as 0.1–0.2 in the deep ocean.
Abstract
The improvement in the climatological behavior of a numerical model as a consequence of the assimilation of surface data is investigated. The model used for this study is a quasigeostrophic (QG) model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height that have been obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations.
Comparisons of the assimilation results with available in situ observations show a significant improvement in the degree of realism of the climatological model behavior, with respect to the model in which no data are assimilated. The remaining discrepancies in the model mean circulation seem to be mainly associated with deficiencies in the mean component of the surface data that are assimilated. On the other hand, the possibility of building into the model more realistic eddy characteristics through the assimilation of the surface eddy field proves very successful in driving components of the mean model circulation that are in relatively good agreement with the available observations. Comparisons with current meter time series during a time period partially overlapping the Geosat mission show that the model is able to “correctly” extrapolate the instantaneous surface eddy signals to depths of approximately 1500 m. The correlation coefficient between current meter and model time series varies from values close to 0.7 in the top 1500 m to values as low as 0.1–0.2 in the deep ocean.
Abstract
In Part I of the present work we studied the transient Rossby wave radiation excited in the far field by a northern boundary forcing. We proposed as a possible mechanism responsible for transient pulses of Rossby waves the sudden growth to finite amplitude and successive amplitude pulsations of a stationary or eastward moving meander of an intense current like the Gulf Stream. In Part I the linear regime was thoroughly analyzed and the major findings were: 1) the transient part of the response initially excited by “switching on” the meander is the one responsible for the Rossby waves radiating and being clearly detectable in the far field; 2) the Rossby radiation is excited only when the meander pulsation frequency is below a critical value.
In this paper, we extend the results of Part I to the highly nonlinear regime addressing the problem of the production of nonlinear, coherent structures through the same boundary forcing mechanisms. The model used is the quasigeostrophic, potential vorticity equation on a βplane with or without topography, in a periodic channel forced by a northern boundary function. The forcing functions are designed to simulate (i) a stationary pulsating meander, (ii) a slowly propagating meander, (iii) a combination of both. This may constitute the simplest possible idealization of meander growth to finite amplitude, its successive steepening and bending with pinching off of eddies as frequently observed in the Gulf Stream system.
Using a boundary forcing idealizing a steady meander with a pulsating amplitude, we excite nonlinear vortices on the βplane without topography. The related regime is called the “modon” regime as, at every pulsation cycle, two circulation patterns are produced, one cyclonic, the other anticyclonic. These vortices are part of a dipole pair and selfadvect eastward or westward depending on their polarity (lowabovehigh or inverse). Through diagnostic tools we show that these vortices are ‘reasonable’ modons, that is highly nonlinear, exact solutions of the equivalent barotropic model on the βplane.
The second type of boundary forcing capitalizes upon a resonance mechanism which was previously studied analytically in the weakly nonlinear regime. In it, an interior vortex is resonantly excited by a meander with constant amplitude which propagates at one of the eigenspeeds of the free linear modes of the channel. This mechanism needs topography to be effective when using quasigeostrophic dynamics. With parabolic relief all along the channel, a resonantly forced nonlinear eddy is excited in the interior which propagates eastward leaving the forcing wave trailing behind it. In the highly nonlinear limit, the excited vortex remains coherent against the effects of dispersion even when the northern boundary forcing is turned off. Through a process of Rossbywave shedding, it adjusts itself to the freemodon solution of the potential vorticity equation with a parabolic topography, as shown by using modon diagnostics.
Abstract
In Part I of the present work we studied the transient Rossby wave radiation excited in the far field by a northern boundary forcing. We proposed as a possible mechanism responsible for transient pulses of Rossby waves the sudden growth to finite amplitude and successive amplitude pulsations of a stationary or eastward moving meander of an intense current like the Gulf Stream. In Part I the linear regime was thoroughly analyzed and the major findings were: 1) the transient part of the response initially excited by “switching on” the meander is the one responsible for the Rossby waves radiating and being clearly detectable in the far field; 2) the Rossby radiation is excited only when the meander pulsation frequency is below a critical value.
In this paper, we extend the results of Part I to the highly nonlinear regime addressing the problem of the production of nonlinear, coherent structures through the same boundary forcing mechanisms. The model used is the quasigeostrophic, potential vorticity equation on a βplane with or without topography, in a periodic channel forced by a northern boundary function. The forcing functions are designed to simulate (i) a stationary pulsating meander, (ii) a slowly propagating meander, (iii) a combination of both. This may constitute the simplest possible idealization of meander growth to finite amplitude, its successive steepening and bending with pinching off of eddies as frequently observed in the Gulf Stream system.
Using a boundary forcing idealizing a steady meander with a pulsating amplitude, we excite nonlinear vortices on the βplane without topography. The related regime is called the “modon” regime as, at every pulsation cycle, two circulation patterns are produced, one cyclonic, the other anticyclonic. These vortices are part of a dipole pair and selfadvect eastward or westward depending on their polarity (lowabovehigh or inverse). Through diagnostic tools we show that these vortices are ‘reasonable’ modons, that is highly nonlinear, exact solutions of the equivalent barotropic model on the βplane.
The second type of boundary forcing capitalizes upon a resonance mechanism which was previously studied analytically in the weakly nonlinear regime. In it, an interior vortex is resonantly excited by a meander with constant amplitude which propagates at one of the eigenspeeds of the free linear modes of the channel. This mechanism needs topography to be effective when using quasigeostrophic dynamics. With parabolic relief all along the channel, a resonantly forced nonlinear eddy is excited in the interior which propagates eastward leaving the forcing wave trailing behind it. In the highly nonlinear limit, the excited vortex remains coherent against the effects of dispersion even when the northern boundary forcing is turned off. Through a process of Rossbywave shedding, it adjusts itself to the freemodon solution of the potential vorticity equation with a parabolic topography, as shown by using modon diagnostics.
Abstract
The dynamical consequences of constraining a numerical model with sea surface height data have been investigated. The model used for this study is a quasigeostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. The assimilations of the surface data is thus equivalent to the prescription of a surface pressure boundary condition. The authors analyzed the mechanisms of the model adjustment and the characteristics of the resultant equilibrium state when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final equilibrium state, two different experiments have been considered: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean plus eddies) has been used. It is shown that the model behavior in the presence of the surface data constraint can be conveniently described in terms of baroclinic Fofonoff modes. The prescribed mean component of the surface data acts as a “surface topography” in this problem. Its presence determines a distortion of the geostrophic contours in the subsurface layers, thus constraining the mean circulation in those layers. The intensity of the mean flow is determined by the inflow/outflow conditions at the open boundaries, as well as by eddy forcing and dissipation.
Abstract
The dynamical consequences of constraining a numerical model with sea surface height data have been investigated. The model used for this study is a quasigeostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. The assimilations of the surface data is thus equivalent to the prescription of a surface pressure boundary condition. The authors analyzed the mechanisms of the model adjustment and the characteristics of the resultant equilibrium state when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final equilibrium state, two different experiments have been considered: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean plus eddies) has been used. It is shown that the model behavior in the presence of the surface data constraint can be conveniently described in terms of baroclinic Fofonoff modes. The prescribed mean component of the surface data acts as a “surface topography” in this problem. Its presence determines a distortion of the geostrophic contours in the subsurface layers, thus constraining the mean circulation in those layers. The intensity of the mean flow is determined by the inflow/outflow conditions at the open boundaries, as well as by eddy forcing and dissipation.
Abstract
A reducedgravity, primitive equation, upperocean GCM is used to study subduction pathways in the Atlantic subtropical and tropical gyres. In order to compare the different responses in the pathways to strong and weak wind stress forcings, Hellerman and Rosenstein (HR) and da Silva (DSV) climatological annualmean and monthly wind stress forcings are used to force the model. It is shown that subtropical–tropical communication is dependent on both the strength and structure of the wind forcing. A comparison between the two experiments shows two results for the North Atlantic: 1) the full communication window between the subtropical and tropical gyres is similar in width despite the difference in the intensity of the winds and 2) the interior exchange window width is substantially larger in the weak forcing experiment (DSV) than the strong forcing experiment (HR), accompanied by a larger transport as well. The South Atlantic exhibits a similar communication between the subtropics and Tropics in both cases. The annualmean of the seasonally varying forcing also supports these results. A twolayer ventilated thermocline model is developed with a zonally varying, even though idealized, wind stress in the North Atlantic, which includes the upward Ekman pumping region absent from the classical ventilated thermocline model. The model shows that the communication window for subduction pathways is a function of the zonal gradient of the Ekman pumping velocity, not the Ekman pumping itself, at outcrop lines and at the boundary between the subtropical and tropical gyres. This solution is validated using three additional GCM experiments. It is shown that the communication windows are primarily explained by the ventilated thermocline model without considering the buoyancy effects. From the GCM experiments, the interior exchange window, which is a part of the communication window and cannot be explained by the ventilated thermocline model, is widened by two factors: 1) eliminating part of the positive Ekman pumping region in the eastern North Atlantic and 2) weakening the Ekman pumping over the whole region. The implications of these results suggest that changes in the wind forcing on the order of the difference in the wind products used here can have a significant effect on the attributes of the communication window and, hence, the thermocline structure at lower latitudes.
Abstract
A reducedgravity, primitive equation, upperocean GCM is used to study subduction pathways in the Atlantic subtropical and tropical gyres. In order to compare the different responses in the pathways to strong and weak wind stress forcings, Hellerman and Rosenstein (HR) and da Silva (DSV) climatological annualmean and monthly wind stress forcings are used to force the model. It is shown that subtropical–tropical communication is dependent on both the strength and structure of the wind forcing. A comparison between the two experiments shows two results for the North Atlantic: 1) the full communication window between the subtropical and tropical gyres is similar in width despite the difference in the intensity of the winds and 2) the interior exchange window width is substantially larger in the weak forcing experiment (DSV) than the strong forcing experiment (HR), accompanied by a larger transport as well. The South Atlantic exhibits a similar communication between the subtropics and Tropics in both cases. The annualmean of the seasonally varying forcing also supports these results. A twolayer ventilated thermocline model is developed with a zonally varying, even though idealized, wind stress in the North Atlantic, which includes the upward Ekman pumping region absent from the classical ventilated thermocline model. The model shows that the communication window for subduction pathways is a function of the zonal gradient of the Ekman pumping velocity, not the Ekman pumping itself, at outcrop lines and at the boundary between the subtropical and tropical gyres. This solution is validated using three additional GCM experiments. It is shown that the communication windows are primarily explained by the ventilated thermocline model without considering the buoyancy effects. From the GCM experiments, the interior exchange window, which is a part of the communication window and cannot be explained by the ventilated thermocline model, is widened by two factors: 1) eliminating part of the positive Ekman pumping region in the eastern North Atlantic and 2) weakening the Ekman pumping over the whole region. The implications of these results suggest that changes in the wind forcing on the order of the difference in the wind products used here can have a significant effect on the attributes of the communication window and, hence, the thermocline structure at lower latitudes.