ALBERT

Description: Changes in the background climate are known to affect El Niño–Southern Oscillation (ENSO) by altering feedbacks that control ENSO’s characteristics. Here, the sensitivity of ENSO variability to the background climate is investigated by utilizing two Community Earth System Model, version 1 (CESM1), simulations in which the solar constant is altered by ±25 W m−2. The resulting stable warm and cold climate mean state simulations differ in terms of ENSO amplitude, frequency, diversity, asymmetry, and seasonality. In the warm run, ENSO reveals a larger amplitude and occurs at higher frequencies relative to the cold and control runs as well as observations. The warm run also features more eastern Pacific El Niños, an increased asymmetry, and a stronger seasonal phase locking. These changes are linked to changes in the mean state via the amplifying and damping feedbacks. In the warm run, a shallower mean thermocline results in a stronger subsurface–surface coupling, whereas the cold run reveals reduced ENSO variability due to a reduced Bjerknes feedback in accordance with a deeper mean thermocline and enhanced surface wind stress. A strong zonal advective and upwelling feedback further contribute to the large ENSO amplitude in the run with a warmer mean state. In the cold run, ENSO events are partly forced by anomalous shortwave radiation. However, in light of the large temperature contrast between the simulations of up to 6 K in the tropical Pacific, the relatively small changes in ENSO variability highlight the robustness of ENSO dynamics under vastly different climate mean states.

Description: Two decades of high-resolution satellite observations and climate modeling studies have indicated strong ocean–atmosphere coupled feedback mediated by ocean mesoscale processes, including semipermanent and meandrous SST fronts, mesoscale eddies, and filaments. The air–sea exchanges in latent heat, sensible heat, momentum, and carbon dioxide associated with this so-called mesoscale air–sea interaction are robust near the major western boundary currents, Southern Ocean fronts, and equatorial and coastal upwelling zones, but they are also ubiquitous over the global oceans wherever ocean mesoscale processes are active. Current theories, informed by rapidly advancing observational and modeling capabilities, have established the importance of mesoscale and frontal-scale air–sea interaction processes for understanding large-scale ocean circulation, biogeochemistry, and weather and climate variability. However, numerous challenges remain to accurately diagnose, observe, and simulate mesoscale air–sea interaction to quantify its impacts on large-scale processes. This article provides a comprehensive review of key aspects pertinent to mesoscale air–sea interaction, synthesizes current understanding with remaining gaps and uncertainties, and provides recommendations on theoretical, observational, and modeling strategies for future air–sea interaction research.

Description: Several years of moored turbulence measurements from xpods at three sites in the equatorial cold tongues of Atlantic and Pacific Oceans yield new insights into proxy estimates of turbulence that specifically target the cold tongues. They also reveal previously unknown wind dependencies of diurnally varying turbulence in the near-critical stratified shear layers beneath the mixed layer and above the core of the Equatorial Undercurrent that we have come to understand as deep cycle (DC) turbulence. Isolated by the mixed layer above, the DC layer is only indirectly linked to surface forcing. Yet, it varies diurnally in concert with daily changes in heating/cooling. Diurnal composites computed from 10-min averaged data at fixed xpod depths show that transitions from daytime to nighttime mixing regimes are increasingly delayed with weakening wind stress t. These transitions are also delayed with respect to depth such that they follow a descent rate of roughly 6 m h-1, independent of t. We hypothesize that this wind-dependent delay is a direct result of wind-dependent diurnal warm layer deepening, which acts as the trigger to DC layer instability by bringing shear from the surface down-ward but at rates much slower than 6 m h-1. This delay in initiation of DC layer instability contributes to a reduction in daily averaged values of turbulence dissipation. Both the absence of descending turbulence in the sheared DC layer prior to arrival of the diurnal warm layer shear and the magnitude of the subsequent descent rate after arrival are roughly predicted by laboratory experiments on entrainment in stratified shear flows.

Description: Numerical weather prediction models operate on grid spacings of a few kilometers, where deep convection begins to become resolvable. Around this scale, the emergence of coherent structures in the planetary boundary layer, often hypothesized to be caused by cold pools, forces the transition from shallow to deep convection. Yet, the kilometer-scale range is typically not resolved by standard surface operational measurement networks. The measurement campaign FESSTVaL aimed at addressing this gap by observing atmospheric variability at the hectometer to kilometer scale, with a particular emphasis on cold pools, wind gusts and coherent patterns in the planetary boundary layer during summer. A unique feature was the distribution of 150 self-developed and low-cost instruments. More specifically, FESSTVaL included dense networks of 80 autonomous cold pool loggers, 19 weather stations and 83 soil sensor systems, all installed in a rural region of 15-km radius in eastern Germany, as well as self-developed weather stations handed out to citizens. Boundary layer and upper air observations were provided by 8 Doppler lidars and 4 microwave radiometers distributed at 3 supersites; water vapor and temperature were also measured by advanced lidar systems and an infrared spectrometer; and rain was observed by a X-band radar. An uncrewed aircraft, multicopters and a small radiometer network carried out additional measurements during a four-week period. In this paper, we present FESSTVaL’s measurement strategy and show first observational results including unprecedented highly-resolved spatio-temporal cold-pool structures, both in the horizontal as well as in the vertical dimension, associated with overpassing convective systems.

Description: As an important external forcing, the effect of the 11-yr solar cycle on the tropical Pacific decadal variability is an interesting question. Here, we systematically investigate the phase-locking of the atmosphere and ocean covariations to the solar cycle in the tropical Pacific and propose a new mechanism to explain these decadal covariations. In both observation/reanalysis datasets and a solar cycle forced sensitivity experiment (named the SOL experiment), the ocean heat content anomalies (OHCa; 300 m) resemble a La Niña–like pattern in the solar cycle ascending phase, and the Walker circulation shifts westward. In the declining phase, the opposite is true. The accumulative solar irradiation directly contributes to this coherent decadal variability via changing the warm water volume and the solar-related heat is redistributed by the ocean dynamic processes. During the 11-yr solar cycle, the Pacific Walker circulation anomalies maintain the OHCa in the western equatorial Pacific and work as negative feedback for the eastern Pacific to help the OHCa phase transition. In addition, oceanic meridional heat transport via the subtropical cells and the propagation of off-equatorial Rossby waves also provide a lagged negative feedback to the OHCa phase transition according to the 11-yr solar cycle. The decadal coupled responses of the tropical Pacific climate system are 2 years more lag in the SOL experiment than in the observation/reanalysis. Significance Statement Here, we propose a new mechanism that the heating effect of the accumulative solar irradiation during the 11-yr solar cycle can be “integrated” into the tropical Pacific OHC and then provide a bottom-up effect on the atmosphere at decadal time scales. The strongly coupled processes in this region amplify the decadal phase-locking of the covariations to the 11-yr solar cycle. Our study demonstrates the role of the 11-yr solar cycle in the tropical Pacific decadal variability and provides a new explanation for the “bottom-up” mechanism of the solar cycle forcing. Our results update the understanding of the tropical Pacific decadal variability and may help to improve climate predictions at decadal time scales.

Description: We investigate the origin of the equatorial Pacific cold sea surface temperature (SST) bias and its link to wind biases, local and remote, in the Kiel Climate Model (KCM). The cold bias is common in climate models participating in the 5 th and 6 th phases of the Coupled Model Intercomparison Project. In the coupled experiments with the KCM, the interannually varying NCEP/CFSR wind stress is prescribed over four spatial domains: globally, over the equatorial Pacific (EP), the northern Pacific (NP) and southern Pacific (SP). The corresponding EP SST bias is reduced by 100%, 52%, 12% and 23%, respectively. Thus, the EP SST bias is mainly attributed to the local wind bias, with small but not negligible contributions from the extratropical regions. Erroneous ocean circulation driven by overly strong winds cause the cold SST bias, while the surface-heat flux counteracts it. Extratropical Pacific SST biases contribute to the EP cold bias via the oceanic subtropical gyres, which is further enhanced by dynamical coupling in the equatorial region. The origin of the wind biases is examined by forcing the atmospheric component of the KCM in a stand-alone mode with observed SSTs and simulated SSTs from the coupled experiments. Wind biases over the EP, NP and SP regions originate in the atmosphere model. The cold EP SST bias substantially enhances the wind biases over all three regions, while the NP and SP SST biases support local amplification of the wind bias. This study suggests that improving surface-wind stress, at and off the equator, is a key to improve mean-state equatorial Pacific SST in climate models.

Description: NORP-SORP Workshop on Polar Fresh Water: Sources, Pathways and Impacts of Freshwater in Northern and Southern Polar Oceans and Seas (SPICE-UP)What: Up to 60 participants at a time and more than twice as many registrants in total from 20 nations and across experience levels met to discuss the current status of research on freshwater in both polar regions, future directions, and synergies between the Arctic and Southern Ocean research communitiesWhen: 19-21 September 2022 Where: Online

Description: This study utilizes observations and a series of idealized experiments to explore whether eastern Pacific (EP)- and central Pacific (CP)-type El Niño–Southern Oscillation (ENSO) events produce surface wind stress responses with distinct spatial structures. We find that the meridionally broader sea surface temperatures (SSTs) during CP events lead to zonal wind stresses that are also meridionally broader than those found during EP-type events, leading to differences in the near-equatorial wind stress curl. These wind spatial structure differences create differences in the associated pre- and post-ENSO event WWV response. For instance, the meridionally narrow winds found during EP events have (i) weaker wind stresses along 58N and 58S, leading to weaker Ekman-induced pre-event WWV changes; and (ii) stronger near-equatorial wind stress curls that lead to a much larger post-ENSO event WWV changes than during CP events. The latter suggests that, in the framework of the recharge oscillator model, the EP events have stronger coupling between sea surface temperatures (SST) and thermocline (WWV), supporting more clearly the phase transition of ENSO events, and therefore, the oscillating nature of ENSO than CP events. The results suggest that the spatial structure of the SST pattern and the related differences in the wind stress curl, are required along with equatorial wind stress to accurately model the WWV changes during EP- and CP-type ENSO events.

Description: Equatorial deep jets (EDJ) are zonal currents along the equator in all three ocean basins that alternate in direction with depth and time. In the Atlantic below the thermocline, they are the dominant variability on interannual timescales. Observations of equatorial deep jets are available but scarce, given the EDJs’ location at depth, their small vertical scale and their long periodicity of several years. In the last few years, Argo floats have added a significant amount of measurements at intermediate depth. In this study we therefore revise estimates of the EDJ scales based on Argo float data. Mostly, we use velocity data at 1000 m depth calculated from float displacement, which yield robust estimates of the Atlantic EDJ period (4.6 years), amplitude distribution, phase distribution, zonal wavelength (146.7°), and meridional structure. We also show that the equatorial amplitude of the EDJs’ first meridional mode Rossby wave component (9.8 cm s −1 ) is larger than that of their Kelvin wave component (2.8 cm s −1 ). Additionally, we present a new estimation of the EDJs’ vertical structure throughout the Atlantic basin, based on an equatorial geostrophic velocity reconstruction from hydrographic Argo float measurements from depths between 400 and 2000 m. Our new estimates from Argo float data provide the first basin-wide assessment of the Atlantic EDJ scales, as well as having smaller uncertainties than estimates from earlier studies.

Description: The physical processes driving the genesis of surface- and subsurface-intensified cyclonic and anticyclonic eddies originating from the coastal current system of the Mauritanian Upwelling Region are investigated using a high-resolution (~1.5 km) configuration of GFDL’s Modular Ocean Model. Estimating an energy budget for the boundary current reveals a baroclinically unstable state during its intensification phase in boreal summer and which is driving eddy generation within the near-coastal region. The mean poleward coastal flow’s interaction with the sloping topography induces enhanced anticyclonic vorticity, with potential vorticity close to zero generated in the bottom boundary layer. Flow separation at sharp topographic bends intensifies the anticyclonic vorticity, and submesoscale structures of low PV coalesce to form anticyclonic vortices. A combination of offshore Ekman transport and horizontal advection determined the amount of SACW in an anticyclonic eddy. A vortex with a relatively dense and low PV core will form an anticyclonic mode-water eddy, which will subduct along isopycnals while propagating offshore and hence be shielded from surface buoyancy forcing. Less contribution of dense SACW promotes the generation of surface anticyclonic eddies as the core is composed of a lighter water mass, which causes the eddy to stay closer to the surface and hence be exposed to surface buoyancy forcing. Simulated cyclonic eddies are formed between the rotational flow of an offshore anticyclonic vortex and a poleward flowing boundary current, with eddy potential energy being the dominant source of eddy kinetic energy. All three types of eddies play a key role in the exchange between the Mauritanian Coastal currents system and the adjacent eastern boundary shadow zone region.

Description: The Earth system is accumulating energy due to human-induced activities. More than 90% of this energy has been stored in the ocean as heat since 1970, with similar to 60% of that in the upper 700 m. Differences in upper-ocean heat content anomaly (OHCA) estimates, however, exist. Here, we use a dataset protocol for 1970-2008-with six instrumental bias adjustments applied to expendable bathythermograph (XBT) data, and mapped by six research groups-to evaluate the spatiotemporal spread in upper OHCA estimates arising from two choices: 1) those arising from instrumental bias adjustments and 2) those arising from mathematical (i.e., mapping) techniques to interpolate and extrapolate data in space and time. We also examined the effect of a common ocean mask, which reveals that exclusion of shallow seas can reduce global OHCA estimates up to 13%. Spread due to mapping method is largest in the Indian Ocean and in the eddy-rich and frontal regions of all basins. Spread due to XBT bias adjustment is largest in the Pacific Ocean within 30 degrees N-30 degrees S. In both mapping and XBT cases, spread is higher for 1990-2004. Statistically different trends among mapping methods are found not only in the poorly observed Southern Ocean but also in the well-observed northwest Atlantic. Our results cannot determine the best mapping or bias adjustment schemes, but they identify where important sensitivities exist, and thus where further understanding will help to refine OHCA estimates. These results highlight the need for further coordinated OHCA studies to evaluate the performance of existing mapping methods along with comprehensive assessment of uncertainty estimates.

Description: Modular Observation Solutions of Earth Systems (MOSES) is a novel observation system that is specifically designed to unravel the impact of distinct, dynamic events on the long-term development of environmental systems. Hydrometeorological extremes such as the recent European droughts or the floods of 2013 caused severe and lasting environmental damage. Modeling studies suggest that abrupt permafrost thaw events accelerate Arctic greenhouse gas emissions. Short-lived ocean eddies seem to comprise a significant share of the marine carbon uptake or release. Although there is increasing evidence that such dynamic events bear the potential for major environmental impacts, our knowledge on the processes they trigger is still very limited. MOSES aims at capturing such events, from their formation to their end, with high spatial and temporal resolution. As such, the observation system extends and complements existing national and international observation networks, which are mostly designed for long-term monitoring. Several German Helmholtz Association centers have developed this research facility as a mobile and modular “system of systems” to record energy, water, greenhouse gas, and nutrient cycles on the land surface, in coastal regions, in the ocean, in polar regions, and in the atmosphere—but especially the interactions between the Earth compartments. During the implementation period (2017–21), the measuring systems were put into operation and test campaigns were performed to establish event-driven campaign routines. With MOSES’s regular operation starting in 2022, the observation system will then be ready for cross-compartment and cross-discipline research on the environmental impacts of dynamic events.

