ALBERT

Description: This review article aims to provide an overview and insight into the most relevant aspects of wind energy development and current state-of-the-art. The industry is in a very mature stage, so it seems to be the right time to take stock of the relevant areas of wind energy use for power generation. For this review, the authors considered the essential aspects of the development of wind energy technology: research, modeling, and prediction of wind speed as an energy source, the technology development of the plants divided into the mechanical and electrical systems and the plant control, and finally the optimal plant operation including the maintenance strategies. The focus is on the development in Europe, with a partial focus on Germany. The authors are employees of the Fraunhofer Institutes, Institute for Energy Economics and Energy Systems Technology and Institute for Wind Energy Systems, who have contributed to the development of this technology for decades.

Description: A simple analytical model of the atmospheric boundary layer (ABL) coupled to sea ice is presented. It describes clear-sky cooling over sea ice during polar night in the presence of leads. The model solutions show that the sea ice concentration and wind speed have a strong impact on the thermal regime over sea ice. Leads cause both a warming of the ABL and an increase of stability over sea ice. The model describes a sharp ABL transition from a weakly stable coupled state to a strongly stable decoupled state when wind speed is decreasing. The threshold value of the transition wind speed is a function of sea ice concentration. The decoupled state is characterized by a large air–surface temperature difference over sea ice, which is further increased by leads. In the coupled regime, air and surface temperatures increase almost linearly with wind speed due to warming by leads and also slower cooling of the ABL. The cooling time scale shows a nonmonotonic dependency on wind speed, being lowest for the threshold value of wind speed and increasing for weak and strong winds. Theoretical solutions agree well with results of a more realistic single-column model and with observations performed at the three Russian “North Pole” drifting stations (NP-35, -37, and -39) and at the Surface Heat Budget of the Arctic Ocean ice camp. Both modeling results and observations show a strong implicit dependency of the net longwave radiative flux at the surface on wind speed.

Description: Heinrich events are among the dominant modes of glacial climate variability. During these events, massive iceberg armadas were released by the Laurentide Ice Sheet, sailed across the Atlantic, and caused large-scale climate changes. We study these events in a fully coupled complex ice sheet–climate model with synchronous coupling between ice sheets and oceans. The ice discharges occur as internal variability of the model with a recurrence period of 5kyr, an event duration of 1–1.5kyr, and a peak discharge rate of about 50mSv, roughly consistent with reconstructions. The climate response shows a two-stage behavior, with freshwater release effects dominating the surge phase and ice-sheet elevation effects dominating in the post-surge phase. As a direct response to the freshwater discharge during the surge phase, the deepwater formation in the North Atlantic decreases and the North Atlantic deepwater cell weakens by 3.5Sv. With the reduced oceanic heat transport, the surface temperatures across the North Atlantic decrease, and the associated reduction in evaporation causes a drying in Europe. The ice discharge lowers the surface elevation in the Hudson Bay area and thus leads to increased precipitation and accelerated ice sheet regrowth in the post-surge phase. Furthermore, the jet stream widens to the north and becomes more zonal. This contributes to a weakening of the subpolar gyre, and a continued cooling over Europe even after the ice discharge. This two-stage behavior can explain previously contradicting model results and understandings of Heinrich Events.

Description: Proxy data and observations suggest that large tropical volcanic eruptions induce a poleward shift of the North Atlantic jet stream in boreal winter. However, there is far from universal agreement in models on this effect and its mechanism, and the possibilities of a corresponding jet shift in the Southern Hemisphere or the summer season have received little attention. Using a hierarchy of simplified atmospheric models, this study examines the impact of stratospheric aerosol on the extratropical circulation over the annual cycle. In particular, the models allow the separation of the dominant shortwave (surface cooling) and longwave (stratospheric warming) impacts of volcanic aerosol. It is found that stratospheric warming shifts the jet poleward in both summer and winter hemispheres. The experiments cannot definitively rule out the role of surface cooling, but provide no evidence that it shifts the jet poleward. Further study with simplified models demonstrates that the response to stratospheric warming is remarkably generic and does not depend critically on the boundary conditions (e.g., the planetary wave forcing) or the atmospheric physics (e.g., the treatment of radiative transfer and moist processes). It does, however, fundamentally involve both zonal-mean and eddy circulation feedbacks. The timescales, seasonality, and structure of the response provide further insight into the mechanism, as well as its connection to modes of intrinsic natural variability. These findings have implications for the interpretation of comprehensive model studies and for post-volcanic prediction

Description: To model tracer spreading in the ocean, Lagrangian simulations in an offline framework are a practical and efficient alternative to solving the advective–diffusive tracer equations online. Differences in both approaches raise the question of whether both methods are comparable. Lagrangian simulations usually use model output averaged in time, and trajectories are not subject to parameterized subgrid diffusion, which is included in the advection–diffusion equations of ocean models. Previous studies focused on diffusivity estimates in idealized models but could show that both methods yield similar results as long as the deformations-scale dynamics are resolved and a sufficient amount of Lagrangian particles is used. This study compares spreading of an Eulerian tracer simulated online and a cloud of Lagrangian particles simulated offline with velocities from the same ocean model. We use a global, eddy-resolving ocean model featuring 1/20° horizontal resolution in the Agulhas region around South Africa. Tracer and particles were released at one time step in the Cape Basin and below the mixed layer and integrated for 3 years. Large-scale diagnostics, like mean pathways of floats and tracer, are almost identical and 1D horizontal distributions show no significant differences. Differences in vertical distributions, seen in a reduced vertical spreading and downward displacement of particles, are due to the combined effect of unresolved subdaily variability of the vertical velocities and the spatial variation of vertical diffusivity. This, in turn, has a small impact on the horizontal spreading behavior. The estimates of eddy diffusivity from particles and tracer yield comparable results of about 4000 m2 s−1 in the Cape Basin.

Description: Satellite observations and output from a high-resolution ocean model are used to investigate how the Loop Current in the Gulf of Mexico affects the Gulf Stream transport through the Florida Straits. We find that the expansion (contraction) of the Loop Current leads to lower (higher) transports through the Straits of Florida. The associated surface velocity anomalies are coherent from the southwestern tip of Florida to Cape Hatteras. A simple continuity-based argument can be used to explain the link between the Loop Current and the downstream Gulf Stream transport: as the Loop Current lengthens (shortens) its path in the Gulf of Mexico, the flow out of the Gulf decreases (increases). Anomalies in the surface velocity field are first seen to the southwest of Florida and within 4 weeks propagate through the Florida Straits up to Cape Hatteras and into the Gulf Stream Extension. In both the observations and the model this propagation can be seen as pulses in the surface velocities. We estimate that the Loop Current variability can be linked to a variability of several Sverdrups (1Sv = 10(6) m(3) s(-1)) through the Florida Straits. The exact timing of the Loop Current variability is largely unpredictable beyond a few weeks and its variability is therefore likely a major contributor to the chaotic/intrinsic variability of the Gulf Stream. However, the time lag between the Loop Current and the flow downstream of the Gulf of Mexico means that if a lengthening/shortening of the Loop Current is observed this introduces some predictability in the downstream flow for a few weeks.

Description: The diffusive layering (DL) form of double-diffusive convection cools the Atlantic Water (AW) as it circulates around the Arctic Ocean. Large DL steps, with heights of homogeneous layers often greater than 10 m, have been found above the AW core in the Eurasian Basin (EB) of the eastern Arctic. Within these DL staircases, heat and salt fluxes are determined by the mechanisms for vertical transport through the high-gradient regions (HGRs) between the homogeneous layers. These HGRs can be thick (up to 5 m and more) and are frequently complex, being composed of multiple small steps or continuous stratification. Microstructure data collected in the EB in 2007 and 2008 are used to estimate heat fluxes through large steps in three ways: using the measured dissipation rate in the large homogeneous layers; utilizing empirical flux laws based on the density ratio and temperature step across HGRs after scaling to account for the presence of multiple small DL interfaces within each HGR; and averaging estimates of heat fluxes computed separately for individual small interfaces (as laminar conductive fluxes), small convective layers (via dissipation rates within small DL layers), and turbulent patches (using dissipation rate and buoyancy) within each HGR. Diapycnal heat fluxes through HGRs evaluated by each method agree with each other and range from ~2 to ~8 W m−2, with an average flux of ~3–4 W m−2. These large fluxes confirm a critical role for the DL instability in cooling and thickening the AW layer as it circulates around the eastern Arctic Ocean.

Description: Coastal trapped waves (CTWs) that propagate poleward along the southwest African shelf potentially leak energy from lower latitudes into the Benguela Upwelling System (BUS). Thus, in addition to local winds, these waves provide an important remote forcing mechanism for the upwelling region. The present study aims at elucidating the nature of CTWs in the northern BUS. To this end, we make use of multisite velocity observations from the Namibian shelf (18°, 20°, 23°S) and examine the alongshore velocity signal for signatures of CTWs by means of wavelet methods. We found that a substantial amount of energy is concentrated within a submonthly to subseasonal frequency band (10–50 days). Based on the coherence and phase spectra of the alongshelf currents, we provide evidence for a predominantly southward phase propagation and establish typical time and length scales of CTWs in the region. It turns out that their properties differ significantly within a few hundred kilometers along the coast. A comparison of the results with theoretical dispersion curves shows that this difference may be explained by variations in the bottom topography. Finally, we investigate the coupling of the alongshore currents with the coastal and equatorial wind stress and highlight regions of potential wave generation.

Description: The turbulent dissipation rate ɛ is a key parameter to many oceanographic processes. Recently gliders have been increasingly used as a carrier for microstructure sensors. Compared to conventional ship-based methods, glider-based microstructure observations allow for long duration measurements under adverse weather conditions, and at lower costs. The incident water velocity U is an input parameter for the calculation of the dissipation rate. Since U can not be measured using the standard glider sensor setup, the parameter is normally computed from a steady-state glider flight model. As ɛ scales with U2 or U4, depending whether it is computed from temperature or shear microstructure, flight model errors can introduce a significant bias. This study is the first to use measurements of in-situ glider flight, obtained with a profiling Doppler velocity log and an electromagnetic current meter, to test and calibrate a flight model, extended to include inertial terms. Compared to a previously suggested flight model, the calibrated model removes a bias of approximately 1 cm s−1 in the incident water velocity, which translates to roughly a factor of 1.2 in estimates of the dissipation rate. The results further indicate that 90% of the estimates of the dissipation rate from the calibrated model are within a factor of 1.1 and 1.2 for measurements derived from microstructure temperature sensors and shear probes, respectively. We further outline the range of applicability of the flight model.

