An Ensemble Kalman Filter multi-tracer assimilation: Determining uncertain ocean model parameters for improved climate-carbon cycle projections

An Ensemble Kalman Filter is applied to assimilate observed tracer fields in various combinations in the Bern3D ocean model. Each tracer combination yields a set of optimal transport parameter values that are used in projections with prescribed CO2 stabilization pathways. The assimilation of temperature and salinity fields yields a too vigorous ventilation of the thermocline and the deep ocean, whereas the inclusion of CFC-11 and radiocarbon improves the representation of physical and biogeochemical tracers and of ventilation time scales. Projected peak uptake rates and cumulative uptake of CO2 by the ocean are around 20% lower for the parameters determined with CFC-11 and radiocarbon as additional target compared to those with salinity and temperature only. Higher surface temperature changes are simulated in the Greenland–Norwegian–Iceland Sea and in the Southern Ocean when CFC-11 is included in the Ensemble Kalman model tuning. These findings highlights the importance of ocean transport calibration for the design of near-term and long-term CO2 emission mitigation strategies and for climate projections.

Amplified inception of European Little Ice Age by sea ice-ocean-atmosphere feedbacks

The inception of the Little Ice Age (~1400–1700 AD) is believed to have been driven by an interplay of external forcing and climate system internal variability. While the hemispheric signal seems to have been dominated by solar irradiance and volcanic eruptions, the understanding of mechanisms shaping the climate on a continental scale is less robust. In an ensemble of transient model simulations and a new type of sensitivity experiments with artificial sea ice growth, the authors identify a sea ice–ocean–atmosphere feedback mechanism that amplifies the Little Ice Age cooling in the North Atlantic–European region and produces the temperature pattern suggested by paleoclimatic reconstructions. Initiated by increasing negative forcing, the Arctic sea ice substantially expands at the beginning of the Little Ice Age. The excess of sea ice is exported to the subpolar North Atlantic, where it melts, thereby weakening convection of the ocean. Consequently, northward ocean heat transport is reduced, reinforcing the expansion of the sea ice and the cooling of the Northern Hemisphere. In the Nordic Seas, sea surface height anomalies cause the oceanic recirculation to strengthen at the expense of the warm Barents Sea inflow, thereby further reinforcing sea ice growth. The absent ocean–atmosphere heat flux in the Barents Sea results in an amplified cooling over Northern Europe. The positive nature of this feedback mechanism enables sea ice to remain in an expanded state for decades up to a century, favoring sustained cold periods over Europe such as the Little Ice Age. Support for the feedback mechanism comes from recent proxy reconstructions around the Nordic Seas.

Upper ocean vertical supply: A neglected primary factor controlling the distribution of neodymium concentrations of open ocean surface waters?

Neodymium (Nd) isotopes are an important geochemical tool to trace the present and past water mass mixing as well as continental inputs. The distribution of Nd concentrations in open ocean surface waters (0�100 m) is generally assumed to be controlled by lateral mixing of Nd from coastal surface currents and by removal through reversible particle scavenging. However, using 228Ra activity as an indicator of coastal water mass influence, surface water Nd concentration data available on key oceanic transects as a whole do not support the above scenario. From a global compilation of available data, we find that more stratified regions are generally associated with low surface Nd concentrations. This implies that upper ocean vertical supply may be an as yet neglected primary factor in determining the basin-scale variations of surface water Nd concentrations. Similar to the mechanism of nutrients supply, it is likely that stratification inhibits vertical supply of Nd from the subsurface thermocline waters and thus the magnitude of Nd flux to the surface layer. Consistently, the estimated required input flux of Nd to the surface layer to maintain the observed concentrations could be nearly two orders of magnitudes larger than riverine/dust flux, and also larger than the model-based estimation on shelf-derived coastal flux. In addition, preliminary results from modeling experiments reveal that the input from shallow boundary sources, riverine input, and release from dust are actually not the primary factors controlling Nd concentrations most notably in the Pacific and Southern Ocean surface waters.

Taking action against ocean acidification: A review of management and policy options

