ALBERT

Description: In this study, the authors discuss two different parameterizations for the effect of mixed layer eddies, one based on ageostrophic linear stability analysis (ALS) and the other one based on a scaling of the potential energy release by eddies (PER). Both parameterizations contradict each other in two aspects. First, they predict different functional relationships between the magnitude of the eddy fluxes and the Richardson number (Ri) related to the background state. Second, they also predict different vertical structure functions for the horizontal eddy fluxes. Numerical simulations for two different configurations and for a large range of different background conditions are used to evaluate the parameterizations. It turns out that PER is better suited to capture the Ri dependency of the magnitude of the eddy fluxes. On the other hand, the vertical structure of the meridional eddy fluxes predicted by ALS is more accurate than that of PER, while the vertical structure of the vertical eddy fluxes is well predicted by both parameterizations. Therefore, this study suggests the use of the magnitude of PER and the vertical structure functions of ALS for an improved parameterization of mixed layer eddy fluxes.

Description: Northern Hemisphere stratospheric variability is investigated with respect to chaotic behavior using time series from three different variables extracted from four different reanalysis products and two numerical model runs with different forcing. The time series show red spectra at all frequencies and the probability distribution functions show persistent deviations from a Gaussian distribution. An exception is given by the numerical model forced with perpetual winter conditions—a case that shows more variability and follows a Gaussian distribution, suggesting that the deviation from Gaussianity found in the observations is due to the transition between summer and winter variability. To search for the presence of a chaotic attractor the correlation dimension and entropy, the Lyapunov spectrum, and the associated Kaplan–Yorke dimension are estimated. A finite value of the dimensions can be computed for each variable and data product, with the correlation dimension ranging between 3.0 and 4.0 and the Kaplan–Yorke dimension between 3.3 and 5.5. The correlation entropy varies between 0.6 and 1.1. The model runs show similar values for the correlation and Lyapunov dimensions for both the seasonally forced run and the perpetual-winter run, suggesting that the structure of a possible chaotic attractor is not determined by the seasonality in the forcing, but must be given by other mechanisms.

Description: Southern Hemisphere (SH) stratospheric variability is investigated with respect to chaotic behavior using time series from three different variables extracted from four different reanalysis products. The results are compared with the same analysis applied to the Northern Hemisphere (NH). The probability density functions (PDFs) for the SH show persistent deviations from a Gaussian distribution. The variability is given by white spectra for low frequencies, a slope of −1 for intermediate frequencies, and −3 slopes for high frequencies. Considering the time series for winter and summer separately, PDFs show a Gaussian distribution and the variability spectra change their slopes, indicating the role of the transition between winter and summer variability in shaping the time series. The correlation (D2) and the Kaplan–Yorke (DKY) dimensions are estimated. A finite value of the dimensions can be computed for each variable and data product, except for the NCEP zonal-mean zonal wind and temperature data, which violate the requirement D2 ≤ DKY, possibly owing to the presence of spurious trends and inconsistencies in the data. The value of D2 ranges between 2.6 and 3.9, while DKY ranges between 3.0 and 4.5. The results show that both D2 and DKY display large variability in their values both for different datasets and for different variables within the same dataset. The variability of the values of D2 and DKY thus leaves open the question about the existence of a low-dimensional attractor or if the finite dimensions of the system are the result of the projection of a larger attractor in a low-dimensional embedding space.

Description: The sensitivities of the Brewer–Dobson circulation (BDC) to different distributions of tropical SST heating are investigated in an idealized aquaplanet model. It is found that an increase in tropical SSTs generally leads to an acceleration of tropical upwelling and an associated reduction in the age of air (AOA) in the polar stratosphere and that the AOA near the subtropical tropopause is correlated with local isentropic mixing of tropospheric air with stratospheric air. The zonal distribution of SST perturbations has a major impact on the vertical and meridional structure of the BDC as compared with other SST characteristics. Zonally localized SST heatings tend to generate a shallow acceleration of the stratospheric residual circulation, enhanced isentropic mixing associated with a weakened stratospheric jet, and a reduction in AOA mostly within the polar vortex. In contrast, SST heatings with a zonally symmetric structure tend to produce a deep strengthening of the stratospheric residual circulation, suppressed isentropic mixing associated with a stronger stratospheric jet, and a decrease of AOA in the entire stratosphere. The shallow versus deep strengthening of the stratospheric residual circulation change has been linked to wave propagation and dissipation in the subtropical lower stratosphere rather than wave generation in the troposphere, and the former can be strongly affected by the vertical position of the subtropical jet. These results suggest that, while the longitudinally localized SST trends under climate change may contribute to the change in the shallow branch of the BDC, the upward shift of the subtropical jet associated with the zonal SST heating can impact the deep branch of the BDC.

Description: The authors test the hypothesis that recent observed trends in surface westerlies in the Southern Hemisphere are directly consequent on observed trends in the timing of stratospheric final warming events. The analysis begins by verifying that final warming events have an impact on tropospheric circulation in a simplified GCM driven by specified equilibrium temperature distributions. Seasonal variations are imposed in the stratosphere only. The model produces qualitatively realistic final warming events whose influence extends down to the surface, much like what has been reported in observational analyses. The authors then go on to study observed trends in surface westerlies composited with respect to the date of final warming events. If the considered hypothesis were correct, these trends would appear to be much weaker when composited with respect to the date of the final warming events. The authors find that this is not the case, and accordingly they conclude that the observed surface changes cannot be attributed simply to this shift toward later final warming events.

Description: Precipitation is highly variable in space and time; hence, rain gauge time series generally exhibit additional random small-scale variability compared to area averages. Therefore, differences between daily precipitation statistics simulated by climate models and gauge observations are generally not only caused by model biases, but also by the corresponding scale gap. Classical bias correction methods, in general, cannot bridge this gap; they do not account for small-scale random variability and may produce artifacts. Here, stochastic model output statistics is proposed as a bias correction framework to explicitly account for random small-scale variability. Daily precipitation simulated by a regional climate model (RCM) is employed to predict the probability distribution of local precipitation. The pairwise correspondence between predictor and predictand required for calibration is ensured by driving the RCM with perfect boundary conditions. Wet day probabilities are described by a logistic regression, and precipitation intensities are described by a mixture model consisting of a gamma distribution for moderate precipitation and a generalized Pareto distribution for extremes. The dependence of the model parameters on simulated precipitation is modeled by a vector generalized linear model. The proposed model effectively corrects systematic biases and correctly represents local-scale random variability for most gauges. Additionally, a simplified model is considered that disregards the separate tail model. This computationally efficient model proves to be a feasible alternative for precipitation up to moderately extreme intensities. The approach sets a new framework for bias correction that combines the advantages of weather generators and RCMs.

Description: There is an ongoing discussion in the community concerning the wave-averaged momentum equations in the hybrid vertically Lagrangian and horizontally Eulerian (VL) framework and, in particular, the form stress term (representing the residual effect of pressure perturbations) which is thought to restrict the handling of higher order waves in terms of a perturbation expansion. The present study shows that the traditional pressure-based form stress term can be transformed into a set of terms that do not contain any pressure quantities but do contain the time derivative of a wave-induced velocity. This wave-induced velocity is referred to as the pseudomomentum in the VL framework, as it is analogous to the generalized pseudomomentum in Andrews and McIntyre. This enables the second expression for the wave-averaged momentum equations in the VL framework (this time for the development of the total transport velocity minus the VL pseudomomentum) to be derived together with the vortex force. The velocity-based expression of the form stress term also contains the residual effect of the turbulent viscosity, which is useful for understanding the dissipation of wave energy leading to transfer of momentum from waves to circulation. It is found that the concept of the virtual wave stress of Longuet-Higgins is applicable to quite general situations: it does not matter whether there is wind forcing or not, the waves can have slow variations, and the viscosity coefficient can vary in the vertical. These results provide a basis for revisiting the surface boundary condition used in numerical circulation models.

Description: Mooring observations and model simulations point to an instability of the Labrador Current (LC) during winter, with enhanced eddy kinetic energy (EKE) at periods between 2 to 5 days, and much less EKE during other seasons. Linear stability analysis using vertical shear and stratification from the model reveals three dominant modes of instability in the LC: - a balanced interior mode with along-flow wavelengths of about 30–45 km, phase velocities of 0.3 m/s, maximal growth rates of 1 d−1 and surface intensified, but deep reaching amplitudes, - a balanced shallow mode with along-flow wavelengths of about 0.3–1.5 km, about three times larger phase speeds and growth rates, but amplitudes confined to the mixed layer (ML), - and an unbalanced symmetric mode with largest growth rates, vanishing phase speeds and along-flow structure, and very small cross-flow wavelengths, also confined to the ML. Both balanced modes are akin to baroclinic instability, but operate at moderate to small Richardson numbers Ri with much larger growth rates as for the quasi-geostrophic limit of Ri ≫ 1. The interior mode is found to be responsible for the instability of the LC during winter. Weak stratification and enhanced vertical shear due to local buoyancy loss and the advection of convective water masses from the interior result in small Ri within the LC, and to three times larger growth rates of the interior mode in March compared to summer and fall conditions. Both the shallow and the symmetric mode are not resolved by the model, but it is suggested that they might also play an important role for the instability in the LC and for lateral mixing.

Description: Sensible and latent heat fluxes were estimated from turbulence measurements gathered during several Atlantic transects of the R/V Polarstern. The inertial dissipation method was used to analyze the data. Resulting bulk transfer coefficients were then applied to the data from the ship’s meteorological system to get continuous time series of the heat fluxes. Combined to the measured downward solar and longwave radiation fluxes allows for an estimate of the total energy budget at the air-sea interface. Comparing these parameterized energy fluxes to ones based on the COARE 3.0 bulk flux algorithm show very strong agreement.

Description: The dynamical origin of the spectral and autocorrelation structure of annular variability in the troposphere is investigated by a deductive approach. Specifically, the structure of the power spectrum and autocorrelation function of the zonal-mean geopotential is analyzed for the case of a quasigeostrophic spherical atmosphere subject to a white noise mechanical forcing applied in a single Hough mode and concentrated at a particular level in the vertical, with vertically uniform Newtonian cooling and Rayleigh drag concentrated at a rigid lower boundary. Analytic expressions for the power spectrum are presented together with expressions for an approximate red noise (i.e., a Lorentzian-shaped) power spectrum. It is found that for an infinitely deep atmosphere the power spectrum can be well approximated by a red noise process for the first few Hough modes (associated with large Rossby heights), provided the distance from the forcing is not larger than about one Rossby height. When a frictional rigid lower boundary is included, however, the approximation is generally bad. The high-frequency part of the power spectrum exhibits near-exponential behavior and the autocorrelation function shows a transition from a rapid decay at short lags to a much slower decay at longer lags, if the thermal and mechanical damping time scales are sufficiently well separated. Since observed annular variability exhibits the same characteristics, the above results lead to the hypothesis that these characteristics may, to some extent, be intrinsic to the linear zonal-mean response problem—although the need for an additional contribution from eddy feedbacks is also implied by the results.

Description: A representation of an equatorial basin mode excited in a shallow water model for a single high order baroclinic vertical normal mode is used as a simple model for the equatorial deep jets. The model is linearized about both a state of rest and a barotropic mean flow corresponding to the observed Atlantic Equatorial Intermediate Current System. We found that the eastward mean flow associated with the North and South Intermediate Counter Currents (NICC and SICC, respectively) effectively shields the Equator from off-equatorial Rossby waves. The westward propagation of these waves is blocked and focusing on the Equator due to beta dispersion is prevented. This leads to less energetic jets along the Equator. On the other hand, the westward barotropic mean flow along the Equator reduces the gradient of absolute vorticity and hence widens the cross-equatorial structure of the basin mode. Increasing lateral viscosity predominantly affects the width of the basin modes’ Kelvin wave component in the presence of the mean flow while the Rossby wave is confined by the flanking NICC and SICC. Independent of the presence of the mean flow, the application of sufficient lateral mixing also hinders the focusing of off-equatorial Rossby waves, which is hence an unlikely feature of a low-frequency basin mode in the real ocean.

Description: A decadal change in the character of ENSO was observed around year 2000 toward weaker-amplitude, higher-frequency events with an increased occurrence of central Pacific El Niños. Here these changes are assessed in terms of the Bjerknes stability index (BJ index), which is a measure of the growth rate of ENSO-related SST anomalies. The individual terms of the index are calculated from ocean reanalysis products separately for the time periods 1980–99 and 2000–10. The spread between the products is large, but they show a robust weakening of the thermocline feedback due to a reduced thermocline slope response to anomalous zonal wind stress as well as a weakened wind stress response to eastern equatorial Pacific SST anomalies. These changes are consistent with changes in the background state of the tropical Pacific: cooler mean SST in the eastern and central equatorial Pacific results in reduced convection there together with a westward shift in the ascending branch of the Walker circulation. This shift leads to a weakening in the relationship between eastern Pacific SST and longitudinally averaged equatorial zonal wind stress. Also, despite a steeper mean thermocline slope in the more recent period, the thermocline slope response to wind stress anomalies weakened due to a smaller zonal wind fetch that results from ENSO-related wind anomalies being more confined to the western basin. As a result, the total BJ index is more negative, corresponding to a more strongly damped system in the past decade compared to the 1980s and 1990s.

Description: In his comment, G. Bürger criticizes the conclusion that inflation of trends by quantile mapping is an adverse effect.He assumes that the argument would be ‘‘based on the belief that long-term trends and along with them future climate signals are to be large scale.’’ His line of argument reverts to the so-called inflated regression. Here it is shown, by referring to previous critiques of inflation and standard literature in statistical modeling as well as weather forecasting, that inflation is built upon a wrong understanding of explained versus unexplained variability and prediction versus simulation. It is argued that a sound regressionbased downscaling can in principle introduce systematic local variability in long-term trends, but inflation systematically deteriorates the representation of trends. Furthermore, it is demonstrated that inflation by construction deteriorates weather forecasts and is not able to correctly simulate small-scale spatiotemporal structure.

Description: The relationship between the Agulhas Current and the Agulhas leakage is not well understood. Here, this is investigated using two basin-scale and two global ocean models, of incrementally increasing resolution. The response of the Agulhas Current is evaluated under a series of sensitivity experiments, in which idealised anomalies, designed to geometrically modulate zonal trade wind stress, are applied across the Indian Ocean basin. The imposed wind stress changes exceed ±2 standard deviations from the annual mean trade winds and, in the case of intensification, are partially representative of recently observed trends. The Agulhas leakage is quantified using complimentary techniques based on Lagrangian virtual floats and Eulerian passive tracer flux. As resolution increases, model behavior converges and the sensitivity of the leakage to Agulhas Current transport anomalies is reduced. In the two eddy-resolving configurations tested, the leakage is insensitive to changes in Agulhas Current transport at 32°S, though substantial eddy kinetic energy anomalies are evident. Consistent with observations, the position of the retroflection remains stable. The decoupling of Agulhas Current variability from the Agulhas leakage suggests that, while correlations between the two may exist, they may not have a clear dynamical basis. It is suggested that present and future Agulhas leakage proxies be considered in the context of potentially transient forcing regimes.