Description: Since 1750, land use change and fossil fuel combustion has led to a 46 % increase in the atmospheric carbon dioxide (CO2) concentrations, causing global warming with substantial societal consequences. The Paris Agreement aims to limiting global temperature increases to well below 2°C above pre-industrial levels. Increasing levels of CO2 and other greenhouse gases (GHGs), such as methane (CH4) and nitrous oxide (N2O), in the atmosphere are the primary cause of climate change. Approximately half of the carbon emissions to the atmosphere is sequestered by ocean and land sinks, leading to ocean acidification but also slowing the rate of global warming. However, there are significant uncertainties in the future global warming scenarios due to uncertainties in the size, nature and stability of these sinks. Quantifying and monitoring the size and timing of natural sinks and the impact of climate change on ecosystems are important information to guide policy-makers’ decisions and strategies on reductions in emissions. Continuous, long-term observations are required to quantify GHG emissions, sinks, and their impacts on Earth systems. The Integrated Carbon Observation System (ICOS) was designed as the European in situ observation and information system to support science and society in their efforts to mitigate climate change. It provides standardized and open data currently from over 140 measurement stations across 12 European countries. The stations observe GHG concentrations in the atmosphere and carbon and GHG fluxes between the atmosphere, land surface and the oceans. This article describes how ICOS fulfills its mission to harmonize these observations, ensure the related long-term financial commitments, provide easy access to well-documented and reproducible high-quality data and related protocols and tools for scientific studies, and deliver information and GHG-related products to stakeholders in society and policy.

Description: In the equatorial Atlantic Ocean, meridional velocity variability exhibits a pronounced peak on intraseasonal timescales whereas zonal velocity dominantly varies on seasonal to interannual timescales. We focus on the intraseasonal meridional velocity variability away from the near-surface layer, its source regions and its pathways into the deep ocean. This deep intraseasonal velocity variability plays a key role in equatorial dynamics as it is an important energy source for the deep equatorial circulation. The results are based on the output of a high-resolution ocean model revealing intraseasonal energy levels along the equator at all depths that are in good agreement with shipboard and moored velocity data. Spectral analyses reveal a pronounced signal of intraseasonal Yanai waves with westward phase velocities and zonal wavelengths longer than 450 km. Different sources and characteristics of these Yanai waves are identified: near the surface between 40°W and 10°W low-baroclinic-mode Yanai waves with periods of around 30 days are exited. These waves have a strong seasonal cycle with a maximum in August. High-frequency Yanai waves (10–20-day period) are excited at the surface east of 10°W. In the region between the North Brazil Current and the Equatorial Undercurrent high-baroclinic-mode Yanai waves with periods between 30 and 40 days are generated. Yanai waves with longer periods (40-80 days) are shed from the Deep Western Boundary Current. The Yanai wave energy is carried along beams east- and downward thus explaining differences in strength, structure and periodicity of the meridional intraseasonal variability in the equatorial Atlantic Ocean.

Description: Marine heatwaves along the coast ofWestern Australia, referred to as Ningaloo Niño, have had dramatic impacts on the ecosystem in the recent decade. A number of local and remote forcing mechanisms have been put forward, however little is known about the depth structure of such temperature extremes. Utilizing an eddy-active global Ocean General Circulation Model, Ningaloo Niño and the corresponding cold Ningaloo Niña events are investigated between 1958-2016, with focus on their depth structure. The relative roles of buoyancy and wind forcing are inferred from sensitivity experiments. Composites reveal a strong symmetry between cold and warm events in their vertical structure and associated large-scale spatial patterns. Temperature anomalies are largest at the surface, where buoyancy forcing is dominant and extend down to 300m depth (or deeper), with wind forcing being the main driver. Large-scale subsurface anomalies arise from a vertical modulation of the thermocline, extending from the western Pacific into the tropical eastern Indian Ocean. The strongest Ningaloo Niños in 2000 and 2011 are unprecedented compound events, where long-lasting high temperatures are accompanied by extreme freshening, which emerges in association with La Niñas, more common and persistent during the negative phase of the Interdecadal Pacific Oscillation. It is shown that Ningaloo Niños during La Nina phases have a distinctively deeper reach and are associated with a strengthening of the Leeuwin Current, while events during El Niño are limited to the surface layer temperatures, likely driven by local atmosphere-ocean feedbacks, without a clear imprint on salinity and velocity.

Description: The Antarctic Slope Front (ASF) is a fundamental feature of the subpolar Southern Ocean that is still poorly observed. In this study we build a statistical climatology of the temperature and salinity fields of the upper 380 m of the Antarctic margin. We use a comprehensive compilation of observational datasets including the profiles gathered by instrumented marine mammals. The mapping method consists first of a decomposition in vertical modes of the combined temperature and salinity profiles. Then the resulting principal components are optimally interpolated on a regular grid and the monthly climatological profiles are reconstructed, providing a physically plausible representation of the ocean. The ASF is located with a contour method and a gradient method applied on the temperature field, two complementary approaches that provide a complete view of the ASF structure. The front extends from the Amundsen Sea to the eastern Weddell Sea and closely tracks the continental shelf break. It is associated with a sharp temperature gradient that is stronger in winter and weaker in summer. The emergence of the front in the Amundsen and Bellingshausen sectors appears to be seasonally variable (slightly more westward in winter than in summer). Investigation of the density gradients across the shelf break indicates a winter slowdown of the baroclinic component of the Antarctic Slope Current at the near surface, in contrast with the seasonal variability of the temperature gradient.

Description: There is a controversy about the nature of multidecadal climate variability in the North Atlantic (NA) region, concerning the roles of ocean circulation and atmosphere–ocean coupling. Here we describe NA multidecadal variability from a version of the Kiel Climate Model, in which both subpolar gyre (SPG)–Atlantic meridional overturning circulation (AMOC) coupling and atmosphere–ocean coupling are essential. The oceanic barotropic and meridional overturning streamfunctions and the sea level pressure are jointly analyzed to derive the leading mode of Atlantic sector variability. This mode accounting for 23.7% of the total combined variance is oscillatory with an irregular periodicity of 25–50 years and an e-folding time of about a decade. SPG and AMOC mutually influence each other and together provide the delayed negative feedback necessary for maintaining the oscillation. An anomalously strong SPG, for example, drives higher surface salinity and density in the NA’s sinking region. In response, oceanic deep convection and AMOC intensify, which, with a time delay of about a decade, reduces SPG strength by enhancing upper-ocean heat content. The weaker gyre leads to lower surface salinity and density in the sinking region, which reduces deep convection and eventually AMOC strength. There is a positive ocean–atmosphere feedback between the sea surface temperature and low-level atmospheric circulation over the southern Greenland area, with related wind stress changes reinforcing SPG changes, thereby maintaining the (damped) multidecadal oscillation against dissipation. Stochastic surface heat flux forcing associated with the North Atlantic Oscillation drives the eigenmode.

Description: The Agulhas Current (AC) creates a sharp temperature gradient with the surrounding ocean, leading to a large turbulent flux of moisture from ocean to atmosphere. We use two simulations of the Weather Research and Forecasting (WRF) model to show the seasonal impact of the warm core of the AC on southern Africa precipitation. In one simulation the sea surface temperature (SST) of the AC is similar to satellite observations, while the second uses satellite SST observations spatially smoothed to reduce the temperature of the core of the AC by ~1.5°C. We show that decreasing the SST of the AC reduces the precipitation of the wettest seasons (austral summer and autumn) inland. Over the ocean, reducing the SST reduces precipitation, low-level wind convergence, SST and SLP Laplacian above the AC in all seasons, consistent with the pressure adjustment mechanism. Moreover, winter precipitation above the Current may be also related to increased latent flux. In summer and autumn, the AC SST reduction is also associated with decreased precipitation further inland (more than 1.5 mm/day), caused by an atmospheric circulation that decreases the horizontal moisture flux from the AC to South Africa. The reduction is also associated with higher geopotential height extending from the surface east and over the AC to the mid-troposphere over southeastern Africa. The westward tilted geopotential height is consistent with the linear response to shallow diabatic heating in midlatitudes. An identical mechanism occurs in spring but is weaker. Winter rainfall response is confined above the AC.

Description: Spatial and temporal variations of nutrient-rich upwelled water across the major eastern boundary upwelling systems are primarily controlled by the surface wind with different, and sometimes contrasting, impacts on coastal upwelling systems driven by alongshore wind and offshore upwelling systems driven by the local wind-stress-curl. Here, concurrently measured wind-fields, satellite-derived Chlorophyll-a concentration along with a state-of-the-art ocean model simulation spanning 2008-2018 are used to investigate the connection between coastal and offshore physical drivers of the Benguela Upwelling System (BUS). Our results indicate that the spatial structure of long-term mean upwelling derived from Ekman theory and the numerical model are fairly consistent across the entire BUS and closely followed by the Chlorophyll-a pattern. The variability of the upwelling from the Ekman theory is proportionally diminished with offshore distance, whereas different and sometimes opposite structures are revealed in the model-derived upwelling. Our result suggests the presence of sub-mesoscale activity (i.e., filaments and eddies) across the entire BUS with a large modulating effect on the wind-stress-curl-driven upwelling off Lüderitz and Walvis Bay. In Kunene and Cape Frio upwelling cells, located in the northern sector of the BUS, the coastal upwelling and open-ocean upwelling frequently alternate each other, whereas they are modulated by the annual cycle and mostly in phase off Walvis Bay. Such a phase relationship appears to be strongly seasonally dependent off Lüderitz and across the southern BUS. Thus, our findings suggest this relationship is far more complex than currently thought and seems to be sensitive to climate changes with short- and far-reaching consequences for this vulnerable marine ecosystem.

Description: Enhanced Southern Ocean ventilation in recent decades has been suggested to be a relevant modulator of the observed changes in ocean heat and carbon uptake. This study focuses on the Southern Ocean midlatitude ventilation changes from the 1960s to the 2010s. A global 1/4° configuration of the NEMO–Louvain-la-Neuve sea ice model, version 2 (LIM2), including the inert tracer CFC-12 (a proxy of ocean ventilation) is forced with the CORE, phase II (CORE-II), and JRA-55 driving ocean (JRA55-do) atmospheric reanalyses. Sensitivity experiments, where the variability of wind stress and/or the buoyancy forcing is suppressed on interannual time scales, are used to unravel the mechanisms driving ventilation changes. Ventilation changes are estimated by comparing CFC-12 interior inventories among the different experiments. All simulations suggest a multidecadal fluctuation of Southern Ocean ventilation, with a decrease until the 1980s–90s and a subsequent increase. This evolution is related to the atmospheric forcing and is caused by the (often counteracting) effects of wind stress and buoyancy forcing. Until the 1980s, increased buoyancy gains caused the ventilation decrease, whereas the subsequent ventilation increase was driven by strengthened wind stress causing deeper mixed layers and a stronger meridional overturning circulation. Wind stress emerges as the main driver of ventilation changes, even though buoyancy forcing modulates its trend and decadal variability. The three Southern Ocean basins take up CFC-12 in distinct density intervals but overall respond similarly to the atmospheric forcing. This study suggests that Southern Ocean ventilation is expected to increase as long as the effect of increasing Southern Hemisphere wind stress overwhelms that of increased stratification.

Description: The cold-water coral Lophelia pertusa builds up bioherms that sustain high biodiversity in the deep ocean worldwide. Photographic monitoring of the polyp activity represents a helpful tool to characterize the health status of the corals and to assess anthropogenic impacts on the microhabitat. Discriminating active polyps from skeletons of white Lophelia pertusa is usually time-consuming and error-prone due to their similarity in color in common RGB camera footage. Acquisition of finer resolved spectral information might increase the contrast between the segments of polyps and skeletons, and therefore could support automated classification and accurate activity estimation of polyps. For recording the needed footage, underwater multispectral imaging systems can be used, but they are often expensive and bulky. Here we present results of a new, light-weight, compact and low-cost deep-sea tunable LED-based underwater multispectral imaging system (TuLUMIS) with eight spectral channels. A brunch of healthy white Lophelia pertusa was observed under controlled conditions in a laboratory tank. Spectral reflectance signatures were extracted from pixels of polyps and skeletons of the observed coral. Results showed that the polyps can be better distinguished from the skeleton by analysis of the eight-dimensional spectral reflectance signatures compared to three-channel RGB data. During a 72-hour monitoring of the coral with a half-hour temporal resolution in the lab, the polyp activity was estimated based on the results of the multispectral pixel classification using a support vector machine (SVM) approach. The computational estimated polyp activity was consistent with that of the manual annotation, which yielded a correlation coefficient of 0.957.

Description: The El Niño-Southern Oscillation (ENSO) is the dominant mode of interannual climate variability on the planet, with far-reaching global impacts. It is therefore key to evaluate ENSO simulations in state-of-the-art numerical models used to study past, present and future climate. Recently, the Pacific Region Panel of the International Climate and Ocean - Variability, Predictability, and Change (CLIVAR) Project, as a part of the World Climate Research Programme (WCRP), led a community-wide effort to evaluate the simulation of ENSO variability, teleconnections and processes in climate models. The new CLIVAR 2020 ENSO metrics package enables model diagnosis, comparison, and evaluation to (1) highlight aspects that need improvement; (2) monitor progress across model generations; (3) help in selecting models that are well suited for particular analyses; (4) reveal links between various model biases, illuminating the impacts of those biases on ENSO and its sensitivity to climate change; and to (5) advance ENSO literacy. By interfacing with existing model evaluation tools, the ENSO metrics package enables rapid analysis of multi-petabyte databases of simulations, such as those generated by the Coupled Model Intercomparison Project phases 5 (CMIP5) and 6 (CMIP6). The CMIP6 models are found to significantly outperform those from CMIP5 for 8 out of 24 ENSO-relevant metrics, with most CMIP6 models showing improved tropical Pacific seasonality and ENSO teleconnections. Only one ENSO metric is significantly degraded in CMIP6, namely the coupling between the ocean surface and subsurface temperature anomalies, while the majority of metrics remain unchanged.

