An Online Ensemble Coupled Data Assimilation Capability for the Community Earth System Model: System Design and Evaluation

The Community Earth System Model (CESM) developed at the National Center of Atmospheric Research (NCAR) has been used worldwide for climate studies. This study extends the efforts of CESM development to include an online (i.e., in-core) ensemble coupled data assimilation system (CESM-ECDA) to enhance CESM’s capability for climate predictability studies and prediction applications. The CESM-ECDA system consists of an online atmospheric data assimilation (ADA) component implemented to both the finite-volume and spectral-element dynamical cores, and an online oceanic data assimilation (ODA) component. In ADA, surface pressures (Ps) are assimilated, while in ODA, gridded sea surface temperature (SST) and ocean temperature and salinity profiles at real Argo locations are assimilated. The system has been evaluated within a perfect twin experiment framework, showing significantly reduced errors of the model atmosphere and ocean states through “observation”-constraints by ADA and ODA. The weakly CDA in which both the online ADA and ODA are conducted during the coupled model integration shows smaller errors of air-sea fluxes than the single ADA and ODA, facilitating the future utilization of cross-covariance between the atmosphere and ocean at the air-sea interface. A three-year CDA reanalysis experiment is also implemented by assimilating Ps, SST and ocean temperature and salinity profiles from the real world spanning the period 1978 to 1980 using 12 ensemble members. Results show that Ps RMSE is smaller than 20CR and SST RMSE is better than ERA-20C and close to CFSR. The success of the online CESM-ECDA system is the first step to implement a high-resolution long-term climate reanalysis once the algorithm efficiency is much improved.

Adiabatic and diabatic signatures of ocean temperature variability

Anthropogenically induced radiative imbalances in the climate system lead to a slow accumulation of heat in the ocean. This warming is often obscured by natural modes of climate variability such as the El Niño-Southern Oscillation (ENSO), which drive substantial ocean temperature changes as a function of depth and latitude. The use of watermass coordinates has been proposed to help isolate forced signals and filter out fast adiabatic processes associated with modes of variability. However, how much natural modes of variability project into these different coordinate systems has not been quantified. Here we apply a rigorous framework to quantify ocean temperature variability using both a quasi-Lagrangian, watermass-based temperature coordinate and Eulerian depth and latitude coordinates in a free-running climate model under pre-industrial conditions. The temperature-based coordinate removes the adiabatic component of ENSO-dominated interannual variability by definition, but a substantial diabatic signal remains. At slower (decadal to centennial) frequencies, variability in the temperature- and depth-based coordinates is comparable. Spectral analysis of temperature tendencies reveals the dominance of advective processes in latitude and depth coordinates while the variability in temperature coordinates is related closely to the surface forcing. Diabatic mixing processes play an important role at slower frequencies where quasi steady-state balances emerge between forcing and mixing in temperature, advection and mixing in depth, and forcing and advection in latitude. While watermass-based analyses highlight diabatic effects by removing adiabatic variability, our work shows that natural variability has a strong diabatic component and cannot be ignored in the analysis of long term trends.

Embedding a One-column Ocean Model (SIT 1.06) in the Community Atmosphere Model 5.3 (CAM5.3; CAM5–SIT v1.0) to Improve Madden–Julian Oscillation Simulation in Boreal Winter

The effect of the air–sea interaction on the Madden–Julian Oscillation (MJO) was investigated using the one-column ocean model Snow–Ice–Thermocline (SIT 1.06) embedded in the Community Atmosphere Model 5.3 (CAM5.3; hereafter CAM5–SIT v1.0). The SIT model with 41 vertical layers was developed to simulate sea surface temperature (SST) and upper-ocean temperature variations with a high vertical resolution that resolves the cool skin and diurnal warm layer and the upper oceanic mixed layer. A series of 30-year sensitivity experiments were conducted in which various model configurations (e.g., coupled versus uncoupled, vertical resolution and depth of the SIT model, coupling domains, and absence of the diurnal cycle) were considered to evaluate the effect of air–sea coupling on MJO simulation. Most of the CAM5–SIT experiments exhibited higher fidelity than the CAM5-alone experiment in characterizing the basic features of the MJO such as spatiotemporal variability and the eastward propagation in boreal winter. The overall MJO simulation performance of CAM5–SIT benefited from (1) better resolving the fine structure of upper-ocean temperature and therefore the air–sea interaction that resulted in more realistic intraseasonal variability in both SST and atmospheric circulation and (2) the adequate thickness and vertical resolution of the oceanic mixed layer. The sensitivity experiments demonstrated the necessity of coupling the tropical eastern Pacific in addition to the tropical Indian Ocean and the tropical western Pacific. Enhanced MJO could be obtained without considering the diurnal cycle in coupling.