Description: Oceanic eddies are an important component in preconditioning the central Labrador Sea (LS) for deep convection and in restratifying the convected water. This study investigates the different sources and impacts of Eddy Kinetic Energy (EKE) and its temporal variability in the LS with the help of a 52-year long hindcast simulation of a 1/20° ocean model. Irminger Rings (IR) are generated in the West Greenland Current (WGC) between 60 and 62°N, mainly affect preconditioning and limit the northward extent of the convection area. The IR exhibit a seasonal cycle and decadal variations linked to the WGC strength, varying with the circulation of the subpolar gyre. The mean and temporal variations of IR generation can be attributed to changes in deep ocean baroclinic and upper ocean barotropic instabilities at comparable magnitudes. The main source of EKE and restratification in the central LS are Convective Eddies (CE). They are generated by baroclinic instabilities near the bottom of the mixed layer during and after convection. The CE have a mid-depth core and reflect the hydrographic properties of the convected water mass with a distinct minimum in potential vorticity. Their seasonal to decadal variability is tightly connected to the local atmospheric forcing and the associated air-sea heat fluxes. A third class of eddies in the LS are the Boundary Current Eddies shed from the Labrador Current (LC). Since they are mostly confined to the vicinity of the LC, these eddies appear to exert only minor influence on preconditioning and restratification.

Description: The Eddy Kinetic Energy (EKE) associated with the Subtropical Countercurrent (STCC) in the western subtropical South Pacific is known to exhibit substantial seasonal and decadal variability. Using an eddy-permitting ocean general circulation model, which is able to reproduce the observed, salient features of the seasonal cycles of shear, stratification, baroclinic production and the associated EKE, we investigate the decadal changes of EKE. We show that the STCC region exhibits, uniquely among the subtropical gyres of the world’s oceans, significant, atmospherically forced, decadal EKE variability. The decadal variations are driven by changing vertical shear between the STCC in the upper 300 m and the South Equatorial Current below, predominantly caused by variations in STCC strength associated with a changing meridional density gradient. In the 1970s, an increased meridional density gradient results in EKE twice as large as in later decades in the model. Utilizing sensitivity experiments, decadal variations in the wind field are shown to be the essential driver. Local wind stress curl anomalies associated with the Interdecadal Pacific Oscillation (IPO) lead to up- and downwelling of the thermocline, inducing strengthening or weakening of the STCC and the associated EKE. Additionally, remote wind stress curl anomalies in the eastern subtropical South Pacific, which are not related to the IPO, generate density anomalies that propagate westward as Rossby waves and can account for up to 30–40 % of the density anomalies in the investigated region.

Description: The North Atlantic Oscillation (NAO) and the Arctic Oscillation (AO) describe the dominant part of the variability in the Northern Hemisphere extratropical troposphere. Due to the strong connection of these patterns with surface climate, recent years have shown an increased interest and an increasing skill in forecasting them. However, it is unclear what the intrinsic limits of short-term predictability for the NAO and AO patterns are. This study compares the variability and predictability of both patterns, using a range of data and index computation methods for the daily NAO/AO indices. Small deviations from Gaussianity are found and characteristic decorrelation time scales of around one week. In the analysis of the Lyapunov spectrum it is found that predictability is not significantly different between the AO and NAO or between reanalysis products. Differences exist however between the indices based on EOF analysis, which exhibit predictability time scales around 12 - 16 days, and the station-based indices, exhibiting a longer predictability of 18 - 20 days. Both of these time scales indicate predictability beyond that currently obtained in ensemble prediction models for short-term predictability. Additional longer-term predictability for these patterns may be gained through local feedbacks and remote forcing mechanisms for particular atmospheric conditions.

Description: The formation of secondary ice in clouds, i.e. ice particles that are created at temperatures above the limit for homogeneous freezing without the direct involvement of a heterogeneous ice nucleus is one of the longest standing puzzles in cloud physics. Here we present comprehensive laboratory investigations on the formation of small ice particles upon the freezing of drizzle-sized cloud droplets levitated in an electrodynamic balance. Four different categories of secondary ice formation (bubble bursting, jetting, cracking, breakup) could be detected and their respective frequencies of occurrence as a function of temperature and droplet size are given. We find that bubble bursting occurs more often than droplet splitting. While we do not observe the shattering of droplets into many large fragments, we find that the average number of small secondary ice particles released during freezing is strongly droplet-size dependent and may well exceed unity for droplets larger than 300 μm in diameter. This leaves droplet fragmentation an important secondary ice process effective at temperatures around -10 °C in clouds where large drizzle droplets are present.

Description: Predicting tropical cyclone (TC) activity becomes more important every year while the understanding of what factors impact them continues to be complicated. El Niño–Southern Oscillation (ENSO) is one of the primary factors impacting the activities in both the Pacific and the Atlantic, but an extensive examination of the fluctuation in this system has yet to be studied in its entirety. This article analyzes the ENSO impacts on the Atlantic tropical cyclone activity during the assessed warm and cold years to show the dominant centennial-scale variation impact. This study looks to plausibly link this variation to the Southern Ocean centennial variability, which is rarely mentioned in any factors affecting the Atlantic tropical cyclone activity. This centennial variability could be used to enhance future work related to predicting tropical cyclones.

Description: Downward wave coupling (DWC) is an important process that characterizes the dynamical coupling between the stratosphere and troposphere via planetary wave reflection. A recent modeling study indicated that natural forcing factors, including sea-surface temperature variability and quasi-biennial oscillation, influence DWC and the associated surface impact in the Northern Hemisphere (NH). In light of this, we further investigate how DWC in the NH is affected by anthropogenic forcings, using a fully coupled chemistry-climate model CESM1 (WACCM). The results indicate that the occurrence of DWC is significantly suppressed in the future, starting later in the seasonal cycle, with more events concentrated in late winter (February-March). The future decrease in DWC events is associated with enhanced wave absorption in the stratosphere due to increased greenhouse gases. The enhanced wave absorption is manifest as more absorbing types of stratospheric sudden warmings, with more events concentrated in early winter. This early winter condition leads to a delay in the development of the upper stratospheric reflecting surface, resulting in a shift in the seasonal cycle of DWC towards late winter. The tropospheric responses to DWC events in the future exhibit different spatial patterns compared to those of the past. In the North Atlantic sector, DWC-induced circulation changes are characterized by a poleward shift and an eastward extension of the tropospheric jet, while in the North Pacific sector, the circulation changes are characterized by a weakening of the tropospheric jet. These responses are consistent with a change in the pattern of DWC-induced synoptic-scale eddy-mean flow interaction.

Description: Decadal variabilities in Indian Ocean subsurface ocean heat content (OHC; 50–300 m) since the 1950s are examined using ocean reanalyses. This study elaborates on how Pacific variability modulates the Indian Ocean on decadal time scales through both oceanic and atmospheric pathways. High correlations between OHC and thermocline depth variations across the entire Indian Ocean Basin suggest that OHC variability is primarily driven by thermocline fluctuations. The spatial pattern of the leading mode of decadal Indian Ocean OHC variability closely matches the regression pattern of OHC on the interdecadal Pacific oscillation (IPO), emphasizing the role of the Pacific Ocean in determining Indian Ocean OHC decadal variability. Further analyses identify different mechanisms by which the Pacific influences the eastern and western Indian Ocean. IPO-related anomalies from the Pacific propagate mainly through oceanic pathways in the Maritime Continent to impact the eastern Indian Ocean. By contrast, in the western Indian Ocean, the IPO induces wind-driven Ekman pumping in the central Indian Ocean via the atmospheric bridge, which in turn modifies conditions in the southwestern Indian Ocean via westward-propagating Rossby waves. To confirm this, a linear Rossby wave model is forced with wind stresses and eastern boundary conditions based on reanalyses. This linear model skillfully reproduces observed sea surface height anomalies and highlights both the oceanic connection in the eastern Indian Ocean and the role of wind-driven Ekman pumping in the west. These findings are also reproduced by OGCM hindcast experiments forced by interannual atmospheric boundary conditions applied only over the Pacific and Indian Oceans, respectively.

Description: Gap flows and the stable boundary layer were studied in northwest Greenland during the aircraft-based Investigation of Katabatic Winds and Polynyas during Summer (IKAPOS) experiment in June 2010. The measurements were performed using the research aircraft POLAR 5 of Alfred Wegener Institute (AWI; Bremerhaven). Besides navigational and basic meteorological instrumentation, the aircraft was equipped with radiation and surface temperature sensors and a turbulence measurement system. In the area of Smith Sound at the southern end of the Nares Strait, a stable, but fully turbulent, boundary layer with strong winds of up to 22 m s−1 was found during conditions of synoptically induced northerly winds through the Nares Strait. Strong surface inversions were present in the lowest 100–200 m. As a consequence of channeling effects, a well-pronounced low-level jet system was documented for each of four flights. The wind maximum is located at 20–50-km distance from the exit of Smith Sound. The 3D boundary layer structure past this gap is studied in detail. The channeling process is consistent with gap flow theory. The flow through the gap and over the surrounding mountains leads to the lowering of isotropic surfaces and the acceleration of the flow. The orographically channeled flow through Smith Sound plays a key role for the formation of the North Water polynya being the largest ice-producing polynya in the Arctic.

Description: This study analyzes the response of the Atlantic meridional overturning circulation (AMOC) to different CO2 concentrations and two ice sheet configurations in simulations with the coupled climate model MPI-ESM. With preindustrial (PI) ice sheets, there are two different AMOC states within the studied CO2 range: one state with a strong and deep upper overturning cell at high CO2 concentrations and one state with a weak and shallow upper cell at low CO2 concentrations. Changes in AMOC variability with decreasing CO2 indicate two stability thresholds. The strong state is stable above the first threshold near 217 ppm, and the weak state is stable below the second threshold near 190 ppm. Between the two thresholds, both states are marginally unstable, and the AMOC oscillates between them on millennial time scales. The weak AMOC state is stable when Antarctic Bottom Water becomes dense and salty enough to replace North Atlantic Deep Water (NADW) in the deep North Atlantic and when the density gain over the North Atlantic becomes too weak to sustain continuous NADW formation. With Last Glacial Maximum (LGM) ice sheets, the density gain over the North Atlantic and the northward salt transport are enhanced with respect to the PI ice sheet case. This enables active NADW formation and a strong AMOC for the entire range of studied CO2 concentrations. The AMOC variability indicates that the simulated AMOC is far away from a stability threshold with LGM ice sheets. The nonlinear relationship among AMOC, CO2, and prescribed ice sheets provides an explanation for the large intermodel spread of AMOC states found in previous coupled LGM simulations.

