Search

Acknowledgements: S.G., C.C., M.K. and M. Stoffel were supported by the Swiss National Science Foundation Sinergia Project CALDERA (CRSII5_183571). S.G. acknowledges A. Harrak (Department of Near and Middle Eastern Civilizations, University of Toronto), F. Espenak (NASA Goddard Space Flight Center), F. Hierink (Institute for Environmental Sciences, University of Geneva) and P. Souyri (Department of East Asian Studies, University of Geneva) for providing advice on the manuscript. F. Lavigne was supported by Institut Universitaire de France (IUF, Academic Institute of France). M.K. received funding from the EUR IPSL – Climate Graduate School project, managed by the ANR within the “Investissements d’avenir” programme under reference ANR-11-IDEX-0004-17-EURE-0006. F. Ludlow received funding from an Irish Research Council Starting Laureate Award (CLICAB project, IRCLA/2017/303). F. Ludlow and Z.Y. also received funding from a European Research Council (ERC) Synergy Grant (4-OCEANS; grant agreement no. 951649) under the European Union’s Horizon 2020 research and innovation programme. M. Sigl received funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 820047). This paper is a product of the Volcanic Impacts on Climate and Society (VICS) working group., Explosive volcanism is a key contributor to climate variability on interannual to centennial timescales1. Understanding the far-field societal impacts of eruption-forced climatic changes requires firm event chronologies and reliable estimates of both the burden and altitude (that is, tropospheric versus stratospheric) of volcanic sulfate aerosol2, 3. However, despite progress in ice-core dating, uncertainties remain in these key factors4. This particularly hinders investigation of the role of large, temporally clustered eruptions during the High Medieval Period (HMP, 1100–1300 ce), which have been implicated in the transition from the warm Medieval Climate Anomaly to the Little Ice Age5. Here we shed new light on explosive volcanism during the HMP, drawing on analysis of contemporary reports of total lunar eclipses, from which we derive a time series of stratospheric turbidity. By combining this new record with aerosol model simulations and tree-ring-based climate proxies, we refine the estimated dates of five notable eruptions and associate each with stratospheric aerosol veils. Five further eruptions, including one responsible for high sulfur deposition over Greenland circa 1182 ce, affected only the troposphere and had muted climatic consequences. Our findings offer support for further investigation of the decadal-scale to centennial-scale climate response to volcanic eruptions.