Cenozoic Antarctic climate evolution based on molecular and isotopic biomarker reconstructions from geological archives in the Ross Sea region

<p>During the Cenozoic Era (the last 65 Ma), Antarctica’s climate has evolved from ice free conditions of the ‘Greenhouse world’, which at its peak (~ 55 Ma) supported near-tropical forests, to the ‘Icehouse’ climate of today with permanent ice sheets, and a very sparse macroflora. This long-term cooling trend is punctuated by a number of major, abrupt, and in some cases, irreversible climate transitions. Reconstructing past changes in vegetation, sea surface temperature, hydroclimate and the carbon cycle require robust geological proxies that in turn can provide insights into climatic thresholds and feedbacks that drove major transitions in the evolution of Antarctica’s ice sheets. Biomarkers allow climate and environmental proxy reconstructions for this region, where other more traditional paleoclimate methods are less suitable. This study has two aims. Firstly to assess the suitability and applicability of biomarkers in Antarctic sediments across a range of depositional settings and ages, and secondly to apply biomarker-based climate proxies to reconstruct environmental and climate conditions during key periods in the development of the Antarctic Ice Sheets. The distribution and abundances of n-alkanes are assessed in Oligocene and Miocene sediments from a terrestrial outcrop locality in the Transantarctic Mountains, and two glaciomarine sediment cores and an ice-distal deep marine core from the western Ross Sea. Comparisons are made with n-alkane distributions in Eocene glacial erratics and sedimentary rocks of the Mesozoic Beacon Supergroup, both likely sources of reworked material. A shift in dominant chain length from n-C₂₉ to n-C₂₇ occurs between the Late Eocene and Early Oligocene, considered a response to a significant climate cooling. Samples from glaciofluvial environments onshore, and subglacial and ice-proximal environments offshore display a reworked n-alkane distribution, characterised by low carbon preference index (CPI), high average chain length (ACL) and high n-C₂₉/n-C₂₇ values. Whereas, samples from lower-energy, more benign lacustrine and ice-distal marine environments predominantly contained contemporary material. Palynomorphs and biomarker proxies based on n-alkanes and glycerol dialkyl glycerol tetraethers (GDGTs) are applied to a Late Oligocene and Early Miocene glaciomarine succession spanning the large transient excursion of the Mi-1 glaciation (~23 Ma) in DSDP Site 270 drill core from the central Ross Sea. While the Late Oligocene is marked by relatively warm conditions, regional cooling initiated a transition into Mi-1. This was likely driven by a combination of decreasing atmospheric CO₂ and an orbital geometry favouring low seasonality and cool summers, leading to an intensification of proto-Antarctic bottom water production as the Ross Sea deepened and cooled. Mi-1 manifests as a regionally cool period, with minimum subsurface temperatures of ~4°C and onshore mean summer temperatures of ~8°C. A negative n-alkane δ¹³C excursion of up to 4.8‰ is interpreted as a vegetation response to cold, restricted growing seasons, with plants driven to lower altitudes and more stunted growth forms. However, ocean temperatures remained too warm for marine-based ice sheets to advance onto the outer continental shelf and over-ride the drill site. The large increase in ice volume associated with this event, implied by global δ¹⁸O records, was probably held on a higher, terrestrial West Antarctica of greater extent than present day. The relative lack of ice rafted debris during Mi-1, suggests the presence of a marginal marine-terminating ice sheet with fringing ice shelves to the south of DSDP site 270, calving icebergs lacking a basal debris layer, similar to those calving from the Ross Ice Shelf today. This extensive ice cover may explain a large decrease in marine n-alkanes at this time restricting marine productivity on the continental shelf. The biomarker data for the Early Miocene in DSDP 270 indicates a relative warming in both terrestrial and marine temperatures following the transient Mi-1 glacial expansion, but an overall baseline cooling of climate between Late Oligocene and the Early Miocene in the Ross Sea embayment. Isoprenoid GDGTs are used to reconstruct a Cenozoic subsurface ocean temperature compilation for the Ross Sea, a key source region of ocean deep water. The ocean temperature TEXL86 calibration and BAYSPAR in standard subsurface mode were considered, through comparison with independent microfossil and sedimentological data, the most appropriate for use in this region. Ocean temperatures cool prior to the Eocene/Oligocene transition and remain cool for the rest of the Cenozoic, with the exception of short periods of relative warmth in the Late Oligocene and Mid-Miocene Climate Optimum, and long-term trends broadly mirror that of the foraminiferal δ¹⁸O record from the deep Pacific. The Δ Ring Index is used to assess non-thermal influences on GDGT distributions, and displays a long term shift from more positive to more negative deviations. This correlates with %GDGT-0, and also relates to a declining trend in the Methane Index, which reflect the contribution of methanogenic and methanotrophic archaea. These changes suggest that these archaea contributed more to the archaeal community in the early to mid Cenozoic, potentially indicating a more anoxic depositional environment in the Ross Sea. The Branched to Isoprenoid Tetraether index (BIT) steadily declines over the Cenozoic, reflecting increasingly hyper-arid conditions onshore, with less active glaciofluvial systems, limited soil development and less ice-free land.</p>