Description: Major sudden stratospheric warmings (SSWs) are extreme events during boreal winter, which not only impact tropospheric weather up to three months but also can influence oceanic variability through wind stress and heat flux anomalies. In the North Atlantic region, SSWs have the potential to modulate deep convection in the Labrador Sea and thereby the strength of the Atlantic meridional overturning circulation. The impact of SSWs on the Northern Hemisphere surface climate is investigated in two coupled climate models: a stratosphere-resolving (high top) and a non-stratosphere-resolving (low top) model. In both configurations, a robust link between SSWs and a negative NAO is detected, which leads to shallower-than-normal North Atlantic mixed layer depth. The frequency of SSWs and the persistence of this link is better captured in the high-top model. Significant differences occur over the Pacific region, where an unrealistically persistent Aleutian low is observed in the low-top configuration. An overrepresentation of SSWs during El Nino conditions in the low-top model is the main cause for this artifact. Our results underline the importance of a proper representation of the stratosphere in a coupled climate model for a consistent surface response in both the atmosphere and the ocean, which, among others, may have implications for oceanic deep convection in the subpolar North Atlantic.

Description: Benthic storms are important for both the energy budget of the ocean and for sediment resuspension and transport. Using 30 years of output from a high-resolution model of the North Atlantic, it is found that most of the benthic storms in the model occur near the western boundary in association with the Gulf Stream and the North Atlantic Current, in regions that are generally co-located with the peak near-bottom eddy kinetic energy. A common feature are meander troughs in the near-surface jets that are accompanied by deep low pressure anomalies spinning up deep cyclones with near-bottom velocities of up to more than 0.5 m/s. A case study of one of these events shows the importance of both baroclinic and barotropic instability of the jet, with energy being extracted from the jet in the upstream part of the meander trough and partly returned to the jet in the downstream part of the meander trough. This motivates examining the 30-year time mean of the energy transfer from the (annual mean) background flow into the eddy kinetic energy. This quantity is shown to be co-located well with the region in which benthic storms and large increases in deep cyclonic relative vorticity occur most frequently, suggesting an important role for mixed barotropic-baroclinic instability driven cyclogenesis in generating benthic storms throughout the model simulation. Regions of largest energy transfer and most frequent benthic storms are found to be the Gulf Stream west of the New England Seamounts and the North Atlantic Current near Flemish Cap.

Description: In this study, we propose a statistical method to validate sea-level reconstructions using geological records known as sea-level indicators (SLIs). SLIs are often the only available data to retrace late-glacial relative sea level (RSL). Determining the RSL from SLI height is not straight forward, the elevation at which an SLI was found usually does not represent the past RSL. In contrast, it has to be related to past RSL by investigating sample’s type, habitat and deposition conditions. For instance, water distribution at which a specific specimen is found today can be related to the indicator's depositional height range. Furthermore, the precision of dating varies between geological samples, and, in case of radiocarbon dating, the age has to be calibrated using a non-linear calibration curve. To avoid an a-priori assumption like normal-distributed uncertainties, we define likelihood functions which take into account the indicative meaning’s available error information and calibration statistics represented by joint probabilities. For this conceptional study, we restrict ourselves to one type of indicators, shallow-water shells, which are usually considered as low-grade samples giving only a lower limit of former sea level, as the depth range in which they live spreads over several tens of meters, and does not follow a normal distribution. The presented method is aimed to serve as a strategy for glacial isostatic adjustment reconstructions, in this case for the German Paleo-Climate Modelling Initiative PalMod (https://www.palmod.de/en) and by extending it to other SLI types.

Description: This study analyzes the response of the Atlantic meridional overturning circulation (AMOC) to different CO2 concentrations and two ice sheet configurations in simulations with the coupled climate model MPI-ESM. With preindustrial (PI) ice sheets, there are two different AMOC states within the studied CO2 range: one state with a strong and deep upper overturning cell at high CO2 concentrations and one state with a weak and shallow upper cell at low CO2 concentrations. Changes in AMOC variability with decreasing CO2 indicate two stability thresholds. The strong state is stable above the first threshold near 217 ppm, and the weak state is stable below the second threshold near 190 ppm. Between the two thresholds, both states are marginally unstable, and the AMOC oscillates between them on millennial time scales. The weak AMOC state is stable when Antarctic Bottom Water becomes dense and salty enough to replace North Atlantic Deep Water (NADW) in the deep North Atlantic and when the density gain over the North Atlantic becomes too weak to sustain continuous NADW formation. With Last Glacial Maximum (LGM) ice sheets, the density gain over the North Atlantic and the northward salt transport are enhanced with respect to the PI ice sheet case. This enables active NADW formation and a strong AMOC for the entire range of studied CO2 concentrations. The AMOC variability indicates that the simulated AMOC is far away from a stability threshold with LGM ice sheets. The nonlinear relationship among AMOC, CO2, and prescribed ice sheets provides an explanation for the large intermodel spread of AMOC states found in previous coupled LGM simulations.

Description: Dynamic vegetation models simulate global vegetation in terms of fractional coverages of a few plant functional types (PFTs). Although these models often share the same concept, they differ with respect to the number and kind of PFTs, complicating the comparability of simulated vegetation distributions. Pollen-based reconstructions are initially only available in form of time-series of individual taxa that are not distinguished in the models. Thus, to evaluate simulated vegetation distributions, the modelling results and pollen-based reconstructions have to be converted into a comparable format. The classical approach is the method of biomisation, but hitherto, PFT-based biomisation methods were only available for individual models. We introduce and evaluate a simple, universally applicable technique to harmonize PFT-distributions by assigning them into nine mega-biomes that follow the definitions commonly used for vegetation reconstructions. The method works well for all state-of the art dynamic vegetation models, independent of the spatial resolution or the complexity of the models. Large biome belts (such as tropical forest) are well represented, but regionally confined biomes (warm-mixed forest, Savanna) are only partly captured. Overall, the PFT-based biomisation is able to keep up with the conventional biomisation approach of forcing biome models (here: BIOME1) with the background climate states. The new method has, however, the advantage that it allows a more direct comparison and evaluation of the vegetation distributions simulated by Earth System Models. Thereby, the new method provides a powerful tool for the evaluation of Earth System Models in general.

Description: We have developed a new module to calculate soil organic carbon (SOC) accumulation in perennially frozen ground in the land surface model JSBACH. Running this offline version of MPI-ESM we have modelled permafrost carbon accumulation and release from the Last Glacial Maximum (LGM) to the Pre-industrial (PI). Our simulated near-surface PI permafrost extent of 16.9Miokm2 is close to observational evidence. Glacial boundary conditions, especially ice sheet coverage, result in profoundly different spatial patterns of glacial permafrost extent. Deglacial warming leads to large-scale changes in soil temperatures, manifested in permafrost disappearance in southerly regions, and permafrost aggregation in formerly glaciated grid cells. In contrast to the large spatial shift in simulated permafrost occurrence, we infer an only moderate increase of total LGM permafrost area (18.3Miokm2) – together with pronounced changes in the depth of seasonal thaw. Reconstructions suggest a larger spread of glacial permafrost towards more southerly regions, but with a highly uncertain extent of non-continuous permafrost. Compared to a control simulation without describing the transport of SOC into perennially frozen ground, the implementation of our newly developed module for simulating permafrost SOC accumulation leads to a doubling of simulated LGM permafrost SOC storage (amounting to a total of ~150PgC). Despite LGM temperatures favouring a larger permafrost extent, simulated cold glacial temperatures – together with low precipitation and low CO2 levels – limit vegetation productivity and therefore prevent a larger glacial SOC build-up in our model. Changes in physical and biogeochemical boundary conditions during deglacial warming lead to an increase in mineral SOC storage towards the Holocene (168PgC at PI), which is below observational estimates (575PgC in continuous and discontinuous permafrost). Additional model experiments clarified the sensitivity of simulated SOC storage to model parameters, affecting long-term soil carbon respiration rates and simulated active layer depths. Rather than a steady increase in carbon release from the LGM to PI as a consequence of deglacial permafrost degradation, our results suggest alternating phases of soil carbon accumulation and loss as an effect of dynamic changes in permafrost extent, active layer depths, soil litter input, and heterotrophic respiration.

Description: The eastern boundary region off Angola encompasses a highly productive ecosystem important for the food security of the coastal population. The fish-stock distribution, however, undergoes large variability on intraseasonal, interannual, and longer time scales. These fluctuations are partly associated with large-scale warm anomalies that are often forced remotely from the equatorial Atlantic and propagate southward, reaching the Benguela upwelling off Namibia. Such warm events, named Benguela Niños, occurred in 1995 and in 2011. Here we present results from an underexplored extensive in situ dataset that was analyzed in the framework of a capacity-strengthening effort. The dataset was acquired within the Nansen Programme executed by the Food and Agriculture Organization of the United Nations and funded by the Norwegian government. It consists of hydrographic and velocity data from the Angolan continental margin acquired biannually during the main downwelling and upwelling seasons over more than 20 years. The mean seasonal changes of the Angola Current from 6° to 17°S are presented. During austral summer the southward Angola Current is concentrated in the upper 150 m. It strengthens from north to south, reaching a velocity maximum just north of the Angola Benguela Front. During austral winter the Angola Current is weaker, but deeper reaching. While the southward strengthening of the Angola Current can be related to the wind forcing, its seasonal variability is most likely explained by coastally trapped waves. On interannual time scales, the hydrographic data reveal remarkable variability in subsurface upper-ocean heat content. In particular, the 2011 Benguela Niño was preceded by a strong subsurface warming of about 2 years’ duration.