Description: This study demonstrates that the generalization that strong anomalous equatorial Pacific westerly (easterly) winds during El Niño (La Niña) events displays strong adjusted warm water volume (WWV) discharges (recharges) is often incorrect. Using ocean model simulations, we categorize the oceanic adjusted responses to strong anomalous equatorial winds into two categories: (i) transitioning (consistent with the above generalization); and (ii) neutral adjusted responses (with negligible WWV re- and discharge) During the 1980-2016 period only 47% of strong anomalous equatorial winds are followed by transitioning adjusted responses, while the remaining are followed by neutral adjusted responses. Moreover, 55% (only 30%) of the strongest winds lead to transitioning adjusted responses during the pre-2000 (post-2000) period in agreement with the previously reported post-2000 decline of WWV lead time to El Niño-Southern Oscillation (ENSO) events. The prominent neutral adjusted WWV response is shown to be largely excited by anomalous wind stress forcing with a weaker curl (on average consistent with a higher ratio of off-equatorial to equatorial wind events) and weaker Rossby wave projection than the transitioning adjusted response. We also identify a prominent ENSO phase asymmetry where strong anomalous equatorial westerly winds (i.e., El Niño events) are roughly 1.6 times more likely to strongly discharge WWV than strong anomalous equatorial easterly winds (i.e., La Niña events) are to strongly recharge WWV. This ENSO phase asymmetry may be added to the list of mechanisms proposed to explain why El Niño events have a stronger tendency to be followed by La Niña events than vice versa.

Description: The extratropical effect of the quasi-biennial oscillation (QBO), known as the Holton-Tan effect, is manifest as a weaker, warmer winter Arctic polar vortex during the east QBO phase. While previous studies have shown that the extratropical QBO signal is caused by the modified propagation of planetary waves in the stratosphere, the mechanism dominating the onset and seasonal development of the Holton-Tan effects remains unclear. Here, the governing wave-mean flow dynamics of the early winter extratropical QBO signal onset and its reversibility is investigated on a synoptic time scale with a finite-amplitude diagnostic using reanalysis and a chemistry-climate model. The extratropical QBO signal onset in October is found to primarily result from modulated stratospheric life cycles of wave pulses entering the stratosphere from the troposphere, rather than from a modulation of their tropospheric wave source. A comprehensive analysis of the wave activity budget during fall, when the stratospheric winter polar vortex starts forming and waves start propagating up into the stratosphere, shows significant differences. During the east QBO phase, the deceleration of the mid-high-latitude stratospheric zonal-mean jet by the upward-propagating wave pulses is less reversible, due to stronger dissipation processes, while during the west phase, a more reversible deceleration of the main polar vortex is found owing to the waves being dissipated at lower latitudes, accompanied by a weak but different response of the tropospheric subtropical jet. From this synoptic wave-event viewpoint, the early season onset of the Holton-Tan effect results from the cumulative effect of the QBO dependent wave-induced deceleration during the life cycle of individual upward wave pulses.

Description: The representation of tropical precipitation is evaluated across three generations of models participating in phases 3, 5, and 6 of the Coupled Model Intercomparison Project (CMIP). Compared to state-of-the-art observations, improvements in tropical precipitation in the CMIP6 models are identified for some metrics, but we find no general improvement in tropical precipitation on different temporal and spatial scales. Our results indicate overall little changes across the CMIP phases for the summer monsoons, the double-ITCZ bias, and the diurnal cycle of tropical precipitation. We find a reduced amount of drizzle events in CMIP6, but tropical precipitation occurs still too frequently. Continuous improvements across the CMIP phases are identified for the number of consecutive dry days, for the representation of modes of variability, namely, the Madden–Julian oscillation and El Niño–Southern Oscillation, and for the trends in dry months in the twentieth century. The observed positive trend in extreme wet months is, however, not captured by any of the CMIP phases, which simulate negative trends for extremely wet months in the twentieth century. The regional biases are larger than a climate change signal one hopes to use the models to identify. Given the pace of climate change as compared to the pace of model improvements to simulate tropical precipitation, we question the past strategy of the development of the present class of global climate models as the mainstay of the scientific response to climate change. We suggest the exploration of alternative approaches such as high-resolution storm-resolving models that can offer better prospects to inform us about how tropical precipitation might change with anthropogenic warming.

Description: Weather and climate variations on subseasonal to decadal time scales can have enormous social, economic, and environmental impacts, making skillful predictions on these time scales a valuable tool for decision-makers. As such, there is a growing interest in the scientific, operational, and applications communities in developing forecasts to improve our foreknowledge of extreme events. On subseasonal to seasonal (S2S) time scales, these include high-impact meteorological events such as tropical cyclones, extratropical storms, floods, droughts, and heat and cold waves. On seasonal to decadal (S2D) time scales, while the focus broadly remains similar (e.g., on precipitation, surface and upper-ocean temperatures, and their effects on the probabilities of high-impact meteorological events), understanding the roles of internal variability and externally forced variability such as anthropogenic warming in forecasts also becomes important. The S2S and S2D communities share common scientific and technical challenges. These include forecast initialization and ensemble generation; initialization shock and drift; understanding the onset of model systematic errors; bias correction, calibration, and forecast quality assessment; model resolution; atmosphere-ocean coupling; sources and expectations for predictability; and linking research, operational forecasting, and end-user needs. In September 2018 a coordinated pair of international conferences, framed by the above challenges, was organized jointly by the World Climate Research Programme (WCRP) and the World Weather Research Programme (WWRP). These conferences surveyed the state of S2S and S2D prediction, ongoing research, and future needs, providing an ideal basis for synthesizing current and emerging developments in these areas that promise to enhance future operational services. This article provides such a synthesis.

Description: There is a zonally oriented teleconnection pattern over the high-latitude Eurasian continent, which is maintained through baroclinic energy conversion. In this study, we investigate the unique features of the maintenance mechanism of this teleconnection. It is found that the baroclinic energy conversion is most efficient in both the mid-troposphere and the lower troposphere, and that the baroclinic energy conversion in the lower troposphere is comparable to that in the mid-troposphere. Further results indicate that the basic state plays a crucial role in the baroclinic energy conversion. For both the mid and lower troposphere, the atmospheric stability is low and the Coriolis parameter is large over high-latitude Eurasia, favoring strong baroclinic energy conversion. Particularly, in the lower troposphere, the atmospheric stability exhibits a clear land-sea contrast, favoring baroclinic energy conversion over the continents rather than the oceans. Furthermore, in the lower troposphere, the in-phase configuration of the meridional wind and temperature anomalies, which results from the strong meridional gradient of mean temperature around the north edge of the Eurasian continent, also significantly contributes to baroclinic energy conversion. This study highlights the role of the basic state of temperature rather than zonal wind in maintaining the high-latitude teleconnection through baroclinic energy conversion.

Description: Subtropical dry zones, located in the Hadley cells' subsidence regions, strongly influence regional climate as well as outgoing longwave radiation. Changes in these dry zones could have significant impact on surface climate as well as on the atmospheric energy budget. This study investigates the behavior of upper-tropospheric dry zones in a changing climate, using the variable upper-tropospheric humidity (UTH), calculated from climate model experiment output as well as from radiances measured with satellite-based sensors. The global UTH distribution shows that dry zones form a belt in the subtropical winter hemisphere. In the summer hemisphere they concentrate over the eastern ocean basins, where the descent regions of the subtropical anticyclones are located. Recent studies with model and satellite data have found tendencies of increasing dryness at the poleward edges of the subtropical subsidence zones. However, UTH calculated from climate simulations with 25 models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) shows these tendencies only for parts of the winter-hemispheric dry belts. In the summer hemisphere, even though differences exist between the simulations, UTH is increasing in most dry zones, particularly in the South and North Pacific Ocean. None of the summer dry zones is expanding in these simulations. Upper-tropospheric dry zones estimated from observational data do not show any robust signs of change since 1979. Apart from a weak drying tendency at the poleward edge of the southern winter-hemispheric dry belt in infrared measurements, nothing indicates that the subtropical dry belts have expanded poleward.

Description: It might be impossible to truly fathom the magnitude of the threat that global-mean sea level rise poses. However, conceptualizing the scale of the solutions required to protect ourselves against global-mean sea level rise, aids in our ability to acknowledge and understand the threat that sea level rise poses. On these grounds, we here discuss a means to protect over 25 million people and important economical regions in northern Europe against sea level rise. We propose the construction of a Northern European Enclosure Dam (NEED) that stretches between France, the United Kingdom and Norway. NEED may seem an overwhelming and unrealistic solution at first. However, our preliminary study suggests that NEED is potentially favorable financially, but also in scale, impacts and challenges compared to that of alternative solutions, such as (managed) migrations and that of country-by-country protection efforts. The mere realization that a solution as considerable as NEED might be a viable and cost-effective protection measure is illustrative of the extraordinary global threat of global-mean sea level rise that we are facing. As such, the concept of constructing NEED showcases the extent of protection efforts that are required if mitigation efforts fail to limit sea level rise.

Description: The comparison of equivalent neutral winds obtained from (a) four WHOI buoys in the subtropics and (b) scatterometer estimates at those locations reveals a root-mean-square (RMS) difference of 0.56-0.76 m/s. To investigate this RMS difference, different buoy wind error sources were examined. These buoys are particularly well suited to examine two important sources of buoy wind errors because: (1) redundant anemometers and a comparison with numerical flow simulations allow us to quantitatively assess flow distortion errors, and (2) one-minute sampling at the buoys allows us to examine the sensitivity of buoy temporal sampling/averaging in the buoy-scatterometer comparisons. The inter-anemometer difference varies as a function of wind direction relative to the buoy wind vane and is consistent with the effects of flow distortion expected based on numerical flow simulations. Comparison between the anemometers and scatterometer winds supports the interpretation that the inter-anemometer disagreement, which can be up to 5% of the wind speed, is due to flow distortion. These insights motivate an empirical correction to the individual anemometer records and subsequent comparison with scatterometer estimates show good agreement.

Description: The North Atlantic (NA) basin-averaged sea surface temperature (NASST) is often used as an index to study climate variability in the NA sector. However, there is still some debate on what drives it. Based on observations and climate models, an analysis of the different influences on the NASST index and its low-pass filtered version, the Atlantic multidecadal oscillation (AMO) index, is provided. In particular, the relationships of the two indices with some of its mechanistic drivers including the Atlantic meridional overturning circulation (AMOC) are investigated. In observations, the NASST index accounts for significant SST variability over the tropical and subpolar NA. The NASST index is shown to lump together SST variability originating from different mechanisms operating on different time scales. The AMO index emphasizes the subpolar SST variability. In the climate models, the SST-anomaly pattern associated with the NASST index is similar. The AMO index, however, only represents pronounced SST variability over the extratropical NA, and this variability is significantly linked to the AMOC. There is a sensitivity of this linkage to the cold NA SST bias observed in many climate models. Models suffering from a large cold bias exhibit a relatively weak linkage between the AMOC and AMO and vice versa. Finally, the basin-averaged SST in its unfiltered form, which has been used to question a strong influence of ocean dynamics on NA SST variability, mixes together multiple types of variability occurring on different time scales and therefore underemphasizes the role of ocean dynamics in the multidecadal variability of NA SSTs.

Description: MILAN was a multidisciplinary, international study examining how the diel variability of sea-surface microlayer biogeochemical properties potentially impacts ocean-atmosphere interaction, in order to improve our understanding of this globally important process. The sea-surface microlayer (SML) at the air-sea interface is 〈 1 mm deep but it is physically, chemically and biologically distinct from the underlying water and the atmosphere above. Wind-driven turbulence and solar radiation are important drivers of SML physical and biogeochemical properties. Given that the SML is involved in all ocean-atmosphere exchanges of mass and energy, its response to solar radiation, especially in relation to how it regulates the air-sea exchange of climate-relevant gases and aerosols, is surprisingly poorly characterised. MILAN (sea-surface MIcroLAyer at Night) was an international, multidisciplinary campaign designed to specifically address this issue. In spring 2017, we deployed diverse sampling platforms (research vessels, radio-controlled catamaran, free-drifting buoy) to study full diel cycles in the coastal North Sea SML and in underlying water, and installed a land-based aerosol sampler. We also carried out concurrent ex situ experiments using several microsensors, a laboratory gas exchange tank, a solar simulator, and a sea spray simulation chamber. In this paper we outline the diversity of approaches employed and some initial results obtained during MILAN. Our observations of diel SML variability, e.g. the influence of changing solar radiation on the quantity and quality of organic material, and diel changes in wind intensity primarily forcing air-sea CO2 exchange, underline the value and the need of multidisciplinary campaigns for integrating SML complexity into the context of air-sea interaction.

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: The structure, transport, and seasonal variability of the West Greenland boundary current system near Cape Farewell are investigated using a high-resolution mooring array deployed from 2014 to 2018. The boundary current system is comprised of three components: the West Greenland Coastal Current, which advects cold and fresh Upper Polar Water (UPW); the West Greenland Current, which transports warm and salty Irminger Water (IW) along the upper slope and UPW at the surface; and the Deep Western Boundary Current, which advects dense overflow waters. Labrador Sea Water (LSW) is prevalent at the seaward side of the array within an offshore recirculation gyre and at the base of the West Greenland Current. The 4-yr mean transport of the full boundary current system is 31.1 ± 7.4 Sv (1 Sv ≡ 106 m3 s-1), with no clear seasonal signal. However, the individual water mass components exhibit seasonal cycles in hydrographic properties and transport. LSW penetrates the boundary current locally, through entrainment/mixing from the adjacent re-circulation gyre, and also enters the current upstream in the Irminger Sea. IW is modified through air–sea interaction during winter along the length of its trajectory around the Irminger Sea, which converts some of the water to LSW. This, together with the seasonal increase in LSW entering the current, results in an anticorrelation in transport between these two water masses. The seasonality in UPW transport can be explained by remote wind forcing and subsequent adjustment via coastal trapped waves. Our results provide the first quantitatively robust observational description of the boundary current in the eastern Labrador Sea.

Description: State of the climate in 2019

Description: The decadal-scale variability in the tropical Pacific has been analyzed herein by means of observations and numerical model simulations. The two leading modes of the sea surface temperature (SST) variability in the central western Pacific are a decadal mode with a period of about 10 yr and the ENSO mode with a dominant period of about 4 yr. The SST anomaly pattern of the decadal mode is ENSO like. The decadal mode, however, explains most variance in the western equatorial Pacific and off the equator. A simulation with an ocean general circulation model (OGCM) forced by reanalysis data is used to explore the origin of the decadal mode. It is found that the variability of the shallow subtropical–tropical overturning cells is an important factor in driving the decadal mode. This is supported by results from a multicentury integration with a coupled ocean–atmosphere general circulation model (CGCM) that realistically simulates tropical Pacific decadal variability. Finally, the sensitivity of the shallow subtropical–tropical overturning cells to greenhouse warming is discussed by analyzing the results of a scenario integration with the same CGCM.