Cenozoic Antarctic climate evolution based on molecular and isotopic biomarker reconstructions from geological archives in the Ross Sea region

<p>During the Cenozoic Era (the last 65 Ma), Antarctica’s climate has evolved from ice free conditions of the ‘Greenhouse world’, which at its peak (~ 55 Ma) supported near-tropical forests, to the ‘Icehouse’ climate of today with permanent ice sheets, and a very sparse macroflora. This long-term cooling trend is punctuated by a number of major, abrupt, and in some cases, irreversible climate transitions. Reconstructing past changes in vegetation, sea surface temperature, hydroclimate and the carbon cycle require robust geological proxies that in turn can provide insights into climatic thresholds and feedbacks that drove major transitions in the evolution of Antarctica’s ice sheets. Biomarkers allow climate and environmental proxy reconstructions for this region, where other more traditional paleoclimate methods are less suitable. This study has two aims. Firstly to assess the suitability and applicability of biomarkers in Antarctic sediments across a range of depositional settings and ages, and secondly to apply biomarker-based climate proxies to reconstruct environmental and climate conditions during key periods in the development of the Antarctic Ice Sheets. The distribution and abundances of n-alkanes are assessed in Oligocene and Miocene sediments from a terrestrial outcrop locality in the Transantarctic Mountains, and two glaciomarine sediment cores and an ice-distal deep marine core from the western Ross Sea. Comparisons are made with n-alkane distributions in Eocene glacial erratics and sedimentary rocks of the Mesozoic Beacon Supergroup, both likely sources of reworked material. A shift in dominant chain length from n-C₂₉ to n-C₂₇ occurs between the Late Eocene and Early Oligocene, considered a response to a significant climate cooling. Samples from glaciofluvial environments onshore, and subglacial and ice-proximal environments offshore display a reworked n-alkane distribution, characterised by low carbon preference index (CPI), high average chain length (ACL) and high n-C₂₉/n-C₂₇ values. Whereas, samples from lower-energy, more benign lacustrine and ice-distal marine environments predominantly contained contemporary material. Palynomorphs and biomarker proxies based on n-alkanes and glycerol dialkyl glycerol tetraethers (GDGTs) are applied to a Late Oligocene and Early Miocene glaciomarine succession spanning the large transient excursion of the Mi-1 glaciation (~23 Ma) in DSDP Site 270 drill core from the central Ross Sea. While the Late Oligocene is marked by relatively warm conditions, regional cooling initiated a transition into Mi-1. This was likely driven by a combination of decreasing atmospheric CO₂ and an orbital geometry favouring low seasonality and cool summers, leading to an intensification of proto-Antarctic bottom water production as the Ross Sea deepened and cooled. Mi-1 manifests as a regionally cool period, with minimum subsurface temperatures of ~4°C and onshore mean summer temperatures of ~8°C. A negative n-alkane δ¹³C excursion of up to 4.8‰ is interpreted as a vegetation response to cold, restricted growing seasons, with plants driven to lower altitudes and more stunted growth forms. However, ocean temperatures remained too warm for marine-based ice sheets to advance onto the outer continental shelf and over-ride the drill site. The large increase in ice volume associated with this event, implied by global δ¹⁸O records, was probably held on a higher, terrestrial West Antarctica of greater extent than present day. The relative lack of ice rafted debris during Mi-1, suggests the presence of a marginal marine-terminating ice sheet with fringing ice shelves to the south of DSDP site 270, calving icebergs lacking a basal debris layer, similar to those calving from the Ross Ice Shelf today. This extensive ice cover may explain a large decrease in marine n-alkanes at this time restricting marine productivity on the continental shelf. The biomarker data for the Early Miocene in DSDP 270 indicates a relative warming in both terrestrial and marine temperatures following the transient Mi-1 glacial expansion, but an overall baseline cooling of climate between Late Oligocene and the Early Miocene in the Ross Sea embayment. Isoprenoid GDGTs are used to reconstruct a Cenozoic subsurface ocean temperature compilation for the Ross Sea, a key source region of ocean deep water. The ocean temperature TEXL86 calibration and BAYSPAR in standard subsurface mode were considered, through comparison with independent microfossil and sedimentological data, the most appropriate for use in this region. Ocean temperatures cool prior to the Eocene/Oligocene transition and remain cool for the rest of the Cenozoic, with the exception of short periods of relative warmth in the Late Oligocene and Mid-Miocene Climate Optimum, and long-term trends broadly mirror that of the foraminiferal δ¹⁸O record from the deep Pacific. The Δ Ring Index is used to assess non-thermal influences on GDGT distributions, and displays a long term shift from more positive to more negative deviations. This correlates with %GDGT-0, and also relates to a declining trend in the Methane Index, which reflect the contribution of methanogenic and methanotrophic archaea. These changes suggest that these archaea contributed more to the archaeal community in the early to mid Cenozoic, potentially indicating a more anoxic depositional environment in the Ross Sea. The Branched to Isoprenoid Tetraether index (BIT) steadily declines over the Cenozoic, reflecting increasingly hyper-arid conditions onshore, with less active glaciofluvial systems, limited soil development and less ice-free land.</p>