Description: Climate reconstructions based on proxy records recovered from marine sediments, such as alkenone records or geochemical parameters measured on foraminifera, play an important role in our understanding of the climate system. They provide information about the state of the ocean ranging back hundreds to millions of years and form the backbone of paleo-oceanography. However, there are many sources of uncertainty associated with the signal recovered from sediment archived proxies. These include seasonal or depth habitat biases in the recorded signal, a frequency dependent reduction in the amplitude of the recorded signal due to bioturbation of the sediment, aliasing of high frequency climate variation onto a nominally annual, decadal or centennial resolution signal, and additional sample processing and measurement error introduced when the proxy signal is recovered. Here we present a forward model for sediment archived proxies that jointly models the above processes, so that the magnitude of their separate and combined effects can be investigated. Applications include the interpretation and analysis of uncertainty in existing proxy records, parameter sensitivity analysis to optimize future studies, and the generation of pseudo-proxy records that can be used to test reconstruction methods. We provide examples, such as the simulation of individual foraminifera records, that demonstrate the usefulness of the forward model for paleoclimate studies. The model is implemented as a user-friendly R package, sedproxy, the use of which we hope will contribute to a better understanding of both the limitations and potential of marine sediment proxies to inform about past climate.

Description: Besides the zonal flow that dominates the seasonal and long-term variability in the equatorial Atlantic, energetic intraseasonal meridional velocity fluctuations are observed in large parts of the water column. We use 15 years of partly full-depth velocity data from an equatorial mooring at 23°W to investigate intraseasonal variability and specifically the downward propagation of intraseasonal energy from the near-surface into the deep ocean. Between 20 and 50 m, intraseasonal variability at 23°W peaks at periods between 30 and 40 days. It is associated with westward-propagating tropical instability waves, which undergo an annual intensification in August. At deeper levels down to about 2000 m considerable intraseasonal energy is still observed. A frequency–vertical mode decomposition reveals that meridional velocity fluctuations are more energetic than the zonal ones for periods 〈 50 days. The energy peak at 30–40 days and at vertical modes 2–5 excludes equatorial Rossby waves and suggests Yanai waves to be associated with the observed intraseasonal energy. Yanai waves that are considered to be generated by tropical instability waves propagate their energy from the near-surface west of 23°W downward and eastward to eventually reach the mooring location. The distribution of intraseasonal energy at the mooring position depends largely on the dominant frequency and the time, depth, and longitude of excitation, while the dominant vertical mode of the Yanai waves plays only a minor role. Observations also show the presence of weaker intraseasonal variability at 23°W below 2000 m that cannot be associated with tropical instability waves.

Description: Multi-year moored velocity observations of the Angola Current near 11°S reveal a weak southward mean flow superimposed by substantial intraseasonal to seasonal variability, including annual and semiannual cycles with distinct baroclinic structures. In the equatorial Atlantic these oscillations are associated with basin-mode resonances of the fourth and second baroclinic modes, respectively. Here, the role of basin-mode resonance and local forcing for the Angola Current seasonality are investigated. A suite of linear shallow-water models for the tropical Atlantic is employed, each model representing a single baroclinic mode forced at a specific period. The annually and semiannually oscillating forcing is given by 1) an idealized zonally uniform zonal forcing restricted to the equatorial band corresponding to a remote equatorial forcing or 2) realistic, spatially-varying Fourier components of wind stress data that include local forcing off Angola, particularly alongshore winds. Model-computed modal amplitudes are scaled to match moored velocity observations from the equatorial Atlantic. The observed annual cycle of alongshore velocity at 11°S is well reproduced by the remote equatorial forcing. Including local forcing slightly improves the agreement between observed and simulated semiannual oscillations at 11°S compared to the purely equatorial forcing. However, the model-computed semiannual cycle lacks amplitude at mid-depth. This could be the result of either underestimating the strength of the second equatorial basin-mode of the fourth baroclinic mode or other processes not accounted for in the shallow-water models. Overall, our findings underline the importance of large-scale linear equatorial wave dynamics for the seasonal variability of the boundary circulation off Angola.

Description: The Lagrangian analysis of sets of particles advected with the flow fields of ocean models are used to study connectivity, i.e. exchange pathways, timescales and volume transports, between distinct oceanic regions. One important factor influencing the dispersion of fluid particles and hence connectivity is the Lagrangian eddy diffusivity, which quantifies the influence of turbulent processes on the rate of particle dispersal. Due to spatial and temporal discretization, turbulence is not fully resolved in modelled velocities, and the concept of eddy diffusivity is used to parametrize the impact of unresolved processes. However, the relations between observational- and model-based Lagrangian eddy diffusivity estimates as well as eddy parameterizations are not clear. This study presents an analysis of the spatially variable near-surface lateral eddy diffusivity estimates obtained from Lagrangian trajectories simulated with 5-day mean velocities from an eddy-resolving ocean model (INALT01) for the Agulhas system. INALT01 features diffusive regimes for dynamically different regions, some of which exhibit strong suppression of eddy mixing by mean flow, and is consistent with the pattern and magnitude of drifter-based eddy diffusivity estimates. Using monthly-mean velocities decreases the estimated diffusivities less than eddy kinetic energy, supporting the idea that large and persistent eddy features dominate eddy diffusivities. For a non-eddying ocean model (ORCA05), Lagrangian eddy diffusivities are greatly reduced, in particular when the Gent and McWilliams parameterization of mesoscale eddies is employed.

Description: Arctic sea ice area (SIA) during late summer and early fall decreased substantially over the last four decades, and its decline accelerated beginning in the early 2000s. Statistical analyses of observations show that enhanced poleward moisture transport from the North Pacific to the Arctic Ocean contributed to the accelerated SIA decrease during the most recent period. As a consequence, specific humidity in the Arctic Pacific sector significantly increased along with an increase of downward longwave radiation beginning in 2002, which led to a significant acceleration in the decline of SIA in the Arctic Pacific sector. The resulting sea ice loss led to increased evaporation in the Arctic Ocean, resulting in a further increase of the specific humidity in mid-to-late fall, thus acting as a positive feedback to the sea ice loss. The overall set of processes is also found in a long control simulation of a coupled climate model.

Description: In the present study, the influence of some major tropical modes of variability on northern hemisphere regional blocking frequency variability during boreal winter is investigated. Reanalysis data and an ensemble experiment with the ECMWF model using relaxation towards the ERA-Interim reanalysis data inside the tropics are used. The tropical modes under investigation are El Niño Southern Oscillation (ENSO), the Madden-Julian Oscillation (MJO) and the upper tropospheric equatorial zonal-mean zonal wind . An early (late) MJO phase refers to the part of the MJO cycle when enhanced (suppressed) precipitation occurs over the western Indian Ocean and suppressed (enhanced) precipitation occurs over the Maritime Continent and the western tropical Pacific. Over the North Pacific sector, it is found that enhanced (suppressed) high latitude blocking occurs in association with El Niño (La Niña) events, late (early) MJO phases and westerly (easterly) . Over central to southern Europe and the east Atlantic, it is found that late MJO phases, as well as a suppressed MJO are leading to enhanced blocking frequency. Furthermore, early (late) MJO phases are followed by blocking anomalies over the western North Atlantic region, similar to those associated with a positive (negative) North Atlantic Oscillation. Over northern Europe, the easterly (westerly) phase of is associated with enhanced (suppressed) blocking. These results are largely confirmed by both the reanalysis and the model experiment.

Description: The Indian Ocean has sustained robust surface warming in recent decades, but the role of multi-decadal variability remains unclear. Using ocean model hindcasts, characteristics of low-frequency Indian Ocean temperature variations are explored. Simulated upper-ocean temperature changes across the Indian Ocean in the hindcast are consistent with those recorded in observational products and ocean reanalyses. Indian Ocean temperatures exhibit strong warming trends since the 1950s limited to the surface and south of 30°S, while extensive subsurface cooling occurs over much of the tropical Indian Ocean. Previous work focused on diagnosing causes of these long-term trends in the Indian Ocean over the second half of the 20th Century. Instead, the temporal evolution of Indian Ocean subsurface heat content is shown here to reveal distinct multi-decadal variations associated with the Pacific Decadal Oscillation and the long-term trends are thus interpreted to result from aliasing of the low-frequency variability. Transmission of the multi-decadal signal occurs via an oceanic pathway through the Indonesian Throughflow and is manifest across the Indian Ocean centered along 12°S as westward propagating Rossby waves modulating thermocline and subsurface heat content variations. Resulting low-frequency changes in the eastern Indian Ocean thermocline depth are associated with decadal variations in the frequency of Indian Ocean Dipole (IOD) events, with positive IOD events unusually common in the 1960s and 1990s with a relatively shallow thermocline. In contrast, the deeper thermocline depth in the 1970s and 1980s is associated with frequent negative IOD and rare positive IOD events. Changes in Pacific wind forcing in recent decades and associated rapid increases in Indian Ocean subsurface heat content can thus affect the basin’s leading mode of variability, with implications for regional climate and vulnerable societies in surrounding countries.

Description: The Argentine margin contains important sedimentological, paleontological and chemical records of regional and local tectonic evolution, sea level, climate evolution and ocean circulation since the opening of the South Atlantic in the Late Jurassic–Early Cretaceous as well as the present-day results of post-depositional chemical and biological alteration. Despite its important location, which underlies the exchange of southern- and northern-sourced water masses, the Argentine margin has not been investigated in detail using scientific drilling techniques, perhaps because the margin has the reputation of being erosional. However, a number of papers published since 2009 have reported new high-resolution and/or multichannel seismic surveys, often combined with multi-beam bathymetric data, which show the common occurrence of layered sediments and prominent sediment drifts on the Argentine and adjacent Uruguayan margins. There has also been significant progress in studying the climatic records in surficial and near-surface sediments recovered in sediment cores from the Argentine margin. Encouraged by these recent results, our 3.5-day IODP (International Ocean Discovery Program) workshop in Buenos Aires (8–11 September 2015) focused on opportunities for scientific drilling on the Atlantic margin of Argentina, which lies beneath a key portion of the global ocean conveyor belt of thermohaline circulation. Significant opportunities exist to study the tectonic evolution, paleoceanography and stratigraphy, sedimentology, and biosphere and geochemistry of this margin.

Description: The 13th European Polar Low Workshop was organized by the European Polar Low Working Group (www.uni-trier.de/index.php?id=20308)and gathered scientists from nine countries focusing on polar mesocyclones in both hemispheres and other mesoscale weather phenomena such as katabatic winds, tip jets, boundary layer fronts, cold air outbreaks, and weather extremes in polar regions. Topics included experimental, climatological, theoretical, modeling, and remote sensing studies. The aim was to bring together scientists and forecasters to present their latest work and recent findings on these topics and to encourage discussions on improving forecasting and understanding of these phenomena.