Description: The current system east of the Grand Banks was intensely observed by World Ocean Circulation Experiment (WOCE) array ACM-6 during 1993–95 with eight moorings, reaching about 500 km out from the shelf edge and covering the water column from about 400-m depth to the bottom. More recently, a reduced array by the Institut für Meerskunde (IfM) at Kiel, Germany, of four moorings was deployed during 1999–2001, focusing on the deep-water flow near the western continental slope. Both sets of moored time series, each about 22 months long, are combined here for a mean current boundary section, and both arrays are analyzed for the variability of currents and transports. A mean hydrographic section is derived from seven ship surveys and is used for geostrophic upper-layer extrapolation and isopycnal subdivision of the mean transports into deep-water classes. The offshore part of the combined section is dominated by the deep-reaching North Atlantic Current (NAC) with currents still at 10 cm s−1 near the bottom and a total northward transport of about 140 Sv (Sv ≡ 106 m3 s−1), with the details depending on the method of surface extrapolation used. The mean flow along the western boundary was southward with the section-mean North Atlantic Deep Water outflow determined to be 12 Sv below the σθ = 27.74 kg m−3 isopycnal. However, east of the deep western boundary current (DWBC), the deep NAC carries a transport of 51 Sv northward below σθ = 27.74 kg m−3, resulting in a large net northward flow in the western part of the basin. From watermass signatures it is concluded that the deep NAC is not a direct recirculation of DWBC water masses. Transport time series for the DWBC variability are derived for both arrays. The variance is concentrated in the period range from 2 weeks to 2 months, but there are also variations at interannual and longer periods, with much of the DWBC variability being related to fluctuations and meandering of the NAC. A significant annual cycle is not recognizable in the combined current and transport time series of both arrays. The moored array results are compared with other evidence on deep outflow and recirculation, including recent models of different types and complexity.

Description: Simulations and seasonal forecasts of tropical Pacific SST and subsurface fields that are based on the global Consortium for Estimating the Circulation and Climate of the Ocean (ECCO) ocean-state estimation procedure are investigated. As compared to similar results from a traditional ENSO simulation and forecast procedure, the hindcast of the constrained ocean state is significantly closer to observed surface and subsurface conditions. The skill of the 12-month lead SST forecast in the equatorial Pacific is comparable in both approaches. The optimization appears to have better skill in the SST anomaly correlations, suggesting that the initial ocean conditions and forcing corrections calculated by the ocean-state estimation do have a positive impact on the predictive skill. However, the optimized forecast skill is currently limited by the low quality of the statistical atmosphere. Progress is expected from optimizing a coupled model over a longer time interval with the coupling statistics being part of the control vector.

Description: In the eastern South Pacific Ocean, at a depth of about 200 m, a salinity minimum is found. This minimum is associated with a particular water mass, the “Shallow Salinity Minimum Water” (SSMW). SSMW outcrops in a fresh tongue (Smin) centered at about 45°S. The Smin appears to emanate from the eastern boundary, against the mean flow. The watermass transformation that creates SSMW and Smin is investigated here. The Smin and SSMW are transformed from saltier and warmer waters originating from the western South Pacific. The freshening and cooling occur when the water is advected eastward at the poleward side of the subtropical gyre. Sources of freshening and cooling are air–sea exchange and advection of water from south of the subtropical gyre. A freshwater and heat budget for the mixed layer reveals that both sources equally contribute to the watermass transformation in the mixed layer. The freshened and cooled mixed layer water is subducted into the gyre interior along the southern rim of the subtropical gyre. Subduction into the zonal flow restricts the transformation of interior properties to diffusion only. A simple advection/diffusion balance reveals diffusion coefficients of order 2000 m2 s−1. The tongue shape of the Smin is explained from a dynamical viewpoint because no relation to a positive precipitation–evaporation balance was found. Freshest Smin values are found to coincide with slowest eastward mixed layer flow that accumulates the largest amounts of freshwater in the mixed layer and creates the fresh tongue at the sea surface. Although the SSMW is the densest and freshest mode of water subducted along the South American coast, the freshening and cooling in the South Pacific affect a whole range of densities (25.0–26.8 kg m−3). The transformed water turns northward with the gyre circulation and contributes to the hydrographic structure of the gyre farther north. Because the South Pacific provides most of the source waters that upwell along the equatorial Pacific, variability in South Pacific hydrography may influence equatorial Pacific hydrography. Because one-half of the transformation is found to be controlled through Ekman transport, variability in wind forcing at the southern rim of the subtropical gyre may be a source for variability of the equatorial Pacific.

Description: An analytical model is developed to study the tidally induced mean circulation in the frontal zone. Four distinct forcing mechanisms are identified, which result in the generation of the counterclockwise Bernoulli cell, the clockwise Ekman cell, the clockwise frontal cell, and the Stokes drift (facing in the direction with the shallow water to the left). The decomposition of the cross-frontal circulation provides a dynamical framework for interpreting and understanding its complex structure. To illustrate the underlying physics, three model configurations are considered pertaining to a homogenous ocean and winter and summer fronts. For a homogeneous ocean, the circulation is dominated by three cells; for the winter front, the offshore Bernoulli cell is strengthened; and for the summer front, two counterrotating cells are found in the vertical direction, associated with the two branches of the front. The dependence of the cell structure on the Ekman, Burger, and other dimensionless numbers is examined.

Description: A multi-model ensemble-based system for seasonal-to-interannual prediction has been developed in a joint European project known as DEMETER (Development of a European Multimodel Ensemble Prediction System for Seasonal to Interannual Prediction). The DEMETER system comprises seven global atmosphere–ocean coupled models, each running from an ensemble of initial conditions. Comprehensive hindcast evaluation demonstrates the enhanced reliability and skill of the multimodel ensemble over a more conventional single-model ensemble approach. In addition, innovative examples of the application of seasonal ensemble forecasts in malaria and crop yield prediction are discussed. The strategy followed in DEMETER deals with important problems such as communication across disciplines, downscaling of climate simulations, and use of probabilistic forecast information in the applications sector, illustrating the economic value of seasonal-to-interannual prediction for society as a whole.

Description: This study focuses on an important aspect of air–sea interaction in models, namely, large-scale, spurious heat fluxes due to false pathways of the Gulf Stream and North Atlantic Current (NAC) in the “storm formation region” south and east of Newfoundland. Although high-resolution eddy-resolving models show some improvement in this respect, results are sensitive to poorly understood, subgrid-scale processes for which there is currently no complete, physically based parameterization. A simple method to correct an ocean general circulation model (OGCM), acting as a practical substitute for a physically based parameterization, is explored: the recently proposed “semiprognostic method,” a technique for adiabatically adjusting flow properties of a hydrostatic OGCM. The authors show that application of the method to an eddy-permitting model of the North Atlantic Ocean yields more realistic flow patterns and watermass characteristics in the Gulf Stream and NAC regions; in particular, spurious surface heat fluxes are reduced. Four simple modifications to the method are proposed, and their benefits are demonstrated. The modifications successfully account for three drawbacks of the original method: reduced geostrophic wave speeds, damped mesoscale eddy activity, and spurious interaction with topography. It is argued that use of a corrected (eddy permitting) OGCM in a coupled modeling system for simulating present climate (as now becomes possible because of increasing computer power) should lead to a more realistic simulation in regions of strong air–sea interaction as compared with that obtained with an uncorrected model. The method is also well suited for the simulation of the uptake and transport of passive tracers, such as anthropogenic carbon dioxide or components of ecosystem models.

Description: Two major water masses dominate the deep layers in the Mariana and Caroline Basins: the Lower Circumpolar Water (LCPW), arriving from the Southern Ocean along the slopes north of the Marshall Islands, and the North Pacific Deep Water (NPDW) reaching the region from the northeastern Pacific Ocean. Hydrographic and moored observations and multibeam echosounding were performed in the East Mariana and the East Caroline Basins to detail watermass distributions and flow paths in the area. The LCPW enters the East Mariana Basin from the east. At about 13°N, however, in the southern part of the basin, a part of this water mass arrives in a southward western boundary flow along the Izu–Ogasawara–Mariana Ridge. Both hydrographic observations and moored current measurements lead to the conclusion that this water not only continues westward to the West Mariana Basin as suggested before, but also provides bottom water to the East Caroline Basin. The critical throughflow regions were identified by multibeam echosounding at the Yap Mariana Junction between the East and West Mariana Basins and at the Caroline Ridge between the East Mariana and East Caroline Basins. The throughflow is steady between the East and West Mariana Basins, whereas more variability is found at the Caroline Ridge. At both locations, throughflow fluctuations are correlated with watermass property variations suggesting layer-thickness changes. The total transport to the two neighboring basins is only about 1 Sverdrup (1Sv ≡ 106 m3 s−1) but has considerable impact on the watermass structure in these basins. Estimates are given for the diapycnal mixing that is required to balance the inflow into the East Caroline Basin. Farther above in the water column, the high-silica tongue of NPDW extends from the east to the far southwestern corner of the East Mariana Basin, with transports being mostly southward across the basin.

Description: Sea surface temperature (SST) observations in the North Atlantic indicate the existence of strong multidecadal variability with a unique spatial structure. It is shown by means of a new global climate model, which does not employ flux adjustments, that the multidecadal SST variability is closely related to variations in the North Atlantic thermohaline circulation (THC). The close correspondence between the North Atlantic SST and THC variabilities allows, in conjunction with the dynamical inertia of the THC, for the prediction of the slowly varying component of the North Atlantic climate system. It is shown additionally that past variations of the North Atlantic THC can be reconstructed from a simple North Atlantic SST index and that future, anthropogenically forced changes in the THC can be easily monitored by observing SSTs. The latter is confirmed by another state-of-the-art global climate model. Finally, the strong multidecadal variability may mask an anthropogenic signal in the North Atlantic for some decades.

Description: A systematic modular approach to investigate the respective roles of the ocean and atmosphere in setting El Niño characteristics in coupled general circulation models is presented. Several state-of-the-art coupled models sharing either the same atmosphere or the same ocean are compared. Major results include 1) the dominant role of the atmosphere model in setting El Niño characteristics (periodicity and base amplitude) and errors (regularity) and 2) the considerable improvement of simulated El Niño power spectratoward lower frequencywhen the atmosphere resolution is significantly increased. Likely reasons for such behavior are briefly discussed. It is argued that this new modular strategy represents a generic approach to identifying the source of both coupled mechanisms and model error and will provide a methodology for guiding model improvement.

Description: The Indonesian Throughflow (ITF) spreading pathways and time scales in the Indian Ocean are investigated using both observational data and two numerical tracer experiments, one being a three-dimensional Lagrangian trajectory experiment and the other a transit-time probability density function (PDF) tracer experiment, in an ocean general circulation model. The model climatology is in agreement with observations and other model results except that speeds of boundary currents are lower. Upon reaching the western boundary within the South Equatorial Current (SEC), the trajectories of the ITF tracers within the thermocline exhibit bifurcation. The Lagrangian trajectory experiment shows that at the western boundary about 38%±5% thermocline ITF water flows southward to join the Agulhas Current, consequently exiting the Indian Ocean, and the rest, about 62%±5%, flows northward to the north of SEC. In boreal summer, ITF water penetrates into the Northern Hemisphere within the Somali Current. The primary spreading pathway of the thermocline ITF water north of SEC is upwelling to the surface layer with subsequent advection southward within the surface Ekman layer toward the southern Indian Ocean subtropics. There it is subducted and advected northward in the upper thermocline to rejoin the SEC. Both the observations and the trajectory experiment suggest that the upwelling occurs mainly along the coast of Somalia during boreal summer and in the open ocean within a cyclonic gyre in the Tropics south of the equator throughout the year. All the ITF water eventually exits the Indian Ocean along the western boundary within the Mozambique Channel and the east coast of Madagascar and, farther south, the Agulhas Current region. The advective spreading time scales, represented by the elapsed time corresponding to the maximum of transit- time PDF, show that in the upper thermocline the ITF crosses the Indian Ocean, from the Makassar Strait to the east coast of the African continent, on a time scale of about 10 yr and reaches the Arabian Sea on a time scale of over 20 yr.

Description: Aspects of the sea level changes in the western Mediterranean Sea are investigated using a numerical tidal model of the Strait of Gibraltar. As a prerequisite, the performance of this model, that is, a two-dimensional, nonlinear, two-layer, boundary-fitted coordinate numerical model based on the hydrostatic approximation on an f plane, is assessed in the simulation of mean and tidal circulation of the Strait of Gibraltar. The model is forced by imposing mean interface and surface displacements as well as M2, S2, O1, and K1 tidal components along the Atlantic and Mediterranean model open boundaries. Model results are compared with observations and with results obtained from a tidal inverse model for the eastern entrance of the Strait of Gibraltar. In general, good agreement is found. A sensitivity study performed by varying different model parameters shows that the model behaves reasonably well in the simulation of the averaged circulation. The model is then used to investigate the climatological sensitivity of the simulated dynamics in the Strait of Gibraltar to changes in the density difference between Atlantic and Mediterranean waters. For this purpose, given a certain density difference between Atlantic and Mediterranean waters, the authors iteratively searched for that sea level drop between the Atlantic and the Mediterranean that fulfills the mass balance of the Mediterranean. It is found that an increase of the density difference leads to an increase of the exchange flow and to an increase of the sea level drop between the two basins. A trend in the sea level drop of O(1 cm yr−1), such as the one observed between 1994 and 1997, is explained by the model as the result of a trend of O(10−4 yr−1) in the relative density difference between the Mediterranean and Atlantic waters. The observed north–south asymmetry in this trend is also captured by the model, and it is found to arise from changes in the along-strait velocity. Results suggest that the dynamics within the Strait of Gibraltar cannot be neglected when sea level changes in the western Mediterranean basin are investigated.

Description: The bifurcation behavior of a conceptual heat–salt oscillator model is analyzed by means of numerical continuation methods. A global (homoclinic) bifurcation acts as an organizing center for the dynamics of the simplified convective model. It originates from a codimension-2 bifurcation in an extended parameter space. Comparison with earlier work by Cessi shows that the intriguing stochastic thermohaline excitability can be understood from the bifurcation structure of the model. It is argued that global bifurcations may play a crucial role in determining long-term variability of the thermohaline circulation.

Description: A method is presented to reconstruct decadal variations of the North Atlantic Oscillation (NAO). The spectral characteristics of the NAO on time scales of decades and longer are of particular interest for the understanding of North Atlantic ocean–atmosphere interactions. The reconstruction is based on a transfer model calibration that uses bandpass-filtered time series. The maximum overlap discrete wavelet transform (MODWT) is applied for decomposing the time series variance into different time scales. A total of 43 proxies, including Greenland ice cores and European tree-ring chronologies, are selected and regionally grouped providing four independent reconstructions for the period 1700–1978. The mean reconstruction agrees well with two recently published reconstructions during most of the time period. However, there are considerable differences in the earliest part before 1750. Running correlations between the reconstructions indicate that time-dependent relations exist among the different NAO reconstructions. The results suggest that the geographical distribution of proxies strongly affects the reconstruction and could explain some of the apparent discrepancies among the reconstructions recently published in literature. In the early eighteenth century, external forcing (solar, volcanic) seems to mask the NAO signature within the proxies