The deep roots of volcanic systems play a key role in the priming, initiation, and duration of eruptions. Causative links between initial magmatic unrest at depth and eruption triggering remain poorly constrained. The 2021 CE eruption at Fagradalsfjall in southwestern Iceland, the first deep-sourced eruption on a spreading-ridge system monitored with modern instrumentation, presents an ideal opportunity for comparing geophysical and petrological data sets to explore processes of deep magma mobilization. We used diffusion chronometry to show that deep magmatic unrest in the roots of volcanic systems can precede apparent geophysical eruption precursors by years, suggesting that early phases of magma accumulation and reorganization can occur in the absence of significant increases in shallow seismicity (

Fluid-rock interactions play a critical role in Earth’s lithosphere and environmental subsurface systems. In the absence of chemical mass transport, mineral-hydration reactions would be accompanied by a solid-volume increase that may induce differential stresses and associated reaction-induced deformation processes, such as dilatant fracturing to increase fluid permeability. However, the magnitudes of stresses that manifest in natural systems remain poorly constrained. We used optical and electron microscopy to show that one of the simplest hydration reactions in nature [MgO + H2O = Mg(OH)2] can induce stresses of several hundred megapascals, with local stresses of as much as ∼1.5 GPa. We demonstrate that these stresses not only cause fracturing but also induce plastic deformation with dislocation densities (1015 m−2) exceeding those typical of tectonically deformed rocks. If these reaction-induced stresses can be transmitted across larger length scales, they may influence the bulk stress state of reacting regions. Moreover, the structural damage induced may be the first step toward catastrophic rock failure, triggering crustal seismicity.

International audience; The geological units on the floor of Jezero crater, Mars, are part of a wider regional stratigraphy of olivine-rich rocks, which extends well beyond the crater. We investigate the petrology of olivine and carbonate-bearing rocks of the Séítah formation in the floor of Jezero. Using multispectral images and x-ray fluorescence data, acquired by the Perseverance rover, we performed a petrographic analysis of the Bastide and Brac outcrops within this unit. We find that these outcrops are composed of igneous rock, moderately altered by aqueous fluid. The igneous rocks are mainly made of coarse-grained olivine, similar to some Martian meteorites. We interpret them as an olivine cumulate, formed by settling and enrichment of olivine through multi-stage cooling of a thick magma body.

A new ice core was drilled in West Antarctica on Skytrain Ice Rise in field season 2018/2019. This 651 m ice core is one of the main targets of the WACSWAIN (WArm Climate Stability of the West Antarctic ice sheet in the last INterglacial) project. A present-day accumulation rate of 13.5 cm w.e. yr−1 was derived. Although the project mainly aims to investigate the last interglacial (115–130 ka), a robust chronology period covering the recent past is needed to constrain the age models for the deepest ice. Additionally, this time period is important for understanding current climatic changes in the West Antarctic region. Here, we present a stratigraphic chronology for the top 184.14 m of the Skytrain ice core based on absolute age tie points interpolated using annual layer counting encompassing the last 2000 years of climate history. Together with a model-based depth–age relationship of the deeper part of the ice core, this will form the ST22 chronology. The chemical composition, dust content, liquid conductivity, water isotope concentration and methane content of the whole core was analysed via continuous flow analysis (CFA) at the British Antarctic Survey. Annual layer counting was performed by manual counting of seasonal variations in mainly the sodium and calcium records. This counted chronology was informed and anchored by absolute age tie points, namely, the tritium peak (1965 CE) and six volcanic eruptions. Methane concentration variations were used to further constrain the counting error. A minimal error of ±1 year at the tie points was derived, accumulating to ± 5 %–10 % of the age in the unconstrained sections between tie points. This level of accuracy enables data interpretation on at least decadal timescales and provides a solid base for the dating of deeper ice, which is the second part of the chronology.

Across a 1000‐km stretch of the River Nile, from the 1st Cataract in southern Egypt to the 4th Cataract in Sudan, many hundreds of drystone walls are located within active channels, on seasonally inundated floodplains or in now‐dry Holocene palaeochannel belts. These walls (or river groynes) functioned as flood and flow control structures and are of a type now commonly in use worldwide. In the Nile Valley, the structures have been subject only to localised investigations, and none have been radiometrically dated. Some were built within living memory to trap nutrient‐rich Nile silts for agriculture, a practice already recorded in the early 19th century C.E. However, others situated within ancient palaeochannel belts indicate construction over much longer time frames. In this paper, we map the distribution of these river groynes using remote sensing and drone survey. We then establish their probable functions and a provisional chronology using ethnoarchaeological investigation and the ground survey, excavation and radiometric dating of the structures in northern Sudan, focusing on the Holocene riverine landscape surrounding the pharaonic settlement of Amara West (c. 1300–1000 B.C.E.). Finally, we consider the historical and economic implications of this form of hydraulic engineering in the Nile Valley over the past three millennia.

Magmatic continental rifts show evidence that discrete rift segments experience episodic intrusive and eruptive events, more commonly termed rifting episodes. However, whether multiple rifting episodes across adjacent rift segments are clustered in time is not well understood. To address this issue, we conduct new radiocarbon dating that constrains the timing of the most recent rifting episode at the Boset magmatic segment of the northern Ethiopian rift, and combine this with historical dating of similar rifting events in the adjacent magmatic segments. New radiocarbon dates of multiple charcoal samples from the base of the most recent fissural lava at the Boset Volcanic Complex indicate that it likely occurred between 1812 and 1919 CE. These dates are similar to those from historical accounts of fissural eruption from the neighbouring Kone (∼1810 CE), and Fantale (∼1770 to 1808 CE) magmatic segments. We conduct new analysis of major and trace element compositions from these historical fissural lavas, as well as from a fresh-looking lava flow from Beru cone near to Kone volcano. The results of the geochemistry from these flows of all three magmatic segments show compositions that vary in the basalt and trachybasalt fields, with sufficient variation to rule out them having erupted from a single dike intrusion episode. This, combined with the scatter in dates from the radiocarbon analysis and historical accounts, along with the location of each eruption in a discrete and spatially offset magmatic segment, favours an interpretation of each magmatic segment experiencing separate rifting episodes but with these being clustered in time. Mechanisms to explain the clustering of rifting episodes are more speculative but could include stress transfer from dike intrusion and deep crustal hydraulic connection in the plumbing system of multiple segments.

Funding: This research was undertaken, in part, thanks to funding from the Canada Excellence Research Chairs Programme and has been financially supported by the Villum Investigator Project IceFlow (Grant No. 16572 to D.D.-J.). Radar development was further supported by funding from the University of Alabama. EGRIP is directed and organised by the Centre for Ice and Climate at the Niels Bohr Institute, University of Copenhagen. S.F. received funding from the German Academic Exchange Service (DAAD): Forschungsstipendien für promovierte Nachwuchswissenschaftlerinnen und -wissenschaftler. M.R.E. was supported by a DFG Emmy Noether grant (grant no. DR 822/3-1). The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation. Publisher PDF

The release of phosphorus (P) from crustal rocks during weathering plays a key role in determining the size of Earth’s biosphere, yet the concentration of P in crustal rocks over time remains controversial. Here, we combine spatial, temporal, and chemical measurements of preserved rocks to reconstruct the lithological and chemical evolution of Earth’s continental crust. We identify a threefold increase in average crustal P concentrations across the Neoproterozoic-Phanerozoic boundary (600 to 400 million years), showing that preferential biomass burial on shelves acted to progressively concentrate P within continental crust. Rapid compositional change was made possible by massive removal of ancient P-poor rock and deposition of young P-rich sediment during an episode of enhanced global erosion. Subsequent weathering of newly P-rich crust led to increased riverine P fluxes to the ocean. Our results suggest that global erosion coupled to sedimentary P-enrichment forged a markedly nutrient-rich crust at the dawn of the Phanerozoic.

The classic work on the morphology of limestone calcite cements done in the 1960s is extended here by utilising growth zones to reconstruct the growth of cement crystals. Only cement composed of fitted polyhedral monocrystals that form by passive crystallisation of calcite on the walls of liquid‐filled, static pores and fissures is considered. Cement can either be initiated by (1) nucleation, when new crystals start but are not attached to their substrate, or (2) seeding, when new crystals are seamlessly connected to and influenced by substrate crystals. After seeding, epitaxial cement growth starts with many sub‐crystals that coalesce distally, followed by layered mantle growth. Junctions between three intercrystalline boundaries in cement aggregates with one interfacial angel = 180° are of two types: the first, enfacial junctions are caused by a pause in the growth of one crystal and the second is caused by movement of all boundaries due to dissolution of adjacent calcite. Growth zone offsetting at some intercrystalline boundaries is caused by dissolution of calcite at boundaries when permeability values are low. The same width to height ratio of mature aggregate crystals is predicted from the shape of the crystal's growth surfaces; dogtooth calcite forms columnar and nail‐head calcite forms tabular‐shaped crystals. Seeding on different sized crystals causes variations in epitaxial growth rate with faster growth on large crystals resulting in a disorganised cement fabric; the variation in epitaxial growth rate is perpetuated into mantle growth. Echinoderm syntaxial crystals dominate many pore cements due to the large size of their seed ossicles, at the same time, syntaxial crystals form on relatively tiny seeds. Texturally mature crystal aggregates with isopachous fabric are initiated from three different substrate to cement arrangements. Calcite cement zones preserve their original positions allowing the investigation of cement's growth and chemical history.

We present an age model for the 651 m deep ice core from Skytrain Ice Rise, situated inland of the Ronne Ice Shelf, Antarctica. The top 2000 years have previously been dated using age markers interpolated through annual layer counting. Below this, we align the Skytrain core to the AICC2012 age model using tie points in the ice and air phase, and we apply the Paleochrono program to obtain the best fit to the tie points and glaciological constraints. In the gas phase, ties are made using methane and, in critical sections, δ18Oair; in the ice phase ties are through 10Be across the Laschamps event and through ice chemistry related to long-range dust transport and deposition. This strategy provides a good outcome to about 108 ka (∼ 605 m). Beyond that there are signs of flow disturbance, with a section of ice probably repeated. Nonetheless values of CH4 and δ18Oair confirm that part of the last interglacial (LIG), from about 117–126 ka (617–627 m), is present and in chronological order. Below this there are clear signs of stratigraphic disturbance, with rapid oscillation of values in both the ice and gas phase at the base of the LIG section, below 628 m. Based on methane values, the warmest part of the LIG and the coldest part of the penultimate glacial are missing from our record. Ice below 631 m appears to be of age > 150 ka., Climate of the Past, 19 (4), ISSN:1814-9324, ISSN:1814-9332

Acknowledgements: We sincerely thank the professional support of the captain and crew of the R/V Maria S. Merian as well as the scientific team on the expedition MSM12/2. We also thank W. Luttmer for technical assistance with the measurement at AWI. Many thanks to G. Mollenhauer, the AWI MICADAS lab, and S. Morton (SUERC) for radiocarbon measurement. Thanks to C. Vogt for XRD analysis. This project was supported by the Deutsche Forschungsgemeinschaft (DFG) through the International Research Training Group IRTG 1904 ArcTrain. Grant NE/I020261/1 from NERC in the ocean acidification program to D.N.S. is also gratefully acknowledged. Furthermore, we acknowledge support by the Open Access Publication Funds of AWI., Funder: Deutsche Forschungsgemeinschaft through "ArcTrain" (GRK1904), Freshwater perturbations are often thought to be associated with abrupt climate changes during the last deglaciation, while many uncertainties remain regarding the exact timing, pathway, mechanism, and influence of meltwater release. Here, we present very well-dated and high-resolution records from the eastern Labrador Sea representing the last 19.000 years, which demonstrate abrupt changes in sea surface characteristics. Four millennial-scale meltwater events have been identified between the last 14.000 and 8.200 years based on independent biomarker proxies and X-ray fluorescence scanning data. These events are characterized by increased sea ice formation and decreased sea surface temperatures which might have occurred within a few decades. We propose these abrupt changes were triggered by meltwater pulsing into the Labrador Sea periodically, resulting from collapse of the Laurentide-Greenland Ice Sheets caused by (sub-)surface ocean warming in the Labrador Sea. Our findings provide more precise information about impact of freshwater forcing on abrupt climate changes, which may help to improve simulations for past and future changes in ocean circulation and climate.

The lavas associated with mantle plumes may sample domains throughout Earth’s mantle and probe its dynamics. However, plume studies are often only able to take snapshots in time, usually of the most recent plume activity, leaving the chemical and geodynamic evolution of major convective upwellings in Earth’s mantle poorly constrained. Here, we report the geodynamically key information of how the lithology and density of a plume change from plume head phase to tail. We use iron stable isotopes and thermodynamic modeling to show that the Galápagos plume has contained small, nearly constant, amounts of dense recycled crust over its 90-million-year history. Despite a temporal evolution in the amount of recycled crust-derived melt in Galápagos-related lavas, we show that this can be explained by plume cooling alone, without associated changes in the plume’s mantle source; results are also consistent with a plume rooted in a lower mantle low-velocity zone also sampling primordial components.

Funder: Centre for Effective Altruism; Grant(s): Long-Term Future Fund Grant (no number) (December 2021, PI: N Wilson), Island nations may have potential long-term survival value for humanity in global catastrophes such as sun-blocking catastrophes from nuclear winter and large magnitude volcanic eruptions. One way to explore this issue further is to understand the impact on islands after the largest historically observed volcanic eruption: that of Mt Tambora in 1815. For each of the 31 large, populated islands selected, we conducted literature searches for relevant historical and palaeoclimate studies. We also analysed results from a reconstruction (EKF400v2), which uses atmospheric-only general circulation model simulations with assimilated observational and proxy data. From the literature review, there was widespread evidence for weather/climate anomalies in 1815-1817 for these islands (29/29 for those with data). But missing data was an issue for other dimensions such as impaired food production (seen in 8 islands out of only 12 with data). Based on the EKF400v2 reconstruction for temperature anomalies (compared to the relatively "non-volcanic" reference period of 1779 to 1808), the islands had lower temperature anomalies in the 1815-1818 period than latitudinally equivalent continental sites (at 100 km and 1000 km inland). This was statistically significant for the great majority of the comparisons for group analyses by hemisphere, oceans, and temperate/tropical zone. When considering just the islands, all but four showed statistically anomalous temperature reductions in the 1816-1817 period (for most p < 0.00001). In the peak impact year of 1816, the lowest anomalies were seen for islands in the Southern Hemisphere (p < 0.0001), the Indian Ocean (p < 0.0001), and in the tropics and subtropics of the Southern Hemisphere (p = 0.0057). In conclusion, the findings of both the literature review and reconstruction simulations suggest climatic impacts of the Tambora eruption for nearly all these 31 large islands, albeit less than for continental sites. Islands with the smallest temperature anomalies were in the Southern Hemisphere, in particular the Indian Ocean and the tropics and subtropics of the Southern Hemisphere.

Funder: Observation and Research Station of Salt Lakes in Tibetan Plateau, Hydromagnesite (Mg5(CO3)4(OH)2 4H2O) is a common hydrated magnesium carbonate mineral found in alkaline lakes on Earth, potentially present on Mars, and is also a key mineral for carbon capture and storage. However, the mechanisms governing its formation in alkaline lakes remain enigmatic. Extensive hydromagnesite formed during the Holocene in the alkaline Dujiali Lake, central Qinghai-Tibetan Plateau, making it an ideal field site to constrain the process of hydromagnesite formation in a modern environmental context. In this study, we report a set of magnesium isotope ratios (26Mg/24Mg expressed as δ26Mg) data from Dujiali Lake on abiotic hydromagnesite (mean=−1.35 ± 0.14‰) modern lake waters (−0.07 to +0.46‰), rivers and groundwater (−0.53 to −1.46‰). These differences in δ26Mg (and also Mg/Ca) are most likely caused by low-Mg carbonate precipitation, a process which fractionates Mg/Ca and δ26Mg. A semi-quantitative box model of the lake chemistry was developed based on carbonate equilibria to investigate the behavior of Mg isotopes during the evolution of the lake chemistry. The modeling results indicate that evaporation concentrates solutes in the lake driving saturation of multiple minerals. Aragonite reaches saturation before hydromagnesite, preferentially removing Ca relative to Mg via aragonite precipitation. This process elevates the Mg/Ca of the lake and the saturation index of hydromagnesite increasing the likelihood of hydromagnesite formation. Given that the Mg/Ca ratio of many alkaline lakes is far below that required for the formation of hydromagnesite, our findings suggest that low-Mg carbonate precipitation may be a common precursor process for abiotic hydromagnesite precipitation in evaporative environments in addition to a high Mg source likely derived from ultramafic rocks.

The alteration of olivine-rich rocks to serpentine minerals, (hydr)oxides, and H2(aq) through serpentinization is long thought to have influenced the distribution of habitable environments on early Mars and the evolution of the early Martian hydrosphere and atmosphere. Nevertheless, the planetary importance of Martian serpentinization has remained a matter of debate; H2 production rates and mechanisms are poorly understood for serpentinization of Fe-rich Martian olivines, and orbital surveys suggest that serpentines are rarely present in the Martian crust. To constrain the process and products of Martian serpentinization, we acquired petrographic, geochemical, mineralogical, and synchrotron-based spectroscopic data on serpentinized Fe-rich olivines from the 1.1 Ga Duluth Complex (Minnesota, USA). These data indicate that serpentinized Fe-rich olivine would have been accompanied by a five-fold increase in H2 production relative to serpentinized terrestrial mantle peridotites. In contrast to previous expectations of mineral products, a vacancy-coupled Fe3+ substitution for Mg2+ in the serpentine octahedral site yields hisingerite (with minor magnetite) as the dominant Fe serpentine mineral at comparatively low temperature and pH, consistent with meteorite mineralogy and in-situ rover data. The widespread occurrence of Fe3+-phyllosilicates in highly magnetized regions of the Martian crust supports the hypothesis that serpentinization was more pervasive on early Mars than currently estimated. Our results provide the Mars 2020 Perseverance Rover mission with new and specific mineralogical and geochemical signatures of serpentinization and H2 production as it approaches the largest olivine-bearing lithological unit identified on the surface of Mars., Leverhulme Trust

The classic Nihewan Fauna, as the representative of early Pleistocene Asian Land Mammal Age Nihewanian, has long been cited for bio- and chronostratigraphic correlation. However, its precise provenance and stratigraphic horizon have remained unsolved till now. The authors successfully extracted the vital information by rediscovering 30 of the original excavation localities. The fossils were catalogued with their provenance through an extensive field survey and comprehensive inspection of the Tianjin Natural History Museum collection. A review of the original description of these fossil localities, using satellite images, and subsequent lithological examination of the Xiashagou strata in the field verified the new findings. The survey produced the first stratigraphic profile calibrated with fossil horizons of the Nihewan formation in the Xiashagou section. Correlated with the published magnetostratigraphical profile of the section, an age of similar to 2.4-1.8 Ma is estimated for the classic Nihewan Fauna.

Sequential analyses of δ13C, δ18O and Δ47 values of calcite and dolomite deposited in millimetre‐sized cavities are reported from the Ronaldsway Member packstones, Isle of Man. The Ronaldsway brachiopods have δ13C values of ca +2.3‰ and δ18O values of ca −7.2‰; carbon is like predicted Carboniferous values, while oxygen values are more negative. The brachiopods show preserved microstructure but have marginal alteration and a streaky cathodoluminescence pattern. Crinoid ossicles have δ13C values of ca +2.3‰ and one with a δ18O value of ca −3.1‰, compatible with Carboniferous marine precipitates; three samples have δ18O values of ca −6.5‰ and are 18O‐depleted. Calcite stages 1 and 2 have δ13C values ca +3.2‰ and δ18O values ca −2.5‰, compatible with Carboniferous sea water. Stage 1 and 2 have non‐luminescent to orange CL zones. Stage 1 and early stage 2 contain red luminescent dolomite micro crystals generated during Mg calcite stabilisation. The Δ47 values for stage 1 and 2 cements indicate temperatures of 86 and 105°C that occurred after the stabilisation of Mg calcite. Stage 3–8 zoned cements preserve their original growth surfaces and their δ13C and δ18O values suggest precipitation during burial and exhumation. The Δ47 values of the brachiopods and crinoids indicate temperatures between 85 and 140°C indicating they were either recrystallised at high temperatures or affected by solid state reordering. To evaluate these alternatives two quantitative models, water–rock reaction and reordering models are performed. The allochems and cements are progressively altered by porewater towards the fluid‐buffered behaviour. The quantitative evaluation of calcite and dolomite solid‐state reordering suggests the elevated clumped isotopic temperatures are produced by interaction with hydrothermal fluids. This study improves understanding by applying previously untried techniques; further Δ47 data and quantifying elemental variations would help further interpretation but the poorly documented post‐depositional history is a drawback.

Ice rises are areas of locally grounded, slow-moving ice adjacent to floating ice shelves. Temperature profiles measured through ice rises contain information regarding changes to their dynamic evolution and external forcings, such as past surface temperatures, past accumulation rates and geothermal heat flux. While previous work has used borehole temperature–depth measurements to infer one or two such parameters, there has been no systematic investigation of parameter sensitivity to the interplay of multiple external forcings and dynamic changes. A one-dimensional vertical heat flow forward model developed here examines how changing forcings affect temperature profiles. Further, using both synthetic data and previous measurements from the Crary Ice Rise in Antarctica, we use our model in a Markov chain Monte Carlo inversion to demonstrate that this method has potential as a useful dating technique that can be implemented at ice rises across Antarctica. However, we also highlight the non-uniqueness of previous ice rise formation dating based on temperature profiles, showing that using nominal values for forcing parameters, without taking into account their realistic uncertainties, can lead to underestimation of dating uncertainty. In particular, geothermal heat flux represents the dominant source of uncertainty in ice rise age estimation. For instance, in Crary Ice Rise higher heat flux values (i.e. about 90 mW m−2) yield grounding timing of 1400 ± 800 years, whereas lower heat flux of around 60 mW m−2 implies earlier ice rise formation and lower uncertainties in the ice rise age estimations (500 ± 250 years). We discuss the utility of this method in choosing future ice drilling sites and conclude that integrating this technique with other indirect dating methods can provide useful constraints on past forcings and changing boundary conditions from in situ temperature–depth measurements.

Crystal clusters are common in both extrusive and plutonic rocks, but the mechanisms by which they form are not well-constrained. Following a consideration of the physics of nucleation, we outline the expected microstructural characteristics of clusters formed by heterogeneous nucleation and those formed by synneusis, together with the ways they might evolve during subsequent grain growth and textural equilibration. By combining analysis of the microstructures in experimental chromite-basalt charges with a detailed microstructural analysis of the UG2 chromitite of the Bushveld layered intrusion using EBSD, we argue that the UG2 chromitite formed by settling and accumulation of single grains and clusters comprising randomly oriented grains produced by the aggregation of previously isolated chromite crystals. Although there is no evidence of epitaxy, at least some of the lowermost chromite grains of the main UG2 chromitite may have nucleated heterogeneously on the silicate grains forming the floor, with subsequent accumulation and sintering of individual grains or clusters. The reduced thickness of chromitites on the steep and overhanging parts of the floor is thus due to the relative difficulty of sticking more grains to the existing layer in these orientations. The absence of any fining-upwards of grains in either the main UG2 chromitite or the associated stringer can be accounted for if both layers were formed by the settling and accumulation of clusters as well as single grains. Comparison with examples of clustered chromite grains in extrusive rocks suggests that aggregation by synneusis is a widespread magmatic process. The ‘chicken-wire’ texture formed by clustered chromite grains commonly found in olivine-rich cumulates is argued to also be formed by gravitational settling, with the possible exception of clusters of chromite grains in relatively thin seams argued to be the result of metasomatism, which may instead have formed by impingement during in situ growth.

Climate during the last glacial period was marked by abrupt instability on millennial timescales that included large swings of temperature in and around Greenland (Daansgard–Oeschger events) and smaller, more gradual changes in Antarctica (AIM events). Less is known about the existence and nature of similar variability during older glacial periods, especially during the early Pleistocene when glacial cycles were dominantly occurring at 41 kyr intervals compared to the much longer and deeper glaciations of the more recent period. Here, we report a continuous millennially resolved record of stable isotopes of planktic and benthic foraminifera at IODP Site U1385 (the “Shackleton Site”) from the southwestern Iberian margin for the last 1.5 million years, which includes the Middle Pleistocene Transition (MPT). Our results demonstrate that millennial climate variability (MCV) was a persistent feature of glacial climate, both before and after the MPT. Prior to 1.2 Ma in the early Pleistocene, the amplitude of MCV was modulated by the 41 kyr obliquity cycle and increased when axial tilt dropped below 23.5∘ and benthic δ18O exceeded ∼3.8 ‰ (corrected to Uvigerina), indicating a threshold response to orbital forcing. Afterwards, MCV became focused mainly on the transitions into and out of glacial states (i.e. inceptions and terminations) and during times of intermediate ice volume. After 1.2 Ma, obliquity continued to play a role in modulating the amplitude of MCV, especially during times of glacial inceptions, which are always associated with declining obliquity. A non-linear role for obliquity is also indicated by the appearance of multiples (82, 123 kyr) and combination tones (28 kyr) of the 41 kyr cycle. Near the end of the MPT (∼0.65 Ma), obliquity modulation of MCV amplitude wanes as quasi-periodic 100 kyr and precession power increase, coinciding with the growth of oversized ice sheets on North America and the appearance of Heinrich layers in North Atlantic sediments. Whereas the planktic δ18O of Site U1385 shows a strong resemblance to Greenland temperature and atmospheric methane (i.e. Northern Hemisphere climate), millennial changes in benthic δ18O closely follow the temperature history of Antarctica for the past 800 kyr. The phasing of millennial planktic and benthic δ18O variation is similar to that observed for MIS 3 throughout much of the record, which has been suggested to mimic the signature of the bipolar seesaw – i.e. an interhemispheric asymmetry between the timing of cooling in Antarctica and warming in Greenland. The Iberian margin isotopic record suggests that bipolar asymmetry was a robust feature of interhemispheric glacial climate variations for at least the past 1.5 Ma despite changing glacial boundary conditions. A strong correlation exists between millennial increases in planktic δ18O (cooling) and decreases in benthic δ13C, indicating that millennial variations in North Atlantic surface temperature are mirrored by changes in deep-water circulation and remineralization of carbon in the abyssal ocean. We find strong evidence that climate variability on millennial and orbital scales is coupled across different timescales and interacts in both directions, which may be important for linking internal climate dynamics and external astronomical forcing.

Funding: This research has been supported by the Leverhulme Trust (grant RPG-2015-246), by a Royal Society Professorship (grant no. RP/R/180003), and by a Marie Curie Career Integration Grant (CIG14-631752). Volcanic eruptions are the dominant cause of natural variability in climate forcing on timescales up to multidecadal. Large volcanic eruptions lead to global-scale climate effects and influence the carbon cycle on long timescales. However, estimating the frequency of eruptions is challenging. Here we assess the frequency at which eruptions with particular deposition fluxes are observed in the EPICA Dome C ice core over the last 200 kyr. Using S isotope analysis we confirm that most of the largest peaks recorded at Dome C are from stratospheric eruptions. The cumulative frequency through 200 kyr is close to linear, suggesting an approximately constant rate of eruptions. There is no evidence for an increase in the rate of events recorded in Antarctica at either of the last two deglaciations. Millennial variability is at the level expected from recording small numbers of eruptions, while multimillennial variability may be partly due to changes in transport efficiency through the Brewer–Dobson circulation. Our record of events with sulfate deposition rates > 20 and >50 mg m−2 contains 678 and 75 eruptions, respectively, over the last 200 kyr. Calibration with data on historic eruptions and analysis of a global Quaternary dataset of terrestrial eruptions indicates that sulfate peaks with deposition rates > 20 and >50 mg m−2 correspond to explosive eruptions of magnitude ≥ 6.5 and ≥7, respectively. The largest recorded eruption deposited just over 300 mg m−2. Publisher PDF

Vent opening hazard models are routinely used as inputs for assessing distal volcanic hazards (lava flows, tephra fallout) in distributed volcanic fields. These vent opening hazard models have traditionally relied on the location of mapped vents; seldom have they taken into account how vents are linked in space and time. We show that inputs needed to appropriately model distal hazards are fundamentally different than thoses required to model near-vent hazards (ground deformation). We provide a computational model to obtain more appropriate eruptive source parameters (ESPs) for distal volcanic hazard sources and show the utility of our code through three examples. The code's strength is that it links events based on the spatio-temporal relationships of vents through heirarchical clustering. The development of the code and its strenghts and weaknesses are discussed. This work challenges previous ideas about ESPs and we hope this work leads to further improvement in hazard assessment methods.

Funder: Isaac Newton Trust in Cambridge, The transition of a passive continental margin into a subduction zone remains a hypothesis because few geological cases have been reported. The North Sulawesi subduction zone is a 5–9 myr system in Southeast Asia that has evolved from a passive continental margin and has long been overlooked by studies of passive to active margin transitions. Here we compare geophysical evidence from the region with our numerical simulation results. We find that the initial subduction of North Sulawesi relies on horizontal stress, where the trench retreat depends on the negative buoyancy of the oceanic lithosphere. Furthermore, less space available for subduction leads to reduced mantle flow caused by subduction and slower trench retreat. These new dynamical constraints indicate that the negative buoyancy of the oceanic plate is a key factor for the trench retreat, even though subduction initiated was induced.

A potential record of Earth’s magnetic field going back 4.2 billion years (Ga) ago is carried by magnetite inclusions in zircon grains from the Jack Hills. This magnetite may be secondary in nature, however, meaning that the magnetic record is much younger than the zircon crystallization age. Here, we use atom probe tomography to show that Pb-bearing nanoclusters in magnetite-bearing Jack Hills zircons formed during two discrete events at 3.4 and

Funding: This study was supported by urgency funding from the U.K.’s Science and Technology Facilities Council (STFC) as part of the project “Curation and Preliminary Examination of the Winchcombe Carbonaceous Chondrite Fall.” Additional work was funded by STFC through grants ST/N000846/1, ST/T002328/1, ST/T506096/1, and ST/W001128/1 (to L.D., M.R.L., and L.J.Ha.); ST/V000675/1 (to K.H.J. and R.H.J.); ST/P005225/1 (to R.T.); ST/S000348/1 (to M.V.G., P.J.W., and M.J.B.); ST/R00143X/1 (to J.B. and L.J.Hi.); ST/S000615/1 (to G.S.C.); ST/V000799/1 (to P.G.); and ST/V000888/1 (to T.E.). A.J.K. and H.C.B. acknowledge funding support from UK Research and Innovation (UKRI) grant MR/T020261/1. P.L.C. acknowledges funding support from UKRI grant MR/S03465X/1. K.H.J. acknowledges funding support from the Royal Society, grant URF\R\201009. L.J.Ha. and M.R.L. acknowledge funding from Natural Environment Research Council (NERC) National Environmental Isotope Facility (NEIF) grant no. 2406.0321. L.D., M.R.L., and L.J.Ha. acknowledge COVID-19 funding support from the University of Glasgow, UK. D.V. was supported in part by NASA cooperative agreement 80NSSC21M0073. P.P. and I.Sy. acknowledge funding from the VEGA agency, project no.1/0421/20. A.J.T.J. acknowledges support from the European Union and the State of Hungary, cofinanced by the European Regional Development Fund in the project of GINOP-2.3.2-15-2016-00009 “ICER.” P.M.S. acknowledges support from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement no. 945298. FRIPON was initiated by funding from ANR (grant N.13-BS05-0009-03), carried out by the Paris Observatory, Muséum National d’Histoire Naturelle, Paris-Saclay University, and Institut Pythéas (LAM-CEREGE). FRIPON data are hosted and processed at Institut Pythéas SIP (Service Informatique Pythéas). The Desert Fireball Network team and Global Fireball Observatory are funded by the Australian Research Council (DP200102073). Direct links between carbonaceous chondrites and their parent bodies in the solar system are rare. The Winchcombe meteorite is the most accurately recorded carbonaceous chondrite fall. Its pre-atmospheric orbit and cosmic-ray exposure age confirm that it arrived on Earth shortly after ejection from a primitive asteroid. Recovered only hours after falling, the composition of the Winchcombe meteorite is largely unmodified by the terrestrial environment. It contains abundant hydrated silicates formed during fluid-rock reactions, and carbon- and nitrogen-bearing organic matter including soluble protein amino acids. The near-pristine hydrogen isotopic composition of the Winchcombe meteorite is comparable to the terrestrial hydrosphere, providing further evidence that volatile-rich carbonaceous asteroids played an important role in the origin of Earth's water. Publisher PDF

Ultra-low velocity zones (ULVZs) are thin anomalous patches on the boundary between the Earth's core and mantle, revealed by their effects on the seismic waves that propagate through them. Here we map a broad ULVZ near the Galápagos hotspot using shear-diffracted waves. Forward modelling assuming a cylindrical shape shows the patch is ~600 km wide, ~20 km high, and its shear velocities are ~25% reduced. The ULVZ is comparable to other broad ULVZs mapped on the core-mantle boundary near Hawaii, Iceland, and Samoa. Strikingly, all four hotspots where the mantle plume appears rooted by these ‘mega-ULVZs’, show similar anomalous isotopic signatures in He, Ne, and W in their ocean island basalts. This correlation suggests mega-ULVZs might be primordial or caused by interaction with the core, and some material from ULVZs is entrained within the plume. For the Galápagos, the connection implies the plume is offset to the west towards the base of the mantle.

Detailed microstructural analysis of three basaltic sills of the Little Minch Sill Complex demonstrates that convection leaves a detectable signature in fully solidified bodies. The presence of dense clusters of equant grains of olivine and clinopyroxene in the central parts of sills can only be accounted for if they formed and were enlarged while suspended in convecting magma, with delayed settling to the sill floor. An associated stratigraphic invariance of plagioclase grain shape is consistent with growth while suspended in convecting magma. These microstructural indicators demonstrate that convection during solidification was vigorous and long-lived in the 135-m-thick picrodolerite-crinanite unit (PCU) of the composite Shiant Isles Main sill and vigorous and likely short-lived in the PCU of the composite Creagan Iar sill. In contrast, convection in the Meall Tuath sill was weak and short-lived: plagioclase grain shape in this sill varies with stratigraphic height, indicative of primarily in situ nucleation and growth at the magma-mush interface, while olivine and clinopyroxene were kept suspended in the overlying convecting magma. The magma in all three sills fractionated during solidification, permitting convection driven by the instability of an upper thermal boundary layer. The comparative vigour and longevity of convection in the Shiant Isles Main sill and the Creagan Iar sill was due to their emplacement above an earlier, still-hot, intrusion, resulting in highly asymmetric cooling.

The geological record of landforms and sediments produced beneath deglaciating ice sheets offers insights into inaccessible glacial processes. Large subglacial valleys formed by meltwater erosion of sediments (tunnel valleys) are widespread in formerly glaciated regions such as the North Sea. Obtaining a better understanding of these features may help with the parameterisation of basal melt rates and the interplay between basal hydrology and ice dynamics in numerical models of past, present, and future ice-sheet configurations. However, the mechanisms and timescales over which tunnel valleys form remain poorly constrained. Here, we present a series of numerical modelling experiments, informed by new observations from high-resolution 3D seismic data (6.25 m bin size, ∼4 m vertical resolution), which test different hypotheses of tunnel valley formation and calculate subglacial water routing, seasonal water discharges, and the rates at which tunnel valleys are eroded beneath deglaciating ice sheets. Networks of smaller or abandoned channels, pervasive slump deposits, and subglacial landforms are imaged inside and at the base of larger tunnel valleys, indicating that these tunnel valleys were carved through the action of migrating smaller channels within tens of kilometres of the ice margin and were later widened by ice-contact erosion. Our model results imply that the drainage of extensive surface meltwater to the ice-sheet bed is the dominant mechanism responsible for tunnel valley formation; this process can drive rapid incision of networks of regularly spaced subglacial tunnel valleys beneath the fringes of retreating ice sheets within hundreds to thousands of years during deglaciation. Combined, our observations and modelling results identify how tunnel valleys form beneath deglaciating mid-latitude ice sheets and have implications for how the subglacial hydrological systems of contemporary ice sheets may respond to sustained climate warming.

Recent thermodynamic and experimental studies have suggested that volatile organic compounds (e.g., methane, formate, and acetate) can be produced and stabilized in subduction zones, potentially playing an important role in the deep carbon cycle. However, field evidence for the high-pressure production and storage of solid organic compounds is missing. Here, we examine forearc serpentinite clasts recovered by drilling mud volcanoes above the Mariana subduction zone. Notable correlations between carbon and iron stable-isotope signatures and fluid-mobile element (B, As and Sb) concentrations provide evidence for the percolation of slab-derived CO 2 -rich aqueous fluids through the forearc mantle. The presence of carbonaceous matter rich in aliphatic moieties within high-temperature clasts (>350°C) demonstrates that molecular hydrogen production associated with forearc serpentinization is an efficient mechanism for the reduction and conversion of slab-derived CO 2 -rich fluids into solid organic compounds. These findings emphasize the need to consider the forearc mantle as an important reservoir of organic carbon on Earth.

The mid-Palaeozoic was a significant interval in the evolution of plants, during which their biogeomorphic influence on terrestrial landscapes increased dramatically. The products of plant-sediment interactions are well known from Carboniferous strata, but earlier evidence dating to the initial rise of arborescent trees in the Devonian is less common. We present new evidence for early plant-sediment interactions from the Late Devonian (Famennian) Harrylock Formation (County Wexford, Ireland), which hosts standing trees that represent Ireland's earliest known fossil forest. The formation records deposition in fluvial and lacustrine environments, the former of which is here shown to have hosted Earth's earliest known log jam as well as early examples of vegetation-induced sedimentary structures (VISS) and a vegetation-stabilized chute channel. Fossilized driftwood preserved in the lacustrine facies contains the earliest evidence for arthropod(?) borings in large vascular plant debris. Together these early examples show that plant-sediment and plant-animal interactions, frequently recorded in Carboniferous strata, were already in existence by the Devonian and may have been previously underreported in other successions of the same age. The sparse Devonian record of plant-related sedimentary phenomena could be explained by a lesser abundance and/or longevity of such features in landscapes prior to subsequent evolutionary innovations in large vascular plants.

The crystal orientation fabric (COF) of ice sheets records the past history of ice sheet deformation and influences present-day ice flow dynamics. Though not widely implemented, coherent ice-penetrating radar is able to detect bulk anisotropic fabric patterns by exploiting the birefringence of ice crystals at radar frequencies, with the assumption that one of the crystallographic axes is aligned in the vertical direction. In this study, we conduct a suite of quad-polarimetric measurements consisting of four orthogonal antenna orientation combinations near the Western Antarctic Ice Sheet (WAIS) Divide ice core site. From these measurements, we are able to quantify the azimuthal fabric asymmetry at this site to a depth of 1400 m at a bulk-averaged resolution of up to 15 m. Our estimates of fabric asymmetry closely match corresponding fabric estimates directly measured from the WAIS Divide ice core. While ice core studies are often unable to determine the absolute fabric orientation due to core rotation during extraction, we are able to identify and conclude that the fabric orientation is depth-invariant to at least 1400 m, equivalent to 6700 years BP (years before 1950) and aligns closely with the modern surface strain direction at WAIS Divide. Our results support the claim that the deformation regime at WAIS Divide has not changed substantially through the majority of the Holocene. Rapid polarimetric determination of bulk fabric asymmetry and orientation compares well with much more laborious sample-based COF measurements from thin ice sections. Because it is the bulk-averaged fabric that ultimately influences ice flow, polarimetric radar methods provide an opportunity for its accurate and widespread mapping and its incorporation into ice flow models.

Sediments deposited by glacial meltwaters (for example, ice‐contact delta deposits) form permeable packages in the subsurface that can act as reservoirs for both water and hydrocarbons. They are also an important source of aggregate for the construction industry. As reservoirs they are challenging to characterize in terms of their structure, flow and storage properties due to their complex depositional history. In this study, ice‐contact deltas of Salpausselkä I and II end moraines in Southern Finland are studied using a combination of geomorphological mapping, sedimentological studies and near surface geophysical methods. Sedimentary logs from isolated outcrops were correlated to ground penetrating radar (GPR) profiles to unravel the internal structure and depositional history of these ice‐contact deltas. Subsequently, electrical resistivity tomography (ERT) and gravity data were analysed to estimate the depth to bedrock and to model porosity distribution within the sediments. Results of the study suggest that the delta deposits have a broad range of porosities (10 to 42%) with lowest values found in the bottomset beds. The most variable porosities are in the subaqueous ice‐contact–fan zone, and consistently high porosities occur in delta foreset/topset facies. Detailed sedimentary logging linked to the GPR data shows heterogeneities such as mud drapes on foresets and kettle holes which are below the resolution of ERT and gravity methods but significantly affect reservoir properties of the deltas. Moreover, oscillation of the ice‐margin may have introduced larger heterogeneities (for example, buried ice marginal ridges, or eskers) into the sedimentary sequence which are atypical for other Gilbert‐type deltas. Finally, subglacially sculpted, highly variable bedrock topography exerts a major control on sediment distribution within the delta making reservoir volume and quality less predictable. This work has implications for present‐day freshwater aquifers and low enthalpy geothermal energy in southern Finland and other deglaciated regions, as well as hydrocarbon exploration of analogous deposits in the subsurface from Pleistocene and pre‐Pleistocene glaciogenic sequences.

Volcanic eruptions can cause significant disruption to society, and numerical models are crucial for forecasting the dispersion of erupted material. Here we assess the skill and limitations of the Met Office's Numerical Atmospheric-dispersion Modelling Environment (NAME) in simulating the dispersion of the sulfur dioxide (SO2) cloud from the 21–22 June 2019 eruption of the Raikoke volcano (48.3∘ N, 153.2∘ E). The eruption emitted around 1.5±0.2 Tg of SO2, which represents the largest volcanic emission of SO2 into the stratosphere since the 2011 Nabro eruption. We simulate the temporal evolution of the volcanic SO2 cloud across the Northern Hemisphere (NH) and compare our model simulations to high-resolution SO2 measurements from the TROPOspheric Monitoring Instrument (TROPOMI) and the Infrared Atmospheric Sounding Interferometer (IASI) satellite SO2 products. We show that NAME accurately simulates the observed location and horizontal extent of the SO2 cloud during the first 2–3 weeks after the eruption but is unable, in its standard configuration, to capture the extent and precise location of the highest magnitude vertical column density (VCD) regions within the observed volcanic cloud. Using the structure–amplitude–location (SAL) score and the fractional skill score (FSS) as metrics for model skill, NAME shows skill in simulating the horizontal extent of the cloud for 12–17 d after the eruption where VCDs of SO2 (in Dobson units, DU) are above 1 DU. For SO2 VCDs above 20 DU, which are predominantly observed as small-scale features within the SO2 cloud, the model shows skill on the order of 2–4 d only. The lower skill for these high-SO2-VCD regions is partly explained by the model-simulated SO2 cloud in NAME being too diffuse compared to TROPOMI retrievals. Reducing the standard horizontal diffusion parameters used in NAME by a factor of 4 results in a slightly increased model skill during the first 5 d of the simulation, but on longer timescales the simulated SO2 cloud remains too diffuse when compared to TROPOMI measurements. The skill of NAME to simulate high SO2 VCDs and the temporal evolution of the NH-mean SO2 mass burden is dominated by the fraction of SO2 mass emitted into the lower stratosphere, which is uncertain for the 2019 Raikoke eruption. When emitting 0.9–1.1 Tg of SO2 into the lower stratosphere (11–18 km) and 0.4–0.7 Tg into the upper troposphere (8–11 km), the NAME simulations show a similar peak in SO2 mass burden to that derived from TROPOMI (1.4–1.6 Tg of SO2) with an average SO2 e-folding time of 14–15 d in the NH. Our work illustrates how the synergy between high-resolution satellite retrievals and dispersion models can identify potential limitations of dispersion models like NAME, which will ultimately help to improve dispersion modelling efforts of volcanic SO2 clouds.

Iron speciation in ancient sedimentary rocks is widely used to reconstruct oceanic redox conditions over geological time, specifically to assess the extent of oxic, euxinic (anoxic containing sulfide), and ferruginous (anoxic containing iron) conditions. We explore how post-depositional sedimentary processes can skew particular geochemical signals in the rock record. One such process is when aqueous sulfide—including that produced in the sediment column—reacts with sedimentary iron, converting non-sulfide, highly reactive iron minerals to iron-sulfide minerals; this can lead to increased preservation of iron as pyrite and an overestimation of seafloor euxinia. We show that sedimentary rocks with higher (>5 wt%) total iron content are more buffered to this effect and thus are a more reliable indicator of true water-column euxinia. When considering this effect in the geological past, we estimate that true euxinia in the mid-Proterozoic may have been as much as fourfold less than previously thought—more in line with other recent paleoredox proxies not based on iron minerals. Marine iron and sulfate concentrations were more equivalent in Proterozoic–Neoproterozoic oceans, suggesting this time period was particularly susceptible to this post-depositional alteration, explaining the extent of euxinia suggested for this geological interval.

Alkaline magmatism develops in the rear arc area of Java, Sunda Arc, at different range of slab depth; ~270-580 km in central section and ~150 km in the east. We collate published geochemical data of volcanic rocks from four alkaline volcanoes (Muria, Lasem, Bawean, and Ringgit-Beser) and perform statistical analysis to evaluate geochemical characteristics of each suite. A set of major and trace elements is scaled and transformed using principal component analysis (PCA) and then followed by implementation of k-means algorithm to cluster the data points based on Euclidian distances. K-means clustering of the dataset suggests that Central Java alkalines are most elevated in K2O and total alkali. The algorithm further suggests that Muria samples can be clustered into two, owing to these components. These two clusters, however, are not well reflected on trace element-based clustering. Lasem volcanics show distinct cluster high in Na2O/K2O and SiO2, while Bawean samples are mixed into both Muria clusters. Ringgit-Beser alkalines show two distinct clusters tied to MgO and enrichment in Ba, Rb, and Sr. Our findings suggest that the potassium and LILE enrichment in these alkaline rocks is independent of slab depth and is most likely regulated by tectonic-related arc segmentation in Java subduction zone.

Funder: Department of Geography, University of Cambridge, Terrestrial archives from the Mediterranean have been crucial to expanding our under-standing of past environmental variability on a range of timescales, however, dating these sequences is often challenging. Tephra deposits can provide crucial age control for detailed environmental reconstructions on sub-centennial timescales. Here, tephra analysis is undertaken for the first time on a sediment core (I-08) from Lake Ioannina, northwest Greece, for the interval spanning 46 to 4 ka BP. Detailed visible and ‘crypto-’tephra analysis identifies deposits associated with explosive volcanism at Italian volcanic sources, including Campi Flegrei, Pantelleria, and the Aeolian Islands. We identify two visible tephra layers, the Campanian Ignimbrite (CI/Y-5; ca. 39.8 ka BP) and Pantelleria Green Tuff (PGT/Y-6; ca. 45.7 ka), as well as the Holocene Vallone del Gabellotto cryptotephra marker (E-1; ca. 8.3 ka BP). Evidence for repeated remobilisation and redeposition of CI tephra material is outlined, and the potential mechanisms and effects of sediment reworking in lake environments are examined. Bayesian modelling, which incorporates the new tephra ages with earlier radiocarbon dates, extends the I-08 core chronology back to ca. 46 ka BP, facilitating direct correlation of the Ioannina sequence to others in the Mediterranean region.

The international ice core community has a target to obtain continuous ice cores stretching back as far as 1.5 Myr. This would provide vital data (including a CO2 profile) allowing us to assess ideas about the cause of the Mid-Pleistocene Transition (MPT). The European Beyond EPICA project and the Australian Million Year Ice Core project each plan to drill such a core in the region known as Little Dome C. Dating the cores will be challenging, and one approach will be to match some of the records obtained with existing marine sediment datasets, informed by similarities in the existing 800 kyr period. Water isotopes in Antarctica have been shown to closely mirror deepwater temperature, estimated from Mg/Ca ratios of benthic foraminifera, in a marine core on the Chatham Rise near to New Zealand. The dust record in ice cores resembles very closely a South Atlantic marine record of iron accumulation rate. By assuming these relationships continue beyond 800 ka, our ice core record could be synchronised to dated marine sediments. This could be supplemented, and allow synchronisation at higher resolution, by the identification of rapid millennial-scale events that are observed both in Antarctic methane records and in emerging records of planktic oxygen isotopes and alkenone sea surface temperature (SST) from the Portuguese Margin. Although published data remain quite sparse, it should also be possible to match 10Be from ice cores to records of geomagnetic palaeo-intensity and authigenic 10Be/9Be in marine sediments. However, there are a number of issues that have to be resolved before the ice core 10Be record can be used. The approach of matching records to a template will be most successful if the new core is in stratigraphic order but should also provide constraints on disordered records if used in combination with absolute radiogenic ages.

Indonesia hosts the largest number of active volcanoes, several of which are renowned for climate-changing historical eruptions. This pedigree might suggest a substantial fraction of global volcanic sulfur emissions from Indonesia and are intrinsically driven by sulfur-rich magmas. However, a paucity of observations has hampered evaluation of these points—many volcanoes have hitherto not been subject to emissions measurements. Here we report new gas measurements from Indonesian volcanoes. The combined SO2 output amounts to 1.15 ± 0.48 Tg/yr. We estimate an additional time-averaged SO2 yield of 0.12-0.54 Tg/yr for explosive eruptions, indicating a total SO2 inventory of 1.27-1.69 Tg/yr for Indonesian. This is comparatively modest—individual volcanoes such as Etna have sustained higher fluxes. To understand this paradox, we compare the geodynamic, petrologic, magma dynamical and shallow magmatic-hydrothermal processes that influence the sulfur transfer to the atmosphere. Results reinforce the idea that sulfur-rich eruptions reflect long-term accumulation of volatiles in the reservoirs.

The evolution of land plants during the Palaeozoic Era transformed Earth’s biosphere 1. Because the Earth's surface and interior are linked by tectonic processes, the linked evolution of the biosphere and sedimentary rocks should be recorded as a near-contemporary shift in the composition of the continental crust. To test this hypothesis, we assessed the isotopic signatures of zircon formed at subduction zones where marine sediments are transported into the mantle 2,3, thereby recording interactions between surface environments and the deep Earth. Using oxygen and lutetium-hafnium isotopes of magmatic zircon that respectively track surface weathering (time-independent) 4 and radiogenic decay (time-dependent) 5, we find a correlation in the composition of continental crust after 430 Myr ago, which is coeval with the onset of enhanced complexity and stability in sedimentary systems related to the evolution of vascular plants. The expansion of terrestrial vegetation brought channelled sand-bed and meandering rivers, muddy floodplains, and thicker soils, lengthening the duration of weathering before final marine deposition 6,7. Collectively, our results suggest that the evolution of vascular plants coupled the degree of weathering and timescales of sediment routing to depositional basins where they were subsequently subducted and melted. The late Palaeozoic isotopic shift of zircon indicates that the greening of the continents was recorded in the deep Earth.

Coastal wetlands are dynamic bio-physical systems in which vegetation affects the movement of water and sediment, which in turn build and maintain the landform and ecosystem. Wetlands are an effective buffer against coastal erosion and flooding, enhance water quality and human health and wellbeing. Numerous field and laboratory experiments have quantified the reduction of waves by coastal ecosystems. Numerical models, however, are only able to capture observed reduction in wave energy when calibration coefficients are obtained by comparison with measured dissipation rates. A deeper understanding of how wave attenuation varies over time, with local flow conditions and ecosystem properties, is still lacking and should be acquired from a greater range of ecosystem types and geographical settings. Few studies have observed the detailed seasonal variations in how coastal wetlands function as wave buffers and how such seasonal variations might be explained. Equally, few studies have focused on the effect of coastal reed beds on wave dynamics. This study addresses both: i) seasonal variability in wave dissipation through reed vegetation and ii) intricate connections between reed vegetation and the physical context (meteorological and topographical) that might explain such variability. We present observations of wind generated wave transformation through two Phragmites australis reed beds in the Razelm-Sinoe Lagoon System, Danube Delta, Romania. We find that seasonal changes in vegetation density and biomass, as well as meteorological conditions, affect observed wave conditions within the first few meters of the reed beds. Our results also show a preferential reduction of higher frequency waves, irrespective of reed stem diameter or density and suggest the potential importance of seasonal vegetation debris to observed wave dissipation. Such complex and non-linear biogeomorphic effects on wave dissipation are not currently well understood or captured in the parameterisation of vegetation-induced wave dissipation. Our study highlights the importance of an accurate and temporally granular quantification of nearshore bathymetry, wetland topography, and vegetation to fully understand, model, and manage bio-physical interactions in coastal wetlands. More specifically, our results point towards the need for spatially and temporally explicit wave decay functions in emergent reed vegetation. This is particularly critical where the accurate evaluation of the flood and erosion risk contribution of any wetland is required as part of nature-based coastal protection solutions.