Description: The regional climate model COSMOin Climate Limited-AreaMode (COSMO-CLM or CCLM) is used with a high resolution of 15km for the entire Arctic for all winters 2002/03–2014/15. The simulations show a high spatial and temporal variability of the recent 2-m air temperature increase in the Arctic. The maximum warming occurs north of Novaya Zemlya in the Kara Sea and Barents Sea between March 2003 and 2012 and is responsible for up to a 208C increase. Land-based observations confirm the increase but do not cover the maximum regions that are located over the ocean and sea ice.Also, the 30-km version of theArctic SystemReanalysis (ASR) is used to verify the CCLM for the overlapping time period 2002/03–2011/12. The differences between CCLM and ASR 2-m air temperatures vary slightly within 18C for the ocean and sea ice area. Thus,ASR captures the extreme warming as well. The monthly 2-m air temperatures of observations and ERA-Interim data show a large variability for the winters 1979–2016. Nevertheless, the air temperature rise since the beginning of the twenty-first century is up to 8 times higher than in the decades before. The sea ice decrease is identified as the likely reason for the warming. The vertical temperature profiles show that the warming has a maximum near the surface, but a 0.58Cyr21 increase is found up to 2 km. CCLM, ASR, and also the coarser resolved ERA-Interim data show that February and March are the months with the highest 2-m air temperature increases, averaged over the ocean and sea ice area north of 708N; for CCLM the warming amounts to an average of almost 58C for 2002/03–2011/12.

Description: The ratio of global mean surface air temperature change to cumulative CO2 emissions, referred to as transient climate response to cumulative CO2 emissions (TCRE), has been shown to be approximately constant on centennial time scales. The mechanisms behind this constancy are not well understood, but previous studies suggest that compensating effects of ocean heat and carbon fluxes, which are governed by the same ocean mixing processes, could be one cause for this approximate constancy. This hypothesis is investigated by forcing different versions of the University of Victoria Earth System Climate Model, which differ in the ocean mixing parameterization, with an idealized scenario of 1% annually increasing atmospheric CO2 until quadrupling of the preindustrial CO2 concentration and constant concentration thereafter. The relationship between surface air warming and cumulative emissions remains close to linear, but the TCRE varies between model versions, spanning the range of 1.2°–2.1°C EgC−1 at the time of CO2 doubling. For all model versions, the TCRE is not constant over time while atmospheric CO2 concentrations increase. It is constant after atmospheric CO2 stabilizes at 1120 ppm, because of compensating changes in temperature sensitivity (temperature change per unit radiative forcing) and cumulative airborne fraction. The TCRE remains approximately constant over time even if temperature sensitivity, determined by ocean heat flux, and cumulative airborne fraction, determined by ocean carbon flux, are taken from different model versions with different ocean mixing settings. This can partially be explained with temperature sensitivity and cumulative airborne fraction following similar trajectories, which suggests ocean heat and carbon fluxes scale approximately linearly with changes in vertical mixing.

Description: The Weddell Sea polynya is a large opening in the open-ocean sea ice cover associated with intense deep convection in the ocean. A necessary condition to form and maintain a polynya is the presence of a strong subsurface heat reservoir. This study investigates the processes that control the stratification and hence the buildup of the subsurface heat reservoir in the Weddell Sea. To do so, a climate model run for 200 years under preindustrial forcing with two eddying resolutions in the ocean (0.25° CM2.5 and 0.10° CM2.6) is investigated. Over the course of the simulation, CM2.6 develops two polynyas in the Weddell Sea, while CM2.5 exhibits quasi-continuous deep convection but no polynyas, exemplifying that deep convection is not a sufficient condition for a polynya to occur. CM2.5 features a weaker subsurface heat reservoir than CM2.6 owing to weak stratification associated with episodes of gravitational instability and enhanced vertical mixing of heat, resulting in an erosion of the reservoir. In contrast, in CM2.6, the water column is more stably stratified, allowing the subsurface heat reservoir to build up. The enhanced stratification in CM2.6 arises from its refined horizontal grid spacing and resolution of topography, which allows, in particular, a better representation of the restratifying effect by transient mesoscale eddies and of the overflows of dense waters along the continental slope.

Description: The current generation of marine biogeochemical modules in Earth system models (ESMs) considers mainly the effect of marine biota on the carbon cycle. We propose to also implement other biologically driven mechanisms in ESMs so that more climate-relevant feedbacks are captured. We classify these mechanisms in three categories according to their functional role in the Earth system: (1) "biogeochemical pumps", which affect the carbon cycling; (2) "biological gas and particle shuttles", which affect the atmospheric composition; and (3) "biogeophysical mechanisms", which affect the thermal, optical, and mechanical properties of the ocean. To resolve mechanisms from all three classes, we find it sufficient to include five functional groups: bulk phyto- and zooplankton, calcifiers, and coastal gas and surface mat producers. We strongly suggest to account for a larger mechanism diversity in ESMs in the future to improve the quality of climate projections.

Description: Seasonal variability in pathways of warm water masses toward the Kangerdlugssuaq Fjord-Glacier system (KF/KG), southeast Greenland, is investigated by backtracking Lagrangian particles seeded at the fjord mouth in a high-resolution regional ocean model simulation in the ice-free and the ice-covered seasons. The waters at KF are a mixture of Atlantic-origin water advected from the Irminger Basin (FF for Faxaflói), the deep waters from the Denmark Strait and the waters from the Arctic Ocean, both represented by the Kögur section (KO). Below 200m depth, the warm water is a mixture of FF and KO water masses, and is warmer in winter than in summer. We find that seasonal differences in pathways double the fraction of FF particles in winter, causing the seasonal warming and salinification. Seasonal temperature variations at the upstream sections (FF and KO) have a negligible impact on temperature variations near the fjord. Successful monitoring of heat flux to the fjord therefore needs to take place close to the fjord, and cannot be inferred from upstream conditions.

Description: By performing two sets of high-resolution atmospheric general circulation model (AGCM) experiments, we find that the atmospheric response to a sea surface temperature (SST) anomaly in the extratropical North Pacific is sensitive to decadal variations of the background SST on which the SST anomaly is superimposed. The response in the first set of experiments, in which the SST anomaly is superimposed on the observed daily SST of 1981-1990, strongly differs from the response in the second experiment, in which the same SST anomaly is superimposed on the observed daily SST of 1991-2000. The atmospheric response over the North Pacific during 1981-1990 is eddy-mediated, equivalent barotropic and concentrated in the east. In contrast, the atmospheric response during 1991-2000 is weaker and strongest in the west. The results are discussed in terms of Rossby wave dynamics, with the proposed primary wave source switching from baroclinic eddy vorticity forcing over the eastern North Pacific in 1981-1990 to mean flow divergence over the western North Pacific in 1991-2000. The wave source changes are linked to the decadal reduction of daily SST variability over the eastern North Pacific and strengthening of the Oyashio Extension front over the western North Pacific. Thus, both daily and frontal aspects of the background SST variability in determining the atmospheric response to extratropical North Pacific SST anomalies are emphasized by our AGCM experiments.

Description: Warm water of open ocean origin on the continental shelf of the Amundsen and Bellingshausen Seas causes the highest basal melt rates reported for Antarctic ice shelves with severe consequences for the ice shelf/ice sheet dynamics. Ice shelves fringing the broad continental shelf in the Weddell and Ross Seas melt at rates orders of magnitude smaller. However, simulations using coupled ice–ocean models forced with the atmospheric output of the HadCM3 SRES-A1B scenario run (CO2 concentration in the atmosphere reaches 700 ppmv by the year 2100 and stays at that level for an additional 100 years) show that the circulation in the southern Weddell Sea changes during the twenty-first century. Derivatives of Circumpolar Deep Water are directed southward underneath the Filchner–Ronne Ice Shelf, warming the cavity and dramatically increasing basal melting. To find out whether the open ocean will always continue to power the melting, the authors extend their simulations, applying twentieth-century atmospheric forcing, both alone and together with prescribed basal mass flux at the end of (or during) the SRES-A1B scenario run. The results identify a tipping point in the southern Weddell Sea: once warm water flushes the ice shelf cavity a positive meltwater feedback enhances the shelf circulation and the onshore transport of open ocean heat. The process is irreversible with a recurrence to twentieth-century atmospheric forcing and can only be halted through prescribing a return to twentieth-century basal melt rates. This finding might have strong implications for the stability of the Antarctic ice sheet.

Description: Kretzschmar et al., in a comment in 2017, use the spread in the output of aerosol–climate models to argue that the models refute the hypothesis (presented in a paper by Stevens in 2015) that for the mid-twentieth-century warming to be consistent with observations, then the present-day aerosol forcing, must be less negative than −1 W m−2. The main point of contention is the nature of the relationship between global SO2 emissions and In contrast to the concave (log-linear) relationship used by Stevens and in earlier studies, whereby becomes progressively less sensitive to SO2 emissions, some models suggest a convex relationship, which would imply a less negative lower bound. The model that best exemplifies this difference, and that is most clearly in conflict with the hypothesis of Stevens, does so because of an implausible aerosol response to the initial rise in anthropogenic aerosol precursor emissions in East and South Asia—already in 1975 this model’s clear-sky reflectance from anthropogenic aerosol over the North Pacific exceeds present-day estimates of the clear-sky reflectance by the total aerosol. The authors perform experiments using a new (observationally constrained) climatology of anthropogenic aerosols to further show that the effects of changing patterns of aerosol and aerosol precursor emissions during the late twentieth century have, for the same global emissions, relatively little effect on These findings suggest that the behavior Kretzschmar et al. identify as being in conflict with the lower bound in Stevens arises from an implausible relationship between SO2 emissions and and thus provides little basis for revising this lower bound.