Description: A surface diurnal warm layer is diagnosed from Seaglider observations, and develops on half the days in the CINDY/DYNAMO Indian Ocean experiment. The diurnal warm layer occurs on days of high solar radiation flux (〉 80 W m−2) and low wind speed (〈 6 m s−1), and preferentially in the inactive stage of the Madden–Julian Oscillation. Its diurnal harmonic has an exponential vertical structure with a depth scale of 4–5 m (dependent on chlorophyll concentration), consistent with forcing by absorption of solar radiation. The effective sea surface temperature (SST) anomaly due to the diurnal warm layer often reaches 0.8°C in the afternoon, with a daily mean of 0.2°C, rectifying the diurnal cycle onto longer time scales. This SST anomaly drives an anomalous flux of 4 W m−2 that cools the ocean. Alternatively, in a climate model where this process is unresolved, this represents an erroneous flux that warms the ocean. A simple model predicts a diurnal warm layer to occur on 30–50% of days across the tropical warm pool. On the remaining days, with low solar radiation and high wind speeds, a residual diurnal cycle is observed by the Seaglider, with a diurnal harmonic of temperature that decreases linearly with depth. As wind speed increases, this already weak temperature gradient decreases further, tending towards isothermal conditions.

Description: Sea surface temperature (SST) anomalies in the eastern equatorial Atlantic are connected to modulations in the strength of the South Atlantic subtropical high-pressure system, referred to as the South Atlantic Anticyclone (SAA). Using ocean and atmosphere reanalysis products we show here that the strength of the SAA from February to May impacts the timing of the cold tongue onset and the intensity of its development in the eastern equatorial Atlantic (EEA) via anomalous tropical wind power. This modulation of the timing and amplitude of the seasonal cold tongue development manifests as anomalous SST events peaking between June and August. The timing and impact of this connection is not completely symmetric for warm and cold events. For cold events, an anomalously strong SAA in February and March leads to positive wind power anomalies from February to June resulting in an early cold tongue onset and subsequent cold SST anomalies in June and July. For warm events the anomalously weak SAA persists until May, generating negative wind power anomalies that lead to a late cold tongue onset as well as a suppression of the cold tongue development and associated warm SST anomalies. Mechanisms by which SAA induced wind power variations south of the equator influence EEA SST are discussed, including ocean adjustment via Rossby and Kelvin wave propagation, meridional advection, and local intraseasonal wind variations

Description: We present a detailed quality assessment of a novel underwater sensor for the measurement of CO2 partial pressure (pCO2) based on surface water field deployments carried out between 2008 and 2011. The commercially available sensor, which is based on membrane equilibration and NDIR spectrometry is small and can be integrated into mobile platforms. It is calibrated in water against a proven flow-through pCO2 instrument within a custom-built calibration setup. The aspect of highest concern with respect to achievable data quality of the sensor is the compensation for signal drift inevitably connected to absorption measurements. We use three means to correct for drift effects: (i) a filter correlation or dual-beam setup, (ii) regular zero gas measurements realized automatically within the sensor and (iii) a zero-based transformation of two sensor calibrations flanking the time of sensor deployment. Three sensors were tested against an underway pCO2 system during two major research cruises providing an in situ temperature range from 7.4 to 30.1°C and pCO2 values between 289 and 445 μatm. The average difference between sensor and reference pCO2 was found to be -0.6 ± 3 μatm with a RMSE of 3.7 μatm.

Description: Ammassalik in southeast Greenland is known for strong wind events that can reach hurricane intensity and cause severe destruction in the local town. Yet, these winds and their impact on the nearby fjord and shelf region have not been studied in detail. Here, data from two meteorological stations and the European Centre for Medium-Range Weather Forecasts Interim Re-Analysis (ERA-Interim) are used to identify and characterize these strong downslope wind events, which are especially pronounced at a major east Greenland fjord, Sermilik Fjord, within Ammassalik. Their local and regional characteristics, their dynamics and their impacts on the regional sea ice cover, and air–sea fluxes are described. Based on a composite of the events it is concluded that wind events last for approximately a day, and seven to eight events occur each winter. Downslope wind events are associated with a deep synoptic-scale cyclone between Iceland and Greenland. During the events, cold dry air is advected down the ice sheet. The downslope flow is accelerated by gravitational acceleration, flow convergence inside the Ammassalik valley, and near the coast by an additional thermal and synoptic-scale pressure gradient acceleration. Wind events are associated with a large buoyancy loss over the Irminger Sea, and it is estimated that they drive one-fifth of the net wintertime loss. Also, the extreme winds drive sea ice out of the fjord and away from the shelf.

Description: In 2013, the vast majority of the monitored climate variables reported here maintained trends established in recent decades. ENSO was in a neutral state during the entire year, remaining mostly on the cool side of neutral with modest impacts on regional weather patterns' around the world. This follows several years dominated by the effects of either La Nina. or El Nino events. According to several independent analyses, 2013 was again among the 10 warmest years on record at the global scale, both at the Earth's surface and through the troposphere. Some regions in the Southern Hemispherehad record or near-record high temperatures for the year. Australia observed its hottest year on record, while Argentina and New Zealand reported their second and third hottest years, respectively. In Antarctica, Amundsen-Scott South Pole Station reported its highest annual temperature since records began in 1957 At the opposite pole, the Arctic observed its seventh warmest year since record's began in the early 20th century. At 20-m depth, record high temperatures were measured at some permafrost stations on the North Slope Of Alaska and in the Brooks Range. In the Northern Hemisphere extratropics, anomalous meridional atmospheric circulation occurred throughout much of the year leading to marked regional extremes of both temperature and precipitation. Cold temperature anomalies during Winter across Eurasia were followed by warm Spring temperature anomalies, which, were linked to a new record Eurasian snow cover extent in the May. Minimum sea ice extent in the :Arctic was the sixth lowest since satellite Observations began in 1979. Including 2013, all seven lowest extents on record have occurred in the past seven years Antarctica, on the other hand, had above average sea ice extent throughout 2013, with 116 days Of new daily high extent records, inclding a new daily maximum sea ice area of 19.57 million km(2) reached on 1 October. ENSO-neutral conditions in the eastern central Pacific Ocean and a negative Pacific decadal oscillation pattern in the North Pacific had the largest impacts on the global sea surface temperature in 2013: The North Pacific reached a historic high temperature in 2013 and on balance the globally-averaged sea surface temperature was among the 10 highest on record. Overall, the salt Content in near-surface ocean waters increased while in intermediate waters it decreased. Global mean sea level continued to rise during 2013, on pace with a trend of 3,2 mm yr(-1) over the past two decades. A portion of this trend (0.5 mm yr(-1)) has been attributed to natural variability associated with the Pacific decadal oscillation as Well as to ongoing contribution from the melting of glaciers and ice sheets and ocean warming. Global tropical cyclone frequency during 2013 was slightly above average with a total of 94 storms, although the North Atlantic Basin had its quietest hurricane season since 1994. In the Western North Pacific Basin, Super Typhoon Haiyan, the deadliest tropical Cyclone of 2013, had 1-minute sustained winds estimated to be 170 kt (OS m s(-1)) on 7 November, the highest Wind speed ever assigned to a tropical cyclone. High storm surge was also associated with Haiyan at it Made landfall over the central Philippines, an area where sea level is currently at historic highs increasing by 200 mm since 1970. In the atmosphere, carbon dioxide, methane, and nitrous oxide all Continued to increase in 2013. As in previous our years, each of these Major greenhouse gases once again reached historic high Concentrations. In the Arctic carbon dioxide and methane:, Increased at the same rate as:the global increase. These increases are likely due to export from lower latitudes rather than a consequence Of increases in Arctic sources such as thawing permafrost., At Mauna Loa, Hawaii, for the first time since measurement began in 1958, the daily average Mixing ratio Of carbon dioxide exceeded 400 ppm on 9 May. The state of these variables, along with dozens of others, and the 2013 climate conditions of regions around the world are discussed in further detail in this 24th, edition of the State of the Climate series

Description: The Denmark Strait Overflow (DSO) supplies about one-third of the North Atlantic Deep Water and is critical to global thermohaline circulation. Knowledge of the pathways of DSO through the Irminger Basin and its transformation there is still incomplete, however. The authors deploy over 10 000 Lagrangian particles at the Denmark Strait in a high-resolution ocean model to study these issues. First, the particle trajectories show that the mean position and potential density of dense waters cascading over the Denmark Strait sill evolve consistently with hydrographic observations. These sill particles transit the Irminger Basin to the Spill Jet section (65.25°N) in 5–7 days and to the Angmagssalik section (63.5°N) in 2–3 weeks. Second, the dense water pathways on the continental shelf are consistent with observations and particles released on the shelf in the strait constitute a significant fraction of the dense water particles recorded at the Angmagssalik section within 60 days (~25%). Some particles circulate on the shelf for several weeks before they spill off the shelf break and join the overflow from the sill. Third, there are two places where the water density following particle trajectories decreases rapidly due to intense mixing: to the southwest of the sill and southwest of the Kangerdlugssuaq Trough on the continental slope. After transformation in these places, the overflow particles exhibit a wide range of densities.

Description: Quantile mapping is routinely applied to correct biases of regional climate model simulations compared to observational data. If the observations are of similar resolution as the regional climate model, quantile mapping is a feasible approach. However, if the observations are of much higher resolution, quantile mapping also attempts to bridge this scale mismatch. Here, it is shown for daily precipitation that such quantile mapping-based downscaling is not feasible but introduces similar problems as inflation of perfect prognosis ("prog") downscaling: the spatial and temporal structure of the corrected time series is misrepresented, the drizzle effect for area means is overcorrected, area-mean extremes are overestimated, and trends are affected. To overcome these problems, stochastic bias correction is required.

Description: Previous attempts to derive the depth-dependent expression of the radiation stress have lead to a debate concerning (i) the applicability of Mellor’s approach to a sloping bottom, (ii) the introduction of the delta function at the mean sea surface in the later papers by Mellor, and (iii) a wave-induced pressure term derived in several recent studies. The authors use an equation system in vertically Lagrangian and horizontally Eulerian (VL) coordinates suitable for a concise treatment of the surface boundary, and obtain an expression for the depth-dependent radiation stress that is consistent with the vertically-integrated expression given by Longuet-Higgins and Stewart. Concerning (i)-(iii) in the above, the difficulty of handling a sloping bottom disappears when wave-averaged momentum equations in the VL coordinates are written for the development of (not the Lagrangian mean velocity but) the Eulerian mean velocity. There is also no delta function at the sea surface in the expression for the depth-dependent radiation stress. The connection between the wave-induced pressure term in the recent studies and the depth-dependent radiation stress term is easily shown by rewriting the pressure-based form stress term in the thickness-weighted-mean (TWM) momentum equations as a velocity-based term which contains the time derivative of the pseudomomentum in the TWM framework.

Description: The upper ocean, including the biologically productive euphotic zone and the mixed layer, has great relevance for studies of physical, biogeochemical, and ecosystem processes and their interaction. Observing this layer with a continuous presence, sampling many of the relevant variables, and with sufficient vertical resolution, has remained a challenge. Here a system is presented which can be deployed on the top of deep-ocean moorings, with a drive mechanism at depths of 150-200m, which mechanically winches a large sensor float and smaller communications float tethered above it to the surface and back down again, typically twice per day for periods up to 1 year. The sensor float can carry several sizeable sensors, and it has enough buoyancy to reach the near surface and for the communications float to pierce the surface even in the presence of strong currents. The system can survive mooring blow-over to 1000m depth. The battery-powered design is made possible by using a balanced energy-conserving principle. Reliability is enhanced with a drive assembly that employs a single rotating part that has no slip rings or rotating seals. The profiling bodies can break the surface to sample the near-surface layer and to establish satellite communication for data relay or reception of new commands. An inductive pass-through mode allows communication with other mooring components throughout the water column beneath the system. A number of successful demonstration deployments have been completed.

Description: In recent years, profiling floats, which form the basis of the successful international Argo observatory, are also being considered as platforms for marine biogeochemical research. This study showcases the utility of floats as a novel tool for combined gas measurements of CO2 partial pressure (pCO2) and O2. These float prototypes were equipped with a small-sized and submersible pCO2 sensor and an optode O2 sensor for high resolution measurements in the surface ocean layer. Four consecutive deployments were carried out during Nov. 2010 and June 2011 near the Cape Verde Ocean Observatory (CVOO) in the eastern tropical North Atlantic. The profiling float performed upcasts every 31 h while measuring pCO2, O2, salinity, temperature and hydrostatic pressure in the upper 200 m of the water column. In order to maintain accuracy, regular pCO2 sensor zeroings at depth and surface, as well as optode measurements in air, were performed for each profile. Through the application of data processing procedures (e.g., time-lag correction) accuracies of float-borne pCO2 measurements were greatly improved (10 – 15 μatm for water column and 5 μatm for surface measurements). O2 measurements yielded an accuracy of 2 μmol kg−1. First results of this pilot study show the possibility of using profiling floats as a platform for detailed and unattended observations of the marine carbon and oxygen cycle dynamics.

Description: The use of a coupled ocean/atmosphere/sea-ice model to hindcast (i.e. historical forecast) recent climate variability is described and illustrated for the cases of the 1976/77 and 1998/99 climate shift events in the Pacific. The initialization is achieved by running the coupled model in partially coupled mode whereby global observed wind stress anomalies are used to drive the ocean/sea-ice component of the coupled model while maintaining the thermodynamic coupling between the ocean/sea-ice and atmosphere components. Here we show that hindcast experiments can successfully capture many features associated with the 1976/77 and 1998/99 climate shifts. For instance, hindcast experiments started from the beginning of 1976 can capture sea surface temperature (SST) warming in the central-eastern equatorial Pacific and the positive phase of the Pacific Decadal Oscillation (PDO) throughout the 9 years following the 1976/77 climate shift, including the deepening of the Aleutian low pressure system. Hindcast experiments started from the beginning of 1998 can also capture part of the anomalous conditions during the 4 years after the 1998/99 climate. We argue that the dynamical adjustment of heat content anomalies that are present in the initial conditions in the tropics is important for the successful hindcast of the two climate shifts.

Description: One of the most remarkable features of the reproductive systems of eubrachyuran crabs is the presence of specialized organs for sperm storage, the seminal receptacles. Descriptions of seminal receptacle morphology, sperm storage time, sperm retention across molts, and the capacity to store multiple ejaculates from different males can help in understanding crab mating strategies as well as in preventing negative effects of male-biased fisheries of heavily harvested species. Metacarcinus edwardsii is the most harvested crab in Chile, but its reproductive biology is largely unstudied. In this study, the morphology of the seminal receptacles of M. edwardsii is characterized from the macroscopic to the microscopic level, during key points in the reproductive cycle. The receptacles of experimentally mated and wild-caught females were included in this analysis. Metacarcinus edwardsii has ventral-type seminal receptacles that are able to retain sperm after molting, and even after extrusion of the eggs. Stratification of multiple ejaculates is clearly observed. In general, the pattern of sperm storage indicates that populations of this species, like those of other cancrid crabs, could have high resilience to the negative effects of the selective harvest of males, principally because females have a great sperm storage capacity.

