scholarly journals Intrinsic Variability of Sea Level from Global Ocean Simulations: Spatiotemporal Scales

2015 ◽  
Vol 28 (10) ◽  
pp. 4279-4292 ◽  
Author(s):  
Guillaume Sérazin ◽  
Thierry Penduff ◽  
Sandy Grégorio ◽  
Bernard Barnier ◽  
Jean-Marc Molines ◽  
...  

Abstract In high-resolution ocean general circulation models (OGCMs), as in process-oriented models, a substantial amount of interannual to decadal variability is generated spontaneously by oceanic nonlinearities: that is, without any variability in the atmospheric forcing at these time scales. The authors investigate the temporal and spatial scales at which this intrinsic oceanic variability has the strongest imprints on sea level anomalies (SLAs) using a ° global OGCM, by comparing a “hindcast” driven by the full range of atmospheric time scales with its counterpart forced by a repeated climatological atmospheric seasonal cycle. Outputs from both simulations are compared within distinct frequency–wavenumber bins. The fully forced hindcast is shown to reproduce the observed distribution and magnitude of low-frequency SLA variability very accurately. The small-scale (L < 6°) SLA variance is, at all time scales, barely sensitive to atmospheric variability and is almost entirely of intrinsic origin. The high-frequency (mesoscale) part and the low-frequency part of this small-scale variability have almost identical geographical distributions, supporting the hypothesis of a nonlinear temporal inverse cascade spontaneously transferring kinetic energy from high to low frequencies. The large-scale (L > 12°) low-frequency variability is mostly related to the atmospheric variability over most of the global ocean, but it is shown to remain largely intrinsic in three eddy-active regions: the Gulf Stream, Kuroshio, and Antarctic Circumpolar Current (ACC). Compared to its ¼° predecessor, the authors’ ° OGCM is shown to yield a stronger intrinsic SLA variability, at both mesoscale and low frequencies.

2011 ◽  
Vol 24 (21) ◽  
pp. 5652-5670 ◽  
Author(s):  
Thierry Penduff ◽  
Mélanie Juza ◽  
Bernard Barnier ◽  
Jan Zika ◽  
William K. Dewar ◽  
...  

Abstract This paper evaluates in a realistic context the local contributions of direct atmospheric forcing and intrinsic oceanic processes on interannual sea level anomalies (SLAs). A ¼° global ocean–sea ice general circulation model, driven over 47 yr by the full range of atmospheric time scales, is quantitatively assessed against altimetry and shown to reproduce most observed features of the interannual SLA variability from 1993 to 2004. Comparing this simulation with a second driven only by the climatological annual cycle reveals that the intrinsic part of the total interannual SLA variance exceeds 40% over half of the open-ocean area and exceeds 80% over one-fifth of it. This intrinsic contribution is particularly strong in eddy-active regions (more than 70%–80% in the Southern Ocean and western boundary current extensions) as predicted by idealized studies, as well as within the 20°–35° latitude bands. The atmosphere directly forces most of the interannual SLA variance at low latitudes and in most midlatitude eastern basins, in particular north of about 40°N in the Pacific. The interannual SLA variance is almost entirely due to intrinsic processes south of the Antarctic Circumpolar Current in the Indian Ocean sector, while half of this variance is forced by the atmosphere north of it. The same simulations were performed and analyzed at 2° resolution as well: switching to this laminar regime yields a comparable forced variability (large-scale distribution and magnitude) but almost suppresses the intrinsic variability. This likely explains why laminar ocean models largely underestimate the interannual SLA variance.


2017 ◽  
Vol 30 (13) ◽  
pp. 4997-5019 ◽  
Author(s):  
Stephan Juricke ◽  
Tim N. Palmer ◽  
Laure Zanna

In global ocean models, the representation of small-scale, high-frequency processes considerably influences the large-scale oceanic circulation and its low-frequency variability. This study investigates the impact of stochastic perturbation schemes based on three different subgrid-scale parameterizations in multidecadal ocean-only simulations with the ocean model NEMO at 1° resolution. The three parameterizations are an enhanced vertical diffusion scheme for unstable stratification, the Gent–McWilliams (GM) scheme, and a turbulent kinetic energy mixing scheme, all commonly used in state-of-the-art ocean models. The focus here is on changes in interannual variability caused by the comparatively high-frequency stochastic perturbations with subseasonal decorrelation time scales. These perturbations lead to significant improvements in the representation of low-frequency variability in the ocean, with the stochastic GM scheme showing the strongest impact. Interannual variability of the Southern Ocean eddy and Eulerian streamfunctions is increased by an order of magnitude and by 20%, respectively. Interannual sea surface height variability is increased by about 20%–25% as well, especially in the Southern Ocean and in the Kuroshio region, consistent with a strong underestimation of interannual variability in the model when compared to reanalysis and altimetry observations. These results suggest that enhancing subgrid-scale variability in ocean models can improve model variability and potentially its response to forcing on much longer time scales, while also providing an estimate of model uncertainty.


