Description: If unstructured meshes are refined to locally represent eddy dynamics in ocean circulation models, a practical question arises on how to vary the resolution and where to deploy the refinement. We propose to use the observed sea surface height variability as the refinement criterion. We explore the utility of this method (i) in a suite of idealized experiments simulating a wind-driven double gyre flow in a stratified circular basin and (ii) in simulations of global ocean circulation performed with FESOM. Two practical approaches of mesh refinement are compared. In the first approach the uniform refinement is confined within the areas where the observed variability exceeds a given threshold. In the second one the refinement varies linearly following the observed variability. The resolution is fixed in time. For the double gyre case it is shown that the variability obtained in a high-resolution reference run can be well captured on variable-resolution meshes if they are refined where the variability is high and additionally upstream the jet separation point. The second approach of mesh refinement proves to be more beneficial in terms of improvement downstream the midlatitude jet. Similarly, in global ocean simulations the mesh refinement based on the observed variability helps the model to simulate high variability at correct locations. The refinement also leads to a reduced bias in the upper-ocean temperature

Description: AGU Chapman Conference: The Agulhas System and Its Role in Changing Ocean Circulation, Climate, and Marine Ecosystems; Stellenbosch, South Africa, 8–12 October 2012

Description: We examine the mean pathways, transit timescales, and transformation of waters flowing from the Pacific and the marginal seas through the Indian Ocean (IO) on their way toward the South Atlantic within a high-resolution ocean/sea-ice model. The model fields are analyzed from a Lagrangian perspective where water volumes are tracked as they enter the IO. The IO contributes 12.6 Sv to Agulhas leakage, which within the model is 14.1 ± 2.2 Sv, the rest originates from the South Atlantic. The Indonesian Through-flow constitutes about half of the IO contribution, is surface bound, cools and salinificates as it leaves the basin within 10–30 years. Waters entering the IO south of Australia are at intermediate depths and maintain their temperature-salinity properties as they exit the basin within 15–35 years. Of these waters, the contribution from Tasman leakage is 1.4 Sv. The rest stem from recirculation from the frontal regions of the Southern Ocean. The marginal seas export 1.0 Sv into the Atlantic within 15–40 years, and the waters cool and freshen on-route. However, the model's simulation of waters from the Gulfs of Aden and Oman are too light and hence overly influenced by upper ocean circulations. In the Cape Basin, Agulhas leakage is well mixed. On-route, temperature-salinity transformations occur predominantly in the Arabian Sea and within the greater Agulhas Current region. Overall, the IO exports at least 7.9 Sv from the Pacific to the Atlantic, thereby quantifying the strength of the upper cell of the global conveyor belt.

Description: Large amounts of methane hydrate locked up within marine sediments are vulnerable to climate change. Changes in bottom water temperatures may lead to their destabilization and the release of methane into the water column or even the atmosphere. In a multimodel approach, the possible impact of destabilizing methane hydrates onto global climate within the next century is evaluated. The focus is set on changing bottom water temperatures to infer the response of the global methane hydrate inventory to future climate change. Present and future bottom water temperatures are evaluated by the combined use of hindcast high-resolution ocean circulation simulations and climate modeling for the next century. The changing global hydrate inventory is computed using the parameterized transfer function recently proposed by Wallmann et al. (2012). We find that the present-day world's total marine methane hydrate inventory is estimated to be 1146Gt of methane carbon. Within the next 100years this global inventory may be reduced by ∼0.03% (releasing ∼473Mt methane from the seafloor). Compared to the present-day annual emissions of anthropogenic methane, the amount of methane released from melting hydrates by 2100 is small and will not have a major impact on the global climate. On a regional scale, ocean bottom warming over the next 100years will result in a relatively large decrease in the methane hydrate deposits, with the Arctic and Blake Ridge region, offshore South Carolina, being most affected.