Ocean acidification has emerged over the last two decades as one of the largest threats to marine organisms and ecosystems. However, most research efforts on ocean acidification have so far neglected management and related policy issues to focus instead on understanding its ecological and biogeochemical implications. This shortfall is addressed here with a systematic, international and critical review of management and policy options. In particular, we investigate the assumption that fighting acidification is mainly, but not only, about reducing CO2 emissions, and explore the leeway that this emerging problem may open in old environmental issues. We review nine types of management responses, initially grouped under four categories: preventing ocean acidification; strengthening ecosystem resilience; adapting human activities; and repairing damages. Connecting and comparing options leads to classifying them, in a qualitative way, according to their potential and feasibility. While reducing CO2 emissions is confirmed as the key action that must be taken against acidification, some of the other options appear to have the potential to buy time, e.g. by relieving the pressure of other stressors, and help marine life face unavoidable acidification. Although the existing legal basis to take action shows few gaps, policy challenges are significant: tackling them will mean succeeding in various areas of environmental management where we failed to a large extent so far.

Mid- to late Holocene Indian Ocean Monsoon variability recorded in four speleothems from Socotra Island, Yemen

Four stalagmites covering the last 7.0 ka were sampled on Socotra, an island in the northern Indian Ocean to investigate the evolution of the northeast Indian Ocean Monsoon (IOM) since the mid Holocene. On Socotra, rain is delivered at the start of the southwest IOM in May–June and at the start of the northeast IOM from September to December. The Haggeher Mountains act as a barrier forcing precipitation brought by the northeast winds to fall preferentially on the eastern side of the island, where the studied caves are located. δ18O and δ13C and Mg/Ca and Sr/Ca signals in the stalagmites reflect precipitation amounts brought by the northeast winds. For stalagmite STM6, this amount effect is amplified by kinetic effects during calcite deposition. Combined interpretation of the stalagmites' signals suggest a weakening of the northeast precipitation between 6.0 and 3.8 ka. After 3.8 ka precipitation intensities remain constant with two superimposed drier periods, between 0 and 0.6 ka and from 2.2 to 3.8 ka. No link can be established with Greenland ice cores and with the summer IOM variability. In contrast to the stable northeast rainy season suggested by the records in this study, speleothem records from western Socotra indicate a wettening of the southwest rainy season on Socotra after 4.4 ka. The local wettening of western Socotra could relate to a more southerly path (more over the Indian Ocean) taken by the southwest winds. Stalagmite STM5, sampled at the fringe between both rain areas displays intermediate δ18O values. After 6.2 ka, similar precipitation changes are seen between eastern Socotra and northern Oman indicating that both regions are affected similarly by the monsoon. Different palaeoclimatologic records from the Arabian Peninsula currently located outside the ITCZ migration pathway display an abrupt drying around 6 ka due to their disconnection from the southwest rain influence. Records that are nowadays still receiving rain by the southwest winds, suggest a more gradual drying reflecting the weakening of the southwest monsoon.

Tracing the origin of Arctic driftwood

Arctic environments, where surface temperatures increase and sea ice cover and permafrost depth decrease, are very sensitive to even slight climatic variations. Placing recent environmental change of the high-northern latitudes in a long-term context is, however, complicated by too short meteorological observations and too few proxy records. Driftwood may represent a unique cross-disciplinary archive at the interface of marine and terrestrial processes. Here, we introduce 1445 driftwood remains from coastal East Greenland and Svalbard. Macroscopy and microscopy were applied for wood anatomical classification; a multi-species subset was used for detecting fungi; and information on boreal vegetation patterns, circumpolar river systems, and ocean current dynamics was reviewed and evaluated. Four conifer (Pinus, Larix, Picea, and Abies) and three deciduous (Populus, Salix, and Betula) genera were differentiated. Species-specific identification also separated Pinus sylvestris and Pinus sibirica, which account for ~40% of all driftwood and predominantly originate from western and central Siberia. Larch and spruce from Siberia or North America represents ~26% and ~18% of all materials, respectively. Fungal colonization caused different levels of driftwood staining and/or decay. Our results demonstrate the importance of combining wood anatomical knowledge with insight on boreal forest composition for successfully tracing the origin of Arctic driftwood. To ultimately reconstruct spatiotemporal variations in ocean currents, and to better quantify postglacial uplift rates, we recommend consideration of dendrochronologically dated material from many more circumpolar sites.

Global ocean storage of anthropogenic carbon