SUMMARY Anisotropy of remanent magnetization and magnetic susceptibility are highly sensitive and important indicators of geological processes which are largely controlled by mineralogical parameters of the ferrimagnetic fraction in rocks. To provide new physical insight into the complex interaction between magnetization structure, shape, and crystallographic relations, we here analyse ‘slice-and-view’ focused-ion-beam (FIB) nano-tomography data with micromagnetic modelling and single crystal hysteresis measurements. The data sets consist of 68 magnetite inclusions in orthopyroxene (Mg60) and 234 magnetite inclusions in plagioclase (An63) were obtained on mineral separates from the Rustenburg Layered Suite of the Bushveld Intrusive Complex, South Africa. Electron backscatter diffraction was used to determine the orientation of the magnetite inclusions relative to the crystallographic directions of their silicate hosts. Hysteresis loops were calculated using the finite-element micromagnetics code MERRILL for each particle in 20 equidistributed field directions and compared with corresponding hysteresis loops measured using a vibrating sample magnetometer (VSM) on silicate mineral separates from the same samples. In plagioclase the ratio of remanent magnetization to saturation magnetization (Mrs/Ms) for both model and measurement agree within 1.0 per cent, whereas the coercivity (Hc) of the average modelled curve is 20 mT lower than the measured value of 60 mT indicating the presence of additional sources of high coercivity in the bulk sample. The VSM hysteresis measurements of the orthopyroxene were dominated by multidomain (MD) magnetite, whereas the FIB location was chosen to avoid MD particles and thus contains only particles with diameters

