Sea Surface Temperature Estimates Research Articles

The bulk parameterizations of turbulent air–sea fluxes in NEMO4: the origin of sea surface temperature differences in a global model study

Abstract. Wind stress and turbulent heat fluxes are the major driving forces that modify the ocean dynamics and thermodynamics. In the Nucleus for European Modelling of the Ocean (NEMO) ocean general circulation model, these turbulent air–sea fluxes (TASFs) can critically impact the simulated ocean characteristics. This paper investigates how the various bulk parameterizations used to calculate turbulent air-sea fluxes in NEMOv4 can lead to substantial differences in the estimation of sea surface temperatures (SSTs). Specifically, we study the contributions of different aspects and assumptions of the bulk parameterizations in driving the SST differences in the NEMO global model configuration at 1/4∘ of horizontal resolution. These aspects include the use of the skin temperature instead of the bulk SST in the computation of turbulent heat flux components and the estimation of wind stress and turbulent heat flux components, which vary in each parameterization due to different bulk transfer coefficients. The analysis of a set of short-term sensitivity experiments where the only change is related to one of the aspects of the bulk parameterizations shows that parameterization-related SST differences are primarily sensitive to wind stress differences and to the implementation of skin temperature in the computation of turbulent heat flux components. In addition, in order to highlight the role of SST–turbulent heat flux negative feedback at play in ocean simulations, we compare the TASF differences obtained using the NEMO ocean model with the estimations by Brodeau et al. (2017), who compared the different bulk parameterizations using prescribed SSTs. Our estimations of turbulent heat flux differences between bulk parameterizations are weaker than those found by Brodeau et al. (2017).

Cenozoic Antarctic Peninsula Temperatures and Glacial Erosion Signals From a Multi‐Proxy Biomarker Study

AbstractTerrestrial climate records for Antarctica, beyond the age limit of ice cores, are restricted to the few unglaciated areas with exposed rock outcrops. Marine sediments on Antarctica's continental shelves contain records of past oceanic and terrestrial environments that can provide important insights into Antarctic climate evolution. The SHALDRIL II (Shallow Drilling on the Antarctic Continental Margin) expedition recovered sedimentary sequences from the eastern side of the Antarctic Peninsula of late Eocene, Oligocene, middle Miocene, and early Pliocene age that provides insights into Cenozoic Antarctic climate and ice sheet development. Here, we use biomarker data to assess atmospheric and oceanic temperatures and glacial reworking from the late Eocene to the early Pliocene. Analyses of hopanes and n‐alkanes indicate increased erosion of mature (thermally altered) soil biomarker components reworked by glacial erosion. Branched glycerol dialkyl glycerol tetraethers from soil bacteria suggest similar air temperatures of 12°C ± 1°C (1σ, n = 46) for months above freezing for Eocene, Oligocene, and Miocene timeslices but much colder (and likely shorter) periods of thaw during the Pliocene (5°C ± 1°C, n = 4) on the Antarctic Peninsula. TEX86‐based (Tetraether index of 86 carbons) sea surface temperature estimates indicate ocean cooling from 7°C ± 3°C (n = 10) in the Miocene to 3°C ± 1°C (n = 3) in the Pliocene, consistent with deep ocean cooling. Resulting temperature records provide useful constraints for ice sheet and climate model simulations seeking to improve understanding of ice sheet response under a range of climate conditions.

Accuracy and Reproducibility of Coral Sr/Ca SIMS Timeseries in Modern and Fossil Corals