Description: The Arctic continental shelf seas hold a globally significant source of freshwater that impacts Arctic Ocean stratification, circulation, and climate. This freshwater can be injected below the surface mixed layer by intense turbulent kinetic energy dissipation events, as resolved by Laptev Sea microstructure observations. The tides provide a major source of energy that can be dissipated and hence drive diapycnal mixing in the Laptev Sea. Multiyear ADCP mooring records from locations across the shelf reveal that semidiurnal tides are dominated by theM2 and S2 constituents, with the largest amplitudes on the outer shelf. Throughoutmost of the shelf, tides are clockwise polarized and sheared by stratification, as characteristic near the M2 critical latitude. Interannual variations of the tidal and shear structures on the inner shelf aremainly determined by the stratification-setting Lena River freshwater plume. In all locations,M2 tides are enhanced under sea ice, and therefore changes in the seasonal ice cover may lead to changes in tides and water column structure. The main conclusions of this study are that (i) tides play a comparatively greater role year-round on the outer shelf relative to the inner shelf; (ii) a sea ice reduction will overall decrease the predictability of the currents, especially on the inner shelf; and (iii) the freshwater distribution directly impacts diapycnal mixing by setting the vertical tidal structure. These combined effects imply that future sea ice loss will increase the variability and vertical mixing of freshwater, particularly on the inner shelf, where the Lena River first enters the Laptev Sea.

Description: A large number of quantities have to be measured and processed to determine the atmospheric-state variables, which are the actual measurands, from aircraft-based measurements. A great part of the dependencies between these quantities depends on the aerodynamic state of the aircraft. Aircraft-based meteorological measurements, hence, require in-flight calibration. Most operators of research aircraft perform some kind of calibration, but the schemes used and the degree they are documented greatly vary. The flight maneuvers and calculation methods required, however, are published in a number of partly overlapping and partly contradictory publications. Some methods are only presented as a minor issue in publications mainly focused on atmospheric processes and are therefore hard to find. For an aircraft user, it is hence challenging to either perform or verify a calibration because of missing comprehensive guidance. This lack was stated on occasion of the in-flight calibration of the German research aircraft Polar5 carried out for the field experiment Investigation of Katabatic Winds and Polynyas during Summer (IKAPOS). In the present paper, a comprehensive review of the existing literature on this field and a practical guide to the wind calibration of a research aircraft to be used for turbulent flux measurements are given.

Description: Accurate projections of stratospheric ozone are required because ozone changes affect exposure to ultraviolet radiation and tropospheric climate. Unweighted multimodel ensemble-mean (uMMM) projections from chemistry–climate models (CCMs) are commonly used to project ozone in the twenty-first century, when ozone-depleting substances are expected to decline and greenhouse gases are expected to rise. Here, the authors address the question of whether Antarctic total column ozone projections in October given by the uMMM of CCM simulations can be improved by using a process-oriented multiple diagnostic ensemble regression (MDER) method. This method is based on the correlation between simulated future ozone and selected key processes relevant for stratospheric ozone under present-day conditions. The regression model is built using an algorithm that selects those process-oriented diagnostics that explain a significant fraction of the spread in the projected ozone among the CCMs. The regression model with observed diagnostics is then used to predict future ozone and associated uncertainty. The precision of the authors’ method is tested in a pseudoreality; that is, the prediction is validated against an independent CCM projection used to replace unavailable future observations. The tests show that MDER has higher precision than uMMM, suggesting an improvement in the estimate of future Antarctic ozone. The authors’ method projects that Antarctic total ozone will return to 1980 values at around 2055 with the 95% prediction interval ranging from 2035 to 2080. This reduces the range of return dates across the ensemble of CCMs by about a decade and suggests that the earliest simulated return dates are unlikely.

Description: Evidence is presented for the notion that some contribution to the recent decadal trends observed in the Southern Hemisphere, including the lack of a strong Southern Ocean surface warming, may have originated from longer-term internal centennial variability originating in the Southern Ocean. The existence of such centennial variability is supported by the instrumental sea surface temperatures (SSTs), a multimillennial reconstruction of Tasmanian summer temperatures from tree rings, and a millennial control integration of the Kiel Climate Model (KCM). The model variability was previously shown to be linked to changes in Weddell Sea deep convection. During phases of deep convection the surface Southern Ocean warms, the abyssal Southern Ocean cools, Antarctic sea ice extent retreats, and the low-level atmospheric circulation over the Southern Ocean weakens. After the halt of deep convection the surface Southern Ocean cools, the abyssal Southern Ocean warms, Antarctic sea ice expands, and the low-level atmospheric circulation over the Southern Ocean intensifies, consistent with what has been observed during the recent decades. A strong sensitivity of the time scale to model formulation is noted. In the KCM, the centennial variability is associated with global-average surface air temperature (SAT) changes of the order of a few tenths of a degree per century. The model results thus suggest that internal centennial variability originating in the Southern Ocean should be considered in addition to other internal variability and external forcing when discussing the climate of the twentieth century and projecting that of the twenty-first century.

Description: The wintertime northern annular mode (NAM) at the surface is known to undergo slow intraseasonal variations in association with stratospheric variability, which leads the surface signal by up to several weeks. The relative contributions, however, of potentially relevant stratosphere–troposphere coupling mechanisms are not yet fully understood. In this study the relative roles of (i) the downward effect of the zonal-mean secondary circulation induced by quasigeostrophic (QG) adjustment to stratospheric wave drag and radiative damping and (ii) wave drag local to the troposphere are estimated. For this purpose, a spectral tendency equation of the QG zonal-mean zonal wind is derived and used, in a first step, to obtain the external mechanical forcing that, in the QG framework, drives exactly the observed stratospheric and tropospheric daily NAM. In a second step, the equation is then integrated in time to reconstruct the daily NAM, but with the forcing restricted to either stratospheric or tropospheric levels, each case leaving a characteristic NAM surface signal. The relative roles of the above-mentioned mechanisms are found to be of similar quantitative importance, but to differ in a qualitative sense. The downward effect of stratospheric QG adjustment is responsible for the initiation of the NAM surface signal, whereas subsequently local tropospheric wave drag actively maintains and persists the signal over several weeks. Furthermore, the downward effect of QG adjustment to stratospheric radiative damping is shown to have only a minor impact, compared to that from stratospheric wave drag. The robustness of these conclusions is demonstrated by a sensitivity study with respect to various model parameters.

Description: The Agulhas Current plays a crucial role in the thermohaline circulation through its leakage into the South Atlantic. Under both past and present climates, the trade winds and westerlies could have the ability to modulate the amount of Indian-Atlantic inflow. Compelling arguments have been put forward suggesting that trade winds alone have little impact on the magnitude of Agulhas leakage. Here, employing three ocean models for robust analysis – a global coarse resolution, a regional eddy-permitting and a nested high-resolution eddy-resolving configuration – and systematically altering the position and intensity of the westerly wind belt in a series of sensitivity experiments, it is shown that the westerlies, in particular their intensity, control the leakage. Leakage responds proportionally to the westerlies intensity up to a certain point. Beyond this, through the adjustment of the large-scale circulation, energetic interactions occur between the Agulhas Return Current and the Antarctic Circumpolar Current that result in a state where leakage no longer increases. This adjustment takes place within 1 to 2 decades. Contrary to previous assertions, our results further show that an equatorward (poleward) shift in westerlies increases (decreases) leakage. This occurs due to the redistribution of momentum input by the winds. It is concluded that the reported present-day leakage increase could therefore reflect an unadjusted oceanic response mainly to the strengthening westerlies over the last few decades.

Description: For the first time in serveral years, the El Nino-Southern Oscillation did not dominate regional climate conditions around the globe. A weak La Ni a dissipated to ENSOneutral conditions by spring, and while El Nino appeared to be emerging during summer, this phase never fully developed as sea surface temperatures in the eastern conditions. Nevertheless, other large-scale climate patterns and extreme weather events impacted various regions during the year. A negative phase of the Arctic Oscillation from mid-January to early February contributed to frigid conditions in parts of northern Africa, eastern Europe, and western Asia. A lack of rain during the 2012 wet season led to the worst drought in at least the past three decades for northeastern Brazil. Central North America also experienced one of its most severe droughts on record. The Caribbean observed a very wet dry season and it was the Sahel's wettest rainy season in 50 years. Overall, the 2012 average temperature across global land and ocean surfaces ranked among the 10 warmest years on record. The global land surface temperature alone was also among the 10 warmest on record. In the upper atmosphere, the average stratospheric temperature was record or near-record cold, depending on the dataset. After a 30-year warming trend from 1970 to 1999 for global sea surface temperatures, the period 2000-12 had little further trend. This may be linked to the prevalence of La Ni a-like conditions during the 21st century. Heat content in the upper 700 m of the ocean remained near record high levels in 2012. Net increases from 2011 to 2012 were observed at 700-m to 2000-m depth and even in the abyssal ocean below. Following sharp decreases in to the effects of La Ni a, sea levels rebounded to reach records highs in 2012. The increased hydrological cycle seen in recent years continued, with more evaporation in drier locations and more precipitation in rainy areas. In a pattern that has held since 2004, salty areas of the ocean surfaces and subsurfaces were anomalously salty on average, while fresher areas were anomalously fresh. Global tropical cyclone activity during 2012 was near average, with a total of 84 storms compared with the 1981-2010 average of 89. Similar to 2010 and 2011, the North Atlantic was the only hurricane basin that experienced above-normal activity. In this basin, Sandy brought devastation to Cuba and parts of the eastern North American seaboard. All other basins experienced either near-or below-normal tropical cyclone activity. Only three tropical cyclones reached Category 5 intensity-all in Bopha became the only storm in the historical record to produce winds greater than 130 kt south of 7 N. It was also the costliest storm to affect the Philippines and killed more than 1000 residents. Minimum Arctic sea ice extent in September and Northern Hemisphere snow cover extent in June both reached new record lows. June snow cover extent is now declining at a faster rate (-17.6% per decade) than September sea ice extent (-13.0% per decade). Permafrost temperatures reached record high values in northernmost Alaska. A new melt extent record occurred on 11-12 July on the Greenland ice sheet; 97% of the ice sheet showed some form of melt, four times greater than the average melt for this time of year. The climate in Antarctica was relatively stable overall. The largest maximum sea ice extent since records begain in 1978 was observed in September 2012. In the stratosphere, warm air led to the second smallest ozone hole in the past two decades. Even so, the springtime ozone layer above Antarctica likely will not return to its early 1980s state until about 2060. Following a slight decline associated with the global 2 emissions from fossil fuel combustion and cement production reached a record 9.5 +/- 0.5 Pg C in 2011 and a new record of 9.7 +/- 0.5 Pg C is estimated for 2012. Atmospheric CO2 concentrations increased by 2.1 ppm in 2012, to 392.6 ppm. In spring 2012, 2 concentration exceeded 400 ppm at 7 of the 13 Arctic observation sites. Globally, other greenhouse gases including methane and nitrous oxide also continued to rise in concentration and the combined effect now represents a 32% increase in radiative forcing over a 1990 baseline. Concentrations of most ozone depleting substances continued to fall.

Description: Mesoscale anticyclonic eddies in the Irminger Sea are observed using a mooring and a glider. Between 2002 and 2009, the mooring observed 53 anticyclones. Using a kinematic model, objective estimates of eddy length scales and velocity structure are made for 16 eddies. Anticyclones had a mean core diameter of 12 km, and their mean peak observed azimuthal speed was 0.1 m s(-1). They had core salinities and potential temperatures of 34.91-34.98 and 4.488-5.34 degrees C, respectively, making them warm and salty features. These properties represent a typical salinity anomaly of 0.03 and a temperature anomaly of 0.28 degrees C from noneddy values. All eddies had small (〈〈 1) Rossby numbers. In 2006, the glider observed two anticyclones having diameters of about 20 km and peak azimuthal speeds of about 0.3 m s(-1). Similar salinity anomalies were detected throughout the Irminger Sea by floats profiling in anticyclones. Two formation regions for the eddies are identified: one to the west of the Reykjanes Ridge and the other off the East Greenland Irminger Current near Cape Farewell close to the mooring. Observations indicate that eddies formed in the former region are larger than eddies observed at the mooring. A clear increase in eddy salinity is observed between 2002 and 2009. The observed breakup of these eddies in winter implies that they are a source of salt for the central gyre. The anticyclones are similar to those found in both the Labrador Sea and Norwegian Sea, making them a ubiquitous feature of the subpolar North Atlantic basins.

Description: In this article we address the causes of the large-scale tropical sea level pressure (SLP) changes during climate change. The analysis we present is based on model simulations, observed trends and the seasonal cycle. In all three cases the regional changes of tropospheric temperature (Ttropos) and SLP are strongly related to each other (considerably stronger than (sea) surface temperature and SLP). This relationship basically follows the Bjerknes Circulation Theorem, with relatively low regional SLP where we have relatively high Ttropos and vice versa. A simple physical model suggests a tropical SLP response to horizontally inhomogeneous warming in the tropical Ttropos, with a sensitivity coefficient of about -1.7 hPa/K. This relationship explains a large fraction of observed and predicted changes in the tropical SLP. It is shown that in climate change model simulations the tropospheric land-sea warming contrast is the most significant structure in the regional Ttropos changes relative to the tropical mean changes. Since the land-sea warming contrast exists in the absent of any atmospheric circulation changes it can be argued that the large-scale response of tropical SLP changes is to first order a response to the tropical land-sea warming contrast. Further, as land-sea warming contrast is mostly available moisture dependent, the models predict a stronger warming and decreasing SLP in the drier regions from South America to Africa and a weaker warming and increasing SLP over the wetter Indo-Pacific warm pool region. This suggests an increase in the potential for deep convection conditions over the Atlantic Sector and a decrease over the Indo-Pacific warm pool region in the future.

Description: Variations in eastern Indian Ocean upper-ocean thermal properties are assessed for the period 1970–2004, with a particular focus on asymmetric features related to opposite phases of Indian Ocean Dipole events, using high-resolution ocean model hindcasts. Sensitivity experiments, where atmospheric forcing variability is restricted to the Indian or Pacific Ocean only, support the interpretation of forcing mechanisms for large-scale asymmetric behavior in eastern Indian Ocean variability. Years are classified according to eastern Indian Ocean subsurface heat content (HC) as proxy of thermocline variations. Years characterized by anomalous low HC feature a zonal gradient in upper-ocean properties near the equator, while high events have a meridional gradient from the tropics into the subtropics. The spatial and temporal characteristics of the seasonal evolution of HC anomalies for the two cases is distinct, as is the relative contribution from Indian Ocean atmospheric forcing versus remote influences from Pacific wind forcing: low events develop rapidly during austral winter/spring in response to Indian Ocean wind forcing associated with an enhanced southeasterly monsoon driving coastal upwelling and a shoaling thermocline in the east; in contrast, formation of anomalous high eastern Indian Ocean HC is more gradual, with anomalies earlier in the year expanding from the Indonesian Throughflow (ITF) region, initiated by remote Pacific wind forcing and transmitted through the ITF via coastal wave dynamics. Implications for seasonal predictions arise with high HC events offering extended lead times for predicting thermocline variations and upper-ocean properties across the eastern Indian Ocean.