The thermochemical structure of lithospheric and asthenospheric mantle exert primary controls on surface topography and volcanic activity. Volcanic rock compositions and mantle seismic velocities provide indirect observations of this structure. Here, we compile and analyze a global database of the distribution and composition of Neogene-Quaternary intraplate volcanic rocks. By integrating this database with seismic tomographic models, we show that intraplate volcanism is concentrated in regions characterized by slow upper mantle shear-wave velocities and by thin lithosphere (i.e., Not applicable

Arthropleura is a genus of giant myriapods that ranged from the early Carboniferous to Early Permian, with some individuals attaining lengths of over 2 metres. While most known fossils of the genus are disarticulated and occur primarily in late Carboniferous (Pennsylvanian) strata, here we report partially articulated Arthropleura remains from the early Carboniferous Stainmore Formation (Serpukhovian; Pendleian) in the Northumberland Basin of northern England. This 76 x 36 cm specimen represents part of an exuvium and is notable because only two comparably articulated giant Arthropleura fossils are previously known. It represents one of the biggest known arthropod fossils and the largest arthropleurid recovered to date, the earliest (Mississippian) body fossil evidence for gigantism in Arthropleura, and the first instance of a giant arthropleurid body fossil within the same regional sedimentary succession as the large arthropod trackway, Diplichnites cuithensis. The remains represent 12-14 anterior Arthropleura tergites, in the form of a partially sand-filled dorsal exoskeleton. The original organism is estimated to have been 55 cm in width and up to 2.63 m in length, weighing c. 50 kg. The specimen is preserved partially in three dimensions within fine sandstone and has been moderately deformed by synsedimentary tectonics. Despite imperfect preservation, the specimen corroborates the hypothesis that Arthropleura had a tough, sclerotized exoskeleton. Sedimentological evidence for a lower delta plain depositional environment supports the contention that Arthropleura preferentially occupied open woody habitats rather than swampy environments, and that it shared such habitats with tetrapods. When viewed in the context of all other global evidence for Arthropleura, the specimen contributes to a dataset that shows the genus had an equatorially restricted palaeogeographic range, achieved gigantism prior to late Palaeozoic peaks in atmospheric oxygen, and was relatively unaffected by climatic events in the late Carboniferous, prior to its extinction in the early Permian.

The year 2021 marked the bicentenary of the birth of James Croll (1821–1890), the self‐educated son of a crofter‐stonemason, whose life was characterised by a dizzying range of occupations and homes, poor health and financial concerns, and yet he became a pioneer of orbital dynamics and ice age climate change with an impressive record of publication. Drawing upon archival information and recently published observations, this paper explores selected aspects of Croll's biography, his scientific connections and controversies, and that area of his life relevant to Quaternary science. He was a 19th century polymath whose multifaceted contributions have been a catalyst for subsequent systems‐based climate science on the grand scale, including the foundations for the seminal work of Milutin Milankovitch on the rhythms of Quaternary environmental change.

Besides providing a unique archaeological assemblage that documents the early emergence of complex behaviour in the human lineage, Border Cave (South Africa) is noteworthy for having yielded hominin remains of at least nine individuals, including the partial cranium Border Cave 1. While the exact provenance of Border Cave 1 is unknown, sequence stratigraphy and ESR dating converge towards an age from about 82 ka to 170 ka. Here we present novel information about the brain, braincase and bony labyrinth of Border Cave 1 and discuss related evolutionary implications. We compare Border Cave 1 to specimens of Early and Middle Pleistocene Homo as well as to fossil and extant Homo sapiens. Virtual segmentation techniques were used to reconstruct the brain and bony labyrinth endocasts, assess the distribution of cranial bone thickness, and identify the vascular and sulcal imprints preserved on the inner surface of the braincase. Our results show that the overall morphology of the brain endocast approximates the globular shape of the modern human brain and differs from the long and low brains seen in Middle Pleistocene fossil hominins. The vascular imprints preserved on the right hemisphere indicate that the middle branch derives from the anterior branch, which is a pattern shared with Neanderthals and modern humans. Bone thickness distribution in the Border Cave 1 cranium resembles the patterns seen in Cro-Magnon 1 and Abri Pataud 1, which both share a diffuse distribution of thickened areas over the frontal region. Finally, the relative size and curvature of the semicircular canals of the bony labyrinth conform to the ancestral configuration shared between Early and Late Pleistocene fossil hominins from Africa and the Levant, as well as modern humans, and distinct from the more derived condition documented within Neanderthals. We discuss the implications of our findings for understanding the biogeography, evolution, and, to some extent, behaviour of fossil Homo sapiens.

During the mid-Pliocene (Zanclean, ca. ∼ 3.9 Ma), parts of the Canadian High Arctic experienced mean annual temperatures that were 14–22°C warmer than today and supported diverse boreal-type forests. The landscapes of this vegetated polar region left behind a fragmented sedimentary record that crops out across several islands in the Canadian Arctic Archipelago as the Beaufort Formation and correlative strata. Paleoecological information from these strata provides a high-fidelity window onto Pliocene environments, and prominent fossil sites yield unparalleled insights into Cenozoic mammal evolution. Significantly, many of the strata reveal evidence for life-sediment interactions in a warm-climate Arctic, most notably in the form of extensive woody debris and phytoclast deposits. This paper presents original field data that refines the sedimentological context of plant debris accumulations from the anactualistic High Arctic forests, most notably at the ‘Fyles Leaf Beds' and ‘Beaver Pond' fossil-bearing sites in the ‘high terrace deposits' of central Ellesmere Island. The former is a remarkably well-preserved, leaf-rich deposit that is part of a complex of facies associations representing lacustrine, fluvio-deltaic and mire deposition above a paleotopographic unconformity. The latter yields tooth-marked woody debris within a peat layer that also contains a rich assemblage of vertebrate and plant fossils including abundant remains from the extinct beaver-group Dipoides. Here we present sedimentological data that provide circumstantial evidence that the woody debris deposit at Beaver Pond could record dam-building in the genus, by comparing the facies motif with new data from known Holocene beaver dam facies in England. Across the Pliocene of the High Arctic region, woody debris accumulations are shown to represent an array of biosedimentary deposits and landforms including mires, driftcretions, woody bedforms, and possible beaver dams, which help to contextualize mammal fossil sites, provide facies models for high-latitude forests, and reveal interactions between life and sedimentation in a vanished world that may be an analogue to that of the near-future.

The Rawnsley Quartzite of South Australia hosts some of the world's most diverse Ediacaran macrofossil assemblages, with many of the constituent taxa interpreted as early representatives of metazoan clades. Globally, a link has been recognized between the taxonomic composition of individual Ediacaran bedding-plane assemblages and specific sedimentary facies. Thorough characterization of fossil-bearing facies is thus of fundamental importance for reconstructing the precise environments and ecosystems in which early animals thrived and radiated, and distinguishing between environmental and evolutionary controls on taxon distribution. This study refines the paleoenvironmental interpretations of the Rawnsley Quartzite (Ediacara Member and upper Rawnsley Quartzite). Our analysis suggests that previously inferred water depths for fossil-bearing facies are overestimations. In the central regions of the outcrop belt, rather than shelf and submarine canyon environments below maximum (storm-weather) wave base, and offshore environments between effective (fair-weather) and maximum wave base, the succession is interpreted to reflect the vertical superposition and lateral juxtaposition of unfossiliferous non-marine environments with fossil-bearing coastal and shoreface settings. Facies comprise: 1, 2) amalgamated channelized and cross-bedded sandstone (major and minor tidally influenced river and estuarine channels, respectively), 3) ripple cross-laminated heterolithic sandstone (intertidal mixed-flat), 4) silty-sandstone (possible lagoon), 5) planar-stratified sandstone (lower shoreface), 6) oscillation-ripple facies (middle shoreface), 7) multi-directed trough- and planar-cross-stratified sandstone (upper shoreface), 8) ripple cross-laminated, planar-stratified rippled sandstone (foreshore), 9) adhered sandstone (backshore), and 10) planar-stratified and cross-stratified sandstone with ripple cross-lamination (distributary channels). Surface trace fossils in the foreshore facies represent the earliest known evidence of mobile organisms in intermittently emergent environments. All facies containing fossils of the Ediacaran macrobiota remain definitively marine. Our revised shoreface and coastal framework creates greater overlap between this classic “White Sea” biotic assemblage and those of younger, relatively depauperate “Nama”-type biotic assemblages located in Namibia. Such overlap lends support to the possibility that the apparent biotic turnover between these assemblages may reflect a genuine evolutionary signal, rather than the environmental exclusion of particular taxa.

Carbonate cave deposits (speleothems) have been used widely for paleoclimate reconstructions; however, few studies have examined the utility of other speleothem-forming minerals for this purpose. Here we demonstrate for the first time that stable isotopes (δ17O, δ18O and δD) of structurally-bound gypsum (CaSO4·2H2O) hydration water (GHW) can be used to infer paleoclimate. Specifically, we used a 63 cm-long gypsum stalactite from Sima Blanca Cave to reconstruct the climate history of SE Spain from ~ 800 BCE to ~ 800 CE. The gypsum stalactite indicates wet conditions in the cave and humid climate from ~ 200 BCE to 100 CE, at the time of the Roman Empire apogee in Hispania. From ~ 100 CE to ~ 600 CE, evaporation in the cave increased in response to regional aridification that peaked at ~ 500–600 CE, roughly coinciding with the transition between the Iberian Roman Humid Period and the Migration Period. Our record agrees with most Mediterranean and Iberian paleoclimate archives, demonstrating that stable isotopes of GHW in subaerial gypsum speleothems are a useful tool for paleoclimate reconstructions.

The Last Interglacial (LIG), a warmer period 130,000–116,000 years before present, is a potential analogue for future climate change. Stronger LIG summertime insolation at high northern latitudes drove Arctic land summer temperatures 4–5 °C higher than in the pre-industrial era. Climate model simulations have previously failed to capture these elevated temperatures, possibly because they were unable to correctly capture LIG sea-ice changes. Here, we show that the latest version of the fully coupled UK Hadley Center climate model (HadGEM3) simulates a more accurate Arctic LIG climate, including elevated temperatures. Improved model physics, including a sophisticated sea-ice melt-pond scheme, result in a complete simulated loss of Arctic sea ice in summer during the LIG, which has yet to be simulated in past generations of models. This ice-free Arctic yields a compelling solution to the long-standing puzzle of what drove LIG Arctic warmth and supports a fast retreat of future Arctic summer sea ice. Arctic climate in the Last Interglacial (LIG)—a warm period 130,000–116,000 years ago—is poorly simulated by modern climate models. A model with improved sea-ice melt-pond physics reproduces LIG Arctic temperatures, suggests an ice-free Arctic during this period and predicts the same by 2035.