Impact of large volcanic events on the marine environment recorded in marine cores

<p>Explosive silicic volcanic eruptions blanket widespread terrestrial and marine areas in ash, and have a profound effect on climate and local ecosystems. Short-term climate effects are caused by the dispersal of ash, but the injection of gas into the stratosphere, with sulphur being particularly important, drives a cooling of the climate that can last several years. These prolonged perturbations have been observed and recorded in recent decades, but despite the importance of the ocean in regulating global atmospheric climate, little is known about how and to what extent the climate signal produced by volcanic eruptions alters the oceanic environment. As the composition of foraminifera tests is highly sensitive to changes in the surrounding environment, a significant sea surface temperature decrease following a large silicic volcanic eruption may be recorded in the tests of live planktic foraminifera, now preserved in marine sediments. This study examines marine cores (and foraminifera within) that contain tephra units from three major volcanic events to determine if changes can be resolved in ocean temperature and/or foraminifera test morphology following large silicic eruptions. The Holocene Taupo, Waimihia and Mamaku tephra units have been identified in a series of marine sediment cores collected from areas with high sedimentation rates off the east coast of North Island, New Zealand. The sources of these eruptions were from two calderas within the Taupo Volcanic Zone, one of the most active and important rhyolitic regions in the world. Sampling of sediment and foraminifera from these cores has been undertaken at 0.5 cm intervals above and below each tephra. This equates to varying sampling resolutions between cores of 5-30 years, with sufficient sampling taken to establish a stratigraphic record of >100 years either side of each tephra unit. A detailed stratigraphy was undertaken on the sediment surrounding all tephra units, including grain size and CaCO₃ analyses, to identify primary and secondary tephra deposits. One core, Tan0810-12 that contained solely Taupo tephra, was selected for foraminiferal analyses to determine changes in ocean temperature and foraminifera test morphology following this eruption. This core was selected based on the results of the stratigraphic analyses that identified the tephra as a primary deposit with minimal bioturbation above the ash layer and a very high sedimentation rate that enabled sub-decadal scale sampling. Scanning electron microscope imaging was employed to identify the presence of surface contaminants on and within the foraminifera tests and allowed observations of test morphology and size. The morphologies of planktic foraminifera species Globigerinoides ruber and Globigerina bulloides showed no obvious change following the Taupo eruption. The Globigerinoides ruber test sizes distinctly decreased for a period after the eruption, while Globigerina bulloides tests slightly increased in size, correlating well with a decrease in sea surface temperature after the eruption as these species prefer warmer and colder temperatures, respectively. This suggests there is potential for test size to be employed as a proxy for temperature change in conjunction with geochemical analyses. Mg/Ca temperature analyses were conducted in situ using laser ablation inductively coupled plasma mass spectrometry. Both species indicated a decrease in sea surface temperatures when comparing results from tests collected below the tephra deposit to those above. Further results indicate ocean temperature may not have recovered for more than 65 years after the eruption. Such a rapid change in the oceanic environment not only has drastic implications for marine ecosystems but also atmospheric climate, and therefore, terrestrial ecosystems. To reduce the margin of error and determine a more exact value of temperature change following the eruption a greater population of foraminifera is needed. Nonetheless, this study highlights the potential of this method in determining how the oceans are impacted by volcanism and how further research is needed to determine the effects of volcanic eruptions on past and future climate.</p>