2018 ◽  
Vol 48 (6) ◽  
pp. 1385-1408 ◽  
Author(s):  
Guillaume Sérazin ◽  
Thierry Penduff ◽  
Bernard Barnier ◽  
Jean-Marc Molines ◽  
Brian K. Arbic ◽  
...  

AbstractA seasonally forced 1/12° global ocean/sea ice simulation is used to characterize the spatiotemporal inverse cascade of kinetic energy (KE). Nonlinear scale interactions associated with relative vorticity advection are evaluated using cross-spectral analysis in the frequency–wavenumber domain from sea level anomaly (SLA) time series. This analysis is applied within four eddy-active midlatitude regions having large intrinsic variability spread over a wide range of scales. Over these four regions, mesoscale surface KE is shown to spontaneously cascade toward larger spatial scales—between the deformation scale and the Rhines scale—and longer time scales (possibly exceeding 10 years). Other nonlinear processes might have to be invoked to explain the longer time scales of intrinsic variability, which have a substantial surface imprint at midlatitudes. The analysis of a fully forced 1/12° hindcast shows that low-frequency and synoptic atmospheric forcing barely affects this inverse KE cascade. The inverse cascade is also at work in a 1/4° simulation, albeit with a weaker intensity, consistent with the weaker intrinsic variability found at this coarser resolution. In the midlatitude North Pacific, the spatiotemporal cascade transfers KE from high-frequency frontal Rossby waves (FRWs), probably generated by baroclinic instability, toward the lower-frequency, westward-propagating mesoscale eddy (WME) field. The WMEs provide local gradients of potential vorticity that support these short Doppler-shifted FRWs. FRWs have periods shorter than 2 months and might be subsampled by altimetric observations, perhaps explaining why the temporal inverse cascade deduced from high-resolution models and mapped altimeter products can be quite different. The nature of the nonlinear interactions between FRWs and WMEs remains unclear but might involve wave turbulence processes.


2019 ◽  
Author(s):  
Donald A. Slater ◽  
Denis Felikson ◽  
Fiamma Straneo ◽  
Heiko Goelzer ◽  
Christopher M. Little ◽  
...  

Abstract. Changes in the ocean are expected to be an important determinant of the Greenland Ice Sheet's future sea level contribution. Yet representing these changes in continental-scale ice sheet models remains challenging due to the small scale of the key physics, and limitations in processing understanding. Here we present the ocean forcing strategy for Greenland Ice Sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change 6th Assessment Report. Beginning from global atmosphere-ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland Ice Sheet models, termed the retreat and submarine melt implementations. The retreat implementation parameterizes glacier retreat as a function of projected submarine melting, is designed to be implementable by all ice sheet models, and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat, and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly-constrained parameterisations and as such, further research on submarine melting, calving and fjord-shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland Ice Sheet models to be systematically and consistently forced by the ocean for the first time, and should therefore result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.


2021 ◽  
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  

<p>The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolution configurations, limiting these studies to individual glaciers or regions over short time scales of decades to a few centuries. To study global and long term interactions, we developed a framework to couple the dynamic ice sheet model PISM with the global ocean general circulation model MOM5 via the ice-shelf cavity module PICO. Since ice-shelf cavities are not resolved by MOM5, but parameterized with the box model PICO, the framework allows the ice sheet and ocean model to be run at resolution of 16 km and 3 degrees, respectively. We present first results from our coupled setup and discuss stability, feedbacks, and interactions of the Antarctic Ice Sheet and the global ocean system on millennial time scales.</p>


2010 ◽  
Vol 23 (15) ◽  
pp. 4243-4254 ◽  
Author(s):  
K. Lorbacher ◽  
J. Dengg ◽  
C. W. Böning ◽  
A. Biastoch

Abstract Some studies of ocean climate model experiments suggest that regional changes in dynamic sea level could provide a valuable indicator of trends in the strength of the Atlantic meridional overturning circulation (MOC). This paper describes the use of a sequence of global ocean–ice model experiments to show that the diagnosed patterns of sea surface height (SSH) anomalies associated with changes in the MOC in the North Atlantic (NA) depend critically on the time scales of interest. Model hindcast simulations for 1958–2004 reproduce the observed pattern of SSH variability with extrema occurring along the Gulf Stream (GS) and in the subpolar gyre (SPG), but they also show that the pattern is primarily related to the wind-driven variability of MOC and gyre circulation on interannual time scales; it is reflected also in the leading EOF of SSH variability over the NA Ocean, as described in previous studies. The pattern, however, is not useful as a “fingerprint” of longer-term changes in the MOC: as shown with a companion experiment, a multidecadal, gradual decline in the MOC [of 5 Sv (1 Sv ≡ 106 m3 s−1) over 5 decades] induces a much broader, basin-scale SSH rise over the mid-to-high-latitude NA, with amplitudes of 20 cm. The detectability of such a trend is low along the GS since low-frequency SSH changes are effectively masked here by strong variability on shorter time scales. More favorable signal-to-noise ratios are found in the SPG and the eastern NA, where a MOC trend of 0.1 Sv yr−1 would leave a significant imprint in SSH already after about 20 years.