The global ocean is a significant sink for anthropogenic carbon (Cant), absorbing roughly a third of human CO2 emitted over the industrial period. Robust estimates of the magnitude and variability of the storage and distribution of Cant in the ocean are therefore important for understanding the human impact on climate. In this synthesis we review observational and model-based estimates of the storage and transport of Cant in the ocean. We pay particular attention to the uncertainties and potential biases inherent in different inference schemes. On a global scale, three data-based estimates of the distribution and inventory of Cant are now available. While the inventories are found to agree within their uncertainty, there are considerable differences in the spatial distribution. We also present a review of the progress made in the application of inverse and data assimilation techniques which combine ocean interior estimates of Cant with numerical ocean circulation models. Such methods are especially useful for estimating the air–sea flux and interior transport of Cant, quantities that are otherwise difficult to observe directly. However, the results are found to be highly dependent on modeled circulation, with the spread due to different ocean models at least as large as that from the different observational methods used to estimate Cant. Our review also highlights the importance of repeat measurements of hydrographic and biogeochemical parameters to estimate the storage of Cant on decadal timescales in the presence of the variability in circulation that is neglected by other approaches. Data-based Cant estimates provide important constraints on forward ocean models, which exhibit both broad similarities and regional errors relative to the observational fields. A compilation of inventories of Cant gives us a "best" estimate of the global ocean inventory of anthropogenic carbon in 2010 of 155 ± 31 PgC (±20% uncertainty). This estimate includes a broad range of values, suggesting that a combination of approaches is necessary in order to achieve a robust quantification of the ocean sink of anthropogenic CO2.

Warm Conveyor Belts in the ERA-Interim Dataset (1979–2010): Part I: Climatology and Potential Vorticity Evolution

A global climatology of warm conveyor belts (WCBs) is presented for the years 1979–2010, based on trajectories calculated with Interim ECMWF Re-Analysis (ERA-Interim) data. WCB trajectories are identified as strongly ascending air parcels (600 hPa in 2 days) near extratropical cyclones. Corroborating earlier studies, WCBs are more frequent during winter than summer and they ascend preferentially in the western ocean basins between 25° and 50° latitude. Before ascending, WCB trajectories typically approach from the subtropics in summer and from more midlatitude regions in winter. Considering humidity, cloud water, and potential temperature along WCBs confirms that they experience strong condensation and integrated latent heating during the ascent (typically >20 K). Liquid and ice water contents along WCBs peak at about 700 and 550 hPa, respectively. The mean potential vorticity (PV) evolution shows typical tropospheric values near 900 hPa, followed by an increase to almost 1 potential vorticity unit (PVU) at 700 hPa, and a decrease to less than 0.5 PVU at 300 hPa. These low PV values in the upper troposphere constitute significant negative anomalies with amplitudes of 1–3 PVU, which can strongly influence the downstream flow. Considering the low-level diabatic PV production, (i) WCBs starting at low latitudes (<40°) are unlikely to attain high PV (due to weak planetary vorticity) although they exhibit the strongest latent heating, and (ii) for those ascending at higher latitudes, a strong vertical heating gradient and high absolute vorticity are both important. This study therefore provides climatological insight into the cloud diabatic formation of significant positive and negative PV anomalies in the extratropical lower and upper troposphere, respectively.

Inter-hemispheric temperature variability over the past millennium

The Earth’s climate system is driven by a complex interplay of internal chaotic dynamics and natural and anthropogenic external forcing. Recent instrumental data have shown a remarkable degree of asynchronicity between Northern Hemisphere and Southern Hemisphere temperature fluctuations, thereby questioning the relative importance of internal versus external drivers of past as well as future climate variability1, 2, 3. However, large-scale temperature reconstructions for the past millennium have focused on the Northern Hemisphere4, 5, limiting empirical assessments of inter-hemispheric variability on multi-decadal to centennial timescales. Here, we introduce a new millennial ensemble reconstruction of annually resolved temperature variations for the Southern Hemisphere based on an unprecedented network of terrestrial and oceanic palaeoclimate proxy records. In conjunction with an independent Northern Hemisphere temperature reconstruction ensemble5, this record reveals an extended cold period (1594–1677) in both hemispheres but no globally coherent warm phase during the pre-industrial (1000–1850) era. The current (post-1974) warm phase is the only period of the past millennium where both hemispheres are likely to have experienced contemporaneous warm extremes. Our analysis of inter-hemispheric temperature variability in an ensemble of climate model simulations for the past millennium suggests that models tend to overemphasize Northern Hemisphere–Southern Hemisphere synchronicity by underestimating the role of internal ocean–atmosphere dynamics, particularly in the ocean-dominated Southern Hemisphere. Our results imply that climate system predictability on decadal to century timescales may be lower than expected based on assessments of external climate forcing and Northern Hemisphere temperature variations5, 6 alone.

