Using the Mellor–Yamada mixing scheme in seasonally ice-covered seas—Corrigendum to: Parameterization of vertical mixing in the Weddell Sea [Ocean Modelling 6 (2004) 83–100]

A coding error in the s-Coordinate Primitive Equation Model (SPEM) has led to misleading statements about the behaviour of the Mellor–Yamada level 2 parameterization of vertical mixing. It has been claimed that the scheme removes static instability only very slowly and preserves statically unstable stratifications for an unrealistic long time. This note corrects this statement by demonstrating that the Mellor–Yamada mixing scheme, if implemented correctly, tends to overestimate rather than underestimate vertical mixing in seasonally ice-covered seas. Similar to other mixing schemes with the same behaviour, this leads to spurious open ocean deep convection, an unrealistic homogenization of the water column, and a significant reduction of sea ice volume.     

Using paleoclimate proxy-data to select optimal realisations in an ensemble of simulations of the climate of the past millennium

We present and describe in detail the advantages and limitations of a technique that combines in an optimal way model results and proxy-data time series in order to obtain states of the climate system consistent with model physics, reconstruction of past radiative forcing and proxy records. To achieve this goal, we select among an ensemble of simulations covering the last millennium performed with a low-resolution 3-D climate model the ones that minimise a cost function. This cost function measures the misfit between model results and proxy records. In the framework of the tests performed here, an ensemble of 30 to 40 simulations appears sufficient to reach reasonable correlations between model results and reconstructions, in configurations for which a small amount of data is available as well as in data-rich areas. Preliminary applications of the technique show that it can be used to provide reconstructions of past large-scale temperature changes, complementary to the ones obtained by statistical methods. Furthermore, as model results include a representation of atmospheric and oceanic circulations, it can be used to provide insights into some amplification mechanisms responsible for past temperature changes. On the other hand, if the number of proxy records is too low, it could not be used to provide reconstructions of past changes at a regional scale.     

Surface temperature control in the North and tropical Pacific during the last glacial maximum

Based on coupled modelling evidence we argue that topographically-induced modifications of the large-scale atmospheric circulation during the last glacial maximum may have led to a reduction of the westerlies, and a slowdown of the Pacific subtropical gyre as well as to an intensification of the Pacific subtropical cell. These oceanic circulation changes generate an eastern North Pacific warming, an associated cooling in the Kuroshio area, as well as a cooling of the tropical oceans, respectively. The tropical cooling pattern resembles a permanent La Niña state which in turn forces atmospheric teleconnection patterns that lead to an enhancement of the subtropical warming by reduced latent and sensible cooling of the ocean. In addition, the radiative cooling due to atmospheric CO₂ and water vapor reductions imposes a cooling tendency in the tropics and subtropics, thereby intensifying the permanent La Niña conditions. The remote North Pacific response results in a warming tendency of the eastern North Pacific which may level off the effect of the local radiative cooling. Hence, a delicate balance between oceanic circulation changes, remotely induced atmospheric flux anomalies as well local radiative cooling is established which controls the tropical and North Pacific temperature anomalies during the last glacial maximum. Furthermore, we discuss how the aftermath of a Heinrich event may have affected glacial temperatures in the Pacific Ocean.     

Modes of climate variability as simulated by a coupled general circulation model. Part I: ENSO-like climate variability and its low-frequency modulation

The El Niño-Southern Oscillation (ENSO) is investigated in a multicentury integration conducted with the coupled general circulation model (CGCM) ECHAM3/LSG. The quasiperiodic interannual oscillations of the simulated equatorial Pacific climate system are due to subsurface temperature anomaly propagation and a positive atmosphere-ocean feedback. The gravest internal wave modes contribute to the generation of these anomalies. The simulated ENSO has a characteristic period of 5–8 years. Due to the coarse resolution of the ocean model the ENSO amplitude is underestimated by a factor of three as compared to observations. The model ENSO is associated with the typical atmospheric teleconnection patterns. Using wavelet statistics two characteristic interdecadal modulations of the ENSO variance are identified. The origins of a 22 and 35?y ENSO modulation as well as the characteristic ENSO response to greenhouse warming simulated by our model are discussed.     

Southern Ocean warming and increased ice shelf basal melting in the twenty-first and twenty-second centuries based on coupled ice-ocean finite-element modelling