AbstractCoral strontium‐to‐calcium ratios (Sr/Ca) provide quantitative estimates of past sea surface temperatures (SST) that allow for the reconstruction of changes in the mean state and climate variations, such as the El Nino‐Southern Oscillation, through time. However, coral Sr/Ca ratios are highly susceptible to diagenesis, which can impart artifacts of 1–2°C that are typically on par with the tropical climate signals of interest. Microscale sampling via Secondary Ion Mass Spectrometry (SIMS) for the sampling of primary skeletal material in altered fossil corals, providing much‐needed checks on fossil coral Sr/Ca‐based paleotemperature estimates. In this study, we employ a set modern and fossil corals from Palmyra Atoll, in the central tropical Pacific, to quantify the accuracy and reproducibility of SIMS Sr/Ca analyses relative to bulk Sr/Ca analyses. In three overlapping modern coral samples, we reproduce bulk Sr/Ca estimates within ±0.3% (1σ). We demonstrate high fidelity between 3‐month smoothed SIMS coral Sr/Ca timeseries and SST (R = −0.5 to −0.8; p < 0.5). For lightly‐altered sections of a young fossil coral from the early‐20th century, SIMS Sr/Ca timeseries reproduce bulk Sr/Ca timeseries, in line with our results from modern corals. Across a moderately‐altered section of the same fossil coral, where diagenesis yields bulk Sr/Ca estimates that are 0.6 mmol too high (roughly equivalent to −6°C artifacts in SST), SIMS Sr/Ca timeseries track instrumental SST timeseries. We conclude that 3–4 SIMS analyses per month of coral growth can provide a much‐needed quantitative check on the accuracy of fossil coral Sr/Ca‐derived estimates of paleotemperature, even in moderately altered samples.

A new global surface temperature reconstruction for the Last Glacial Maximum

Abstract. We present a new reconstruction of surface air temperature and sea surface temperature for the Last Glacial Maximum. The method blends model fields and sparse proxy-based point estimates through a data assimilation approach. Our reconstruction updates that of Annan and Hargreaves (2013), using the full range of general circulation model (GCM) simulations which contributed to three generations of the PMIP database, three major compilations of gridded sea surface temperature (SST) and surface air temperature (SAT) estimates from proxy data, and an improved methodology based on an ensemble Kalman filter. Our reconstruction has a global annual mean surface air temperature anomaly of -4.5±0.9 ∘C relative to the pre-industrial climate. This is slightly colder than the previous estimate of Annan and Hargreaves (2013), with an upwards revision on the uncertainty due to different methodological assumptions. It is, however, substantially less cold than the recent reconstruction of Tierney et al. (2020). We show that the main reason for this discrepancy is in the choice of prior. We recommend the use of the multi-model ensemble of opportunity as potentially offering a credible prior, but it is important that the range of models included in the PMIP ensembles represent the main sources of uncertainty as realistically and comprehensively as practicable if they are to be used in this way.

Coupled Coccolith-Based Temperature and Productivity High-Resolution Reconstructions in the Eastern Equatorial Pacific During the Last Deglaciation and the Holocene

We present a new high-resolution reconstruction of annual sea-surface temperatures (SSTa) and net primary productivity (NPP) using novel coccolithophore-based models developed for the Eastern Equatorial Pacific (EEP). We combined published coccolithophore census counts from core-tops in the Eastern Pacific with 32 new samples from the Equatorial region, to derive a new statistical model to reconstruct SSTa. Results show that the addition of the new EEP samples improves existing coccolithophore-based SST-calibrations, and allow reconstructing SSTa in the EEP with higher confidence. We also merged the relative abundance of deep-photic species Florisphaera profunda in the same surface sediment samples with existing calibration datasets for tropical regions, to reconstruct annual NPP. Both temperature and productivity calibrations were successfully applied to fossil coccolith data from Ocean Drilling Project Site 1240, in the EEP. The coccolith-based SSTa estimates show a cooling during the Last Glacial Maximum (LGM) and the Younger Dryas, and warming at the start of the Holocene. This pattern differs in the timing and magnitude of the temperature changes from other available SST-reconstructions based on biogeochemical and faunal proxies. We discuss these discrepancies to be the result of different proxy sensitivities to insolation forcing, seasonal bias, and/or preservation artifacts. Reconstructed annual NPP shows a general decreasing trend from the late last glacial period to recent times, which we relate to the weakening of wind-driven equatorial upwelling towards the Holocene. We also calculated carbon export using our SSTa and NPP reconstructions, and compared to other geochemical-based reconstructions for the same location. Our coupled SSTa-NPP reconstruction provides key data to more fully assess the evolution of primary and export productivity as well as organic carbon burial in the EEP, with implications for its role in global biogeochemical cycles across glacial terminations.