Impact of large volcanic events on the marine environment recorded in marine cores

<p>Explosive silicic volcanic eruptions blanket widespread terrestrial and marine areas in ash, and have a profound effect on climate and local ecosystems. Short-term climate effects are caused by the dispersal of ash, but the injection of gas into the stratosphere, with sulphur being particularly important, drives a cooling of the climate that can last several years. These prolonged perturbations have been observed and recorded in recent decades, but despite the importance of the ocean in regulating global atmospheric climate, little is known about how and to what extent the climate signal produced by volcanic eruptions alters the oceanic environment. As the composition of foraminifera tests is highly sensitive to changes in the surrounding environment, a significant sea surface temperature decrease following a large silicic volcanic eruption may be recorded in the tests of live planktic foraminifera, now preserved in marine sediments. This study examines marine cores (and foraminifera within) that contain tephra units from three major volcanic events to determine if changes can be resolved in ocean temperature and/or foraminifera test morphology following large silicic eruptions. The Holocene Taupo, Waimihia and Mamaku tephra units have been identified in a series of marine sediment cores collected from areas with high sedimentation rates off the east coast of North Island, New Zealand. The sources of these eruptions were from two calderas within the Taupo Volcanic Zone, one of the most active and important rhyolitic regions in the world. Sampling of sediment and foraminifera from these cores has been undertaken at 0.5 cm intervals above and below each tephra. This equates to varying sampling resolutions between cores of 5-30 years, with sufficient sampling taken to establish a stratigraphic record of >100 years either side of each tephra unit. A detailed stratigraphy was undertaken on the sediment surrounding all tephra units, including grain size and CaCO₃ analyses, to identify primary and secondary tephra deposits. One core, Tan0810-12 that contained solely Taupo tephra, was selected for foraminiferal analyses to determine changes in ocean temperature and foraminifera test morphology following this eruption. This core was selected based on the results of the stratigraphic analyses that identified the tephra as a primary deposit with minimal bioturbation above the ash layer and a very high sedimentation rate that enabled sub-decadal scale sampling. Scanning electron microscope imaging was employed to identify the presence of surface contaminants on and within the foraminifera tests and allowed observations of test morphology and size. The morphologies of planktic foraminifera species Globigerinoides ruber and Globigerina bulloides showed no obvious change following the Taupo eruption. The Globigerinoides ruber test sizes distinctly decreased for a period after the eruption, while Globigerina bulloides tests slightly increased in size, correlating well with a decrease in sea surface temperature after the eruption as these species prefer warmer and colder temperatures, respectively. This suggests there is potential for test size to be employed as a proxy for temperature change in conjunction with geochemical analyses. Mg/Ca temperature analyses were conducted in situ using laser ablation inductively coupled plasma mass spectrometry. Both species indicated a decrease in sea surface temperatures when comparing results from tests collected below the tephra deposit to those above. Further results indicate ocean temperature may not have recovered for more than 65 years after the eruption. Such a rapid change in the oceanic environment not only has drastic implications for marine ecosystems but also atmospheric climate, and therefore, terrestrial ecosystems. To reduce the margin of error and determine a more exact value of temperature change following the eruption a greater population of foraminifera is needed. Nonetheless, this study highlights the potential of this method in determining how the oceans are impacted by volcanism and how further research is needed to determine the effects of volcanic eruptions on past and future climate.</p>