We utilise a global finite-element sea ice–ocean model (FESOM), focused on the Antarctic marginal seas, to analyse projections of ice shelf basal melting in a warmer climate. Ice shelf–ocean interaction is described using a three-equation system with a diagnostic computation of temperature and salinity at the ice–ocean interface. A tetrahedral mesh with a minimumhorizontal resolution of 4 km and hybrid vertical coordinates is used. Ice shelf draft, cavity geometry, and global ocean bathymetry have been derived from the RTopo-1 data set. The model is forced with the atmospheric output from two climate models: (1) the Hadley Centre Climate Model (HadCM3) and (2) Max Planck Institute’s ECHAM5/MPI-OM coupled climate model. Results from experiments forced with their twentieth century output are used to evaluate the modelled present-day ocean state. Sea ice coverage is largely realistic in both simulations; modelled ice shelf basal melt rates compare well with observations in both cases, but are consistently smaller for ECHAM5/MPI-OM. Projections for future ice shelf basal melting are computed using atmospheric output for the Intergovernmental Panel on Climate Change (IPCC) scenarios E1 and A1B. In simulations forced with ECHAM5 data, trends in ice shelf basal melting are small. In contrast, decreasing convection along the Antarctic coast in HadCM3 scenarios leads to a decreasing salinity on the continental shelf and to intrusions of warm deep water of open ocean origin. In the case of the Filchner–Ronne Ice Shelf (FRIS), this water reaches deep into the cavity, so that basal melting increases by a factor of 4 to 6 compared to the present value of about 90 Gt/year. By the middle of the twenty-second century, FRIS becomes the dominant contributor to total ice shelf basal mass loss in these simulations. Our results indicate that the surface freshwater fluxes on the continental shelves may be crucial for the future of especially the large cold water ice shelves in the Southern Ocean.     

The use of tracers to assess drill-mud penetration depth into sandstone cores during deep drilling: method development and application

For the utilization of deep saline aquifers in the frame of geotechnical use, such as geological sequestration of CO₂, H₂ or energy storage, a baseline characterization of pristine reservoir rock cores is required to monitor changes in the indigenous microbial communities and pore fluids, and to study alterations in rock characteristics resulting from interaction with geological storage technologies. However, drilling procedures and technical fluids, particularly drill mud, are sources of core contamination. To measure the penetration of drill mud into the cores, three tracers (fluorescein, microspheres, and 4′,6-diamidino-2-phenylindole stained bacteria) were tested under laboratory conditions. The flow of drill mud into core samples was induced by applying uniaxial pressure differentials to the core, and the penetration depth was microscopically determined for each tracer. Fluorescein was extracted from the rock samples and quantified fluorometrically. The results indicate that all tested tracers are suitable for tracking drill-mud penetration. The actual penetration depth seems to be related to differences in mineral composition and texture as well as microfractures. Among all tested tracers, fluorescein labelling is the simplest, cheapest and most accurate method for analyzing the contamination of rock cores by technical fluids. The application of this tracer was successfully applied during two deep drilling campaigns at the CO₂ storage pilot site in Ketzin, Germany. The results highlight that the use of tracers is indispensable to ensure the quality of core samples for microbiological and biogeochemical analysis.     

Impact of diurnal atmosphere–ocean coupling on tropical climate simulations using a coupled GCM

The impacts of diurnal atmosphere–ocean (air–sea) coupling on tropical climate simulations are investigated using the SNU coupled GCM. To investigate the effect of the atmospheric and oceanic diurnal cycles on a climate simulation, a 1-day air–sea coupling interval experiment is compared to a 2-h coupling experiment. As previous studies have suggested, cold temperature biases over equatorial western Pacific regions are significantly reduced when diurnal air–sea coupling strategy is implemented. This warming is initiated by diurnal rectification and amplified further by the air–sea coupled feedbacks. In addition to its effect on the mean climatology, the diurnal coupling has also a distinctive impact on the amplitude of the El Nino-Southern Oscillation (ENSO). It is demonstrated that a weakening of the ENSO magnitude is caused by reduced (increased) surface net heat fluxes into the ocean during El Nino (La Nina) events. Primarily, decreased (increased) incoming shortwave radiation during El Nino (La Nina) due to cloud shading is responsible for the net heat fluxes associated with ENSO.