Description: This article investigates the dynamics and temporal evolution of the Atlantic Multidecadal Variability (AMV) in a coupled climate model. The model contains a correction to the North Atlantic flow field to improve the path of the North Atlantic Current, thereby alleviating the surface cold bias, a common problem with climate models, and offering a unique opportunity to study the AMV in a model. Changes in greenhouse gas forcing or aerosol loading are not considered. A striking feature of our results is the contrast between the western and eastern sides of the subpolar gyre in the model. On the western side, heat supply from the ocean plays a major role, with most of this heat being given up to the atmosphere in the warm phase, largely symmetrically about the time of the AMV maximum. By contrast, on the eastern side, the ocean gains heat from the atmosphere, with relatively little role for ocean heat supply in the years before the AMV maximum. Thereafter, the balance changes with heat now being removed from the eastern side by the ocean leading to a reducing ocean heat content, behavior we associate with the establishment of an intergyre gyre at the time of the AMV maximum. In the warm phase, melting sea-ice leads to a freshening of surface waters northeast of Greenland which travel southward into the Irminger and Labrador Sea, shutting down convection and terminating the AMV warm phase.

Description: For decades oceanographers have understood the Atlantic meridional overturning circulation (AMOC) to be primarily driven by changes in the production of deep-water formation in the subpolar and subarctic North Atlantic. Indeed, current Intergovernmental Panel on Climate Change (IPCC) projections of an AMOC slowdown in the twenty-first century based on climate models are attributed to the inhibition of deep convection in the North Atlantic. However, observational evidence for this linkage has been elusive: there has been no clear demonstration of AMOC variability in response to changes in deep-water formation. The motivation for understanding this linkage is compelling, since the overturning circulation has been shown to sequester heat and anthropogenic carbon in the deep ocean. Furthermore, AMOC variability is expected to impact this sequestration as well as have consequences for regional and global climates through its effect on the poleward transport of warm water. Motivated by the need for a mechanistic understanding of the AMOC, an international community has assembled an observing system, Overturning in the Subpolar North Atlantic Program (OSNAP), to provide a continuous record of the transbasin fluxes of heat, mass, and freshwater, and to link that record to convective activity and water mass transformation at high latitudes. OSNAP, in conjunction with the Rapid Climate Change–Meridional Overturning Circulation and Heatflux Array (RAPID–MOCHA) at 26°N and other observational elements, will provide a comprehensive measure of the three-dimensional AMOC and an understanding of what drives its variability. The OSNAP observing system was fully deployed in the summer of 2014, and the first OSNAP data products are expected in the fall of 2017.

Description: Well-known problems trouble coupled general circulation models of the eastern Atlantic and Pacific Ocean basins. Model climates are significantly more symmetric about the equator than is observed. Model sea surface temperatures are biased warm south and southeast of the equator, and the atmosphere is too rainy within a band south of the equator. Near-coastal eastern equatorial SSTs are too warm, producing a zonal SST gradient in the Atlantic opposite in sign to that observed. The U.S. Climate Variability and Predictability Program (CLIVAR) Eastern Tropical Ocean Synthesis Working Group (WG) has pursued an updated assessment of coupled model SST biases, focusing on the surface energy balance components, on regional error sources from clouds, deep convection, winds, and ocean eddies; on the sensitivity to model resolution; and on remote impacts. Motivated by the assessment, the WG makes the following recommendations: 1) encourage identification of the specific parameterizations contributing to the biases in individual models, as these can be model dependent; 2) restrict multimodel intercomparisons to specific processes; 3) encourage development of high-resolution coupled models with a concurrent emphasis on parameterization development of finer-scale ocean and atmosphere features, including low clouds; 4) encourage further availability of all surface flux components from buoys, for longer continuous time periods, in persistently cloudy regions; and 5) focus on the eastern basin coastal oceanic upwelling regions, where further opportunities for observational–modeling synergism exist.

Description: Downward wave coupling occurs when an upward propagating planetary wave from the troposphere decelerates the flow in the upper stratosphere, and forms a downward reflecting surface that redirects waves back to the troposphere. To test this mechanism and potential factors influencing the downward wave coupling, three 145-year sensitivity simulations with NCAR’s Community Earth System Model (CESM-WACCM), a state-of-the-art high-top chemistry-climate model, are analyzed. The results show that the QBO and SST variability significantly impact downward wave coupling. Without the QBO, the occurrence of downward wave coupling is significantly suppressed. In contrast, stronger and more persistent downward wave coupling occurs when SST variability is excluded. The above influence on the occurrence of downward wave coupling is mostly due to a direct influence of the QBO and SST variability on stratospheric planetary wave source and propagation. The strengths of the tropospheric circulation and surface responses to a given downward wave coupling event, however, behave differently. The surface anomaly is significantly weaker (stronger) in the experiment with fixed SSTs (without QBO), even though the statistical signal of downward coupling is strongest (weakest) in this experiment. This apparent mismatch is explained by the differences in the strength of the synoptic-scale eddy-mean flow feedback and the possible contribution of SST anomalies in the North Atlantic during DWC event. The weaker synoptic-scale eddy-mean flow feedback, and the absence of the positive NAO-related SST-tripole pattern in the fixed SST experiment are consistent with a weaker tropospheric response in this experiment. The results highlight the importance of synoptic-scale eddies in setting the tropospheric response to downward wave coupling.

Description: There is evidence that the strengthened stratospheric westerlies arising from the Antarctic ozone hole–induced cooling cause a polar mesospheric warming and a subsequent cooling in the lower thermosphere. While previous studies focus on the role of nonresolved (gravity) wave drag filtering, here the role of resolved (planetary) wave drag and radiative forcing on the Antarctic mesosphere and lower thermosphere (MLT) is explored in detail. Using simulations with NCAR’s Community Earth System Model, version 1 (Whole Atmosphere Community Climate Model) [CESM1(WACCM)], it is found that in late spring and early summer the anomalous polar mesospheric warming induced by easterly nonresolved wave drag is dampened by anomalous dynamical cooling induced by westerly resolved wave drag. This resolved wave drag is attributed to planetary-scale wave (k = 1–3) activity, which is generated in situ as a result of increased instability of the summer mesospheric easterly jet induced by the ozone hole. On the other hand, the anomalous cooling in the polar lower thermosphere induced by westerly nonresolved wave drag is enhanced by anomalous dynamical cooling due to westerly resolved wave drag. In addition, radiative effects from increased greenhouse gases during the ozone hole period contribute partially to the cooling in the polar lower thermosphere. The polar MLT temperature response to the Antarctic ozone hole is, through thermal wind balance, accompanied by the downward migration of anomalous zonal-mean wind from the lower thermosphere to the stratopause. The results highlight that a proper accounting of both dynamical and radiative effects is required in order to correctly attribute the causes of the polar MLT response to the Antarctic ozone hole.

Description: The transport of dissolved oxygen (O2) from the surface ocean into the interior is a critical process sustaining aerobic life in mesopelagic ecosystems, but its rates and sensitivity to climate variations are poorly understood. Using a circulation model constrained to historical variability by assimilation of observations, we show that the North Pacific thermocline effectively takes up O2 primarily by expanding the area through which O2-rich mixed layer water is detrained into the thermocline. The outcrop area during the critical winter season varies in concert with the Pacific Decadal Oscillation (PDO). When the central North Pacific Ocean is in a cold phase, the winter outcrop window for the Central Mode Water class (CMW; a neutral density range of γ = 25.6 - 26.6) expands southward allowing more O2-rich surface water to enter the ocean’s interior. An increase in volume flux of water to the CMW density class is partly compensated by a reduced supply to the shallower densities of Subtropical Mode Water (γ = 24.0 - 25.5). The thermocline has become better oxygenated since the 1980s due partly to strong O2 uptake. Positive O2 anomalies appear first near the outcrop and subsequently downstream in the subtropical gyre. In contrast to the O2 variations within the ventilated thermocline, observed O2 in Intermediate Water (density range of γ = 26.7 – 27.2) shows a declining trend over the past half-century, a trend not explained by the open ocean water mass formation rate.

Description: This study investigates the interaction of the Quasi-Biennial Oscillation (QBO) and the El Niño-Southern Oscillation (ENSO) in the troposphere separately for the North Pacific and North Atlantic region. Three 145-year model simulations with NCAR’s Community Earth Sytem Model (CESM-WACCM) are analyzed where only natural and no anthropogenic forcings are considered. These long simulations allow us to obtain statistically reliable results from an exceptional large number of cases for each combination of the QBO (westerly and easterly) and ENSO phases (El Niño and La Niña). Two different analysis methods were applied to investigate where nonlinearity might play a role in QBO-ENSO interactions. The analyses reveal that the stratospheric equatorial QBO anomalies extend down to the troposphere over the North Pacific during Northern hemisphere winter only during La Niña and not during El Niño events. The Aleutian low is deepened during QBO westerly (QBOW) as compared to QBO easterly (QBOE) conditions, and the North Pacific subtropical jet is shifted northward during La Niña. In the North Atlantic, the interaction of QBOW with La Niña conditions (QBOE with El Niño) results in a positive (negative) North Atlantic Oscillation (NAO) pattern. For both regions, nonlinear interactions between the QBO and ENSO might play a role. The results provide potential to enhance the skill of tropospheric seasonal predictions in the North Atlantic and North Pacific region.

Description: The Equatorial Deep Jets (EDJs) are an ubiquitous feature of the equatorial oceans; in the Atlantic Ocean, they are the dominant mode of interannual variability of the zonal flow at intermediate depth. On the basis of more than 10 years of moored observations of zonal velocity at 23°W, the vertically propagating EDJs are best described as superimposed oscillations of the 13th to the 23th baroclinic modes with a dominant oscillation period for all modes of 1650 days. This period is close to the resonance period of the respective gravest equatorial basin mode for the dominant vertical modes 16 and 17. It is argued that since the equatorial basin mode is composed of linear equatorial waves, a linear reduced gravity model can be employed for each baroclinic mode, driven by spatially homogeneous zonal forcing oscillating with the EDJ period. The fit of the model solutions to observations at 23°W yields a basin wide reconstruction of the EDJs and the associated vertical structure of their forcing. From the resulting vertical profile of mean power input and vertical energy flux on the equator, it follows that the EDJs are locally maintained over a considerable depth range, from 500-2500 m, with the maximum power input and vertical energy flux at 1300 m. The strong dissipation closely ties the apparent vertical propagation of energy to the vertical distribution of power input and, together with the EDJs’ prevailing downward phase propagation, require the phase of the forcing of the EDJs to propagate downward.