Relative Merits of Optimal Estimation and Non-Linear Retrievals of Sea-Surface Temperature from MODIS

We compared the results of an Optimal Estimation (OE) based approach for the retrieval of the skin sea surface temperature (SSTskin) with those of the traditional non-linear sea surface temperature (NLSST) algorithm. The retrievals were from radiance measurements in two infrared channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) on the NASA satellite Aqua. The OE used a reduced state vector of SST and total column water vapor (TCWV). The SST and atmospheric profiles of temperature and humidity from ERA5 provided prior knowledge, and we made reasonable assumptions about the variance of these fields. An atmospheric radiative transfer model was used as the forward model to simulate the MODIS measurements. The performances of the retrieval approaches were assessed by comparison with in situ measurements. We found that the OESST reduces the satellite–in situ bias, but mostly for retrievals with an already small bias between in situ and the prior SST. The OE approach generally fails to improve the SST retrieval when that difference is large. In such cases, the NLSST often provides a better estimate of the SST than the OE. The OESST also underperforms NLSST in areas that include large horizontal SST gradients.

Formula omitted] Expands the linear dynamic range of the alkenone sea surface temperature proxy

The alkenone paleothermometer, U37K′, has been used to estimate past sea surface temperatures (SST) for nearly forty years. Despite its popularity, U37K′ is conventionally limited to water temperatures between 5 and 28.5 °C. At both low and high temperature ranges [0–8; >24 °C], U37K′ displays nonlinear behavior and significant scattering, leading to major uncertainties in reconstructed SSTs. With data from 171 surface sediments across the global ocean, we demonstrate that U38MEK′, an unsaturation index based on the ratio of di and tri unsaturated methyl alkenones with 38 carbon atoms (C38ME), extends the linear dynamic range of alkenone SST estimates from 0 to 29.5 °C SST. Accurate measurements of U38MEK′ are made possible by recent advances in chromatographic techniques that provide excellent resolution of C38 alkenones, as well as improved separation of all alkenone homologues. Our newly calibrated U38MEK′ permits more accurate SST reconstruction in warmer ocean waters and climates. U38MEK′ also carries less risk than U37K′ of biased SST estimates in regions (e.g. estuaries, closed seas) where both open ocean and coastal alkenone producers (e.g., Isochrysis galbana) are present because the latter do not produce C38 methyl ketones but do produce C37 methyl ketones. The efficacy of the U38MEK′ paleothermometer is shown by downcore analyses of early Pliocene age sediments from ODP sites 806 and 1143, as well as Holocene age sediments from the Black Sea. The U38ETK′ paleothermometer is also evaluated and found to be unsuitable for paleoceanographic applications.

Uncertainty in satellite sea surface temperature with respect to air temperature, dust level, wind speed and solar position

Sea surface temperature (SST) plays a significant role in studying climate change and oceanic environmental processes. SST has always been measured with buoys and thermometers in the field and by remote sensing-based techniques. The use of remotely sensed SST is significantly more common for mapping SST in wide areas and at high temporal resolution. Remotely sensed SST, however, may reveal uncertainties due to factors such as water depth and turbidity. This is in addition to the atmospheric effects, such as dust, fog, wind, solar angles, and humidity. Thus, in this work, SST estimation by the Moderate Resolution Imaging Spectroradiometer (MODIS) is compared with the measured in-situ SST in the Arabian Gulf with respect to the aforementioned atmospheric and oceanic variables. Results show that SST estimation is affected by air temperature, resulting in a higher SST bias in the summer than in winter. Further, the SST estimation is affected by water depth (i.e. bathymetry). As the water depth decreases below 30 m, the SST bias increases due to the sea bottom effect. SST estimates are also influenced by the instantaneous changes in dust levels, particularly when the aerosol optical thickness exceeds 0.3. However, when SST is averaged over each month, this impact is not prominent because the dust events are occurring occasionally, so they do not impact the average. On the other hand, wind speed and solar angles do not contribute to SST residuals , as the bias is caused predominantly by air temperature and shallow water depth.