A 200 Year Sulfate Record from Sixteen Antarctic Ice Cores and Associations With Southern Ocean Sea-Ice Extent

Chemistry data from 16, 50-115 m deep, sub-annually dated ice cores are used to investigate spatial and temporal concentration variability of sea-salt (ss) SO42- and excess (xs) SO42- over West Antarctica and the South Pole for the last 200 years. Low-elevation ice-core sites in western West Antarctica contain higher concentrations Of SO42- as a result of cyclogenesis over the Ross Ice Shelf and proximity to the Ross Sea Polynya. Linear correlation analysis of 15 West Antarctic ice-core SO42- time series demonstrates that at several sites concentrations Of ssSO(4)(2-) are higher when sea-ice (SIE) extent is greater, and the inverse for XSS04. Concentrations Of XSS04 from the South Pole site (East Antarctica) are associated with SIE from the Weddell region, and West Antarctic XSSO42- concentrations are associated with SIE from the Bellingshausen-Amundsen-Ross region. The only notable rise of the last 200 years in xsSO(4)(2-), around 1940, is not related to SIE fluctuations and is most likely a result of increased xsSO(4)(2-) production in the mid-low latitudes and/or an increase in transport efficiency from the mid-low latitudes to central West Antarctica. These high-resolution records show that the source types and source areas Of ssSO(4)(2-) and xsSO(4)(2-) delivered to eastern and western West Antarctica and the South Pole differ from site to site but can best be resolved using records from spatial ice-core arrays such as the International Trans-Antarctic Scientific Expedition (ITASE).

State of the Antarctic and Southern Ocean Climate System

This paper reviews developments in our understanding of the state of the Antarctic and Southern Ocean climate and its relation to the global climate system over the last few millennia. Climate over this and earlier periods has not been stable, as evidenced by the occurrence of abrupt changes in atmospheric circulation and temperature recorded in Antarctic ice core proxies for past climate. Two of the most prominent abrupt climate change events are characterized by intensification of the circumpolar westerlies (also known as the Southern Annular Mode) between similar to 6000 and 5000 years ago and since 1200-1000 years ago. Following the last of these is a period of major trans-Antarctic reorganization of atmospheric circulation and temperature between A. D. 1700 and 1850. The two earlier Antarctic abrupt climate change events appear linked to but predate by several centuries even more abrupt climate change in the North Atlantic, and the end of the more recent event is coincident with reorganization of atmospheric circulation in the North Pacific. Improved understanding of such events and of the associations between abrupt climate change events recorded in both hemispheres is critical to predicting the impact and timing of future abrupt climate change events potentially forced by anthropogenic changes in greenhouse gases and aerosols. Special attention is given to the climate of the past 200 years, which was recorded by a network of recently available shallow firn cores, and to that of the past 50 years, which was monitored by the continuous instrumental record. Significant regional climate changes have taken place in the Antarctic during the past 50 years. Atmospheric temperatures have increased markedly over the Antarctic Peninsula, linked to nearby ocean warming and intensification of the circumpolar westerlies. Glaciers are retreating on the peninsula, in Patagonia, on the sub-Antarctic islands, and in West Antarctica adjacent to the peninsula. The penetration of marine air masses has become more pronounced over parts of West Antarctica. Above the surface, the Antarctic troposphere has warmed during winter while the stratosphere has cooled year-round. The upper kilometer of the circumpolar Southern Ocean has warmed, Antarctic Bottom Water across a wide sector off East Antarctica has freshened, and the densest bottom water in the Weddell Sea has warmed. In contrast to these regional climate changes, over most of Antarctica, near-surface temperature and snowfall have not increased significantly during at least the past 50 years, and proxy data suggest that the atmospheric circulation over the interior has remained in a similar state for at least the past 200 years. Furthermore, the total sea ice cover around Antarctica has exhibited no significant overall change since reliable satellite monitoring began in the late 1970s, despite large but compensating regional changes. The inhomogeneity of Antarctic climate in space and time implies that recent Antarctic climate changes are due on the one hand to a combination of strong multidecadal variability and anthropogenic effects and, as demonstrated by the paleoclimate record, on the other hand to multidecadal to millennial scale and longer natural variability forced through changes in orbital insolation, greenhouse gases, solar variability, ice dynamics, and aerosols. Model projections suggest that over the 21st century the Antarctic interior will warm by 3.4 degrees +/- 1 degrees C, and sea ice extent will decrease by similar to 30%. Ice sheet models are not yet adequate enough to answer pressing questins about the effect of projected warming on mass balance and sea level. Considering the potentially major impacts of a warming climate on Antarctica, vigorous efforts are needed to better understand all aspects of the highly coupled Antarctic climate system as well as its influence on the Earth's climate and oceans.