Description: A German national project coordinates research on improving a global decadal climate prediction system for future operational use. MiKlip, an eight-year German national research project on decadal climate prediction, is organized around a global prediction system comprising the climate model MPI-ESM together with an initialization procedure and a model evaluation system. This paper summarizes the lessons learned from MiKlip so far; some are purely scientific, others concern strategies and structures of research that targets future operational use. Three prediction-system generations have been constructed, characterized by alternative initialization strategies; the later generations show a marked improvement in hindcast skill for surface temperature. Hindcast skill is also identified for multi-year-mean European summer surface temperatures, extra-tropical cyclone tracks, the Quasi-Biennial Oscillation, and ocean carbon uptake, among others. Regionalization maintains or slightly enhances the skill in European surface temperature inherited from the global model and also displays hindcast skill for wind-energy output. A new volcano code package permits rapid modification of the predictions in response to a future eruption. MiKlip has demonstrated the efficacy of subjecting a single global prediction system to a major research effort. The benefits of this strategy include the rapid cycling through the prediction-system generations, the development of a sophisticated evaluation package usable by all MiKlip researchers, and regional applications of the global predictions. Open research questions include the optimal balance between model resolution and ensemble size, the appropriate method for constructing a prediction ensemble, and the decision between full-field and anomaly initialization. Operational use of the MiKlip system is targeted for the end of the current decade, with a recommended generational cycle of two to three years.

Description: Latent heat fluxes (LHF) play an essential role in the global energy budget and are thus important for understanding the climate system. Satellite-based remote sensing permits a large-scale determination of LHF, which, amongst others, are based on near-surface specific humidity qa. However, the qa random retrieval error (Etot) remains unknown. Here, a novel approach is presented to quantify the error contributions to pixel-level qa of the Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite (HOAPS, version 3.2) dataset. The methodology makes use of multiple triple collocation (MTC) analysis between 1995-2008 over the global ice-free oceans. Apart from satellite records, these datasets include selected ship records extracted from the Seewetteramt Hamburg (SWA) archive and the International Comprehensive Ocean-Atmosphere Data Set (ICOADS), serving as the in-situ ground reference. The MTC approach permits the derivation of Etot as the sum of model uncertainty EM and sensor noise EN, while random uncertainties due to in-situ measurement errors (Eins) and collocation (EC) are isolated concurrently. Results show an Etot average of 1.1 ± 0.3 g kg-1, whereas the mean EC (Eins) is in the order of 0.5 ± 0.1 g kg-1 (0.5 ± 0.3 g kg-1). Regional analyses indicate a maximum of Etot exceeding 1.5 g kg-1 within humidity regimes of 12-17 g kg-1, associated with the single-parameter, multilinear qa retrieval applied in HOAPS. Multi-dimensional bias analysis reveals that global maxima are located off the Arabian Peninsula.

Description: Seasonal variability of the tropical Atlantic circulation is dominated by the annual cycle, but semi-annual variability is also pronounced, despite weak forcing at that period. Here we use multi-year, full-depth velocity measurements from the central equatorial Atlantic to analyze the vertical structure of annual and semi-annual variations of zonal velocity. A baroclinic modal decomposition finds that the annual cycle is dominated by the 4th mode and the semi-annual cycle by the 2nd mode. Similar local behavior is found in a high-resolution general circulation model. This simulation reveals that the annual and semi-annual cycles of the respective dominant baroclinic modes are associated with characteristic basin-wide structures. Using an idealized linear reduced-gravity model to simulate the dynamics of individual baroclinic modes, it is shown that the observed circulation variability can be explained by resonant equatorial basin modes. Corollary simulations of the reduced-gravity model with varying basin geometry (i.e. square basin versus realistic coastlines) or forcing (i.e. spatially uniform versus spatially variable wind) show a structural robustness of the simulated basin modes. A main focus of this study is the seasonal variability of the Equatorial Undercurrent (EUC) as identified in recent observational studies. Main characteristics of the observed EUC including seasonal variability of transport, core depth, and maximum core velocity can be explained by the linear superposition of the dominant equatorial basin modes as obtained from the reduced-gravity model.

Description: In Ammassalik, in southeast Greenland, downslope winds can reach hurricane intensity and represent a hazard for the local population and environment. They advect cold air down the ice sheet and over the Irminger Sea, where they drive large ocean–atmosphere heat fluxes over an important ocean convection region. Earlier studies have found them to be associated with a strong katabatic acceleration over the steep coastal slopes, flow convergence inside the valley of Ammassalik, and—in one instance—mountain wave breaking. Yet, for the general occurrence of strong downslope wind events, the importance of mesoscale processes is largely unknown. Here, two wind events—one weak and one strong—are simulated with the atmospheric Weather Research and Forecasting (WRF) Model with different model and topography resolutions, ranging from 1.67 to 60 km. For both events, but especially for the strong one, it is found that lower resolutions underestimate the wind speed because they misrepresent the steepness of the topography and do not account for the underlying wave dynamics. If a 5-km model instead of a 60-km model resolution in Ammassalik is used, the flow associated with the strong wind event is faster by up to 20 m s−1. The effects extend far downstream over the Irminger Sea, resulting in a diverging spatial distribution and temporal evolution of the heat fluxes. Local differences in the heat fluxes amount to 20%, with potential implications for ocean convection.

Description: We perform eddy-resolving and high-vertical-resolution numerical simulations of the circulation in an idealized equatorial Atlantic Ocean in order to explore the formation of the deep equatorial circulation (DEC) in this basin. Unlike in previous studies, the deep equatorial intraseasonal variability (DEIV) that is believed to be the source of the DEC is generated internally by instabilities of the upper ocean currents. Two main simulations are discussed: Solution 1, configured with a rectangular basin and with wind forcing that is zonally and temporally uniform; and Solution 2, with realistic coastlines and with an annual cycle of wind forcing varying zonally. Somewhat surprisingly, Solution 1 produces the more realistic DEC: The large-vertical-scale currents (Equatorial Intermediate Currents or EICs) are found over a large zonal portion of the basin, and the small-vertical-scale equatorial currents (Equatorial Deep Jets or EDJs) form low-frequency, quasi-resonant, baroclinic equatorial basin modes with phase propagating mostly downward, consistent with observations. We demonstrate that both types of currents arise from the rectification of DEIV, consistent with previous theories. We also find that the EDJs contribute to maintaining the EICs, suggesting that the nonlinear energy transfer is more complex than previously thought. In Solution 2, the DEC is unrealistically weak and less spatially coherent than in the first simulation probably because of its weaker DEIV. Using intermediate solutions, we find that the main reason for this weaker DEIV is the use of realistic coastlines in Solution 2. It remains to be determined, what needs to be modified or included to obtain a realistic DEC in the more realistic configuration.

Description: A yet unexplained drift of (some) oxygen optodes during storage/transport and thus significant deviations from factory/laboratory calibrations have been a major handicap for autonomous oxygen observations. Optode drift appears to be systematic and is predominantly a slope effect due to reduced oxygen sensitivity. A small contribution comes from a reduced luminophore lifetime, which causes a small positive offset. A reliable in situ reference is essential to correct such a drift. Traditionally, this called for a ship-based reference cast, which poses some challenges for opportunistic float deployments. This study presents an easily implemented alternative using near-surface/in-air measurements of an Aanderaa optode on a 10-cm stalk and compares it to the more traditional approaches (factory, laboratory, and in situ deployment calibration). In-air samples show a systematic bias depending on the water saturation, which is likely caused by occasional submersions of the standard-height stalk optode. Linear regression of measured in-air supersaturation against in-water supersaturation (using ancillary meteorological data to define the saturation level) robustly removes this bias and thus provides a precise (0.2%) and accurate (1%) in situ correction that is available throughout the entire instrument’s lifetime.

Description: Predictability on seasonal time scales over the North Atlantic–Europe region is assessed using a seasonal prediction system based on an initialized version of the Max Planck Institute Earth System Model (MPI-ESM). For this region, two of the dominant predictors on seasonal time scales are El Niño–Southern Oscillation (ENSO) and sudden stratospheric warming (SSW) events. Multiple studies have shown a potential for improved North Atlantic predictability for either predictor. Their respective influences are however difficult to disentangle, since the stratosphere is itself impacted by ENSO. Both El Niño and SSW events correspond to a negative signature of the North Atlantic Oscillation (NAO), which has a major influence on European weather. This study explores the impact on Europe by separating the stratospheric pathway of the El Niño teleconnection. In the seasonal prediction system, the evolution of El Niño events is well captured for lead times of up to 6 months, and stratospheric variability is reproduced with a realistic frequency of SSW events. The model reproduces the El Niño teleconnection through the stratosphere, involving a deepened Aleutian low connected to a warm anomaly in the northern winter stratosphere. The stratospheric anomaly signal then propagates downward into the troposphere through the winter season. Predictability of 500-hPa geopotential height over Europe at lead times of up to 4 months is shown to be increased only for El Niño events that exhibit SSW events, and it is shown that the characteristic negative NAO signal is only obtained for winters also containing major SSW events for both the model and the reanalysis data.

Description: Classical theory concerning theEliassen–Palmrelation is extended in this study to allowfor a unified treatment of midlatitude inertia–gravity waves (MIGWs), midlatitude Rossby waves (MRWs), and equatorial waves (EQWs). A conservation equation for what the authors call the impulse-bolus (IB) pseudomomentum is useful, because it is applicable to ageostrophic waves, and the associated three-dimensional flux is parallel to the direction of the group velocity of MRWs. The equation has previously been derived in an isentropic coordinate system or a shallow-water model. The authors make an explicit comparison of prognostic equations for the IB pseudomomentum vector and the classical energy-based (CE) pseudomomentum vector, assuming inviscid linear waves in a sufficiently weak mean flow, to provide a basis for the former quantity to be used in an Eulerian time-mean (EM) framework. The authors investigate what makes the three-dimensional fluxes in the IB and CE pseudomomentum equations look in different directions. It is found that the two fluxes are linked by a gauge transformation, previously unmentioned, associated with a divergence-form wave-induced pressure L. The quantity L vanishes for MIGWs and becomes nonzero for MRWs and EQWs, and it may be estimated using the virial theorem. Concerning the effect of waves on the mean flow, L represents an additional effect in the pressure gradient term of both (the three-dimensional versions of) the transformed EM momentum equations and the merged form of the EMmomentumequations, the latter of which is associated with the nonacceleration theorem.