MIS 5e sea-level history along the Pacific coast of North America

Abstract. The primary last interglacial, marine isotope substage (MIS) 5e records on the Pacific coast of North America, from Washington (USA) to Baja California Sur (Mexico), are found in the deposits of erosional marine terraces. Warmer coasts along the southern Golfo de California host both erosional marine terraces and constructional coral reef terraces. Because the northern part of the region is tectonically active, MIS 5e terrace elevations vary considerably, from a few meters above sea level to as much as 70 m above sea level. The primary paleo-sea-level indicator is the shoreline angle, the junction of the wave-cut platform with the former sea cliff, which forms very close to mean sea level. Most areas on the Pacific coast of North America have experienced uplift since MIS 5e time, but the rate of uplift varies substantially as a function of tectonic setting. Chronology in most places is based on uranium-series ages of the solitary coral Balanophyllia elegans (erosional terraces) or the colonial corals Porites and Pocillopora (constructional reefs). In areas lacking corals, correlation to MIS 5e often can be accomplished using amino acid ratios of fossil mollusks, compared to similar ratios in mollusks that also host dated corals. Uranium-series (U-series) analyses of corals that have experienced largely closed-system histories range from ∼124 to ∼118 ka, in good agreement with ages from MIS 5e reef terraces elsewhere in the world. There is no geomorphic, stratigraphic, or geochronological evidence for more than one high-sea stand during MIS 5e on the Pacific coast of North America. However, in areas of low uplift rate, the outer parts of MIS 5e terraces apparently were re-occupied by the high-sea stand at ∼100 ka (MIS 5c), evident from mixes of coral ages and mixes of molluscan faunas with differing thermal aspects. This sequence of events took place because glacial isostatic adjustment processes acting on North America resulted in regional high-sea stands at ∼100 and ∼80 ka that were higher than is the case in far-field regions, distant from large continental ice sheets. During MIS 5e time, sea surface temperatures (SSTs) off the Pacific coast of North America were higher than is the case at present, evident from extralimital southern species of mollusks found in dated deposits. Apparently, no wholesale shifts in faunal provinces took place, but in MIS 5e time, some species of bivalves and gastropods lived hundreds of kilometers north of their present northern limits, in good agreement with SST estimates derived from foraminiferal records and alkenone-based reconstructions in deep-sea cores. Because many areas of the Pacific coast of North America have been active tectonically for much or all of the Quaternary, many earlier interglacial periods are recorded as uplifted, higher-elevation terraces. In addition, from southern Oregon to northern Baja California, there are U-series-dated corals from marine terraces that formed at ∼80 ka, during MIS 5a. In contrast to MIS 5e, these terrace deposits host molluscan faunas that contain extralimital northern species, indicating cooler SST at the end of MIS 5. Here I present a review and standardized database of MIS 5e sea-level indicators along the Pacific coast of North America and the corresponding dated samples. The database is available in Muhs et al. (2021b; https://doi.org/10.5281/zenodo.5903285).

Sea ice changes in the southwest Pacific sector of the Southern Ocean during the last 140 000 years

Abstract. Sea ice expansion in the Southern Ocean is believed to have contributed to glacial–interglacial atmospheric CO2 variability by inhibiting air–sea gas exchange and influencing the ocean's meridional overturning circulation. However, limited data on past sea ice coverage over the last 140 ka (a complete glacial cycle) have hindered our ability to link sea ice expansion to oceanic processes that affect atmospheric CO2 concentration. Assessments of past sea ice coverage using diatom assemblages have primarily focused on the Last Glacial Maximum (∼21 ka) to Holocene, with few quantitative reconstructions extending to the onset of glacial Termination II (∼135 ka). Here we provide new estimates of winter sea ice concentrations (WSIC) and summer sea surface temperatures (SSST) for a full glacial–interglacial cycle from the southwestern Pacific sector of the Southern Ocean using the modern analog technique (MAT) on fossil diatom assemblages from deep-sea core TAN1302-96. We examine how the timing of changes in sea ice coverage relates to ocean circulation changes and previously proposed mechanisms of early glacial CO2 drawdown. We then place SSST estimates within the context of regional SSST records to better understand how these surface temperature changes may be influencing oceanic CO2 uptake. We find that winter sea ice was absent over the core site during the early glacial period until MIS 4 (∼65 ka), suggesting that sea ice may not have been a major contributor to early glacial CO2 drawdown. Sea ice expansion throughout the glacial–interglacial cycle, however, appears to coincide with observed regional reductions in Antarctic Intermediate Water production and subduction, suggesting that sea ice may have influenced intermediate ocean circulation changes. We observe an early glacial (MIS 5d) weakening of meridional SST gradients between 42 and 59∘ S throughout the region, which may have contributed to early reductions in atmospheric CO2 concentrations through its impact on air–sea gas exchange.

A Simple Monte Carlo Approach to Estimate the Uncertainties of SST and δ18Osw Inferred From Coral Proxies

AbstractThe skeletons of massive corals inhabiting tropical and subtropical oceans grow over centuries. The chemical and isotopic compositions of their skeletons reflect changes in the marine environment and provide an archive of sea surface temperature (SST) and salinity at a high temporal resolution spanning decades to centuries. Stable oxygen isotope ratios in coral skeletons (δ18Oc) are an excellent proxy for SST and the stable oxygen isotopic composition in seawater (δ18Osw). Variations in δ18Osw serve as a hydrological proxy and are routinely estimated by subtracting the SST component from δ18Oc. The SST component of δ18Oc is inferred from coral Sr/Ca (SST proxy). Both Sr/Ca and δ18Oc are converted to SST units based on Sr/Ca‐ and (δ18Oc–δ18Osw)‐SST calibrations. However, the exact slope values of the Sr/Ca‐SST and (δ18Oc–δ18Osw)‐SST calibrations are challenging to assess, which hampers the reliable estimation of SST, δ18Osw, and their uncertainties. Here, we present a simple Monte Carlo approach to estimate variations in SST and δ18Osw together with their confidence intervals. Our approach improves on previous error propagation methods by including intercolonial differences of calibration slopes in addition to analytical errors. In addition, Monte Carlo experiments using adjusted errors suggest that uncertainties in the calibration slopes have a minor influence on estimates of SST and δ18Osw.

A comparison of orbital-resolution, Late Pleistocene Alkenone and foraminiferal assemblage-based sea surface temperature reconstructions from the Southwest Pacific

Global and regional reconstructions of past climate conditions often incorporate sea surface temperature (SST) estimates from multiple proxies because not every paleotemperature proxy is applicable in all geographic locations. This practice of assimilating estimates from different proxies in global or regional temperature syntheses makes the implicit assumption that estimates derived from different proxies can be meaningfully intercompared. However, evidence to support the validity of this assumption is limited. Using paleotemperature data from sediments collected from ODP Site 1125 in the Southwest Pacific, we conduct a ∼1 Myr, orbital-scale SST proxy comparison of a recently published alkenone-derived SST record with a previously published foraminiferal assemblage-based SST record. These alkenone and foraminiferal assemblage SST datasets show strong structural similarity and yield remarkably similar estimates for basic climate metrics, including mean, median, standard deviation, and range. Statistical analysis indicates that the correlation between the two SST records is highly significant. In the spectral domain, the records share the same dominant 100 kyr beat, are coherent and in phase with each other at this frequency, and have the same coherence and phase relationship with benthic foraminiferal δ18O. Results from this work demonstrate that these two proxies would yield very similar estimates for the paleoclimate metrics most commonly used in empirical paleoclimate reconstructions that seek to document the evolution of climate over this interval. However, significant disparities between SST estimates derived from the two proxies exist for some time periods, particularly during glacial and interglacial extrema. This comparison suggests that treating estimates from these proxies as equivalent in studies that focus on short time windows (e.g. a few thousand to tens-of-thousands of years), particularly in investigations that seek to characterize glacial or interglacial extrema, could be potentially problematic. However, the sensitive location of Site 1125, just north of the Subtropical Front, likely accentuates the difference between temperature estimates from these proxies, which may be attenuated in other oceanographic settings. We attribute the discrepancies between the two SST records to two main causes: seasonal leakages of cold water across the Subtropical Front during glacial extrema and differential seasonality of maximum alkenone and foraminiferal production.

Validation of satellite-based sea surface temperature products against in situ observations off the western coast of Sumatra

This study validated the sea surface temperature (SST) datasets from the Group for High-Resolution SST Multi Product Ensemble (GMPE), National Oceanic and Atmospheric Administration (NOAA) Optimal Interpolation (OI) SST version 2 and 2.1 (OIv2 and OIv2.1), and Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) in the area off the western coast of Sumatra against in situ observations. Furthermore, the root mean square differences (RMSDs) of OIv2, OIv2.1, and ECCO2 were investigated with respect to GMPE, whose small RMSD < 0.2 K against in situ observations confirmed its suitability as a reference. Although OIv2 showed a large RMSD (1–1.5 K) with a significant negative bias, OIv2.1 (RMSD < 0.4 K) improved remarkably. In the average SST distributions for December 2017, the differences among the 4 datasets were significant in the areas off the western coast of Sumatra, along the southern coast of Java, and in the Indonesian inland sea. These results were consistent with the ensemble spread distribution obtained with GMPE. The large RMSDs of OIv2 corresponded to high clouds, and it was suggested that the change in the satellites used for SST estimation contributed to the improvement in OIv2.1.