Description: The deep circulation and related transports of the southern Labrador Sea are determined from direct current observations from ship surveys and a moored current-meter array. The measurements covered a time span from summer 1997 to 1999 and show a well-defined deep boundary current extending approximately out to the 3300-m depth contour and weak reverse currents farther offshore. The flow has a strong barotropic component, and significant baroclinic flow is only found in the shallow Labrador Current at the shelf break and associated with a deep core of Denmark Strait Overflow Water. The total deep-water transport below σΘ = 27.74 kg m−3 was 26 ± 5 Sv (Sv ≡ 106 m3 s−1) comprising Labrador Sea Water (LSW), Gibbs Fracture Zone Water (GFZW), and Denmark Strait Overflow Water (DSOW). Intraseasonal variability of the flow and transport was high, ranging from 15 to 35 Sv, and the annual means differed by 17%. A seasonal cycle is confined to the shallow Labrador Current; in its deeper part, where the mean flow is still strong, no obvious seasonality could be detected. The transport of the interior anticyclonic recirculation was estimated from lowered acoustic Doppler current profiler stations and geostrophy, yielding about 9 Sv. Thus, the net deep-water outflow from the Labrador Sea was about 17 Sv. The baroclinic transport of GFZW and DSOW referenced to the depth of the isopycnal σΘ = 27.80 kg m−3 is only about one-third of the total transport in these layers. Longer-term variations of the total transports are not represented well by the baroclinic contribution.

Description: On seasonal time scales, ENSO prediction has become feasible in an operational framework in recent years. On decadal to multidecadal time scales, the variability of the oceanic circulation is assumed to provide a potential for climate prediction. To investigate the decadal predictability of the coupled atmosphere–ocean general circulation model (AOGCM) European Centre-Hamburg model version 5/Max Planck Institute Ocean Model (ECHAM5/MPI-OM), a 500-yr-long control integration and “perfect model” predictability experiments are analyzed. The results show that the sea surface temperatures (SSTs) of the North Atlantic, Nordic Seas, and Southern Ocean exhibit predictability on multidecadal time scales. Over the ocean, the predictability of surface air temperature (SAT) is very similar to that of SST. Over land, there is little evidence of decadal predictability of SAT except for some small maritime-influenced regions of Europe. The AOGCM produces predictable signals in lower-tropospheric temperature and precipitation over the North Atlantic, but not in sea level pressure.

Description: The huge warming of the Arctic that started in the early 1920s and lasted for almost two decades is one of the most spectacular climate events of the twentieth century. During the peak period 1930–40, the annually averaged temperature anomaly for the area 60°–90°N amounted to some 1.7°C. Whether this event is an example of an internal climate mode or is externally forced, such as by enhanced solar effects, is presently under debate. This study suggests that natural variability is a likely cause, with reduced sea ice cover being crucial for the warming. A robust sea ice–air temperature relationship was demonstrated by a set of four simulations with the atmospheric ECHAM model forced with observed SST and sea ice concentrations. An analysis of the spatial characteristics of the observed early twentieth-century surface air temperature anomaly revealed that it was associated with similar sea ice variations. Further investigation of the variability of Arctic surface temperature and sea ice cover was performed by analyzing data from a coupled ocean–atmosphere model. By analyzing climate anomalies in the model that are similar to those that occurred in the early twentieth century, it was found that the simulated temperature increase in the Arctic was related to enhanced wind-driven oceanic inflow into the Barents Sea with an associated sea ice retreat. The magnitude of the inflow is linked to the strength of westerlies into the Barents Sea. This study proposes a mechanism sustaining the enhanced westerly winds by a cyclonic atmospheric circulation in the Barents Sea region created by a strong surface heat flux over the ice-free areas. Observational data suggest a similar series of events during the early twentieth-century Arctic warming, including increasing westerly winds between Spitsbergen and Norway, reduced sea ice, and enhanced cyclonic circulation over the Barents Sea. At the same time, the North Atlantic Oscillation was weakening.

Description: Oceanic ecosystems altered by interdecadal climate variability may provide a feedback to the physical climate by phytoplankton affecting heat fluxes into the upper ocean and dimethylsulfide fluxes into the atmosphere

Description: Intrinsic oscillations of stable geophysical surface frontal currents of the unsteady, nonlinear, reduced-gravity shallow-water equations on an f plane are investigated analytically and numerically. For frictional (Rayleigh) currents characterized by linear horizontal velocity components and parabolic cross sections, the primitive equations are reduced to a set of coupled nonlinear ordinary differential equations. In the inviscid case, two periodic analytical solutions of the nonlinear problem describing 1) the inertially reversing horizontal displacement of a surface frontal current having a fixed parabolic cross section and 2) the cross-front pulsation of a coastal current emerging from a motionless surface frontal layer are presented. In a linear and in a weakly nonlinear context, analytical expressions for field oscillations and their frequency shift relative to the inertial frequency are presented. For the fully nonlinear problem, solutions referring to a surface frontal coastal current are obtained analytically and numerically. These solutions show that the currents oscillate always superinertially, the frequency and the amplitude of their oscillations depending on the magnitude of the initial disturbance and on the squared current Rossby number. In a linear framework, it is shown that disturbances superimposed on the surface frontal current are standing waves within the bounded region, the frequencies of which are inertial/superinertial for the first mode/higher modes. In the same frame, a zeroth mode, which could be interpreted as the superposition of an inertial wave on a background vorticity field, would formally yield subinertial frequencies. For surface frontal currents affected by Rayleigh friction, it is shown that the magnitude of the mean current decays according to a power law and that the oscillations decay faster, because this decay follows an exponential law. Implications of the intrinsic oscillations and of their rapid dissipation for the near-inertial motion in an active ambient ocean are discussed.

Description: Aspects of the dynamics of warm-core eddies evolving in a deep ocean are investigated using the results of laboratory experiments and numerical simulations. The vortices, produced experimentally in a system brought to solid body rotation by rapidly lifting a bottomless cylinder containing freshwater immersed in a salty ambient fluid, show clearly the presence of inertial oscillations: deepenings and contractions, shoalings and expansions, alternate during an exact inertial period. These pulsations, though predicted analytically and simulated numerically, had never been measured before for surface eddies having aspect ratios, as well as Rossby and Burger numbers, typical of geophysical warm-core eddies. The spatial structure of the vortex radial and tangential velocity components is analyzed using the experimental results and numerical simulations carried out by means of a layered, nonlinear, reduced-gravity frontal model. It is found that, while the dependence of the vortex radial velocity on the vortex radius evolves toward linearity as time elapses, different spatial structures seem to be possible for the vortex tangential velocity dependence. This behavior, which strongly differs from the “pulson” dynamics, is instead consistent with recently found analytical solutions of the nonlinear, reduced-gravity shallow-water equations describing the dynamics of warm-core eddies on an f plane.

Description: Two large-scale free-drifting isobaric-floats experiments, “SOFARGOS”/Marine Science and Technology Programme, phase 2 (MAST2) and Mass Transfer and Ecosystem Response (MATER)/MAST3, undertaken in 1994–95 in the northwestern Mediterranean Sea and in 1997–98 in the Algerian Basin, respectively, have revealed for the first time that Western Mediterranean Deep Water, newly formed by deep convection in the Gulf of Lion (the so-called Medoc site), can be advected several hundreds of kilometers away from the formation area by anticyclonic submesoscale coherent vortices (SCVs). This behavior implies that SCVs participate actively in the large-scale thermohaline circulation and deep ventilation of the western Mediterranean Sea. These SCVs are characterized by small radius (5 km), very low potential vorticity, high aspect ratio (0.1), and extended lifetime (〉0.5 yr).

Description: Bulk properties of the Denmark Strait overflow (DSO) plume observed in velocity and hydrography surveys undertaken in 1997 and 1998 are described. Despite the presence of considerable short-term variability, it is found that the pathway and evolution of the plume density anomaly are remarkably steady. Bottom stress measurements show that the pathway of the plume core matches well with a rate of descent controlled by friction. The estimated entrainment rate diagnosed from the rate of plume dilution with distance shows a marked increase in entrainment at approximately 125 km from the sill, leading to a net dilution consistent with previous reports of a doubling of overflow transport measured by current meter arrays. The entrainment rate increase is likely related to the increased topographic slopes in the region, compounded by a decrease in interface stratification as the plume is diluted and enters a denser background.

Description: A new shipboard current profiler, a 75-kHz ocean surveyor, was operationally used during two research cruises in the tropical Atlantic and the subpolar North Atlantic, respectively. Here, a report is presented on the first experience with this instrument in two very different current regimes, in the Tropics with large vertical shears, and in the subpolar regime with mainly barotropic flow. The ocean surveyor continuously measured currents in the upper ocean from near the surface to about 500–700-m depth. The measurement range showed a dependence on the regional and temporal variations of scattering particles and on the intensity of swell and wind waves. Statistical comparisons are performed with on-station lowered acoustic Doppler current profiler (LADCP) profiles and underway measurements by classic shipboard acoustic Doppler current profiler (ADCP) measurements. Accuracy estimates for hourly averaged ocean surveyor currents result in errors of about 1 cm s–1 for on-station data and of 2–4 cm s–1 for underway measurements, depending on the regional abundance of scatterers and on the weather conditions encountered.

Description: A new mechanism is proposed that explains two key features of the observed El Niño–Southern Oscillation (ENSO) phenomenon—its irregularity and decadal amplitude changes. Using a low-order ENSO model, the authors show that the nonlinearities in the tropical heat budget can lead to bursting behavior characterized by decadal occurrences of strong El Niño events. La Niña events are not affected, a feature that is also seen in ENSO observations. One key result of this analysis is that decadal variability in the Tropics can be generated without invoking extratropical processes or stochastic forcing. The El Niño bursting behavior simulated by the low-order ENSO model can be understood in terms of the concept of homoclinic and heteroclinic connections. It is shown that this new model for ENSO amplitude modulations and irregularity, although difficult to prove, might explain some features of ENSO dynamics seen in more complex climate models and the observations.

Description: Processes that influence the volume and heat transport across the Greenland–Scotland Ridge system are investigated in a numerical model with ° horizontal resolution. The focus is on the sensitivity of cross-ridge transports and the reaction of the subpolar North Atlantic Ocean circulation to changes in wind stress and buoyancy forcing on seasonal to interannual timescales. A general relation between changes in wind stress or cross-ridge density contrasts and the overturning transport of Greenland–Iceland–Norwegian Seas source water is established from a series of idealized experiments. The relation is used subsequently to interpret changes in an experiment over the years 1992–97 with realistic forcing. On seasonal and interannual timescales there is a clear correlation between heat flux and wind stress curl variability. The realistic model suggests a steady decrease in the strength of the cyclonic subpolar gyre of the North Atlantic with a corresponding decrease in heat transport during the 1990s

Description: The role of mean and stochastic freshwater forcing on the generation of millennial-scale climate variability in the North Atlantic is studied using a low-order coupled atmosphere–ocean–sea ice model. It is shown that millennial-scale oscillations can be excited stochastically, when the North Atlantic Ocean is fresh enough. This finding is used in order to interpret the aftermath of massive iceberg surges (Heinrich events) in the glacial North Atlantic, which are characterized by an excitation of Dansgaard–Oeschger events. Based on model results, it is hypothesized that Heinrich events trigger Dansgaard–Oeschger cycles and that furthermore the occurrence of Heinrich events is dependent on the accumulated climatic effect of a series of Dansgaard–Oeschger events. This scenario leads to a coupled ocean–ice sheet oscillation that shares many similarities with the Bond cycle. Further sensitivity experiments reveal that the timescale of the oscillations can be decomposed into stochastic, linear, and nonlinear deterministic components. A schematic bifurcation diagram is used to compare theoretical results with paleoclimatic data.

Description: Synoptic-scale variability in the air–sea turbulent fluxes in the areas of midlatitudinal western boundary currents is analyzed. In the Gulf Stream area, ocean–atmosphere fluxes on synoptic time- and space scales are clearly coordinated with the propagating synoptic-scale atmospheric transients. The statistical analysis of 6-hourly resolution sea level pressure and surface turbulent fluxes from the NCEP–NCAR reanalysis for the period from 1948 to 2000 in the area of strong sea surface temperature gradients in the Gulf Stream gives strong proof for the association between the propagating cyclones and synoptic patterns of surface turbulent fluxes. It is shown that sea–air interaction in this area is controlled by the sharpness of surface temperature gradients in the ocean and by the intensity of the advection of the air masses in different parts of cyclones during the cold-air and warm-air outbreaks. A simple parameter based on the joint consideration of the characteristics of sea surface temperature and sea level pressure fields is used to characterize the synoptic variability of air–sea turbulent fluxes. The effectiveness of the relationship between surface temperature and surface pressure on one side and air–sea flux anomalies on the other vary from year to year in phase with variability in the frequencies of deep atmospheric cyclones in the Gulf Stream area. The limits of applicability of the approach, its sensitivity to higher-resolution sea surface temperature data, and the possibility of its further applications are discussed.

Description: Recent observational studies have shown that the centers of action of interannual variability of the North Atlantic Oscillation (NAO) were located farther eastward during winters of the period 1978–97 compared to previous decades of the twentieth century. In this study, which focuses on the winter season (December–March), new diagnostics characterizing this shift are presented. Further, the importance of this shift for NAO-related interannual climate variability in the North Atlantic region is discussed. It is shown that an NAO-related eastward shift in variability can be found for a wide range of different parameters like the number of deep cyclones, near-surface air temperature, and turbulent surface heat flux throughout the North Atlantic region. By using a near-surface air temperature dataset that is homogenous with respect to the kind of observations used, it is shown that the eastward shift is not an artifact of changes in observational practices that took place around the late 1970s. Finally, an EOF-based Monte Carlo test is developed to quantify the probability of changes in the spatial structure of interannual NAO variability for a relatively short (20 yr) time series given multivariate “white noise.” It is estimated that the likelihood for differences in the spatial structure of the NAO between two independent 20-yr periods, which are similar (as measured by the angle and pattern correlation between two NAO patterns) to the observed differences, to occur just by chance is about 18%. From the above results it is argued that care has to be taken when conclusions about long-term properties of NAO-related climate variability are being drawn from relatively short recent observational data (e.g., 1978–97).