2006 ◽  
Vol 23 (4) ◽  
pp. 619-629 ◽  
Author(s):  
Rui M. Ponte

Abstract For a dynamical interpretation of sea level records, estimates are needed of the isostatic, or so-called inverted barometer, signals (ηib) associated with the ocean response to atmospheric loading. Seasonal and longer-period ηib signals are evaluated over the global ocean for the period 1958–2000 using monthly sea level pressure fields from two different atmospheric reanalyses. Variability and linear trends in ηib agree well for the two reanalyses in most regions but less so over the Southern Ocean, where uncertainties in ηib seem to be largest. The standard deviation of ηib ranges from <1 cm in equatorial regions to >7 cm in the regions of the Aleutian and Iceland lows and parts of the Southern and Arctic Oceans. When compared to a global tide gauge dataset, both seasonal and interannual ηib signals are found to contribute importantly to the sea level variance in many mid- and high-latitude records, with seasonal signals important as well in tropical records from India and Southeast Asia. For these records, subtracting ηib from the data can lead to changes in variance of 40% or more. Over the period of study, linear trends in ηib are mostly negative at low and midlatitudes and can cause negative biases in tide gauge estimates of global mean sea level rise that are comparable in magnitude to the effects of postglacial rebound. In agreement with previous findings, ηib signals are found to introduce anomalous behavior in local records (e.g., substantially weaker upward trends in the Mediterranean), and their removal can also reduce formal trend uncertainties. Accounting for ηib effects can be even more important when analyzing relatively short (decadal) records, such as those available from satellite altimetry.


2011 ◽  
Vol 41 (2) ◽  
pp. 365-377 ◽  
Author(s):  
Thomas Kilpatrick ◽  
Niklas Schneider ◽  
Emanuele Di Lorenzo

Abstract The generation of variance by anomalous advection of a passive tracer in the thermocline is investigated using the example of density-compensated temperature and salinity anomalies, or spiciness. A coupled Markov model is developed in which wind stress curl forces the large-scale baroclinic ocean pressure that in turn controls the anomalous geostrophic advection of spiciness. The “double integration” of white noise atmospheric forcing by this Markov model results in a frequency (ω) spectrum of large-scale spiciness proportional to ω−4, so that spiciness variability is concentrated at low frequencies. An eddy-permitting regional model hindcast of the northeast Pacific (1950–2007) confirms that time series of large-scale spiciness variability are exceptionally smooth, with frequency spectra ∝ ω−4 for frequencies greater than 0.2 cpy. At shorter spatial scales (wavelengths less than ∼500 km), the spiciness frequency spectrum is whitened by mesoscale eddies, but this eddy-forced variability can be filtered out by spatially averaging. Large-scale and long-term measurements are needed to observe the variance of spiciness or any other passive tracer subject to anomalous advection in the thermocline.


2021 ◽  
Vol 14 (4) ◽  
pp. 2011-2028
Author(s):  
Qing Li ◽  
Luke Van Roekel

Abstract. A multiscale modeling approach for studying the ocean surface turbulent mixing is explored by coupling an ocean general circulation model (GCM) MPAS-Ocean with the Parallelized Large Eddy Simulation Model (PALM). The coupling approach is similar to the superparameterization approach that has been used to represent the effects of deep convection in atmospheric GCMs. However, the focus of this multiscale modeling approach is on the small-scale turbulent mixing and their interactions with the larger-scale processes in the ocean, so that a more flexible coupling strategy is used. To reduce the computational cost, a customized version of PALM is ported on the general-purpose graphics processing unit (GPU) with OpenACC, achieving 10–16 times overall speedup as compared to running on a single CPU. Even with the GPU-acceleration technique, a superparameterization-like approach to represent the ocean surface turbulent mixing in GCMs using embedded high fidelity and three-dimensional large eddy simulations (LESs) over the global ocean is still computationally intensive and infeasible for long simulations. However, running PALM regionally on selected MPAS-Ocean grid cells is shown to be a promising approach moving forward. The flexible coupling between MPAS-Ocean and PALM allows further exploration of the interactions between the ocean surface turbulent mixing and larger-scale processes, as well as future development and improvement of ocean surface turbulent mixing parameterizations for GCMs.


Sign in / Sign up

Export Citation Format

Share Document