Description: European Polar Low Working Group This workshop summarized the current state of PL research in the Arctic and Antarctic. A couple of related projects are in the planning phase or already funded. The creation of a PL database for the Norwegian Sea in the frame of the Sea Surface Temperature and Altimeter Synergy (STARS) project (http://projects.met.no stars) will provide a valuable resource for future research and, potentially, predictability improvements. The maintenance of this database and the creation of similar databases for other polar areas including satellite and NWP data are highly recommended. There is also a need for free and timely access to satellite data, in particular to SAR data to fill the gap caused by the mission end of Envisat. With the increasing resolution of climate models, mesoscale processes such as polar MCs will have to be considered in international research programs such as the World Climate Research Programme (WCRP) Polar Climate Predictability Initiative and the World Weather Research Programme (WWRP) Polar Predictability Project.

Description: The mechanisms involved in setting the annual cycle of the Florida Current transport are revisited using an adjoint model approach. Adjoint sensitivities of the Florida Current transport to wind stress reproduce a realistic seasonal cycle with an amplitude of ~1.2 Sv (1 Sv ≡ 106 m3 s−1). The annual cycle is predominantly determined by wind stress forcing and related coastal upwelling (downwelling) north of the Florida Strait along the shelf off the North American coast. Fast barotropic waves propagate these anomalies southward and reach the Florida Strait within a month, causing an amplitude of ~1 Sv. Long baroclinic planetary Rossby waves originating from the interior are responsible for an amplitude of ~0.8 Sv but have a different phase. The sensitivities corresponding to the first baroclinic mode propagate westward and are highly influenced by topography. Considerable sensitivities are only found west of the Mid-Atlantic Ridge, with maximum values at the western shelf edge. The second baroclinic mode also has an impact on the Florida Current variability, but only when a mean flow is present. A second-mode wave train propagates southwestward from the ocean bottom on the western side of the Mid-Atlantic Ridge between ~36° and 46°N and at Flemish Cap, where the mean flow interacts with topography, to the surface. Other processes such as baroclinic waves along the shelf and local forcing within the Florida Strait are of minor importance.

Description: The residual effect of surface gravity waves on mean flows in the upper ocean is investigated using thickness weighted mean (TWM) theory applied in a vertically Lagrangian and horizontally Eulerian coordinate system. Depth-dependent equations for the conservation of volume, momentum, and energy are derived. These equations allow for (i) finite amplitude fluid motions, (ii) the horizontal divergence of currents and (iii) a concise treatment of both the kinematic and viscous boundary conditions at the sea surface. Under the assumptions of steady and monochromatic waves and a uniform turbulent viscosity, the TWM momentum equations are used to illustrate the pressure- and viscosity-induced momentum fluxes through the surface that are implicit in previous studies of the wave-induced modification of the classical Ekman spiral problem. The TWM approach clarifies, in particular, the surface momentum flux associated with the so-called virtual wave stress of Longuet-Higgins. Overall the TWM framework can be regarded as an alternative to the three-dimensional Lagrangian mean framework of Pierson. Moreover the TWM framework can be used to include the residual effect of surface waves in large-scale circulation models. In specific models that carry the TWM velocity appropriate for advecting tracers as their velocity variable, the turbulent viscosity term should be modified so that the viscosity acts only on the Eulerian mean velocity.

Description: Many aspects of barnacle body form are known to be developmentally plastic. Perhaps the most striking examples of such plasticity occur in their feeding legs and unusually long penises, the sizes and shapes of which can change dramatically and adaptively with changes in conspecific density and local water flow conditions. However, whether variation in overall appendage form is mirrored by structural responses in cuticle and muscle is not known. In order to determine how structural variation underlies phenotypic plasticity in barnacle appendages, we examined barnacles occurring at low and high population densities from one wave-protected and one wave-exposed site. We used histological sectioning and fluorescence microscopy of feeding legs and penises to compare cuticle thickness, muscle thickness, and muscle organization, and artificial penis inflation to compare penis extensibility. We observed striking differences in cuticle thickness, muscle thickness, and muscle organization between sites that differed in water velocity, but we found no clear differences associated with variation in conspecific density. Penis extensibility also did not differ consistently between sites. These results are consistent with an adaptive explanation for much of the remarkable and complex variation in barnacle feeding leg and penis morphology among sites that differ in water velocity.

Description: The tropical Atlantic wind response to El Niño forcing is robust, with weakened northeast trade winds north of the equator and strengthened southeast trade winds along and south of the equator. However, the relationship between sea surface temperature (SST) anomalies in the eastern equatorial Pacific and Atlantic is inconsistent, with El Niño events followed sometimes by warm and other times by cold boreal summer anomalies in the Atlantic cold tongue region. Using observational data and a hindcast simulation of the Nucleus for European Modeling of the Ocean (NEMO) global model at 0.5° resolution (NEMO-ORCA05), this inconsistent SST relationship is shown to be at least partly attributable to a delayed negative feedback in the tropical Atlantic that is active in years with a warm or neutral response in the eastern equatorial Atlantic. In these years, the boreal spring warming in the northern tropical Atlantic that is a typical response to El Niño is pronounced, setting up a strong meridional SST gradient. This leads to a negative wind stress curl anomaly to the north of the equator that generates downwelling Rossby waves. When these waves reach the western boundary, they are reflected into downwelling equatorial Kelvin waves that reach the cold tongue region in late boreal summer to counteract the initial cooling that is due to the boreal winter wind stress response to El Niño. In contrast, this initial cooling persists or is amplified in years in which the boreal spring northern tropical Atlantic warming is weak or absent either because of a positive North Atlantic Oscillation (NAO) phase or an early termination of the Pacific El Niño event.

Description: Initial-value predictability measures the degree to which the initial state can influence predictions. In this paper, the initial-value predictability of six atmosphere–ocean general circulation models in the North Pacific and North Atlantic is quantified and contrasted by analyzing long control integrations with time invariant external conditions. Through the application of analog and multivariate linear regression methodologies, average predictability properties are estimated for forecasts initiated from every state on the control trajectories. For basinwide measures of predictability, the influence of the initial state tends to last for roughly a decade in both basins, but this limit varies widely among the models, especially in the North Atlantic. Within each basin, predictability varies regionally by as much as a factor of 10 for a given model, and the locations of highest predictability are different for each model. Model-to-model variations in predictability are also seen in the behavior of prominent intrinsic basin modes. Predictability is primarily determined by the mean of forecast distributions rather than the spread about the mean. Horizontal propagation plays a large role in the evolution of these signals and is therefore a key factor in differentiating the predictability of the various models.

Description: Sudden stratospheric warmings are prominent examples of dynamical wave–mean flow interactions in the Arctic stratosphere during Northern Hemisphere winter. They are characterized by a strong temperature increase on time scales of a few days and a strongly disturbed stratospheric vortex. This work investigates a wide class of supervised learning methods with respect to their ability to classify stratospheric warmings, using temperature anomalies from the Arctic stratosphere and atmospheric forcings such as ENSO, the quasi-biennial oscillation (QBO), and the solar cycle. It is demonstrated that one representative of the supervised learning methods family, namely nonlinear neural networks, is able to reliably classify stratospheric warmings. Within this framework, one can estimate temporal onset, duration, and intensity of stratospheric warming events independently of a particular pressure level. In contrast to classification methods based on the zonal-mean zonal wind, the approach herein distinguishes major, minor, and final warmings. Instead of a binary measure, it provides continuous conditional probabilities for each warming event representing the amount of deviation from an undisturbed polar vortex. Additionally, the statistical importance of the atmospheric factors is estimated. It is shown how marginalized probability distributions can give insights into the interrelationships between external factors. This approach is applied to 40-yr and interim ECMWF (ERA-40/ERA-Interim) and NCEP–NCAR reanalysis data for the period from 1958 through 2010.

Description: To analyze the probability density distributions of surface turbulent heat fluxes, the authors apply the twoparametric modified Fisher–Tippett (MFT) distribution to the sensible and latent turbulent heat fluxes recomputed from 6-hourly NCEP–NCAR reanalysis state variables for the period from 1948 to 2008. They derived the mean climatology and seasonal cycle of the location and scale parameters of the MFT distribution. Analysis of the parameters of probability distributions identified the areas where similar surface turbulent fluxes are determined by the very different shape of probability density functions. Estimated extreme turbulent heat fluxes amount to 1500–2000 W m22 (for the 99th percentile) and can exceed 2000 W m22 for higher percentiles in the subpolar latitudes and western boundary current regions. Analysis of linear trends and interannual variability in the mean and extreme fluxes shows that the strongest trends in extreme fluxes (more than 15 W m22 decade21) in the western boundary current regions are associated with the changes in the shape of distribution. In many regions changes in extreme fluxes may be different from those for the mean fluxes at interannual and decadal time scales. The correlation between interannual variability of themean and extreme fluxes is relatively low in the tropics, the SouthernOcean, and the Kuroshio Extension region.Analysis of probability distributions in turbulent fluxes has also been used in assessing the impact of sampling errors in theVoluntaryObserving Ship (VOS)-based surface flux climatologies, allowed for the estimation of the impact of sampling in extreme fluxes. Although sampling does not have a visible systematic effect onmean fluxes, sampling uncertainties result in the underestimation of extreme flux values exceeding 100 W m22 in poorly sampled regions.

Description: Antarctic Intermediate Water (AAIW) is a dominant Southern Hemisphere water mass that spreads from its formation regions just north of the Antarctic Circumpolar Current (ACC) to at least 20°S in all oceans. This study uses an isopycnal climatology constructed from Argo conductivity–temperature–depth (CTD) profile data to define the current state of the AAIW salinity minimum (its core) and thence compute anomalies of AAIW core pressure, potential temperature, salinity, and potential density since the mid-1970s from ship-based CTD profiles. The results are used to calculate maps of temporal property trends at the AAIW core, where statistically significant strong circumpolar shoaling (30–50 dbar decade−1), warming (0.05°–0.15°C decade−1), and density reductions [up to −0.03 (kg m−3) decade−1] are found. These trends are strongest just north of the ACC in the southeast Pacific and Atlantic Oceans and decrease equatorward. Salinity trends are generally small, with their sign varying regionally. Bottle data are used to extend the AAIW core potential temperature anomaly analysis back to 1925 in the Atlantic and to ~1960 elsewhere. The modern warm AAIW core conditions appear largely unprecedented in the historical record: biennially and zonally binned median AAIW core potential temperatures within each ocean basin are, with the notable exception of the subtropical South Atlantic in the 1950s–70s, 0.2–1°C colder than modern values. Zonally averaged sea surface temperature anomalies around the AAIW formation latitudes in each ocean and sectoral southern annular mode indices are used to put the AAIW core property trends and variations into context.

Description: Anthropogenic emissions of greenhouse gases could lead to undesirable effects on oceans in coming centuries. Drawing on recommendations published by the German Advisory Council on Global Change, levels of unacceptable global marine change (so-called guardrails) are defined in terms of global mean temperature, sea level rise, and ocean acidification. A global-mean climate model [the Aggregated Carbon Cycle, Atmospheric Chemistry and Climate Model (ACC2)] is coupled with an economic module [taken from the Dynamic Integrated Climate–Economy Model (DICE)] to conduct a cost-effectiveness analysis to derive CO2 emission pathways that both minimize abatement costs and are compatible with these guardrails. Additionally, the “tolerable windows approach” is used to calculate a range of CO2 emissions paths that obey the guardrails as well as a restriction on mitigation rate. Prospects of meeting the global mean temperature change guardrail (2° and 0.2°C decade−1 relative to preindustrial) depend strongly on assumed values for climate sensitivity: at climate sensitivities 〉3°C the guardrail cannot be attained under any CO2 emissions reduction strategy without mitigation of non-CO2 greenhouse gases. The ocean acidification guardrail (0.2 unit pH decline relative to preindustrial) is less restrictive than the absolute temperature guardrail at climate sensitivities 〉2.5°C but becomes more constraining at lower climate sensitivities. The sea level rise and rate of rise guardrails (1 m and 5 cm decade−1) are substantially less stringent for ice sheet sensitivities derived in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report, but they may already be committed to violation if ice sheet sensitivities consistent with semiempirical sea level rise projections are assumed.

Description: Historical hydrographic data (1940s–2010) show a distinct cross-slope difference of the lower halocline water (LHW) over the Laptev Sea continental margins. Over the slope, the LHW is on average warmer and saltier by 0.2°C and 0.5 psu, respectively, relative to the off-slope LHW. The LHW temperature time series constructed from the on-slope historical records are related to the temperature of the Atlantic Water (AW) boundary current transporting warm water from the North Atlantic Ocean. In contrast, the on-slope LHW salinity is linked to the sea ice and wind forcing over the potential upstream source region in the Barents and northern Kara Seas, as also indicated by hydrodynamic model results. Over the Laptev Sea continental margin, saltier LHW favors weaker salinity stratification that, in turn, contributes to enhanced vertical mixing with underlying AW.

Description: The equatorial deep jets (EDJ) are a striking feature of the equatorial ocean circulation. In the Atlantic Ocean, the EDJ are associated with a vertical scale of between 300 and 700 m, a time scale of roughly 4.5 years and upward energy propagation to the surface. It has been found that the meridional width of the EDJ is roughly 1.5 times larger than expected based on their vertical scale. Here we use a shallow water model for a high order baroclinic vertical normal mode to argue that mixing of momentum along isopycnals can explain the enhanced width. A lateral eddy viscosity of 300 m2 s−1 10 is found to be sufficient to account for the width implied by observations.

Description: The suggestion is advanced that the remarkably low static stability of Antarctic surface waters may arise from a feedback loop involving global deep-water temperatures. If deep-water temperatures are too warm, this promotes Antarctic convection, thereby strengthening the inflow of Antarctic Bottom Water into the ocean interior and cooling the deep ocean. If deep waters are too cold, this promotes Antarctic stratification allowing the deep ocean to warm because of the input of North Atlantic Deep Water. A steady-state deep-water temperature is achieved such that the Antarctic surface can barely undergo convection. A two-box model is used to illustrate this feedback loop in its simplest expression and to develop basic concepts, such as the bounds on the operation of this loop. The model illustrates the possible dominating influence of Antarctic upwelling rate and Antarctic freshwater balance on global deep-water temperatures.

Description: The importance of decadal climate variability (DCV) research is being increasingly recognized, including by the World Climate Research Program (WCRP) and the Intergovernmental Panel on Climate Change (IPCC). An improved understanding of DCV is very important because stakeholders and policymakers want to know the likely climate trajectory for the coming decades for applications to water resources, agriculture, energy, and infrastructure development. Responding to this demand, many climate modeling groups in the United States, Europe, Japan, and elsewhere are gearing up to assess the potential for decadal climate predictions. The magnitudes of regional DCV often exceed those associated with the trends resulting from anthropogenic changes. Therefore, differentiating between the two is also very important for planning, implementation, and national and international treaties.

Description: Simple idealized layered models and primitive equation models show that the meridional gradient of the zonally averaged pressure has no direct relation with the meridional flow. This demonstrates a contradiction in an often-used parameterization in zonally averaged models. The failure of this parameterization reflects the inconsistency between the model of Stommel and Arons and the box model of Stommel, as previously pointed out by Straub. A new closure is proposed. The ocean is divided in two dynamically different regimes: a narrow western boundary layer and an interior ocean; zonally averaged quantities over these regions are considered. In the averaged equations three unknowns appear: the interior zonal pressure difference Delta p(i), the zonal pressure difference Delta p(b) of the boundary layer, and the zonal velocity us at the interface between the two regions. Here Delta p(i) is parameterized using a frictionless vorticity balance, Delta p(b), by the difference of the mean pressure in the interior and western boundary, and u(delta) by the mean zonal velocity of the western boundary layer. Zonally resolved models, a layer model, and a primitive equation model validate the new parameterization by comparing with the respective zonally averaged counterparts. It turns out that the zonally averaged models reproduce well the buoyancy distribution and the meridional flow in the zonally resolved model versions with respect to the mean and time changes.

Description: Continuous estimates of the oceanic meridional heat transport in the Atlantic are derived from the Rapid Climate Change–Meridional Overturning Circulation (MOC) and Heatflux Array (RAPID–MOCHA) observing system deployed along 26.5°N, for the period from April 2004 to October 2007. The basinwide meridional heat transport (MHT) is derived by combining temperature transports (relative to a common reference) from 1) the Gulf Stream in the Straits of Florida; 2) the western boundary region offshore of Abaco, Bahamas; 3) the Ekman layer [derived from Quick Scatterometer (QuikSCAT) wind stresses]; and 4) the interior ocean monitored by “endpoint” dynamic height moorings. The interior eddy heat transport arising from spatial covariance of the velocity and temperature fields is estimated independently from repeat hydrographic and expendable bathythermograph (XBT) sections and can also be approximated by the array. The results for the 3.5 yr of data thus far available show a mean MHT of 1.33 ± 0.40 PW for 10-day-averaged estimates, on which time scale a basinwide mass balance can be reasonably assumed. The associated MOC strength and variability is 18.5 ± 4.9 Sv (1 Sv ≡ 106 m3 s−1). The continuous heat transport estimates range from a minimum of 0.2 to a maximum of 2.5 PW, with approximately half of the variance caused by Ekman transport changes and half caused by changes in the geostrophic circulation. The data suggest a seasonal cycle of the MHT with a maximum in summer (July–September) and minimum in late winter (March–April), with an annual range of 0.6 PW. A breakdown of the MHT into “overturning” and “gyre” components shows that the overturning component carries 88% of the total heat transport. The overall uncertainty of the annual mean MHT for the 3.5-yr record is 0.14 PW or about 10% of the mean value.