Cloud-Free Sea-Surface-Temperature Image Reconstruction From Anomaly Inpainting Network

Sea-surface temperature (SST) images obtained by satellites contain noise and missing SSTs due to cloud covers. We propose a method for reconstructing denoised cloud-free SST images via deep-learning-based image inpainting. For denoizing, we use data-assimilation images to train a reconstruction network by considering the physical correctness of SSTs. For reconstruction stability, we introduce <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">anomaly inpainting network</i> , which does not directly complete missing SSTs but estimates the difference between the unobserved SSTs and the average SSTs. SSTs do not fluctuate much over a few days; thus, we can use recent average SSTs as a rough estimation of SSTs and can assume that the SST difference will be within a specific range. We conducted experiments to evaluate our method with satellite SST images and <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in situ</i> SST data. The results indicate that our method with anomaly inpainting network qualitatively and quantitatively outperformed conventional SST image inpainting methods.

Impact of Cloud Reduction on MODIS Thermal Infrared Sea Surface Temperature Retrieval

Considering the significant influence of clouds on thermal infrared (TIR) data sea surface temperature (SST) retrieval, this study focuses on the reduction of cloud influence on Moderate Resolution Imaging Spectroradiometer (MODIS) longwave SST retrieval. First, the quality level (QL) 3 MODIS SSTs are classified into cloudy or clear-sky pixels based on MODIS cloud mask products. The cloudy SST pixels, flagged as confident cloudy or probably cloudy in the cloud masks, are further classified into three groups according to their associated daytime cloud-top and optical properties. Taking the SSTs measured by 11 buoys over one year as the reference data, those three groups of cloudy SSTs are significantly underestimated, with biases of −33.45 °C, −6.35 °C, and −4.72 °C. The QL 3 clear-sky SSTs are identified as probably and confident clear in the cloud masks and are not influenced by clouds, with a bias of nearly 0 °C. Then, three support vector regression (SVR) models are individually proposed for the three groups of cloudy SSTs. Taking the cloud-top and optical parameters as inputs of the proposed SVR models, we can obtain SSTs with a bias of 0 °C and a root-mean-square error (RMSE) of less than 1.6 °C for the three groups of cloudy SSTs. The method proposed in this study shows the potential for TIR SST estimation under some cloudy conditions using satellite remote sensing cloud products. Considering the RMSE of 0.4 °C in operational sea surface temperature and sea ice analysis, further study is needed before actual applying the method proposed in this study.