Sea Level Pressure Variability Over the Southern Indian Ocean Inferred from a Glaciochemical Record in Princess Elizabeth Land, East Antarctica

A 250-year, high-resolution, multivariate ice core record from LGB65 (70degrees50'07"S, 77degrees04'29"E; 1850 m asl), Princess Elizabeth Land (PEL), is used to investigate sea level pressure (SLP) variability over the southern Indian Ocean (SIO). Empirical orthogonal function (EOF) analysis reveals that the first EOF (EOF1) of the glaciochemical record from LGB65 represents most of the variability in sea salt throughout the 250-year record. EOF1 is negatively correlated (95% confidence level and higher) to instrumental mean sea level pressure (MSLP) at Kerguelen and New Amsterdam islands, SIO. On the basis of comparison with NCEP/NCAR reanalysis, strong correlations were found between sea-salt variations and a quasi-stationary low that lies to the north of Prydz Bay, SIO. Comparison with a 250-year-long summer transpolar index (STPI) inferred from sub-Antarctic tree ring records reveals strong coherency. Decadal-scale SLP variability over SIO suggests shifting of the polar vortex. Prominent decadal-scale deepening of the southern Indian Ocean low (SIOL) exists circa 1790, 1810, 1835, 1860, 1880, 1900, and 1940 A. D., continuously after the 1970s, and prominent weakening circa 1750, 1795, 1825, 1850, 1870, 1890, 1910, and 1955 A. D. The LGB65 sea-salt record is characterized by significant decadal-scale variability with a strong similar to21-year periodic structure (99.9% confidence level). The relationship between LGB65 sea salt and solar irradiance changes shows that this periodicity is possibly the solar Hale cycle ( 22 years).

Impact of solar versus volcanic activity variations on tropospheric temperatures and precipitation during the Dalton Minimum

The aim of this work is to elucidate the impact of changes in solar irradiance and energetic particles versus volcanic eruptions on tropospheric global climate during the Dalton Minimum (DM, AD 1780–1840). Separate variations in the (i) solar irradiance in the UV-C with wavelengths λ < 250 nm, (ii) irradiance at wavelengths λ > 250 nm, (iii) in energetic particle spectrum, and (iv) volcanic aerosol forcing were analyzed separately, and (v) in combination, by means of small ensemble calculations using a coupled atmosphere–ocean chemistry–climate model. Global and hemispheric mean surface temperatures show a significant dependence on solar irradiance at λ > 250 nm. Also, powerful volcanic eruptions in 1809, 1815, 1831 and 1835 significantly decreased global mean temperature by up to 0.5 K for 2–3 years after the eruption. However, while the volcanic effect is clearly discernible in the Southern Hemispheric mean temperature, it is less significant in the Northern Hemisphere, partly because the two largest volcanic eruptions occurred in the SH tropics and during seasons when the aerosols were mainly transported southward, partly because of the higher northern internal variability. In the simulation including all forcings, temperatures are in reasonable agreement with the tree ring-based temperature anomalies of the Northern Hemisphere. Interestingly, the model suggests that solar irradiance changes at λ < 250 nm and in energetic particle spectra have only an insignificant impact on the climate during the Dalton Minimum. This downscales the importance of top–down processes (stemming from changes at λ < 250 nm) relative to bottom–up processes (from λ > 250 nm). Reduction of irradiance at λ > 250 nm leads to a significant (up to 2%) decrease in the ocean heat content (OHC) between 0 and 300 m in depth, whereas the changes in irradiance at λ < 250 nm or in energetic particles have virtually no effect. Also, volcanic aerosol yields a very strong response, reducing the OHC of the upper ocean by up to 1.5%. In the simulation with all forcings, the OHC of the uppermost levels recovers after 8–15 years after volcanic eruption, while the solar signal and the different volcanic eruptions dominate the OHC changes in the deeper ocean and prevent its recovery during the DM. Finally, the simulations suggest that the volcanic eruptions during the DM had a significant impact on the precipitation patterns caused by a widening of the Hadley cell and a shift in the intertropical convergence zone.

Air-Sea Interactions During the Passage of a Winter Storm Over the Gulf Stream: A Three-Dimensional Coupled Atmosphere-Ocean Model Study

A three-dimensional, regional coupled atmosphere-ocean model with full physics is developed to study air-sea interactions during winter storms off the U. S. east coast. Because of the scarcity of open ocean observations, models such as this offer valuable opportunities to investigate how oceanic forcing drives atmospheric circulation and vice versa. The study presented here considers conditions of strong atmospheric forcing (high wind speeds) and strong oceanic forcing (significant sea surface temperature (SST) gradients). A simulated atmospheric cyclone evolves in a manner consistent with Eta reanalysis, and the simulated air-sea heat and momentum exchanges strongly affect the circulations in both the atmosphere and the ocean. For the simulated cyclone of 19-20 January 1998, maximum ocean-to-atmosphere heat fluxes first appear over the Gulf Stream in the South Atlantic Bight, and this results in rapid deepening of the cyclone off the Carolina coast. As the cyclone moves eastward, the heat flux maximum shifts into the region near Cape Hatteras and later northeast of Hatteras, where it enhances the wind locally. The oceanic response to the atmospheric forcing is closely related to the wind direction. Southerly and southwesterly winds tend to strengthen surface currents in the Gulf Stream, whereas northeasterly winds weaken the surface currents in the Gulf Stream and generate southwestward flows on the shelf. The oceanic feedback to the atmosphere moderates the cyclone strength. Compared with a simulation in which the oceanic model always passes the initial SST to the atmospheric model, the coupled simulation in which the oceanic model passes the evolving SST to the atmospheric model produces higher ocean-to-atmosphere heat flux near Gulf Stream meander troughs. This is due to wind-driven lateral shifts of the stream, which in turn enhance the local northeasterly winds. Away from the Gulf Stream the coupled simulation produces surface winds that are 5 similar to 10% weaker. Differences in the surface ocean currents between these two experiments are significant on the shelf and in the open ocean.

Regulation of Phytoplankton Carbon to Chlorophyll Ratio By Light, Nutrients and Temperature in the Equatorial Pacific Ocean: A Basin-Scale Model

The complex effects of light, nutrients and temperature lead to a variable carbon to chlorophyll (C:Chl) ratio in phytoplankton cells. Using field data collected in the Equatorial Pacific, we derived a new dynamic model with a non-steady C:Chl ratio as a function of irradiance, nitrate, iron, and temperature. The dynamic model is implemented into a basin-scale ocean circulation-biogeochemistry model and tested in the Equatorial Pacific Ocean. The model reproduces well the general features of phytoplankton dynamics in this region. For instance, the simulated deep chlorophyll maximum (DCM) is much deeper in the western warm pool (similar to 100 m) than in the Eastern Equatorial Pacific (similar to 50 m). The model also shows the ability to reproduce chlorophyll, including not only the zonal, meridional and vertical variations, but also the interannual variability. This modeling study demonstrates that combination of nitrate and iron regulates the spatial and temporal variations in the phytoplankton C:Chl ratio in the Equatorial Pacific. Sensitivity simulations suggest that nitrate is mainly responsible for the high C:Chl ratio in the western warm pool while iron is responsible for the frontal features in the C:Chl ratio between the warm pool and the upwelling region. In addition, iron plays a dominant role in regulating the spatial and temporal variations of the C:Chl ratio in the Central and Eastern Equatorial Pacific. While temperature has a relatively small effect on the C:Chl ratio, light is primarily responsible for the vertical decrease of phytoplankton C:Chl ratio in the euphotic zone.