Description: Aim: The lives of juvenile leatherback turtles are amongst the most enigmatic of all marine mega-vertebrates. For these cryptic organisms, ocean models provide important insights into their dispersion from natal sites. Here, corroborated by fisheries bycatch data, we simulate spatio-temporal variation in hatchling dispersion patterns over five decades from the World's largest leatherback turtle nesting region. Location: Equatorial Central West Africa (3.5°N to −6°S) spanning the Gulf of Guinea in the North, Gabon and the Republic/Democratic Republic of the Congo in the South. Results: Due to dynamic oceanic conditions at these equatorial latitudes, dispersion scenarios differed significantly: (1) along the north to south gradient of the study region, (2) seasonally and (3) between years. From rookeries to the north of the equator, simulated hatchling retention rates within the Gulf of Guinea were very high (〉99%) after 6 months of drift, whilst south of the equator, retention rates were as low as c. 6% with the majority of simulated hatchlings dispersing west into the South Atlantic Ocean with the South Equatorial Current. Seasonal dispersion variability was driven by wind changes arising from the yearly north/southward migration of the intertropical convergence zone resulting in the increasing westerly dispersion of hatchlings throughout the hatching season. Annual variability in wind stress drove a long-term trend for decreased retention within the Gulf of Guinea and increased westerly dispersion into habitats in the South Atlantic Ocean. Main conclusions: Shifts in dispersion habitats arising from spatio-temporal oceanic variability expose hatchlings to different environments and threats that will influence important life history attributes such as juvenile growth/survival rates; anticipated to impact the population dynamics and size/age structure of populations into adulthood. The impacts of local and dynamic oceanic conditions thus require careful considerations, such as subregional management, when managing marine populations of conservation concern.

Description: The North Atlantic Current (NAC) is subject to variability on multiannual to decadal time scales, influencing the transport of volume, heat, and freshwater from the subtropical to the eastern subpolar North Atlantic (NA). Current observational time series are either too short or too episodic to study the processes involved. Here we compare the observed continuous NAC transport time series at the western flank of the Mid-Atlantic Ridge (MAR) and repeat hydrographic measurements at the OVIDE line in the eastern Atlantic with the NAC transport and circulation in the high-resolution (1/20°) ocean model configuration VIKING20 (1960–2008). The modeled baroclinic NAC transport relative to 3400 m (24.5 ± 7.1 Sv) at the MAR is only slightly lower than the observed baroclinic mean of 27.4 ± 4.7 Sv from 1993 to 2008, and extends further north by about 0.5°. In the eastern Atlantic, the western NAC (WNAC) carries the bulk of the transport in the model, while transport estimates based on hydrographic measurements from five repeated sections point to a preference for the eastern NAC (ENAC). The model is able to simulate the main features of the subpolar NA, providing confidence to use the model output to analyze the influence of the North Atlantic Oscillation (NAO). Model based velocity composites reveal an enhanced NAC transport across the MAR of up to 6.7 Sv during positive NAO phases. Most of that signal (5.4 Sv) is added to the ENAC transport, while the transport of the WNAC was independent of the NAO.

Description: The North Brazil Current (NBC) constitutes a bottleneck for the mean northward return flow of the Atlantic Meridional Overturning Circulation (AMOC) in the tropical South Atlantic. Previous studies suggested a link between interannual to multidecadal NBC and AMOC transport variability and proposed to use NBC observations as an index for the AMOC. Here we use a set of hindcast, sensitivity, and perturbation experiments performed within a hierarchy of ocean general circulation models to show that decadal to multidecadal buoyancy-forced changes in the basin-scale AMOC transport indeed manifest themselves in the NBC. The relation is, however, masked by a strong interannual to decadal wind-driven gyre variability of the NBC. While questioning the NBC transport as a direct index for the AMOC, the results support its potential merit for an AMOC monitoring system, provided that the wind-driven circulation variability is properly accounted for.

Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.

Description: Mesoscale dynamics of the Agulhas Current system determine the exchange between the Indian and Atlantic oceans, thereby influencing the global overturning circulation. Using a series of ocean model experiments compared to observations, we show that the representation of mesoscale eddies in the Agulhas ring path improves with increasing resolution of submesoscale flows. Simulated submesoscale dynamics are validated with time‐mean horizontal‐wavenumber spectra from satellite sea surface temperature measurements and mesoscale dynamics with spectra from sea surface height. While the Agulhas ring path in a nonsubmesoscale‐resolving (1/20)° configuration is associated with too less power spectral densities on all scales and too steep spectral slopes, the representation of the mesoscale dynamics improves when the diffusion and the dissipation of the model are reduced and some small‐scale features are resolved. Realistic power spectral densities over all scales are achieved when additionally the horizontal resolution is increased to (1/60)° and a larger portion of the submesoscale spectrum is resolved. Results of an eddy detection algorithm applied to the model outputs as well as to a gridded sea surface height satellite product show that in particular strong cyclones are much better represented when submesoscale flows are resolved by the model. The validation of the submesoscale dynamics with sea surface temperature spectra provides guidance for the choice of advection schemes and explicit diffusion and dissipation as well as for further subgrid‐scale parameterizations. For the Agulhas ring path, the use of upstream biased advection schemes without explicit diffusion and dissipation is found to be associated with realistically simulated submesoscales.