Description: A new type of ocean general circulation model with simplified physics is described and tested for various simple wind-driven circulation problems. The model consists of the vorticity balance of the depth-averaged flow and a hierarchy of equations for “vertical moments” of density and baroclinic velocity. The first vertical density moment is the (vertically integrated) potential energy, which is used to describe the predominant link between the barotropic and the baroclinic oceanic flow in the presence of sloping topography. Tendency equations for the vertical moments of density and baroclinic velocity and an appropriate truncation of the coupled hierarchy of moments are derived that, together with the barotropic vorticity balance, yield a closed set of equations describing the barotropic–baroclinic interaction (BARBI) model of the oceanic circulation. Idealized companion experiments with a numerical implementation of the BARBI model and a primitive equation model indicate that wave propagation properties and baroclinic adjustments are correctly represented in BARBI in midlatitudes as well as in equatorial latitudes. Furthermore, a set of experiments with a realistic application to the Atlantic/Southern Ocean system reproduces important aspects that have been previously reported by studies of gyre circulations and circumpolar currents using full primitive equation models

Description: A simple stochastic atmosphere model is coupled to a realistic model of the North Atlantic Ocean. A north–south SST dipole, with its zero line centered along the subpolar front, influences the atmosphere model, which in turn forces the ocean model by surface fluxes related to the North Atlantic Oscillation. The coupled system exhibits a damped decadal oscillation associated with the adjustment through the ocean model to the changing surface forcing. The oscillation consists of a fast wind-driven, positive feedback of the ocean and a delayed negative feedback orchestrated by overturning circulation anomalies. The positive feedback turns out to be necessary to distinguish the coupled oscillation from that in a model without any influence from the ocean to the atmosphere. Using a novel diagnosing technique, it is possible to rule out the importance of baroclinic wave processes for determining the period of the oscillation, and to show the important role played by anomalous geostrophic advection in sustaining the oscillation.

Description: This study investigates and accounts for the influence of various ice cloud parameters on the retrieval of the surface solar radiation budget (SSRB) from reflected flux at the top of the atmosphere (TOA). The optical properties of ice clouds depend on ice crystal shape, size distribution, water content, and the vertical profiles of geometric and microphysical structure. As a result, the relationship between the SSRB and TOA-reflected flux for an ice cloud atmosphere is more complex and differs from that for water cloud and cloudless atmospheres. The sensitivities of the relationship between the SSRB and TOA-reflected flux are examined with respect to various ice cloud parameters. Uncertainties in the retrieval of the SSRB due to inadequate knowledge of various ice cloud parameters are evaluated thoroughly. The uncertainty study is concerned with both pure ice clouds and multiphase clouds (ice cloud above water cloud). According to the magnitudes of errors in the SSRB retrieval caused by different input variables, parameterized correction terms were introduced. If the input variables are known accurately, errors in the retrieval of the SSRB under a wide range of ice cloud conditions are expected to diminish substantially, to less than 10 W m−2 for 91% of the simulated ice cloud cases. In comparison, the same accuracy may be attained for only 19% of the retrievals for the same ice cloud cases using the retrieval algorithm designed for non-ice-cloud conditions.

Description: North Atlantic synoptic-scale processes are analyzed by bandpassing 6-hourly NCEP–NCAR reanalysis data (1958–98) for several synoptic ranges corresponding to ultrahigh-frequency variability (0.5–2 days), synoptic-scale variability (2–6 days), slow synoptic processes (6–12 days), and low-frequency variability (12–30 days). Climatological patterns of the intensity of synoptic processes are not collocated for different ranges of variability, especially in the lower troposphere. Intensities of synoptic processes demonstrate opposite trends between the North American coast and in the northeast Atlantic. Although north of 40°N the intensity of ultrahigh-frequency variability and synoptic-scale processes show similar interannual variability, further analysis indicates that secular changes, and decadal-scale and interannual variability in the intensities of synoptic processes may not be necessarily consistent for different synoptic timescales. Magnitudes of winter ultrahigh-frequency variability are highly correlated with the intensity of synoptic-scale processes in the 1960s and early 1970s. However, they show little agreement with each other during the last two decades, pointing to the remarkable change in atmospheric variability over the North Atlantic in late 1970s. North Atlantic ultrahigh-frequency variability in winter is highly correlated with surface temperature gradient anomalies in the Atlantic–American sector. These gradients are computed from the merged fields of SST and surface temperature over the continent. They demonstrate a dipolelike pattern associated with the North American coast on one hand, with the subpolar SST front and continental Canada on the other. High-frequency variability and its synoptic counterpart demonstrate different relationships with the North Atlantic Oscillation. Reliability of these results and their sensitivity to the filtering procedures are addressed by comparison to radiosonde data and application of alternative filters.

Description: Turbulent fluxes of momentum and sensible heat were estimated from sonic anemometer measurements gathered over the Labrador Sea during a winter cruise of the R/V Knorr. The inertial dissipation method was used to calculate turbulent fluxes of momentum. The resulting drag coefficients agree well with earlier findings. Sensible heat fluxes were computed using both cross-correlation and inertial dissipation techniques. There is good agreement between results from both methods, although there is more scatter in the correlation fluxes than the dissipation fluxes. The inertial dissipation method gives reasonable results even under conditions of high wind speeds and low air temperatures, which combined with the relatively warm sea surface temperatures lead to sensible heat fluxes of several hundred watts per square meter. Sensible heat fluxes obtained from the sonic anemometer measurements agree well with bulk turbulent fluxes according to the formulation of Isemer and Hasse.

Description: Lowered acoustic Doppler current profilers (LADCPs) have matured from an experimental instrument to an operational hydrographic tool to study ocean dynamics. The data processing, however, is still in a rather primitive state. First, a method to estimate bottom-track velocities using the standard water profile data was developed. Then inverse solutions are presented that enhance the standard data processing by adding external constraints such as bottom-referenced velocity profiles. Depending on the depth of the profile and the ADCP range the inclusion of bottom-track data can reduce the local velocity errors by a significant factor. The least squares framework also allows for simplified error analysis of the LADCP system and some of the trade-offs are discussed.

Description: Zonal transports of North Atlantic Deep Water (NADW) in the South Atlantic are determined. For this purpose the circulation of intermediate and deep water masses is established on the basis of hydrographic sections from the World Ocean Circulation Experiment (WOCE) and some pre-WOCE sections, using temperature, salinity, nutrients, and anthropogenic tracers. Multiple linear regression is applied to infer missing parameters in the bottle dataset. A linear box-inverse model is used for a set of closed boxes given by sections and continental boundaries. After performing a detailed analysis of water mass distribution, 11 layers are prescribed. Neutral density surfaces are selected as layer interfaces, thus improving the description of water mass distribution in the transition between the subtropical and subpolar latitudes. Constraints for the inverse model include integral meridional salt and phosphorus transports, overall salt and silica conservation, and transports from moored current meter observations. Inferred transport numbers for the mean meridional thermohaline overturning are given. Persistent zonal NADW transport bands are found in the western South Atlantic, in particular eastward flow of relatively new NADW between 20° and 25°S and westward flow of older NADW to the north of this latitude range. The axis of the eastward transport band corresponds to the core of property distributions in this region, suggesting Wüstian flow. Part of the eastward flow appears to cross the Mid-Atlantic Ridge at the Rio de Janeiro Fracture Zone. Results are compared qualitatively with deep float observations and results from general circulation models

Description: Fifteen profiling floats were injected into the deep boundary current off Labrador. They were ballasted to drift in the core depth of Labrador Sea Water (LSW) at 1500-m depth and were deployed in two groups during March and July/August 1997. Initially, for about three months, the floats were drifting within the boundary current, and the flow vectors were used to determine the mean horizontal structure of the Deep Labrador Current, which was found to be about 100 km wide with an average core speed of 18 cm s−1. North of Flemish Cap the boundary current encounters complicated topography around “Orphan Knoll,” and there the LSW outflow splits up into different routes. One obvious LSW path is eastward through the Charlie Gibbs Fracture Zone and another route is a narrow recirculation toward the central Labrador Sea. A surprising result was that none of the floats were able to follow the boundary current southward to the Grand Banks area and exit into the subtropics. Trajectories and temperature profiles of the eastward drifting floats indicate the importance of the North Atlantic Current for dispersing the floats, even at the level of LSW.

Description: The analysis of high-resolution oceanographic data referring to velocity measurements carried out by means of a vessel-mounted acoustic Doppler current profiler on 12 November 2000 in the equatorial Atlantic, at 44°W between 4.5° and 6°N, reveals the presence of three large-amplitude internal solitary waves superimposed on the velocity field associated with the North Equatorial Countercurrent (NECC). These waves were found in the deep ocean, more than 500 km off the continental shelf and far from regions of topographic variations. They propagated toward the north-northeast, strongly inclined with respect to the main axis of the NECC and perpendicular to the Brazilian shelf, as well as to the North Brazil Current, and were characterized by maximum horizontal velocities of about 2 m s−1 and maximum vertical velocities of about 20 cm s−1. The large magnitudes of the measured velocities indicate that the observed waves represent disturbances evolving in a strongly stratified ocean. The distance separating the waves (about 70 km) indicates that the observed features cannot be considered as elements of a single train of internal solitary waves. The waves consist, instead, of truly disconnected, pulselike intense solitary disturbances. This behavior, which strongly differs from that typically observed for trains of tidally generated internal solitary waves, indicates that different mechanisms were possibly involved in their generation and/or evolution.

Description: The ventilation of the permanent thermocline of the Southern Hemisphere gyres is quantified using climatological and synoptic observational data. Ventilation is estimated with three independent methods: the kinematic method provides subduction rates from the vertical and horizontal fluxes through the base of the mixed layer, the water age uses in situ age distribution of thermocline waters, and the annual-mean water mass formation through air–sea interaction is calculated. All three independent estimates agree within their error bars, which are admittedly large. The subduction rates are mainly controlled through their vertical and lateral components with only minor transient eddy contributions. The vertical transfer, derived from Ekman pumping, ventilates over most of the areas of the subtropical gyres, while lateral transfer occurs mainly along the Subtropical and Subantarctic Fronts, where it injects mode and intermediate waters. For the permanent thermocline the overall ventilation of the South Atlantic is about 21 Sv (Sv ≡ 106 m3 s−1). Of this, lateral transfer contributes 10 Sv, mainly in the Brazil–Malvinas confluence zone and to the northeast of Drake Passage. The effective vertical transfer at the bottom of the mixed layer is only two-thirds of the Ekman pumping due to strong northward forcing of the mixed layer itself. The Indian Ocean is ventilated at a rate of 35 Sv with equal lateral and vertical contributions. The South Pacific's overall ventilation is 44 Sv of which the lateral input contributes little more than half. West of 130°W, the South Pacific is ventilated through Ekman pumping and with only minor lateral transfer. In the east lateral transfer dominates between 10° and 20°S and along the Subantarctic Front in a narrow density range. Combining overall transports with earlier estimates for the Northern Hemisphere gives a ventilation of the World Ocean's permanent thermocline of about 160 Sv. Analysis of atmospheric reanalysis air–sea flux data reveals an overall increase in the formation of thermocline waters for all three Southern Hemisphere oceans.

Description: The quasi-decadal salinity fluctuations in the upper 300 m of the Labrador Sea are investigated by partitioning all available salinity station data since 1948 by region and bottom depth. There are major freshwater anomalies in the early 1970s (the Great Salinity Anomaly), mid-1980s, and early 1990s. These vary in amplitude throughout the region, being least on the shelf and greatest over the slope region near the Labrador Current. The Labrador Sea cannot be considered a simple conduit for freshwater anomalies originating in the East Greenland Current. There is evidence that local processes modulate the anomaly. The freshwater anomalies in the Labrador Current are approximately twice as large as those in the East Greenland Current. The Baffin Island Current flowing southward through the western Davis Strait is the only local source of freshwater with sufficient volume to account for this increase. The propagation speed, 2–3 cm s−1, of the anomaly along the Labrador Sea margin is much less than the advection speed indicating a highly damped system. The connection of the North Atlantic Oscillation (NAO) with these quasi-decadal salinity fluctuations is most obvious in the Labrador Sea interior, where increased surface buoyancy flux during positive NAO drives deep convective mixing and thus terminates the fresh surface anomalies. Less clear are the processes by which NAO-forced changes of lateral freshwater flux modulate the salinity along the margin. The authors propose a feedback mechanism where, during years of low wind speed, freshwater accumulates offshore of the slope front in the surface layer. The increased upper-layer buoyancy prohibits further mixing, and low salinities persist.

Description: In 1997, a unique hydrographic and chlorofluorocarbon (CFC: component CFC-11) dataset was obtained in the subpolar North Atlantic. To estimate the synopticity of the 1997 data, the recent temporal evolution of the CFC and Labrador Sea Water (LSW) thickness fields are examined. In the western Atlantic north of 50°N, the LSW thickness decreased considerably from 1994–97, while the mean CFC concentrations did not change much. South of 50°N and in the eastern Atlantic, the CFC concentration increased with little or no change in the LSW thickness. On shorter timescales, local anomalies due to the presence of eddies are observed, but for space scales larger than the eddies the dataset can be treated as being synoptic over the 1997 observation period. The spreading of LSW in the subpolar North Atlantic is described in detail using gridded CFC and LSW thickness fields combined with Profiling Autonomous Lagrangian Circulation Explorer (PALACE) float trajectories. The gridded fields are also used to calculate the CFC-11 inventory in the LSW from 40° to 65°N, and from 10° to 60°W. In total, 2300 ± 250 tons of CFC-11 (equivalent to 16.6 million moles) were brought into the LSW by deep convection. In 1997, 28% of the inventory was still found in the Labrador Sea west of 45°W and 31% of the inventory was located in the eastern Atlantic. The CFC inventory in the LSW was used to estimate the lower limits of LSW formation rates. At a constant formation rate, a value of 4.4–5.6 Sv (Sv ≡ 106 m3 s−1) is obtained. If the denser modes of LSW are ventilated only in periods with intense convection, the minimum formation rate of LSW in 1988–94 is 8.1–10.8 Sv, and 1.8–2.4 Sv in 1995–97

Description: Time series of hydrographic and transient tracer (3H and 3He) observations from the central Labrador Sea collected between 1991 and 1996 are presented to document the complex changes in the tracer fields as a result of variations in convective activity during the 1990s. Between 1991 and 1993, as atmospheric forcing intensified, convection penetrated to progressively increasing depths, reaching 2300 m in the winter of 1993. Over that period the potential temperature (θ)/salinity (S) properties of Labrador Sea Water stayed nearly constant as surface cooling and downward mixing of freshwater was balanced by excavating and upward mixing of the warmer and saltier Northeast Atlantic Deep Water. It is shown that the net change in heat content of the water column (150–2500 m) between 1991 and 1993 was negligible compared to the estimated mean heat loss over that period (110 W m−2), implying that the lateral convergence of heat into the central Labrador Sea nearly balances the atmospheric cooling on a surprisingly short timescale. Interestingly, the 3H–3He age of Labrador Sea Water increased during this period of intensifying convection. Starting in 1995, winters were milder and convection was restricted to the upper 800 m. Between 1994 and 1996, the evolution of 3H–3He age is similar to that of a stagnant water body. In contrast, the increase in θ and S over that period implies exchange of tracers with the boundaries via both an eddy-induced overturning circulation and along-isopycnal stirring by eddies [with an exchange coefficient of O(500 m2 s−1)]. The authors construct a freshwater budget for the Labrador Sea and quantitatively demonstrate that sea ice meltwater is the dominant cause of the large annual cycle of salinity in the Labrador Sea, both on the shelf and the interior. It is shown that the transport of freshwater by eddies into the central Labrador Sea (140 cm between March and September) can readily account for the observed seasonal freshening. Finally, the authors discuss the role of the eddy-induced overturning circulation with regard to transport and dispersal of the newly ventilated Labrador Sea Water to the boundary current system and compare its strength (2–3 Sv) to the diagnosed buoyancy-forced formation rate of Labrador Sea Water.