Description: The water mass structure of the Arctic Ocean is remarkable, for its intermediate (depth range ~150–900 m) layer is filled with warm (temperature 〉0°C) and salty water of Atlantic origin (usually called the Atlantic Water, AW). This water is carried into and through the Arctic Ocean by the pan-Arctic boundary current, which moves cyclonically along the basins’ margins (Fig. 1). This system provides the largest input of water, heat, and salt into the Arctic Ocean; the total quantity of heat is substantial, enough to melt the Arctic sea ice cover several times over. By utilizing an extensive archive of recently collected observational data, this study provides a cohesive picture of recent large-scale changes in the AW layer of the Arctic Ocean. These recent observations show the warm pulse of AW that entered the Arctic Ocean in the early 1990s finally reached the Canada Basin during the 2000s. The second warm pulse that entered the Arctic Ocean in the mid-2000s has moved through the Eurasian Basin and is en route downstream. One of the most intriguing results of these observations is the realization of the possibility of uptake of anomalous AW heat by overlying layers, with possible implications for an already-reduced Arctic ice cover.

Description: Vertical mixing in the bottom boundary layer and pycnocline of the Laptev Sea is evaluated from a rapidly sampled 12-h time series of microstructure temperature, conductivity, and shear observations collected under 100% sea ice during October 2008. The bottom boundary turbulent kinetic energy dissipation was observed to be enhanced (ϵ ∼ 10−4 W m−3) beyond background levels (ϵ ∼ 10−6 W m−3), extending up to 10 m above the seabed when simulated tidal currents were directed on slope. Upward heat fluxes into the halocline-class waters along the Laptev Sea seabed peaked at ∼4–8 W m−2, averaging out to ∼2 W m−2 over the 12-h sampling period. In the Laptev Sea pycnocline, an isolated 2-h episode of intense dissipation (ϵ ∼ 10−3 W m−3) and vertical diffusivities was observed that was not due to a localized wind event. Observations from an acoustic Doppler current meter moored in the central Laptev Sea near the M2 critical latitude are consistent with a previous model in which mixing episodes are driven by an enhancement of the pycnocline shear resulting from the alignment of the rotating pycnocline shear vector with the under-ice stress vector. Upward cross-pycnocline heat fluxes from the Arctic halocline peaked at ∼54 W m−2, resulting in a 12-h average of ∼12 W m−2. These results highlight the intermittent nature of Arctic shelf sea mixing processes and how these processes can impact the transformation of Arctic Ocean water masses. The observations also clearly demonstrate that absence or presence of sea ice profoundly affects the availability of near-inertial kinetic energy to drive vertical mixing on the Arctic shelves.

Description: The eddy correlation technique is rapidly becoming an established method for resolving dissolved oxygen fluxes in natural aquatic systems. This direct and noninvasive determination of oxygen fluxes close to the sediment by simultaneously measuring the velocity and the dissolved oxygen fluctuations has considerable advantages compared to traditional methods. This paper describes the measurement principle and analyzes the spatial and temporal scales of those fluctuations as a function of turbulence levels. The magnitudes and spectral structure of the expected fluctuations provide the required sensor specifications and define practical boundary conditions for the eddy correlation instrumentation and its deployment. In addition, data analysis and spectral corrections are proposed for the usual nonideal conditions, such as the time shift between the sensor pair and the limited frequency response of the oxygen sensor. The consistency of the eddy correlation measurements in a riverine reservoir has been confirmed—observing a night–day transition from oxygen respiration to net oxygen production, ranging from −20 to +5 mmol m−2 day−1—by comparing two physically independent, eddy correlation instruments deployed side by side. The natural variability of the fluctuations calls for at least ∼1 h of flux data record to achieve a relative accuracy of better than ∼20%. Although various aspects still need improvement, eddy correlation is seen as a promising and soon-to-be widely applied method in natural waters.

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

Description: Changes in the ventilation of the oxygen minimum zone (OMZ) of the tropical North Atlantic are studied using oceanographic data from 18 research cruises carried out between 28.5° and 23°W during 1999–2008 as well as historical data referring to the period 1972–85. In the core of the OMZ at about 400-m depth, a highly significant oxygen decrease of about 15 μmol kg−1 is found between the two periods. During the same time interval, the salinity at the oxygen minimum increased by about 0.1. Above the core of the OMZ, within the central water layer, oxygen decreased too, but salinity changed only slightly or even decreased. The scatter in the local oxygen–salinity relations decreased from the earlier to the later period suggesting a reduced filamentation due to mesoscale eddies and/or zonal jets acting on the background gradients. Here it is suggested that latitudinally alternating zonal jets with observed amplitudes of a few centimeters per second in the depth range of the OMZ contribute to the ventilation of the OMZ. A conceptual model of the ventilation of the OMZ is used to corroborate the hypothesis that changes in the strength of zonal jets affect mean oxygen levels in the OMZ. According to the model, a weakening of zonal jets, which is in general agreement with observed hydrographic evidences, is associated with a reduction of the mean oxygen levels that could significantly contribute to the observed deoxygenation of the North Atlantic OMZ.

Description: The Atlantic meridional overturning circulation (AMOC) makes the strongest oceanic contribution to the meridional redistribution of heat. Here, an observation-based, forty-eight-month-long time series of the vertical structure and strength of the AMOC at 26.5°N is presented. From April 2004 to April 2008 the AMOC had a mean strength of 18.7 ±2.1 Sv with fluctuations of 4.8 Sv rms. The best guess of the peak-to-peak amplitude of the AMOC seasonal cycle is 6.7 Sv, with a maximum strength in autumn and a minimum in spring. While seasonality in the AMOC was commonly thought to be dominated by the northward Ekman transport, this study reveals that fluctuations of the geostrophic mid-ocean and Gulf Stream transports of 2.2 Sv and 1.7 Sv rms, respectively, are substantially larger than those of the Ekman component (1.2 Sv rms). A simple model based on linear dynamics suggests that the seasonal cycle is dominated by wind stress curl forcing at the eastern boundary of the Atlantic. Seasonal geostrophic AMOC anomalies might represent an important and previously underestimated component of meridional transport and storage of heat in the subtropical North Atlantic. There is evidence that the seasonal cycle observed here is representative of much longer intervals. Previously, hydrographic snapshot estimates between 1957 and 2004 had suggested a long-term decline of the AMOC by 8 Sv. This study suggests that aliasing of seasonal AMOC anomalies might have accounted for a large part of the inferred slowdown.

Description: Analysis of modern and historical observations demonstrates that the temperature of the intermediate-depth (150–900 m) Atlantic water (AW) of the Arctic Ocean has increased in recent decades. The AW warming has been uneven in time; a local 1°C maximum was observed in the mid-1990s, followed by an intervening minimum and an additional warming that culminated in 2007 with temperatures higher than in the 1990s by 0.24°C. Relative to climatology from all data prior to 1999, the most extreme 2007 temperature anomalies of up to 1°C and higher were observed in the Eurasian and Makarov Basins. The AW warming was associated with a substantial (up to 75–90 m) shoaling of the upper AW boundary in the central Arctic Ocean and weakening of the Eurasian Basin upper-ocean stratification. Taken together, these observations suggest that the changes in the Eurasian Basin facilitated greater upward transfer of AW heat to the ocean surface layer. Available limited observations and results from a 1D ocean column model support this surmised upward spread of AW heat through the Eurasian Basin halocline. Experiments with a 3D coupled ice–ocean model in turn suggest a loss of 28–35 cm of ice thickness after 50 yr in response to the 0.5 W m−2 increase in AW ocean heat flux suggested by the 1D model. This amount of thinning is comparable to the 29 cm of ice thickness loss due to local atmospheric thermodynamic forcing estimated from observations of fast-ice thickness decline. The implication is that AW warming helped precondition the polar ice cap for the extreme ice loss observed in recent years.

Description: The 20th century Northern Hemisphere surface climate exhibits a long-term warming trend, largely caused by anthropogenic forcing, and natural decadal climate variability superimposed on it. This study addresses the possible origin and strength of internal decadal climate variability in the Northern Hemisphere during the recent decades. We present results from a set of climate model simulations that suggest natural internal multidecadal climate variability in the North Atlantic-Arctic Sector could have considerably contributed to the Northern Hemisphere surface warming since 1980. Although covering only a few percent of the earth’s surface, the Arctic may have provided the largest share in this. It is hypothesized that a stronger Meridional Overturning Circulation in the Atlantic and the associated increase in northward heat transport enhanced the heat loss from the ocean to the atmosphere in the North Atlantic region, and especially in the North Atlantic portion of the Arctic due to anomalously strong sea ice melt. The model results stress the potential importance of natural internal multidecadal variability originating in the North Atlantic-Arctic Sector in generating inter-decadal climate changes not only on a regional, but possibly also on a hemispheric and even global scale.

Description: The Vasco—Cirene field experiment, in January—February 2007, targeted the Seychelles—Chagos thermocline ridge (SCTR) region, with the main purpose of investigating Madden—Julian Oscillation (MJO)-related SST events. The Validation of the Aeroclipper System under Convective Occurrences (Vasco) experiment (Duvel et al. 2009) and Cirene cruise were designed to provide complementary views of air—sea interaction in the SCTR region. While meteorological balloons were deployed from the Seychelles as a part of Vasco, the Research Vessel (R/V) Suroît was cruising the SCTR region as a part of Cirene. more: The Vasco—Cirene program explores how strong air—sea interactions promoted by the shallow thermocline and high sea surface temperature in the Seychelles—Chagos thermocline ridge results in marked variability at synoptic, intraseasonal, and interannual time scales. The Cirene oceanographic cruise collected oceanic, atmospheric, and air—sea flux observations in this region in January—February 2007. The contemporaneous Vasco field experiment complemented these measurements with balloon deployments from the Seychelles. Cirene also contributed to the development of the Indian Ocean observing system via deployment of a mooring and 12 Argo profilers. Unusual conditions prevailed in the Indian Ocean during January and February 2007, following the Indian Ocean dipole climate anomaly of late 2006. Cirene measurements show that the Seychelles—Chagos thermocline ridge had higher-than-usual heat content with subsurface anomalies up to 7°C. The ocean surface was warmer and fresher than average, and unusual eastward currents prevailed down to 800 m. These anomalous conditions had a major impact on tuna fishing in early 2007. Our dataset also sampled the genesis and maturation of Tropical Cyclone Dora, including high surface temperatures and a strong diurnal cycle before the cyclone, followed by a 1.5°C cooling over 10 days. Balloonborne instruments sampled the surface and boundary layer dynamics of Dora. We observed small-scale structures like dry-air layers in the atmosphere and diurnal warm layers in the near-surface ocean. The Cirene data will quantify the impact of these finescale features on the upper-ocean heat budget and atmospheric deep convection.

Description: The western Pacific subtropical high (WPSH) is closely related to Asian climate. Previous examination of changes in the WPSH found a westward extension since the late 1970s, which has contributed to the inter-decadal transition of East Asian climate. The reason for the westward extension is unknown, however. The present study suggests that this significant change of WPSH is partly due to the atmosphere's response to the observed Indian Ocean-western Pacific (IWP) warming. Coordinated by a European Union's Sixth Framework Programme, Understanding the Dynamics of the Coupled Climate System (DYNAMITE), five AGCMs were forced by identical idealized sea surface temperature patterns representative of the IWP warming and cooling. The results of these numerical experiments suggest that the negative heating in the central and eastern tropical Pacific and increased convective heating in the equatorial Indian Ocean/ Maritime Continent associated with IWP warming are in favor of the westward extension of WPSH. The SST changes in IWP influences the Walker circulation, with a subsequent reduction of convections in the tropical central and eastern Pacific, which then forces an ENSO/Gill-type response that modulates the WPSH. The monsoon diabatic heating mechanism proposed by Rodwell and Hoskins plays a secondary reinforcing role in the westward extension of WPSH. The low-level equatorial flank of WPSH is interpreted as a Kelvin response to monsoon condensational heating, while the intensified poleward flow along the western flank of WPSH is in accord with Sverdrup vorticity balance. The IWP warming has led to an expansion of the South Asian high in the upper troposphere, as seen in the reanalysis.