Surface Temperature trends in the Mediterranean Sea from MODIS data during years 2003–2019

Sea Surface Temperature is a variable recognized as an Essential Climate Variable (ECV) by the Global Climate Observation System (GCOS), due to its determinant influence on climate dynamics, from micro scale to global levels The aim of this paper is to estimate Sea Surface Temperature trends in the Mediterranean Sea during years 2003–2019 by using the MODIS Level 3 SST Thermal IR 8 Day 4km V2019.0. Results show an SST increase of 0.040 ± 0.001 °C/yr. The seasonal maximum trend is associated to summer 0.070 ± 0.001 °C/yr, followed by winter, (0.040 ± 0.001) °C/yr, autumn 0.030 ± 0.001 °C/yr and spring, 0.020 ± 0.001 °C/yr. The total period analyzed has been divided into ten-year time spans, showing a stable increase of 0.055 °C/yr in average, from period 2005–2014 onwards. In absolute SST values terms, the parameter range is of 0.85°C from year 2005 to 2019. We have also analyzed the spatial variability of the parameter by dividing the Mediterranean Sea into nine sub-basins: Alboran Sea, Balearic Sea (Iberian Sea), Mediterranean Western Basin, Mediterranean Eastern Basin, Ligurian Sea, Tyrrhenian Sea, Ionian Sea, Adriatic Sea and Aegean Sea. The results show warming trends from 0.02 °C/yr for the Alboran Sea to 0.07°C/yr for the Adriatic Sea. Results have been validated by using data from a local observational buoy system, obtaining a coefficient of determination of 0.97.