Genesis and evolution of magmas according to data on hot heterogeneous accretion of the Earth

The obtained numerous proofs of hot heterogeneous accretion of the Earth lead to a fundamentally new solution to the problems of genesis and evolution of magmas. According to these data, the Earth’s core was formed earlier than the silicate mantle as a result of the agglutination of iron particles of the protoplanetary disk under the influence of magnetic forces, because with a small body size, these forces were billions of times more powerful thangravitational ones. The accretion of the silicate mantle created a global magmatic ocean under the influence of impact heat release. Its bottom part crystallized and fractionated as a result of the pressure increase of the formed upper parts. Cumulates formed the ultrabasic mantle, and residual melts formed the magmatic ocean. The increase in ocean temperature and depth caused the evolution of bottom residual melts from acidic to ultrabasic, the appearance of corresponding layers in the ocean, and the reverse geothermal gradient in the mantle. As a result of the cooling and crystallization of the ocean from top to bottom after 3.8 billion years ago early Precambrian crystal complexes, acidic crust, and the lithosphere of ancient platforms were formed. The separation of residual melts from various layers caused the evolution of magmatism on them from acidic to akaline-ultramafic and kimberlite. Heating of the mantle by a high-temperature core led to the appearance of a direct geothermal gradient at the end of the Proterozoic, convection in the mantle, and modern geodynamic environments. In them, magmas are formed by the frictional and decompression melting of the differentiates of the magmatic ocean.

Evolution of mean ocean temperature in Marine Isotope Stage 4

Deglaciations are characterized by relatively fast and near-synchronous changes in ice sheet volume, ocean temperature, and atmospheric greenhouse gas concentrations, but glacial inception occurs more gradually. Understanding the evolution of ice sheet, ocean, and atmosphere conditions from interglacial to glacial maximum provides insight into the interplay of these components of the climate system. Using noble gas measurements in ancient ice samples, we reconstruct mean ocean temperature (MOT) from 74 to 59.7 ka, covering the Marine Isotope Stage (MIS) 5a–4 boundary, MIS 4, and part of the MIS 4–3 transition. Comparing this MOT reconstruction to previously published MOT reconstructions from the last and penultimate deglaciation, we find that the majority of the last interglacial–glacial ocean cooling must have occurred within MIS 5. MOT reached equally cold conditions in MIS 4 as in MIS 2 (−2.7 ± 0.3 ∘C relative to the Holocene, −0.1 ± 0.3 ∘C relative to MIS 2). Using a carbon cycle model to quantify the CO2 solubility pump, we show that ocean cooling can explain most of the CO2 drawdown (32 ± 4 of 40 ppm) across MIS 5. Comparing MOT to contemporaneous records of benthic δ18O, we find that ocean cooling can also explain the majority of the δ18O increase across MIS 5 (0.7 ‰ of 1.3 ‰). The timing of ocean warming and cooling in the record and the comparison to coeval Antarctic isotope data suggest an intimate link between ocean heat content, Southern Hemisphere high-latitude climate, and ocean circulation on orbital and millennial timescales.