Description: Comparisons are made between a time series of meteorological surface layer observational data taken on board the R/V Knorr, and model analysis data from the European Centre for Medium-Range Weather Forecasting (ECMWF) and the National Centers for Environmental Prediction (NCEP). The observational data were gathered during a winter cruise of the R/V Knorr, from 6 February to 13 March 1997, as part of the Labrador Sea Deep Convection Experiment. The surface layer observations generally compare well with both model representations of the wintertime atmosphere. The biases that exist are mainly related to discrepancies in the sea surface temperature or the relative humidity of the analyses. The surface layer observations are used to generate bulk estimates of the surface momentum flux, and the surface sensible and latent heat fluxes. These are then compared with the model-generated turbulent surface fluxes. The ECMWF surface sensible and latent heat flux time series compare reasonably well, with overestimates of only 13% and 10%, respectively. In contrast, the NCEP model overestimates the bulk fluxes by 51% and 27%, respectively. The differences between the bulk estimates and those of the two models are due to different surface heat flux algorithms. It is shown that the roughness length formula used in the NCEP reanalysis project is inappropriate for moderate to high wind speeds. Its failings are acute for situations of large air–sea temperature difference and high wind speed, that is, for areas of high sensible heat fluxes such as the Labrador Sea, the Norwegian Sea, the Gulf Stream, and the Kuroshio. The new operational NCEP bulk algorithm is found to be more appropriate for such areas. It is concluded that surface turbulent flux fields from the ECMWF are within the bounds of observational uncertainty and therefore suitable for driving ocean models. This is in contrast to the surface flux fields from the NCEP reanalysis project, where the application of a more suitable algorithm to the model surface-layer meteorological data is recommended

Description: The so-called equatorial stacked jets are analyzed with ship-board observations and moored time series from the Atlantic Ocean. The features are identified and isolated by comparing vertical wavenumber spectra at the equator with those a few degrees from the equator. Mode-filtering gives clear views of the jets in meridional sections, the typical extent being ±1° in latitude. The vertical structure can be well described (explaining 82% of the variance) by N−1-stretched cosines, with a Gaussian amplitude tapering in the vertical. The stretched wavelengths are somewhat variable. Fitting jets of a fixed (stretched) wavelength to four moored sensors in the depth range 1300–1900 m, allows one to track the vertical phase of the jets with an rms error of 30°–45°. The resulting fit from a 20-month moored time series shows long periods of unchanging jet conditions and intermittent times of high variability. There is no significant vertical propagation on these timescales nor a seasonal reversal. Using a composite from many different experiments, interannual variability is visible, however. A possible mechanism for the stacked jets is inertial instability, resulting from background meridional shears at the equator. A condition is that the Ertel potential vorticity becomes zero somewhere, due to meridional asymmetries in the zonal flows. The ship-board observations show that this may be approximately fulfilled by the instantaneous zonal low-mode flows at various depths, resulting from an excess of zonal momentum south of the equator most of the time. Inertial instability should act to redistribute this zonal momentum, and our mooring data show indeed persistent northward momentum flux, but not at the depth levels expected. The momentum transport might suggest that the jets can also flux or mix other properties across the equator.

Description: The interannual variability of the tropical Indian Ocean sea surface temperature (SST) is studied with observational data and a hierarchy of coupled general circulation models (CGCMs). Special attention is given to the question whether an oscillatory dipole mode exists in the tropical Indian Ocean region with centers east and west of 80°E. Our observational analyses indicate that dipole-like variability can be explained as an oscillatory mode only in the context of ENSO (El Nino/Southern Oscillation). A dipole-like structure in the SST anomalies independent of ENSO was found also. Our series of coupled model experiments shows that ocean dynamics is not important to this type of dipole-like SST variability. It is forced by surface heat flux anomalies that are integrated by the thermal inertia ofthe oceanic mixed layer, which reddens the SST spectrum.

Description: Empirical orthogonal function (EOF) analyses (rotated or not) are widely used in climate research. In recent years there have been several studies in which EOF analyses were used to highlight potential physical mechanisms associated with climate variability. For example, several SST modes were identified such as the “Tropical Atlantic Dipole,” the “Tropical Indian Ocean Dipole,” and different SLP modes in the Northern Hemisphere winter. In this note it is emphasized that caution should be used when trying to interpret these statistically derived modes and their significance. Indeed, from a synthetic example it is shown that patterns derived from EOF analyses can be misleading at times and associated with very little climate physics.

Description: Experiments with a suite of North Atlantic general circulation models are used to examine the sources of eddy kinetic energy (EKE) in the Labrador Sea. A high-resolution model version (112°) quantitatively reproduces the observed signature. A particular feature of the EKE in the Labrador Sea is its pronounced seasonal cycle, with a maximum intensity in early winter, as already found in earlier studies based on altimeter data. In contrast to a previously advanced hypothesis, the seasonally varying eddy field is not related to a forcing by high-frequency wind variations but can be explained by a seasonally modulated instability of the West Greenland Current (WGC). The main source of EKE in the Labrador Sea is an energy transfer due to Reynolds interaction work (barotropic instability) in a confined region near Cape Desolation where the WGC adjusts to a change in the topographic slope: Geostrophic contours tend to converge upstream of Cape Desolation, such that the topographically guided WGC narrows as well and becomes barotropically unstable. The eddies spawned from the WGC instability area, dominating the EKE in the interior Labrador Sea, are predominantly anticyclonic with warm and saline cores in the upper kilometer of the water column, while the few cyclones originating as well from the instability area show a more depth-independent structure. Companion experiments with a ⅓° model exhibit the strength of the WGC, influenced by either changes in the wind stress or heat flux forcing, as a leading factor determining seasonal to interannual changes of EKE in the Labrador Sea

Description: Aspects of the decay of stable frontal warm-core eddies in the deep ocean are investigated using a new numerical layered “frontal” model that solves the nonlinear, reduced-gravity, shallow-water equations for a horizontally inhomogeneous, viscous fluid on an f plane. After a discussion on aspects of the numerical techniques implemented to allow for the eddy expansions and contractions at the sea surface, for the first time the capability of a numerical model of reproducing the evolution of analytical nonstationary frontal vortices is explored. This step is necessary, as far as different phenomena related to the dynamics of these oceanic features are to be studied numerically. In fact the comparison between numerical and analytical inviscid solutions allows for a quantification of the numerical dissipation affecting the simulated solutions. This dissipation is found to be very small in this numerical model: The simulated lifetimes are larger than those of most of the frontal eddies observed in the World Ocean. On this basis, the eddy decay due to interfacial (linear and quadratic) friction, harmonic horizontal momentum diffusion, as well as linear ambient-water entrainment is investigated. It is found that interfacial friction represents a much more efficient mechanism than horizontal diffusion and water entrainment in inducing the eddy decay as well as in damping the eddy pulsations. It is thus suggested that internal wave radiation due to vortex pulsation can represent a relevant mechanism for the dissipation of the vortex energy in a stratified ambient ocean only episodically. Finally, a critical discussion about the appropriateness of the different approximations assumed in the investigation is presented. In particular, the appropriateness of the reduced-gravity assumption is discussed. Results are consistent with those obtained analytically in the frame of the frontal-geostrophic theory: Although the effect of an active ambient layer on the vortex dynamics is found to be virtually absent only for unrealistically large water depths, it appears that the reduced-gravity model describes warm-core eddies acceptably for values of the ratio between maximum vortex thickness and total water depth typical for Gulf Stream rings.

Description: Zonally symmetric fluctuations of the midlatitude westerly winds characterize the primary mode of atmospheric variability in the Southern Hemisphere during all seasons. This is true not only in observations but also in an unforced 15 000-yr integration of a coarse-resolution (R15) coupled ocean–atmosphere model. Here it is documented how this mode of atmospheric variability, known as the Southern Annular Mode (SAM), generates ocean circulation and sea ice variations in the model integration on interannual to centennial timescales that are tightly in phase with the SAM. The positive phase of the SAM is associated with an intensification of the surface westerlies over the circumpolar ocean (around 60°S), and a weakening of the surface westerlies farther north. This induces Ekman drift to the north at all longitudes of the circumpolar ocean, and Ekman drift to the south at around 30°S. Through mass continuity, the Ekman drift generates anomalous upwelling along the margins of the Antarctic continent, and downwelling around 45°S. The anomalous flow diverging from the Antarctic continent also increases the vertical tilt of the isopycnals in the Southern Ocean, so that a more intense circumpolar current is also closely associated with positive SAM. In addition, the anomalous divergent flow advects sea ice farther north, resulting in an increase in sea ice coverage. Finally, positive SAM drives increases in poleward heat transport at about 30°S, while decreases occur in the circumpolar region. Ocean and sea ice anomalies of the opposite sign occur when the SAM is negative. The ocean and sea ice fluctuations associated with the SAM constitute a significant fraction of simulated ocean variability poleward of 30°S year-round. The robustness of the mechanisms relating the SAM to oceanic variability suggests that the SAM is likely an important source of large-scale variability in the real Southern Hemisphere ocean.

Description: Two configurations of a primitive-equation model of the North Atlantic are analyzed with respect to the simulated cycling of energy, mass, and heat in the upper ocean. One model is eddy-permitting (1/3° horizontal resolution), the other one is eddy-resolving (1/9° resolution), with both models using identical topographies and identical forcing fields at the surface and lateral boundaries. Besides showing some improvement in the simulated mean circulation and heat budgets, the eddy-resolving model reaches good agreement with satellite altimeter measurements of sea surface height variability. An unexpected finding of the model intercomparison is that simulated winter mixed layer depths in mid and high latitudes turn out to be systematically shallower by some 50 to 500 m in the higher resolution run, thereby agreeing better with observations than the 1/3° model results. This model improvement is related to enhanced levels of baroclinic instability leading to a decrease in potential energy and an associated increase in stratification. In the high-resolution model, shear-induced tilting of lateral density gradients generates stratification within the mixed layer itself, at a rate sufficient to set off an average surface heat loss of 5 W m–2 in mid and high latitudes. Although this is small compared to present uncertainties in surface heat fluxes, the resulting reduction in mixed layer depths may be important for an accurate simulation of water mass formation, air–sea gas exchange, and marine biological production. With traditional formulations of mixed layer physics assuming that properties are set by purely vertical mixing, and parameterizations of lateral subgrid-scale mixing often being tapered to zero in the mixed layer, present mixing schemes would have to be modified in order to account for eddy-induced generation of stratification in the surface mixed layer in noneddy-resolving ocean models.

Description: Time series of hydrographic and transient tracer (3H and 3He) observations from the central Labrador Sea collected between 1991 and 1996 are presented to document the complex changes in the tracer fields as a result of variations in convective activity during the 1990s. Between 1991 and 1993, as atmospheric forcing intensified, convection penetrated to progressively increasing depths, reaching 2300 m in the winter of 1993. Over that period the potential temperature (θ)/salinity (S) properties of Labrador Sea Water stayed nearly constant as surface cooling and downward mixing of freshwater was balanced by excavating and upward mixing of the warmer and saltier Northeast Atlantic Deep Water. It is shown that the net change in heat content of the water column (150–2500 m) between 1991 and 1993 was negligible compared to the estimated mean heat loss over that period (110 W m−2), implying that the lateral convergence of heat into the central Labrador Sea nearly balances the atmospheric cooling on a surprisingly short timescale. Interestingly, the 3H–3He age of Labrador Sea Water increased during this period of intensifying convection. Starting in 1995, winters were milder and convection was restricted to the upper 800 m. Between 1994 and 1996, the evolution of 3H–3He age is similar to that of a stagnant water body. In contrast, the increase in θ and S over that period implies exchange of tracers with the boundaries via both an eddy-induced overturning circulation and along-isopycnal stirring by eddies [with an exchange coefficient of O(500 m2 s−1)]. The authors construct a freshwater budget for the Labrador Sea and quantitatively demonstrate that sea ice meltwater is the dominant cause of the large annual cycle of salinity in the Labrador Sea, both on the shelf and the interior. It is shown that the transport of freshwater by eddies into the central Labrador Sea (140 cm between March and September) can readily account for the observed seasonal freshening. Finally, the authors discuss the role of the eddy-induced overturning circulation with regard to transport and dispersal of the newly ventilated Labrador Sea Water to the boundary current system and compare its strength (2–3 Sv) to the diagnosed buoyancy-forced formation rate of Labrador Sea Water.

Description: Vertical profiles of horizontal currents and hydrographic measurements from three cruises along 80.5°E from the coast of Sri Lanka to 6°S between December 1990 and September 1994 are used to investigate the scales of the Indian Ocean deep jets as well as internal wave parameters and dissipation at the equator. The deep jets at 80.5°E have a vertical wavelength of 660 sm (stretched meters) and amplitudes exceeding 10 cm s−1 in zonal velocity. They are observed throughout the water column and their flow direction reverses at 2° off the equator. The vertical positions of the jets differ among the cruises and are consistent with a flow reversal between the data collected in winter and summer. During September 1994, the jets were less pronounced. Due to the meridional distribution of their zonal velocity and the phase relationship between zonal velocity and vertical displacement, the jets are best described as nondispersive first-mode equatorial Rossby waves. The hydrographic data revealed thick layers of low stratification with vertical scales of 15–55 m in the upper 2000 m of the water column. They are found primarily within 1° of the equator and there is some evidence of correlation between the vertical position as well as the extent and the high strain zones of the deep jets. At vertical wavenumbers larger than those of the deep jets, shear and strain levels are five times larger than at off-equatorial locations and the compliant internal wave range (“roll-off range”) begins at a smaller wavenumber (kc ≈ 0.02 cpsm). An estimate of the average dissipation rate within the deep jets yielded = 7.5 × 10−10 W kg−1 between 500- and 2000-m depth. The elevated finescale internal wave field appears to be the main cause for the existence of the low stratification layers.