Description: [in “State of the Climate in 2014” : Special Supplement to the Bulletin of the American Meteorological Society Vol. 96, No. 7, July 2015]

Description: Atmospheric deposition contributes potentially significant amounts of the nutrients iron, nitrogen and phosphorus (via mineral dust and anthropogenic aerosols) to the oligotrophic tropical North Atlantic Ocean. Transport pathways, deposition processes and source strengths contributing to this atmospheric flux are all highly variable in space and time. Atmospheric sampling was conducted during 28 research cruises through the Eastern Tropical North Atlantic (ETNA) over a 12 year period and a substantial dataset of measured concentrations of nutrients and trace metals in aerosol and rainfall over the region was acquired. This database was used to quantify (on a spatial- and seasonal-basis) the atmospheric input of ammonium, nitrate, soluble phosphorus and soluble and total iron, aluminium and manganese to the ETNA. The magnitude of atmospheric input varies strongly across the region, with high rainfall rates associated with the Inter-tropical Convergence Zone contributing to high wet deposition fluxes in the south, particularly for soluble species. Dry deposition fluxes of species associated with mineral dust exhibited strong seasonality, with highest fluxes associated with winter-time low-level transport of Saharan dust. Overall (wet plus dry) atmospheric inputs of soluble and total trace metals were used to estimate their soluble fractions. These also varied with season and were generally lower in the dry north than in the wet south. The ratio of ammonium plus nitrate to soluble iron in deposition to the ETNA was lower than the N:Fe requirement for algal growth in all cases, indicating the importance of the atmosphere as a source of excess iron.

Description: We investigated the effect of hydrostatic pressure of up to 6000 dbar on Aanderaa and Sea-Bird oxygen optodes both in the laboratory and in the field. The overall pressure response is a reduction in the O2 reading by 3 – 4 % per 1000 dbar which is closely linear with pressure and increases with temperature. Closer inspection reveals two superimposed processes with opposite effect: an O2-independent pressure response on the luminophore which increases optode O2 readings and an O2-dependent change in luminescence quenching which decreases optode O2 readings. The latter process dominates and is mainly due to a shift in the equilibrium between sensing membrane and sea water under elevated pressures. If only the dominant O2-dependent process is considered, Aanderaa and Sea-Bird optodes differ in their pressure response. Compensation of the O2-independent process, however, yields a uniform O2 dependence for Aanderaa optodes with standard foil and fast-response foil as well as Sea-Bird optodes. A new scheme to calculate optode O2 from raw data is proposed to account for the two processes. The overall uncertainty of the optode pressure correction amounts to 0.3 % per 1000 dbar, mainly due to variability between sensors.

Description: The upper tropical Atlantic Ocean has markedly warmed since the 1960s. It has been shown that this warming was not due to local heat fluxes, and that the trade winds that drive the coastal and equatorial upwelling have intensified rather than weakened. Remote forcing might thus have played an important role. Here model experiments are used to investigate the contribution from an increased inflow of warm Indian Ocean water through Agulhas leakage. A high-resolution hindcast experiment with interannually varying forcing for the time period 1948 to 2007, in which Agulhas leakage increases by about 45% from the 1960s to the early 2000s, reproduces the observed warming trend. To tease out the role of Agulhas leakage, a sensitivity experiment designed to only increase Agulhas leakage is used. Compared to a control simulation it shows a pronounced warming in the upper tropical Atlantic Ocean. A Lagrangian trajectory analysis confirms that a significant portion of Agulhas leakage water reaches the upper 300m of the tropical Atlantic Ocean within two decades, and that the tropical Atlantic warming in the sensitivity experiment is mainly due to water of Agulhas origin. Therefore, it is suggested that the increased trade winds since the 1960s favor upwelling of warmer subsurface waters, which in parts originate from the Agulhas, leading to higher SSTs in the tropics

Description: Variations in the global tropospheric zonal mean zonal wind ([U]) during boreal winter are investigated using Rotated Empirical Orthogonal Functions applied to monthly means. The first two modes correspond to the Northern and Southern Annular Mode and modes 3 and 4 represent variability in the tropics. One is related to El Niño Southern Oscillation and the other has variability that is highly correlated with the time series of [U] at 150 hPa between 5°N and 5°S ([U150]E) and is related to activity of the Madden-Julian Oscillation. The extratropical response to [U150]E is investigated using linear regressions of 500 hPa geopotential height onto the [U150]E time series. We make use of reanalysis data and of the ensemble mean output from a relaxation experiment using the European Center for Medium Range Weather Forecasts model in which the tropical atmosphere is relaxed towards reanalysis data. The regression analysis reveals that a shift of the Aleutian low and a wave train across the North Atlantic are associated with [U150]E. We find that the subtropical waveguides and the link between the North Pacific and North Atlantic are stronger during the easterly phase of [U150]E. The wave train over the North Atlantic is associated with Rossby wave sources over the subtropical North Pacific and North America. Finally, we show that a linear combination of both [U150]E and the Quasi Biennial Oscillation in the lower stratosphere can explain the circulation anomalies of the anomalously cold European winter of 1962/63 when both were in an extreme easterly phase.

Description: To investigate the influence of atmospheric model resolution on the representation of daily precipitation extremes, ensemble simulations with the atmospheric general circulation model ECHAM5 at different horizontal (T213 to T31) and vertical (L31 to L19) resolutions and forced with observed sea surface temperatures and sea ice concentrations have been carried out for 01/1982 - 09/2010. All results have been compared with the highest resolution, which has been validated against observations. Resolution affects both the representation of physical processes and the averaging of precipitation across grid boxes. The latter, in particular, smoothes out localized extreme events. These effects have been disentangled by averaging precipitation simulated at the highest resolution to the corresponding coarser grid. Extremes are represented by seasonal maxima, modeled by the generalized extreme value distribution. Effects of averaging and representation of physical processes vary with region and season. In the tropical summer hemisphere, extreme precipitation is reduced by up to 30% due to the averaging effect, and a further 65% owing to a coarser representation of physical processes. Towards mid- to high latitudes, the latter effect reduces to 20%; in the winter hemisphere it vanishes towards the poles. A strong drop is found between T106 and T63 in the convection dominated tropics. At the lowest resolution, northern hemisphere winter precipitation extremes, mainly caused by large scale weather systems, are in general represented reasonably well. Coarser vertical resolution causes an equatorward shift of maximum extreme precipitation in the tropics. The impact of vertical resolution on mean precipitation is less pronounced; for horizontal resolution it is negligible.

Description: Many aspects of the coupling between the ocean and atmosphere at the mesoscale (on the order of 20–100 km) remain unknown. While recent observations from the Southern Ocean revealed that circular fronts associated with oceanic mesoscale eddies leave a distinct imprint on the overlying wind, cloud coverage, and rain, the mechanisms responsible for explaining these atmospheric changes are not well established. Here the atmospheric response above mesoscale ocean eddies is investigated utilizing a newly developed coupled atmosphere–ocean regional model [Consortium for Small-Scale Modeling–Regional Ocean Modelling System (COSMO-ROMS)] configured at a horizontal resolution of ~10 km for the South Atlantic and run for a 3-month period during austral winter of 2004. The model-simulated changes in surface wind, cloud fraction, and rain above the oceanic eddies are very consistent with the relationships inferred from satellite observations for the same region and time. From diagnosing the model’s momentum balance, it is shown that the atmospheric imprint of the oceanic eddies are driven by the modification of vertical mixing in the atmospheric boundary layer, rather than secondary flows driven by horizontal pressure gradients. This is largely due to the very limited ability of the atmosphere to adjust its temperature over the time scale it takes for an air parcel to pass over these mesoscale oceanic features. This results in locally enhanced vertical gradients between the ocean surface and the overlying air and thus a rapid change in turbulent mixing in the atmospheric boundary layer and an associated change in the vertical momentum flux.

Description: This study examines the role of processes transporting tracers across the Polar Front (PF) in the depth interval between the surface and major topographic sills, which this study refers to as the “PF core.” A preindustrial control simulation of an eddying climate model coupled to a biogeochemical model [GFDL Climate Model, version 2.6 (CM2.6)– simplified version of the Biogeochemistry with Light Iron Nutrients and Gas (miniBLING) 0.1° ocean model] is used to investigate the transport of heat, carbon, oxygen, and phosphate across the PF core, with a particular focus on the role of mesoscale eddies. The authors find that the total transport across the PF core results from a ubiquitous Ekman transport that drives the upwelled tracers to the north and a localized opposing eddy transport that induces tracer leakages to the south at major topographic obstacles. In the Ekman layer, the southward eddy transport only partially compensates the northward Ekman transport, while below the Ekman layer, the southward eddy transport dominates the total transport but remains much smaller in magnitude than the near-surface northward transport. Most of the southward branch of the total transport is achieved below the PF core, mainly through geostrophic currents. This study finds that the eddy-diffusive transport reinforces the southward eddy-advective transport for carbon and heat, and opposes it for oxygen and phosphate. Eddy-advective transport is likely to be the leading-order component of eddy-induced transport for all four tracers. However, eddy-diffusive transport may provide a significant contribution to the southward eddy heat transport due to strong along-isopycnal temperature gradients.

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: Eutrophication combined with climate change has caused ephemeral filamentous macroalgae to increase and drifts of seaweed cover large areas of some Baltic Sea sites during summer. In ongoing projects, these mass occurrences of drifting filamentous macroalgae are being harvested to mitigate eutrophication, with preliminary results indicating considerable nutrient reduction potential. In the present study, an energy assessment was made of biogas production from the retrieved biomass for a Baltic Sea pilot case. Use of different indicators revealed a positive energy balance. The energy requirements corresponded to about 30%–40% of the energy content in the end products. The net energy gain was 530–800 MJ primary energy per ton wet weight of algae for small-scale and large-scale scenarios, where 6 000 and 13 000 tonnes dwt were harvested, respectively. However, the exergy efficiency differed from the energy efficiency, emphasising the importance of taking energy quality into consideration when evaluating energy systems. An uncertainty analysis indicated parametric uncertainty of about 25%–40%, which we consider to be acceptable given the generally high sensitivity of the indicators to changes in input data, allocation method, and system design. Overall, our evaluation indicated that biogas production may be a viable handling strategy for retrieved biomass, while harvesting other types of macroalgae than red filamentous species considered here may render a better energy balance due to higher methane yields.

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.