The Yorktown Formation: Improved Stratigraphy, Chronology, and Paleoclimate Interpretations from the U.S. Mid-Atlantic Coastal Plain

The Yorktown Formation records paleoclimate conditions along the mid-Atlantic Coastal Plain during the mid-Piacenzian Warm Period (3.264 to 3.025 Ma), a climate interval of the Pliocene in some ways analogous to near future climate projections. To gain insight into potential near future changes, we investigated Yorktown Formation outcrops and cores in southeastern Virginia, refining the stratigraphic framework. We analyzed 485 samples for alkenone-based sea surface temperature (SST) and productivity estimates from the Holland and Dory cores, an outcrop at Morgarts Beach, Virginia, and the lectostratotype of the Yorktown Formation at Rushmere, Virginia, and analyzed planktonic foraminferal assemblage data from the type section. Using the structure of the SST record, we improved the chronology of the Yorktown Formation by establishing the maximum age ranges of the Rushmere (3.3–3.2 Ma) and Morgarts Beach (3.2–3.15 Ma) Members. SST values for these members average ~26 °C, corroborating existing sclerochronological data. Increasing planktonic foraminifer abundance, productivity, and species diversity parallel increasing SST over the MIS M2/M1 transition. These records constitute the greatest temporal concentration of paleoecological estimates within the Yorktown Formation, aiding our understanding of western North Atlantic temperature patterns, seasonality and ocean circulation during this interval. We provide a chronologic framework for future studies analyzing ecological responses to profound climate change.

Vertical decoupling in Late Ordovician anoxia due to reorganization of ocean circulation