Description: Seasonal reconstructions of the Southern Hemisphere annular mode (SAM) index are derived to extend the record before the reanalysis period, using station sea level pressure (SLP) data as predictors. Two reconstructions using different predictands are obtained: one [Jones and Widmann (JW)] based on the first principal component (PC) of extratropical SLP and the other (Fogt) on the index of Marshall. A regional-based SAM index (Visbeck) is also considered.These predictands agree well post-1979; correlations decline in all seasons except austral summer for the full series starting in 1958. Predictand agreement is strongest in spring and summer; hence agreement between the reconstructions is highest in these seasons. The less zonally symmetric SAM structure in winter and spring influences the strength of the SAM signal over land areas, hence the number of stations included in the reconstructions. Reconstructions from 1865 were, therefore, derived in summer and autumn and from 1905 in winter and spring. This paper examines the skill of each reconstruction by comparison with observations and reanalysis data. Some of the individual peaks in the reconstructions, such as the most recent in austral summer, represent a full hemispheric SAM pattern, while others are caused by regional SLP anomalies over the locations of the predictors. The JW and Fogt reconstructions are of similar quality in summer and autumn, while in winter and spring the Marshall index is better reconstructed by Fogt than the PC index is by JW. In spring and autumn the SAM shows considerable variability prior to recent decades.

Description: Extratropical cyclones and how they may change in a warmer climate have been investigated in detail with a high-resolution version of the ECHAM5 global climate model. A spectral resolution of T213 (63 km) is used for two 32-yr periods at the end of the twentieth and twenty-first centuries and integrated for the Intergovernmental Panel on Climate Change (IPCC) A1B scenario. Extremes of pressure, vorticity, wind, and precipitation associated with the cyclones are investigated and compared with a lower-resolution simulation.Comparison with observations of extreme wind speeds indicates that the model reproduces realistic values. This study also investigates the ability of the model to simulate extratropical cyclones by computing composites of intense storms and contrasting them with the same composites from the 40-yr ECMWF Re-Analysis (ERA-40). Composites of the time evolution of intense cyclones are reproduced with great fidelity; in particular the evolution of central surface pressure is almost exactly replicated, but vorticity, maximum wind speed, and precipitation are higher in the model. Spatial composites also show that the distributions of pressure, winds, and precipitation at different stages of the cyclone life cycle compare well with those from ERA-40, as does the vertical structure. For the twenty-first century, changes in the distribution of storms are very similar to those of previous study. There is a small reduction in the number of cyclones but no significant changes in the extremes of wind and vorticity in both hemispheres. There are larger regional changes in agreement with previous studies. The largest changes are in the total precipitation, where a significant increase is seen. Cumulative precipitation along the tracks of the cyclones increases by some 11% per track, or about twice the increase in global precipitation, while the extreme precipitation is close to the globally averaged increase in column water vapor (some 27%). Regionally, changes in extreme precipitation are even higher because of changes in the storm tracks.

Description: A new, non-flux-corrected, global climate model is introduced, the Kiel Climate Model (KCM), which will be used to study internal climate variability from interannual to millennial time scales and climate predictability of the first and second kind. The version described here is a coarse-resolution version that will be employed in extended-range integrations of several millennia. KCM's performance in the tropical Pacific with respect to mean state, annual cycle, and El Nino-Southern Oscillation (ENSO) is described. Additionally, the tropical Pacific response to global warming is studied.Overall, climate drift in a multicentury control integration is small. However, KCM exhibits an equatorial cold bias at the surface of the order 1 degrees C, while strong warm biases of several degrees are simulated in the eastern tropical Pacific on both sides off the equator, with maxima near the coasts. The annual and semiannual cycles are realistically simulated in the eastern and western equatorial Pacific, respectively. ENSO performance compares favorably to observations with respect to both amplitude and period. An ensemble of eight greenhouse warming simulations was performed, in which the CO2 concentration was increased by 1% yr(-1) until doubling was reached, and stabilized thereafter. Warming of equatorial Pacific sea surface temperature (SST) is, to first order, zonally symmetric and leads to a sharpening of the thermocline. ENSO variability increases because of global warming: during the 30-yr period after CO2 doubling, the ensemble mean standard deviation of Nino-3 SST anomalies is increased by 26% relative to the control, and power in the ENSO band is almost doubled. The increased variability is due to both a strengthened (22%) thermocline feedback and an enhanced (52%) atmospheric sensitivity to SST; both are associated with changes in the basic state. Although variability increases in the mean, there is a large spread among ensemble members and hence a finite probability that in the "model world" no change in ENSO would be observed.

Description: The temporal evolution of the strength of the Atlantic Meridional Overturning Circulation (AMOC) in the subtropical North Atlantic is affected by both remotely forced, basin-scale meridionally coherent, climate-relevant transport anomalies, such as changes in high-latitude deep water formation rates, and locally forced transport anomalies, such as eddies or Rossby waves, possibly associated with small meridional coherence scales, which can be considered as noise. The focus of this paper is on the extent to which local eddies and Rossby waves when impinging on the western boundary of the Atlantic affect the temporal variability of the AMOC at 26.5 degrees N. Continuous estimates of the AMOC at this latitude have been made since April 2004 by combining the Florida Current, Ekman, and midocean transports with the latter obtained from continuous density measurements between the coasts of the Bahamas and Morocco, representing, respectively, the western and eastern boundaries of the Atlantic at this latitude.Within 100 km of the western boundary there is a threefold decrease in sea surface height variability toward the boundary, observed in both dynamic heights from in situ density measurements and altimetric heights. As a consequence, the basinwide zonally integrated upper midocean transport shallower than 1000 m-as observed continuously between April 2004 and October 2006-varies by only 3.0 Sv (1 Sv = 10(6) m(3) s(-1)) RMS. Instead, upper midocean transports integrated from western boundary stations 16, 40, and 500 km offshore to the eastern boundary vary by 3.6, 6.0, and 10.7 Sv RMS, respectively. The reduction in eddy energy toward the western boundary is reproduced in a nonlinear reduced-gravity model suggesting that boundary-trapped waves may account for the observed decline in variability in the coastal zone because they provide a mechanism for the fast equatorward export of transport anomalies associated with eddies impinging on the western boundary. An analytical model of linear Rossby waves suggests a simple scaling for the reduction in thermocline thickness variability toward the boundary. Physically, the reduction in amplitude is understood as along-boundary pressure gradients accelerating the fluid and rapidly propagating pressure anomalies along the boundary. The results suggest that the local eddy field does not dominate upper midocean transport or AMOC variability at 26.5 degrees N on interannual to decadal time scales.

Description: Several recent general circulation model studies discuss the predictability of the Indian Ocean dipole (IOD) mode, suggesting that it is predictable because of coupled ocean-atmosphere interactions in the Indian Ocean. However, it is not clear from these studies how much of the predictability is due to the response to El Nino. It is shown in this note that a simple statistical model that treats the Indian Ocean as a red noise process forced by tropical Pacific SST shows forecast skills comparable to those of recent general circulation model studies. The results also indicate that some of the eastern tropical Indian Ocean SST predictability in recent studies may indeed be beyond the skill of the simple model proposed in this note, indicating that dynamics in the Indian Ocean may have caused this improved predictability in this region. The model further indicates that the IOD index may be the least predictable index of Indian Ocean SST variability. The model is proposed as a null hypothesis for Indian Ocean SST predictions.

Description: A characteristic feature of global warming is the land-sea contrast, with stronger warming over land than over oceans. Recent studies find that this land-sea contrast also exists in equilibrium global change scenarios, and it is caused by differences in the availability of surface moisture over land and oceans. In this study it is illustrated that this land-sea contrast exists also on interannual time scales and that the ocean-land interaction is strongly asymmetric. The land surface temperature is more sensitive to the oceans than the oceans are to the land surface temperature, which is related to the processes causing the land-sea contrast in global warming scenarios. It suggests that the ocean's natural variability and change is leading to variability and change with enhanced magnitudes over the continents, causing much of the longer-time-scale (decadal) global-scale continental climate variability. Model simulations illustrate that continental warming due to anthropogenic forcing (e. g., the warming at the end of the last century or future climate change scenarios) is mostly (80%-90%) indirectly forced by the contemporaneous ocean warming, not directly by local radiative forcing.

Description: Atmospheric pressure observations from the Southern Hemisphere are used to estimate monthly and annually averaged indexes of the southern annular mode (SAM) back to 1884. This analysis groups all relevant observations in the following four regions: one for Antarctica and three in the subtropical zone. Continuous surface pressure observations are available at a number of locations in the subtropical regions since the end of the nineteenth century. However, year-round observations in the subpolar region near the Antarctic continent began only during the 1940-60 period. The shorter Antarctic records seriously compromise the length of a traditionally estimated SAM index. To improve the situation "proxy'' estimates of Antarctic sea level pressure anomalies are provided based on the concept of atmospheric mass conservation poleward of 208S. This allows deriving a longer SAM index back to 1884. Several aspects of the new record, its statistical properties, seasonal trends, and the regional pressure anomaly correlations, are presented.

Description: Statistical analysis of observations (including atmospheric reanalysis and forced ocean model simulations) is used to address two questions: First, does an analogous mechanism to that of El Niño–Southern Oscillation (ENSO) exist in the equatorial Atlantic or Indian Ocean? Second, does the intrinsic variability in these basins matter for ENSO predictability? These questions are addressed by assessing the existence and strength of the Bjerknes and delayed negative feedbacks in each tropical basin, and by fitting conceptual recharge oscillator models, both with and without interactions among the basins. In the equatorial Atlantic the Bjerknes and delayed negative feedbacks exist, although weaker than in the Pacific. Equatorial Atlantic variability is well described by the recharge oscillator model, with an oscillatory mixed ocean dynamics–sea surface temperature (SST) mode present in boreal spring and summer. The dynamics of the tropical Indian Ocean, however, appear to be quite different: no recharge–discharge mechanism is found. Although a positive Bjerknes-like feedback from July to September is found, the role of heat content seems secondary. Results also show that Indian Ocean interaction with ENSO tends to damp the ENSO oscillation and is responsible for a frequency shift to shorter periods. However, the retrospective forecast skill of the conceptual model is hardly improved by explicitly including Indian Ocean SST. The interaction between ENSO and the equatorial Atlantic variability is weaker. However, a feedback from the Atlantic on ENSO appears to exist, which slightly improves the retrospective forecast skill of the conceptual model.

Description: The Agulhas Current system has been analyzed in a nested high-resolution ocean model and compared to observations. The model shows good performance in the western boundary current structure and the transports off the South African coast. This includes the simulation of the northward-flowing Agulhas Undercurrent. It is demonstrated that fluctuations of the Agulhas Current and Undercurrent around 50–70 days are due to Natal pulses and Mozambique eddies propagating downstream. A sensitivity experiment that excludes those upstream perturbations significantly reduces the variability as well as the mean transport of the undercurrent. Although the model simulates undercurrents in the Mozambique Channel and east of Madagascar, there is no direct connection between those and the Agulhas Undercurrent. Virtual float releases demonstrate that topography is effectively blocking the flow toward the north.

Description: For the first time, an analytical theory and a very high-resolution, frontal numerical model, both based on the unsteady, nonlinear, reduced-gravity shallow water equations on a beta plane, have been used to investigate aspects of the migration of homogeneous surface, frontal warm-core eddies on a beta plane. Under the assumption that, initially, such vortices are surface circular anticyclones of paraboloidal shape and having both radial and azimuthal velocities that are linearly dependent on the radial coordinate (i.e., circular pulsons of the first order), approximate analytical expressions are found that describe the nonstationary trajectories of their centers of mass for an initial stage as well as for a mature stage of their westward migration. In particular, near-inertial oscillations are evident in the initial migration stage, whose amplitude linearly increases with time, as a result of the unbalanced vortex initial state on a beta plane. Such an initial amplification of the vortex oscillations is actually found in the first stage of the evolution of warm-core frontal eddies simulated numerically by means of a frontal numerical model initialized using the shape and velocity fields of circular pulsons of the first order. In the numerical simulations, this stage is followed by an adjusted, complex nonstationary state characterized by a noticeable asymmetry in the meridional component of the vortex's horizontal pressure gradient, which develops to compensate for the variations of the Coriolis parameter with latitude. Accordingly, the location of the simulated vortex's maximum depth is always found poleward of the location of the simulated vortex's center of mass. Moreover, during the adjusted stage, near-inertial oscillations emerge that largely deviate from the exactly inertial ones characterizing analytical circular pulsons: a superinertial and a subinertial oscillation in fact appear, and their frequency difference is found to be an increasing function of latitude. A comparison between vortex westward drifts simulated numerically at different latitudes for different vortex radii and pulsation strengths and the corresponding drifts obtained using existing formulas shows that, initially, the simulated vortex drifts correspond to the fastest predicted ones in many realistic cases. As time elapses, however, the development of a beta-adjusted vortex structure, together with the effects of numerical dissipation, tend to slow down the simulated vortex drift.

Description: Wind-induced near-inertial energy has been believed to be an important source for generating the ocean mixing required to maintain the global meridional overturning circulation. In the present study, the near-inertial energy budget in a realistic (1)/(12)degrees model of the North Atlantic Ocean driven by synoptically varying wind forcing is examined. The authors find that nearly 70% of the wind-induced near-inertial energy at the sea surface is lost to turbulent mixing within the top 200 m and, hence, is not available to generate diapycnal mixing at greater depth. Assuming this result can be extended to the global ocean, it is estimated that the wind-induced near-inertial energy available for ocean mixing at depth is, at most, 0.1 TW. This confirms a recent suggestion that the role of wind-induced near-inertial energy in sustaining the global overturning circulation might have been overemphasized.

Description: A new field of study, "decadal prediction," is emerging in climate science. Decadal prediction lies between seasonal/interannual forecasting and longer-term climate change projections, and focuses on time-evolving regional climate conditions over the next 10-30 yr. Numerous assessments of climate information user needs have identified this time scale as being important to infrastructure planners, water resource managers, and many others. It is central to the information portfolio required to adapt effectively to and through climatic changes. At least three factors influence time-evolving regional climate at the decadal time scale: 1) climate change commitment (further warming as the coupled climate system comes into adjustment with increases of greenhouse gases that have already occurred), 2) external forcing, particularly from future increases of greenhouse gases and recovery of the ozone hole, and 3) internally generated variability. Some decadal prediction skill has been demonstrated to arise from the first two of these factors, and there is evidence that initialized coupled climate models can capture mechanisms of internally generated decadal climate variations, thus increasing predictive skill globally and particularly regionally. Several methods have been proposed for initializing global coupled climate models for decadal predictions, all of which involve global time-evolving threedimensional ocean data, including temperature and salinity. An experimental framework to address decadal predictability/prediction is described in this paper and has been incorporated into the coordinated Coupled Model Intercomparison Model, phase 5 (CMIP5) experiments, some of which will be assessed for the IPCC Fifth Assessment Report (AR5). These experiments will likely guide work in this emerging field over the next 5 yr.