Description: A unique open-ocean upwelling exists in the tropical South Indian Ocean (SIO), a result of the negative wind curl between the southeasterly trades and equatorial westerlies, raising the thermocline in the west. Analysis of in situ measurements and a model-assimilated dataset reveals a strong influence of subsurface thermocline variability on sea surface temperature (SST) in this upwelling zone. El Niño–Southern Oscillation (ENSO) is found to be the dominant forcing for the SIO thermocline variability, with SST variability off Sumatra, Indonesia, also making a significant contribution. When either an El Niño or Sumatra cooling event takes place, anomalous easterlies appear in the equatorial Indian Ocean, forcing a westward-propagating downwelling Rossby wave in the SIO. In phase with this dynamic Rossby wave, there is a pronounced copropagation of SST. Moreover, a positive precipitation anomaly is found over, or just to the south of, the Rossby wave–induced positive SST anomaly, resulting in a cyclonic circulation in the surface wind field that appears to feedback onto the SST anomaly. Finally, this downwelling Rossby wave also increases tropical cyclone activity in the SIO through its SST effect. This coupled Rossby wave thus offers potential predictability for SST and tropical cyclones in the western SIO. These results suggest that models that allow for the existence of upwelling and Rossby wave dynamics will have better seasonal forecasts than ones that use a slab ocean mixed layer. The lagged-correlation analysis shows that SST anomalies off Java, Indonesia, tend to precede those off Sumatra by a season, a time lead that may further increase the Indian Ocean predictability.

Description: The Eurofloat experiment was a joint initiative to examine the large-scale spreading of Mediterranean Water (MW) and Labrador Sea Water in the northeast North Atlantic. RAFOS float data from the southern (MW) portion of the Eurofloat experiment have been examined in conjunction with historical float data in order to calculate quasi-Eulerian means in an effort to separate and quantify the constituents of the spreading of the MW tongue east of the Mid-Atlantic Ridge. While recent studies focussed chiefly on the role of meddies in the shaping of the MW tongue, this analysis also examines the tongue's second constituent, that is, the “background” (non-meddy advective and diffusive) flow. The results suggest the existence of two regimes approximately to the north and south of the 36°N parallel (i.e., the latitude of the Gulf of Cadiz), which are distinguished by different types of dominant spreading mechanisms for MW. To the south of the Gulf of Cadiz, the background flow shows an incoherent and weak mean, whereas the mean velocity of the salt enhanced meddies is strong and to the southwest. In contrast, to the north of 36°N the mean velocity of the meddies seems to be less pronounced and the background flow is shown to be a major component in the northwestward spreading of the MW tongue. The two regimes are separated by the Azores Current, which previously has been hypothesized to act as a dynamic barrier to the southward advective spreading of the background regime, which the meddies are able to penetrate because of their high kinetic energy. Overall, the meddies are calculated to contribute to approximately half of the total salinity anomaly flux.

Description: A method for combining ground-based passive microwave radiometer retrievals of integrated liquid water (LWP), radar reflectivity profiles (Z), and statistics of a cloud model is proposed for deriving cloud liquid water profiles (LWC). A dynamic cloud model is used to determine Z–LWC relations and their errors as functions of height above cloud base. The cloud model is also used to develop an LWP algorithm based on simulations of brightness temperatures of a 20–30-GHz radiometer. For the retrieval of LWC, the radar determined Z profile, the passive microwave retrieved LWP, and a model climatology are combined by an inverse error covariance weighting method. Model studies indicate that LWC retrievals with this method result in rms errors that are about 10%–20% smaller in comparison to a conventional LWC algorithm, which constrains the LWC profile exactly to the measured LWP. According to the new algorithm, errors in the range of 30%–60% are to be anticipated when profiling LWC. The algorithm is applied to a time series measurement of a stratocumulus layer at GKSS in Geesthacht, Germany. The GKSS 95-GHz cloud radar, a 20–30-GHz microwave radiometer, and a laser ceilometer were collocated within a 5-m radius and operated continuously during the measurement period. The laser ceilometer was used to confirm the presence of drizzle-sized drops.

Description: The authors derive a string function that describes the propagation of large-scale, potentially large amplitude, baroclinic energy anomalies in a two-layer ocean with variable topography and rotation parameter. The generality of the two-layer results allows results for the 1-layer, 1.5-layer, inverted 1.5-layer, lens, and dome models to be produced as limiting-cases. The string function is a scalar field that acts as a streamfunction for the propagation velocity. In the linear case the string function is simply c2o/f, where co is the background baroclinic shallow water wave speed, and typically describes propagation poleward on the eastern boundaries, westward (with some topographic steering) over the middle ocean, and equatorward on the western boundaries. In the more general nonlinear case, the string function is locally distorted by the anomaly. In the fully nonlinear examples of a lens or dome, there is no rest or background string function; the string function is generated entirely by the disturbance and propagation is due to asymmetric distribution of the anomalous mass over the string function contours. It is shown that conventional beta/topographic propagation results (e.g., beta drift of eddies, the Nof speed of cold domes) can be obtained as limiting cases of the string function. The string function provides, however, more general propagation velocities that are also usually simpler to derive. The first baroclinic mode string function for the global oceans is calculated from hydrographic data. The westward propagation speeds in the ocean basins as derived from the meridional gradient of the string function are typically two to five times faster than those expected from standard theory and agree well with the propagation speeds observed for long baroclinic Rossby waves in the TOPEX/Poseidon data.

Description: A general circulation ocean model has been used to study the formation and propagation mechanisms of North Atlantic Oscillation (NAO)-generated temperature anomalies along the pathway of the North Atlantic Current (NAC). The NAO-like wind forcing generates temperature anomalies in the upper 440 m that propagate along the pathway of the NAC in general agreement with the observations. The analysis of individual components of the ocean heat budget reveals that the anomalies are primarily generated by the wind stress anomaly-induced oceanic heat transport divergence. After their generation they are advected with the mean current. Surface heat flux anomalies account for only one-third of the total temperature changes. Along the pathway of the NAC temperature anomalies of opposite signs are formed in the first and second halves of the pathway, a pattern called here the North Atlantic dipole (NAD). The response of the ocean depends fundamentally on Rt = (L/υ)/τ, the ratio between the time it takes for anomalies to propagate along the NAC [(L/υ) 10 years] compared to the forcing period τ. The authors find that for NAO periods shorter than 4 years (Rt 〉 1) the response in the subpolar region is mainly determined by the local forcing. For NAO periods longer than 32 years (Rt 〈 1); however, the SST anomalies in the northeastern part of the NAD become controlled by ocean advection. In the subpolar region maximal amplitudes of the temperature response are found for intermediate (decadal) periods (Rt 1) where the propagation of temperature anomalies constructively interferes with the local forcing. A comparison of the NAO-generated propagating temperature anomalies with those found in observations will be discussed.

Description: A densely spaced hydrographic survey of the northern Irminger Basin together with satellite-tracked near-surface drifters confirm the intense mesoscale variability within and above the Denmark Strait overflow. In particular, the drifters show distinct cyclonic vortices over the downslope edge of the outflow plume. Growing perturbations such as these can be attributed to the baroclinic instability of a density current. A primitive equation model with periodic boundaries is used to simulate the destabilization of an idealized dense filament on a continental slope that resembles the northeastern Irminger Basin. Unstable waves evolve rapidly if the initial temperature profile is perturbed with a sinusoidal anomaly that exceeds a certain cutoff wavelength. As the waves grow to large amplitudes isolated eddies of both signs develop. Anticyclones form initially within the dense filament and are rich in overflow water. In contrast, cyclones form initially with their center in the ambient water but wrap outflow water around their center, thus containing a mixture of both water types. The nonlinear advection of waters that were originally located within the front between both water masses contributes most significantly to the stronger intensification of the cyclones in comparison with anticyclones. The frontal waters carry positive relative vorticity into the center of the cyclone. The process bears therefore some resemblance to atmospheric frontal cyclogenesis. After saturation there is a bottom jet of overflow water that is confined by counterrotating eddies: anticyclones upslope and cyclones downslope of the overflow core. The parameter dependence of the maximum growth rate is studied, and the implications of eddy-induced mixing for the water mass modification is discussed.

Description: Observations from the WOCE PCM-1 moored current meter array east of Taiwan for the period September 1994 to May 1996 are used to derive estimates of the Kuroshio transport at the entrance to the East China Sea. Three different methods of calculating the Kuroshio transport are employed and compared. These methods include 1) a “direct” method that uses conventional interpolation of the measured currents and extrapolation to the surface and bottom to estimate the current structure, 2) a “dynamic height” method in which moored temperature measurements from moorings on opposite sides of the channel are used to estimate dynamic height differences across the current and spatially averaged baroclinic transport profiles, and 3) an “adjusted geostrophic” method in which all moored temperature measurements within the array are used to estimate a relative geostrophic velocity field that is referenced and adjusted by the available direct current measurements. The first two methods are largely independent and are shown to produce very similar transport results. The latter two methods are particularly useful in situations where direct current measurements may have marginal resolution for accurate transport estimates. These methods should be generally applicable in other settings and illustrate the benefits of including a dynamic height measuring capability as a backup for conventional direct transport calculations. The mean transport of the Kuroshio over the 20-month duration of the experiment ranges from 20.7 to 22.1 Sv (1 Sv ≡ 106 m3 s−1) for the three methods, or within 1.3 Sv of each other. The overall mean transport for the Kuroshio is estimated to be 21.5 Sv with an uncertainty of 2.5 Sv. All methods show a similar range of variability of ±10 Sv with dominant timescales of several months. Fluctuations in the transport are shown to have a robust vertical structure, with over 90% of the transport variance explained by a single vertical mode. The moored transports are used to determine the relationship between Kuroshio transport and sea-level difference between Taiwan and the southern Ryukyu Islands, allowing for long-term monitoring of the Kuroshio inflow to the East China Sea.

Description: Transient eddies in the atmosphere induce a poleward transport of heat and moisture. A moist static energy budget of the surface layer is determined from the NCEP reanalysis data to evaluate the impact of the storm track. It is found that the transient eddies induce a cooling and drying of the surface layer with a monthly mean maximum of 60 W m−2. The cooling in the midlatitudes extends zonally over the entire basin. The impact of this cooling and drying on surface heat fluxes, sea surface temperature (SST), water mass transformation, and vertical structure of the Pacific is investigated using an ocean model coupled to an atmospheric mixed layer model. The cooling by atmospheric storms is represented by adding an eddy-induced transfer velocity to the mean velocity in an atmospheric mixed layer model. This is based on a parameterization of tracer transport by eddies in the ocean. When the atmospheric mixed layer model is coupled to an ocean model, realistic SSTs are simulated. The SST is up to 3 K lower due to the cooling by storms. The additional cooling leads to enhanced transformation rates of water masses in the midlatitudes. The enhanced shallow overturning cells affect even tropical regions. Together with realistic SST and deep winter mixed layer depths, this leads to formation of homogeneous water masses in the upper North Pacific, in accordance to observations.

Description: The dynamics of the Rhine outflow plume in the proximity of the river mouth is investigated by using remote sensing data and numerical simulations. The remote sensing data consist of 41 synthetic aperture radar (SAR) images acquired by the First and Second European Remote Sensing satellites ERS-1 and ERS-2 over the outflow region of the river Rhine. Most of them show sea surface signatures of oceanic phenomena, for example, surface current and wind variations, ship wakes, and oil slicks. In particular, in 36 of these images pronounced frontal features are visible as narrow zones of mainly enhanced, sometimes enhanced/reduced radar backscatter that can be associated with the Rhine surface front. Within the area enclosed by the frontal line, large zones characterized by a lower radar backscatter than in the outer area are often visible. The analysis of the ERS SAR images suggests that the form and the location of the frontal features are mainly linked to the semidiurnal tidal phase in the outflow region, although their variability suggests also that they weakly depend on river discharge, residual currents, and neap-spring tidal cycle. In order to test this observational hypothesis, the results obtained from the analysis of the ERS SAR images are compared with the results obtained from the numerical simulation of the hydrodynamics of the Rhine outflow region carried out using a two-layer, frontal model, which is based on the nonlinear, hydrostatic shallow-water equations on an f plane. The model is forced by prescribing tidal and residual currents and river discharge at the open boundaries. Several simulations are performed by varying the values of these forcing parameters. The numerical results corroborate the observational conjecture: It is found that the form and the location of the simulated interface outcropping lines in the proximity of the river mouth are mainly determined by the semidiurnal tidal phase in the outflow region and that river discharge, residual currents, and neap-spring tidal cycle contribute only secondarily to their determination. Inserting the simulated surface velocity field into a simple radar-imaging model that relates the modulation of the backscattered radar power to the surface velocity convergence in radar look direction, narrow, elongated bands of enhanced radar backscatter emerge near the model frontal line while patches of low radar backscatter appear within the simulated Rhine plume area. The consistency of the model results with the results obtained from the analysis of the SAR images enables one to infer a mean spatial and temporal evolution of the Rhine outflow plume over a semidiurnal tidal cycle from the analysis of spaceborne SAR images acquired during different tidal cycles over the Rhine outflow area and suggests the possibility of using numerical modeling, in conjunction with the analysis of spaceborne measurements, for monitoring the oceanic variability in the Rhine outflow area

Description: Meridional transports of mass, heat, nutrients, and carbon across coast-to-coast WOCE and pre-WOCE sections between 11°S and 45°S in the South Atlantic are calculated using an inverse model. Usually salt preservation is used as a condition in the inverse model, and only in the case of heat transport the condition of zero total mass transport is taken instead. Other constraints include silica conservation, prescribed southward fluxes of salt and phosphate, and transports in the southward Brazil Current and in the northward Antarctic Bottom Water flow obtained from WOCE moored current meter arrays. The constraints set the underdetermined system of linear equations of the inverse model whose solutions depend on weights, scales, and matrix ranks. The discussion emphasizes the sensitivity of the fluxes to changes in the model input. The transports given in the following are obtained as the means of “reasonable” solutions at 30°S. The error numbers in parentheses include uncertainties due to wind stress and temporal variability, the numbers without parentheses do not contain these terms:0.53 ± 0.03 (0.09) Tg s−1 mass to the south, 0.29 ± 0.05 (0.24) PW heat to the north, 15 ± 120 (500) kmol s−1 oxygen to the south, 121 ± 22 (75) kmol s−1 nitrate to the south, 64 ± 110 (300) silica to the north, and 1997 ± 215 (600) kmol s−1 dissolved inorganic carbon to the south. The above errors in transports are obviously dominated by uncertainties in wind stress and temporal variability. The divergence in meridional heat and mass transport is consistent with integral surface flux changes between corresponding zonal bands. The mass compensation of southward flowing North Atlantic Deep Water occurs to a greater extent in the warm surface waters than in the Antarctic Intermediate Water below. If one follows the arguments of earlier authors on the relation between meridional fluxes and the significance of the two possible pathways for the global thermohaline circulation, the warm water path south of Africa seems to be somewhat more important than the cold water path through Drake Passage.