Geochemical redox proxies indicate that seafloor anoxia occurred during the latest Ordovician glacial maximum, coincident with the second pulse of the Late Ordovician mass extinction. However, expanded anoxia in a glacial climate strikingly contrasts with the warming-associated Mesozoic anoxic events and raises questions as to both the causal mechanism of ocean deoxygenation and its relationship with extinction. Here we firstly report iodine-to-calcium ratio (I/Ca) data that document increased upper-ocean oxygenation despite the concurrent expansion of seafloor anoxia. We then resolve these apparently conflicting observations as well as their relationship to global climate by means of a series of Earth system model simulations. Applying available Late Ordovician (Hirnantian) sea-surface temperature estimates from oxygen isotope studies as constraints, alongside our I/Ca data, leads us to identify a scenario in which Hirnantian glacial conditions permit both the spread of seafloor anoxia and increased upper-ocean oxygenation. Our simulated mechanism of a reorganization of global ocean circulation, with reduced importance of northern-sourced waters and a poorer ventilated and deoxygenated deep ocean has parallels with Pleistocene state transitions in Atlantic meridional overturning (despite a very different continental configuration) and suggests that no simple and predictable relationship between past climate state and oxygenation may exist.

A revised mid-Pliocene composite section centered on the M2 glacial event for ODP Site 846

Abstract. The composite section from ODP Site 846 has provided key data sets for Pliocene stable isotope and paleoclimatic time series. We document here apparent outliers in previously published data sets for stable isotopes and alkenone-derived sea surface temperature (SST) estimates in the Pliocene interval containing the M2 glaciation (ca. 3.290–3.3 Ma) by tying high-resolution core measurements to a continuous downhole conductivity log. We generate a revised sequence of new stable isotopic and alkenone measurements across the M2 event that correlate well to the revised splices of color reflectance and gamma ray attenuation porosity evaluator data from Site 846, and to a new composite section produced at equatorial Pacific ODP Site 850. A new composite splice for Site 846 is proposed, along with composite isotope and alkenone time series that should be integrated into revised Pliocene paleoclimatic stacks.

Correcting Observational Biases in Sea Surface Temperature Observations Removes Anomalous Warmth during World War II

AbstractMost historical sea surface temperature (SST) estimates indicate warmer World War II SSTs than expected from forcing and internal climate variability. If real, this World War II warm anomaly (WW2WA) has important implications for decadal variability, but the WW2WA may also arise from incomplete corrections of biases associated with bucket and engine room intake (ERI) measurements. To better assess the origins of the WW2WA, we develop five different historical SST estimates (reconstructions R1–R5). Using uncorrected SST measurements from the International Comprehensive Ocean–Atmosphere Data Set (ICOADS) version 3.0 (R1) gives a WW2WA of 0.41°C. In contrast, using only buckets (R2) or ERI observations (R3) gives WW2WAs of 0.18° and 0.08°C, respectively, implying that uncorrected biases are the primary source of the WW2WA. We then use an extended linear-mixed-effect method to quantify systematic differences between subsets of SSTs and develop groupwise SST adjustments based on differences between pairs of nearby SST measurements. Using all measurements after applying groupwise adjustments (R4) gives a WW2WA of 0.13°C [95% confidence interval (c.i.): 0.01°–0.26°C] and indicates that U.S. and U.K. naval observations are the primary cause of the WW2WA. Finally, nighttime bucket SSTs are found to be warmer than their daytime counterparts during WW2, prompting a daytime-only reconstruction using groupwise adjustments (R5) that has a WW2WA of 0.09°C (95% c.i.: −0.01° to 0.18°C). R5 is consistent with the range of internal variability found in either the CMIP5 (95% c.i.: −0.10° to 0.10°C) or CMIP6 ensembles (95% c.i.: −0.11° to 0.10°C). These results support the hypothesis that the WW2WA is an artifact of observational biases, although further data and metadata analyses will be important for confirmation.