Description: Recent observations show dramatic changes of the Arctic atmosphere–ice–ocean system. Here the authors demonstrate, through the analysis of a vast collection of previously unsynthesized observational data, that over the twentieth century the central Arctic Ocean became increasingly saltier with a rate of freshwater loss of 239 ± 270 km3 decade−1. In contrast, long-term (1920–2003) freshwater content (FWC) trends over the Siberian shelf show a general freshening tendency with a rate of 29 ± 50 km3 decade−1. These FWC trends are modulated by strong multidecadal variability with sustained and widespread patterns. Associated with this variability, the FWC record shows two periods in the 1920s–30s and in recent decades when the central Arctic Ocean was saltier, and two periods in the earlier century and in the 1940s–70s when it was fresher. The current analysis of potential causes for the recent central Arctic Ocean salinification suggests that the FWC anomalies generated on Arctic shelves (including anomalies resulting from river discharge inputs) and those caused by net atmospheric precipitation were too small to trigger long-term FWC variations in the central Arctic Ocean; to the contrary, they tend to moderate the observed long-term central-basin FWC changes. Variability of the intermediate Atlantic Water did not have apparent impact on changes of the upper–Arctic Ocean water masses. The authors’ estimates suggest that ice production and sustained draining of freshwater from the Arctic Ocean in response to winds are the key contributors to the salinification of the upper Arctic Ocean over recent decades. Strength of the export of Arctic ice and water controls the supply of Arctic freshwater to subpolar basins while the intensity of the Arctic Ocean FWC anomalies is of less importance. Observational data demonstrate striking coherent long-term variations of the key Arctic climate parameters and strong coupling of long-term changes in the Arctic–North Atlantic climate system. Finally, since the high-latitude freshwater plays a crucial role in establishing and regulating global thermohaline circulation, the long-term variations of the freshwater content discussed here should be considered when assessing climate change and variability.

Description: Data from an array of six moorings deployed east of Abaco, Bahamas, along 26.5°N during March 2004–May 2005 are analyzed. These moorings formed the western boundary array of a transbasin observing system designed to continuously monitor the meridional overturning circulation and meridional heat flux in the subtropical North Atlantic, under the framework of the joint U.K.–U.S. Rapid Climate Change (RAPID)–Meridional Overturning Circulation (MOC) Program. Important features of the western boundary circulation include the southward-flowing deep western boundary current (DWBC) below 1000 m and the northward-flowing “Antilles” Current in the upper 1000 m. Transports in the western boundary layer are estimated from direct current meter observations and from dynamic height moorings that measure the spatially integrated geostrophic flow between moorings. The results of these methods are combined to estimate the time-varying transports in the upper and deep ocean over the width of the western boundary layer to a distance of 500 km offshore of the Bahamas escarpment. The net southward transport of the DWBC across this region, inclusive of northward deep recirculation, is −26.5 Sv (Sv ≡ 106 m3 s−1), which is divided nearly equally between upper (−13.9 Sv) and lower (−12.6 Sv) North Atlantic Deep Water (NADW). In the top 1000 m, 6.0 Sv flows northward in a thermocline-intensified jet near the western boundary. These transports are found to agree well with historical current meter data in the region collected between 1986 and 1997. Variability in both shallow and deep components of the circulation is large, with transports above 1000 m varying between −15 and +25 Sv and deep transports varying between −60 and +3 Sv. Much of this transport variability, associated with barotropic fluctuations, occurs on relatively short time scales of several days to a few weeks. Upon removal of the barotropic fluctuations, slower baroclinic transport variations are revealed, including a temporary stoppage of the lower NADW transport in the DWBC during November 2004.

Description: A tangent linear adjoint for a low-resolution dynamical model of the atmosphere is used to derive the optimal forcing perturbations for all state variables such that after a specified lead time the model response has a given projection, in terms of an energy norm, on the pattern associated with the 51-yr trend in the Northern Hemisphere winter tropospheric circulation, 1948/49–1998/99. A feature of the derived forcing sensitivity is a Rossby wave–like feature that emanates from the western tropical Pacific and is associated with the deepening of the Aleutian low, whereas an annular pattern in the forcing sensitivity in the uppermost model level is shown to be associated with the pattern of the trend over the Euro-Atlantic/Asian sectors, including the upward trend in the North Atlantic Oscillation index. The authors argue that the Rossby wave–type feature is consistent with studies that have argued a role for the upward trend in tropical sea surface temperature during the 51-yr period. On the other hand, the authors interpret the annular pattern in the forcing sensitivity as being consistent with studies that have argued that the trend over the Euro-Atlantic sector was associated with influences from the stratosphere. In particular, a nonlinear model driven by the optimal forcing perturbation applied only to the top model level is successful at reproducing the trend pattern with the correct amplitude in the Euro-Atlantic sector, but implies a trend over the North Pacific toward a weaker Aleutian low, contrary to what was observed but similar to the spatial pattern associated with the northern annular mode. These results show that the adjoint approach can shed light on previous apparently different interpretations of the trend. The study also presents a successful application of a tangent linear adjoint model to a climate problem.

Description: Recent observations show dramatic changes of the Arctic atmosphere–ice–ocean system. Here the authors demonstrate, through the analysis of a vast collection of previously unsynthesized observational data, that over the twentieth century the central Arctic Ocean became increasingly saltier with a rate of freshwater loss of 239 ± 270 km3 decade−1. In contrast, long-term (1920–2003) freshwater content (FWC) trends over the Siberian shelf show a general freshening tendency with a rate of 29 ± 50 km3 decade−1. These FWC trends are modulated by strong multidecadal variability with sustained and widespread patterns. Associated with this variability, the FWC record shows two periods in the 1920s–30s and in recent decades when the central Arctic Ocean was saltier, and two periods in the earlier century and in the 1940s–70s when it was fresher. The current analysis of potential causes for the recent central Arctic Ocean salinification suggests that the FWC anomalies generated on Arctic shelves (including anomalies resulting from river discharge inputs) and those caused by net atmospheric precipitation were too small to trigger long-term FWC variations in the central Arctic Ocean; to the contrary, they tend to moderate the observed long-term central-basin FWC changes. Variability of the intermediate Atlantic Water did not have apparent impact on changes of the upper–Arctic Ocean water masses. The authors’ estimates suggest that ice production and sustained draining of freshwater from the Arctic Ocean in response to winds are the key contributors to the salinification of the upper Arctic Ocean over recent decades. Strength of the export of Arctic ice and water controls the supply of Arctic freshwater to subpolar basins while the intensity of the Arctic Ocean FWC anomalies is of less importance. Observational data demonstrate striking coherent long-term variations of the key Arctic climate parameters and strong coupling of long-term changes in the Arctic–North Atlantic climate system. Finally, since the high-latitude freshwater plays a crucial role in establishing and regulating global thermohaline circulation, the long-term variations of the freshwater content discussed here should be considered when assessing climate change and variability.

Description: The causes and characteristics of interannual–decadal variability of the meridional overturning circulation (MOC) in the North Atlantic are investigated with a suite of basin-scale ocean models [the Family of Linked Atlantic Model Experiments (FLAME)] and global ocean–ice models (ORCA), varying in resolution from medium to eddy resolving (½°–1/12°), using various forcing configurations built on bulk formulations invoking atmospheric reanalysis products. Comparison of the model hindcasts indicates similar MOC variability characteristics on time scales up to a decade; both model architectures also simulate an upward trend in MOC strength between the early 1970s and mid-1990s. The causes of the MOC changes are examined by perturbation experiments aimed selectively at the response to individual forcing components. The solutions emphasize an inherently linear character of the midlatitude MOC variability by demonstrating that the anomalies of a (non–eddy resolving) hindcast simulation can be understood as a superposition of decadal and longer-term signals originating from thermohaline forcing variability, and a higher-frequency wind-driven variability. The thermohaline MOC signal is linked to the variability in subarctic deep-water formation, and rapidly progressing to the tropical Atlantic. However, throughout the subtropical and midlatitude North Atlantic, this signal is effectively masked by stronger MOC variability related to wind forcing and, especially north of 30°–35°N, by internally induced (eddy) fluctuations.

Description: The relative impact of the subtropical North and South Pacific Oceans on the tropical Pacific climate mean state and variability is estimated using an ocean–atmosphere–sea ice coupled general circulation model. Tailored experiments are performed in which the model is forced by idealized sea surface temperature anomalies (SSTAs) in the subtropics of both hemispheres. The main results of this study suggest that subtropical South Pacific climate variations play a dominant role in tropical Pacific decadal variability and in the decadal modulation of El Niño–Southern Oscillation (ENSO). In response to a 2°C warming in the subtropical South Pacific, the equatorial Pacific SST increases by about 0.6°C, approximately 65% larger than the change in the North Pacific experiment. The subtropics affect equatorial SST mainly through atmosphere–mixed layer interactions in the South Pacific experiments; the response is mostly accomplished within a decade. The “oceanic tunnel” dominates in the North Pacific experiments; the response takes at least 100 yr to be accomplished. Similar sensitivity experiments conducted with the stand-alone atmosphere model showed that both air–sea interactions and ocean dynamics are crucial in shaping the tropical climate response. The statistics of ENSO exhibit significant changes in amplitude and frequency in response to a warming/cooling of the subtropical South Pacific: a 2°C warming (cooling) of subtropical South Pacific SST reduces (increases) the interannual standard deviation by about 30% (20%) and shortens (lengthens) the ENSO period. The simulated changes in the equatorial zonal SST gradient are the main contributor to the modulation of ENSO variability. The simulated intensification (weakening) of the annual cycle in response to an enhanced warming (cooling) in subtropical South Pacific partly explains the shifts in frequency, but may also lead to a weaker (stronger) ENSO. The subtropical North Pacific thermal forcing did not change the statistical properties of ENSO as strongly.

Description: The manner in which monthly mean sea surface temperature anomalies (SSTAs) show enhanced variance at the annual period in the extratropics (an annual peak in the variance spectrum) is illustrated by observations and model simulations. A mechanism, related to the reemergence of winter SST anomalies, is proposed to explain the annual peak in SST spectrum. The idea is supported by the analysis of a hierarchy of models, including Intergovernmental Panel on Climate Change model simulations. The results of the model experiments further suggest that the annual peak is either weak or absent if decadal SST variability is forced by local air–sea interaction. However, if ocean subsurface temperature variability forces decadal SST variability, the annual peak is much stronger. Strong annual peaks may therefore be seen as an indication of ocean-forced decadal SST variability in the extratropics.

Description: The shallow subtropical–tropical cells (STC) of the Atlantic Ocean have been studied from the output fields of a 50-yr run of the German partner of the Estimating the Circulation and Climate of the Ocean (GECCO) consortium assimilation model. Comparison of GECCO with time-mean observational estimates of density and meridional currents at 10°S and 10°N, which represent the boundaries between the tropics and subtropics in GECCO, shows good agreement in transports of major currents. The variability of the GECCO wind stress in the interior at 10°S and 10°N remains consistent with the NCEP forcing, although temporary changes can be large. On pentadal and longer time scales, an STC loop response is found between the poleward Ekman divergence and STC-layer convergence at 10°S and 10°N via the Equatorial Undercurrent (EUC) at 23°W, where the divergence leads the EUC and the convergence, suggesting a “pulling” mechanism via equatorial upwelling. The divergence is also associated with changes in the eastern equatorial upper-ocean heat content. Within the STC layer, partial compensation of the western boundary current (WBC) and the interior occurs at 10°S and 10°N. For the meridional overturning circulation (MOC) at 10°S it is found that more than one-half of the variability in the upper limb can be explained by the WBC. The explained MOC variance can be increased to 85% by including the geostrophic (Sverdrup) part of the wind-driven transports.

Description: Observed sea surface temperatures (SSTs) in the North Atlantic from 1958 through 2000, as well as data from an ocean model simulation driven with the atmospheric variability observed during the same period, are examined using multichannel singular spectrum analysis. The two leading oscillatory modes are associated with a multidecadal and a quasi-decadal period. The former is connected to a basinwide uniform SST pattern and changes in the deep North Atlantic meridional overturning circulation. The quasi-decadal mode involves a tripolar SST anomaly pattern forced by atmospheric variability with a spatial structure resembling that of the North Atlantic Oscillation (NAO). The upper ocean’s dynamical response to this NAO variability provides an instantaneous positive feedback to the SST pattern, while a delayed negative feedback is due to shallow overturning circulation anomalies.

Description: Shortly after the advent of the first imaging passive microwave sensor on board a research satellite an anomalous climate feature was observed within the Weddell Sea. During the years 1974–1976, a 250 × 103 km2 area within the seasonal sea ice cover was virtually free of winter sea ice. This feature, the Weddell Polynya, was created as sea ice formation was inhibited by ocean convection that injected relatively warm deep water into the surface layer. Though smaller, less persistent polynyas associated with topographically induced upwelling at Maud Rise frequently form in the area, there has not been a reoccurrence of the Weddell Polynya since 1976. Archived observations of the surface layer salinity within the Weddell gyre suggest that the Weddell Polynya may have been induced by a prolonged period of negative Southern Annular Mode (SAM). During negative SAM the Weddell Sea experiences colder and drier atmospheric conditions, making for a saltier surface layer with reduced pycnocline stability. This condition enables Maud Rise upwelling to trigger sustained deep-reaching convection associated with the polynya. Since the late 1970s SAM has been close to neutral or in a positive state, resulting in warmer, wetter conditions over the Weddell Sea, forestalling repeat of the Weddell Polynya. A contributing factor to the Weddell Polynya initiation may have been a La Niña condition, which is associated with increased winter sea ice formation in the polynya area. If the surface layer is made sufficiently salty due to a prolonged negative SAM period, perhaps aided by La Niña, then Maud Rise upwelling meets with positive feedback, triggering convection, and a winter persistent Weddell Polynya.

Description: A generalization of the transformed Eulerian and temporal residual means is presented. The new formulation uses rotational fluxes of buoyancy, and the full hierarchy of statistical density moments, to reduce the cross-isopycnal eddy flux to the physically relevant component associated with the averaged water mass properties. The resulting eddy-induced diapycnal diffusivity vanishes for adiabatic, statistically steady flow, and is related to either the growth or decay of mesoscale density variance and/or the covariance between small-scale forcing (mixing) and density fluctuations, such as that associated with the irreversible removal of density variance by dissipation. The relationship between the new formulation and previous approaches is described and is illustrated using results from an eddying channel model. The formalism is quite general and applies to all kinds of averaging and to any tracer (not just density).

Description: Using the same approach as in Part I, here it is shown how sampling problems in voluntary observing ship (VOS) data affect conclusions about interannual variations and secular changes of surface heat fluxes. The largest uncertainties in linear trend estimates are found in relatively poorly sampled regions like the high-latitude North Atlantic and North Pacific as well as the Southern Ocean, where trends can locally show opposite signs when computed from the regularly sampled and undersampled data. Spatial patterns of shorter-period interannual variability, quantified through the EOF analysis, also show remarkable differences between the regularly sampled and undersampled flux datasets in the Labrador Sea and northwest Pacific. In particular, it is shown that in the Labrador Sea region, in contrast to regularly sampled NCEP–NCAR reanalysis fluxes, VOS-like sampled NCEP–NCAR reanalysis fluxes neither show significant interannual variability nor significant trends. These regions, although quite localized covering small parts of the globe, play a crucial role for the coupled atmosphere–ocean system. In the Labrador Sea, for instance, interannual and decadal-scale changes of the surface net heat fluxes are known to affect oceanic convection and, thus, the meridional overturning circulation of the Atlantic Ocean. From a discussion of current atmospheric data assimilation systems it is argued that in poorly sampled regions reanalysis products are superior to VOS-based products for studying interannual and interdecadal variations of atmosphere–ocean interaction. In well-sampled regions, on the other hand, conclusions about surface heat flux variations are relatively insensitive to the choice of the flux products used (VOS versus reanalysis data). The results are confirmed for two different datasets, that is, ECMWF 40-yr Re-Analysis (ERA-40) data and seasonal integrations with a recent version of the ECMWF model in which no actual data were assimilated.

Description: On 19 October 2000, Hurricane Michael merged with an approaching baroclinic trough over the western North Atlantic Ocean south of Nova Scotia. As the hurricane moved over cooler sea surface temperatures (SSTs; less than 25°C), it intensified to category-2 intensity on the Saffir–Simpson hurricane scale [maximum sustained wind speeds of 44 m s−1 (85 kt)] while tapping energy from the baroclinic environment. The large “hybrid” storm made landfall on the south coast of Newfoundland with maximum sustained winds of 39 m s−1 (75 kt) causing moderate damage to coastal communities east of landfall. Hurricane Michael presented significant challenges to weather forecasters. The fundamental issue was determining which of two cyclones (a newly formed baroclinic low south of Nova Scotia or the hurricane) would become the dominant circulation center during the early stages of the extratropical transition (ET) process. Second, it was difficult to predict the intensity of the storm at landfall owing to competing factors: 1) decreasing SSTs conducive to weakening and 2) the approaching negatively tilted upper-level trough, favoring intensification. Numerical hindcast simulations using the limited-area Mesoscale Compressible Community model with synthetic vortex insertion (cyclone bogus) prior to the ET of Hurricane Michael led to a more realistic evolution of wind and pressure compared to running the model without vortex insertion. Specifically, the mesoscale model correctly simulates the hurricane as the dominant circulation center early in the transition process, versus the baroclinic low to its north, which was the favored development in the runs not employing vortex insertion. A suite of experiments is conducted to establish the sensitivity of the ET to various initial conditions, lateral driving fields, domain sizes, and model parameters. The resulting storm tracks and intensities fall within the range of the operational guidance, lending support to the possibility of improving numerical forecasts using synthetic vortex insertion prior to ET in such a model.

Description: In this paper, a version of the European Centre for Medium-Range Weather Forecasts (ECMWF) operational model is used to (i) diagnose the diabatic heating associated with the winter North Atlantic Oscillation (NAO) and (ii) assess the role of this heating in the dynamics of the NAO in the model. Over the North Atlantic sector, the NAO-related diabatic heating is dominated above the planetary boundary layer by the latent heat release associated with precipitation, and within the boundary layer by vertical diffusion associated with sensible heat flux from the ocean. An association between La Niña–El Niño–type conditions in the tropical Pacific and the positive/negative NAO is found in model runs using initial conditions and sea surface temperature (SST) lower boundary conditions from the period 1982–2001, but not in a companion set of model runs for the period 1962–81. Model experiments are then described in which the NAO-related diabatic heating diagnosed from the 1982–2001 control run is applied as a constant forcing in the model temperature equation using both 1982–2001 and 1962–81 model setups. To assess the local feedback from the diabatic heating, the specified forcing is first restricted to the North Atlantic sector alone. In this case, the model response (in an ensemble mean sense) is suggestive of a weak negative feedback, but exhibits more baroclinic structure and has its centers of action shifted compared to those of the NAO. On the other hand, forcing with only the tropical Pacific part of the diabatic heating leads to a robust model response in both the 1982–2001 and 1962–81 model setups. The model response projects on to the NAO with the same sign as that used to diagnose the forcing, arguing that the link between the tropical Pacific and the NAO is real in the 1982–2001 control run. The missing link in the corresponding run for 1962–81 is a result of a change in the tropical forcing between the two periods, and not the extratropical flow regime.

Description: The depth of winter convection in the central Labrador Sea is strongly influenced by the prevailing stratification in late summer. For this late summer stratification salinity is as important as temperature, and in the upper water layers salinity even dominates. To analyze the source of the spring and summer freshening in the central region, seasonal freshwater cycles have been constructed for the interior Labrador Sea, the West Greenland Current, and the Labrador Current. It is shown that none of the local freshwater sources is responsible for the spring–summer freshening in the interior, which appears to occur in two separate events in April to May and July to September. Comparing the timing and volume estimates of the seasonal freshwater cycles of the boundary currents with the central Labrador Sea helps in understanding the origin of the interior freshwater signals. The first smaller pulse cannot be attributed clearly to either of the boundary currents. The second one is about three times stronger and supplies 60% of the seasonal summer freshwater. Transport estimates and calculated mixing properties provide evidence that its source is the West Greenland Current. The finding implies a connection also on interannual time scales between Labrador Sea surface salinity and freshwater sources in the West Greenland Current and farther upstream in the East Greenland Current. The freshwater input from the West Greenland Current thus also is the likely pathway for the known modulation of Labrador Sea Water mass formation by freshwater export from the Arctic (via the East Greenland Current), which implies some predictability on longer time scales.

Description: Multichannel singular spectrum analysis (MSSA) of surface zonal wind, sea surface temperature (SST), 20° isotherm depth, and surface zonal current observations (between 1990 and 2004) identifies three coupled ocean–atmosphere modes of variability in the tropical Pacific: the El Niño–Southern Oscillation (ENSO), the annual cycle, and a mode with a 14–18-month period, which is referred to as sub-ENSO in this study. The sub-ENSO mode accounts for the near 18-month (near annual) variability prior to (following) the 1997/98 El Niño event. It was strongest during this El Niño event, with SST anomalies exceeding 1°C. Sub-ENSO peak SST anomalies are ENSO-like in structure and are associated with eastward propagating heat content variations. However, the SST anomalies are preceded by and in near quadrature with relatively strong remotely forced westward propagating zonal current variations, suggesting the sub-ENSO mode arises from the zonal-advective feedback. The sub-ENSO mode is found to exist also in an intermediate complexity model (ICM) of the tropical Pacific. A heat budget analysis of the model’s sub-ENSO mode shows it indeed arises from the zonal-advective feedback. In the model, both ENSO and sub-ENSO modes coexist, but there is a weak nonlinear interaction between them. Experiments also show that the observed changes in sub-ENSO’s characteristics may be explained by changes in the relative importance of zonal and vertical advection SST tendencies.

Description: A coordinated set of global coupled climate model [atmosphere–ocean general circulation model (AOGCM)] experiments for twentieth- and twenty-first-century climate, as well as several climate change commitment and other experiments, was run by 16 modeling groups from 11 countries with 23 models for assessment in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). Since the assessment was completed, output from another model has been added to the dataset, so the participation is now 17 groups from 12 countries with 24 models. This effort, as well as the subsequent analysis phase, was organized by the World Climate Research Programme (WCRP) Climate Variability and Predictability (CLIVAR) Working Group on Coupled Models (WGCM) Climate Simulation Panel, and constitutes the third phase of the Coupled Model Intercomparison Project (CMIP3). The dataset is called the WCRP CMIP3 multimodel dataset, and represents the largest and most comprehensive international global coupled climate model experiment and multimodel analysis effort ever attempted. As of March 2007, the Program for Climate Model Diagnostics and Intercomparison (PCMDI) has collected, archived, and served roughly 32 TB of model data. With oversight from the panel, the multimodel data were made openly available from PCMDI for analysis and academic applications. Over 171 TB of data had been downloaded among the more than 1000 registered users to date. Over 200 journal articles, based in part on the dataset, have been published so far. Though initially aimed at the IPCC AR4, this unique and valuable resource will continue to be maintained for at least the next several years. Never before has such an extensive set of climate model simulations been made available to the international climate science community for study. The ready access to the multimodel dataset opens up these types of model analyses to researchers, including students, who previously could not obtain state-of-the-art climate model output, and thus represents a new era in climate change research. As a direct consequence, these ongoing studies are increasing the body of knowledge regarding our understanding of how the climate system currently works, and how it may change in the future.

Description: Sampling uncertainties in the voluntary observing ship (VOS)-based global ocean–atmosphere flux fields were estimated using the NCEP–NCAR reanalysis and ECMWF 40-yr Re-Analysis (ERA-40) as well as seasonal forecasts without data assimilation. Air–sea fluxes were computed from 6-hourly reanalyzed individual variables using state-of-the-art bulk formulas. Individual variables and computed fluxes were subsampled to simulate VOS-like sampling density. Random simulation of the number of VOS observations and simulation of the number of observations with contemporaneous sampling allowed for estimation of random and total sampling uncertainties respectively. Although reanalyses are dependent on VOS, constituting an important part of data assimilation input, it is assumed that the reanalysis fields adequately reproduce synoptic variability at the sea surface. Sampling errors were quantified by comparison of the regularly sampled (i.e., 6 hourly) and subsampled monthly fields of surface variables and fluxes. In poorly sampled regions random sampling errors amount to 2.5°–3°C for air temperature, 3 m s−1 for the wind speed, 2–2.5 g kg−1 for specific humidity, and 15%–20% of the total cloud cover. The highest random sampling errors in surface fluxes were found for the sensible and latent heat flux and range from 30 to 80 W m−2. Total sampling errors in poorly sampled areas may be higher than random ones by 60%. In poorly sampled subpolar latitudes of the Northern Hemisphere and throughout much of the Southern Ocean the total sampling uncertainty in the net heat flux can amount to 80–100 W m−2. The highest values of the uncertainties associated with the interpolation/extrapolation into unsampled grid boxes are found in subpolar latitudes of both hemispheres for the turbulent fluxes, where they can be comparable with the sampling errors. Simple dependencies of the sampling errors on the number of samples and the magnitude of synoptic variability were derived. Sampling errors estimated from different reanalyses and from seasonal forecasts yield qualitatively comparable spatial patterns, in which the actual values of uncertainties are controlled by the magnitudes of synoptic variability. Finally, estimates of sampling uncertainties are compared with the other errors in air–sea fluxes and the reliability of the estimates obtained is discussed.

Description: Output from an eddy-resolving model of the North Atlantic Ocean is used to estimate values for the thickness diffusivity κ appropriate to the Gent and McWilliams parameterization. The effect of different choices of rotational eddy fluxes on the estimated κ is discussed. Using the raw fluxes (no rotational flux removed), large negative values (exceeding −5000 m2 s−1) of κ are diagnosed locally, particularly in the Gulf Stream region and in the equatorial Atlantic. Removing a rotational flux based either on the suggestion of Marshall and Shutts or the more general theory of Medvedev and Greatbatch leads to a reduction of the negative values, but they are still present. The regions where κ 〈 0 correspond to regions where eddies are acting to increase, rather than decrease (as in baroclinic instability) the mean available potential energy. In the subtropical gyre, κ ranges between 500 and 2000 m2 s−1, rapidly decreasing to zero below the thermocline in all cases. Rotational fluxes and κ are also estimated using an optimization technique. In this case, |κ| can be reduced or increased by construction, but the regions where κ 〈 0 are still present and the optimized rotational fluxes also remain similar to a priori values given by the theoretical considerations. A previously neglected component (ν) of the bolus velocity is associated with the horizontal flux of buoyancy along, rather than across, the mean buoyancy contours. The ν component of the bolus velocity is interpreted as a streamfunction for eddy-induced advection, rather than diffusion, of mean isopycnal layer thickness, showing up when the lateral eddy fluxes cannot be described by isotropic diffusion only. All estimates show a similar large-scale pattern for ν, implying westward advection of isopycnal thickness over much of the subtropical gyre. Comparing ν with a mean streamfunction shows that it is about 10% of the mean in midlatitudes and even larger than the mean in the Tropics.

Description: Combined simultaneous satellite observations are used to evaluate the performance of parameterizations of the microphysical and optical properties of cirrus clouds used for radiative flux computations in climate models. Atmospheric and cirrus properties retrieved from Television and Infrared Observation Satellite (TIROS-N) Operational Vertical Sounder (TOVS) observations are given as input to the radiative transfer model developed for the Met Office climate model to simulate radiative fluxes at the top of the atmosphere (TOA). Simulated cirrus shortwave (SW) albedos are then compared to those retrieved from collocated Scanner for Radiation Budget (ScaRaB) observations. For the retrieval, special care has been given to angular direction models. Three parameterizations of cirrus ice crystal optical properties are represented in the Met Office radiative transfer model. These parameterizations are based on different physical approximations and different hypotheses on crystal habit. One parameterization assumes pristine ice crystals and two ice crystal aggregates. By relating the cirrus ice water path (IWP) retrieved from the effective infrared emissivity to the cirrus SW albedo, differences between the parameterizations are amplified. This study shows that pristine crystals seem to be plausible only for cirrus with IWP less than 30 g m−2. For larger IWP, ice crystal aggregates lead to cirrus SW albedos in better agreement with the observations. The data also indicate that climate models should allow the cirrus effective ice crystal diameter (De) to increase with IWP, especially in the range up to 30 g m−2. For cirrus with IWP less than 20 g m−2, this would lead to SW albedos that are about 0.02 higher than the ones of a constant De of 55 μm.

Description: A model of the subpolar North Atlantic Ocean is used to study different aspects of ventilation and water mass transformation during a year with moderate convection intensity in the Labrador Sea. The model realistically describes the salient features of the observed hydrographic structure and current system, including boundary currents and recirculations. Ventilation and transformation rates are defined and compared. The transformation rate of Labrador Sea Water (LSW), defined in analogy to several observational studies, is 6.3 Sv (Sv ≡ 106 m3 s−1) in the model. Using an idealized ventilation tracer, mimicking analyses based on chlorofluorocarbon inventories, an LSW ventilation rate of 10 Sv is found. Differences between both rates are particularly significant for those water masses that are partially transformed into denser water masses during winter. The main export route of the ventilated LSW is the deep Labrador Current (LC). Backward calculation of particle trajectories demonstrates that about one-half of the LSW leaving the Labrador Sea within the deep LC originates in the mixed layer during that same year. Near the offshore flank of the deep LC at about 55°W, the transformation of LSW begins in January and is at a maximum in February/March. While the export of transformed LSW out of the central Labrador Sea continues for several months, LSW generated near the boundary current is exported more rapidly, with maximum transport rates during March/April within the deep LC.