Conceptual model of fracture-limited sea cliff erosion: Erosion of the seaward tilted flyschs of Socoa,Basque Country,France

Sea cliff morphology and erosion rates are modulated by several factors, including rock control that reflects both lithology and rock structure. Erosion is anticipated to preferentially exploit ‘fractures’, broadly meant as any discontinuity in an otherwise continuous medium, where the rock mass is weakest. Unpicking the direct control of such fractures on the spatial and temporal pattern of erosion remains, however, challenging. To analyse how such fractures control erosion, we monitored the evolution of a 400 m-long stretch of highly structured sedimentary cliffs in Socoa, Basque Country, France. The rock is known as the Socoa flysch formation. This formation combines decimetre-thick turbidites composed of repeat triplets of medium to strong calcareous sandstone, laminated siltstones and argillaceous marls. The sequence plunges at 45° into the sea with a shore-parallel strike. The cliffs are cross-cut by two normal and reverse fault families, with 10–100 m alongshore spacing, with primary and secondary strata-bound fractures perpendicular to the bedding, which combined delimit the cliff rock mass into discrete blocks that are exploited by the erosion process. Erosion, and sometimes plucking, of such beds and blocks on the cliff face was monitored using ground-based structure-from-motion (SfM) photogrammetry, over the course of 5.7 years between 2011 and 2017. To compare with longer time change, cliff-top retreat rate was assessed using SfM-orthorectified archive aerial photographs spanning 1954–2008. We show that the 13,250 m² cliff face released 4500 blocks exceeding 1.45 × 10⁻³ m³, removing a total volume of 170 m³. This equates to an average cliff erosion rate of 3.4 mm/year, which is slightly slower than the 54-year-long local cliff-top retreat (10.8 ± 1.8 mm/year). The vertical distribution of erosion reflects the height of sea water inundation, where the maximum erosion intensity occurs ca. 2 m above high spring-tide water level. Alongshore, the distribution of rockfall scars is concentrated along bed edges bounding cross-cutting faults; the extent of block detachment is controlled by secondary tectonic joints, which may extend through several beds locally sharing similar mechanical strength; and rockfall depth is always a multiple of bed thickness. Over the longer term, we explain block detachment and resultant cliff collapse as a cycle. Erosion nucleates on readily exploitable fractures but elsewhere, the sea only meets defect-free medium-strong to strong rock slabs offering few morphological features for exploitation. Structurally delimited blocks are quarried, and with sufficient time, carve semi-elliptic scars reaching progressively deeper strata to be eroded. Lateral propagation of erosion is directed along mechanical weaknesses in the bedding, and large episodic collapses affect the overhanging slabs via sliding on the weak marl beds. Collapse geometry is confined to one or several triplets of turbidite beds, but never reaches deeper into the cliff than the eroded depth at the foot. We contend that this fracture-limited model of sea-cliff erosion, inferred from the Socoa site dynamics and its peculiar sets of fractures, applies more broadly to other fractured cliff contexts, albeit with site-specific geometries. The initiation of erosion, the propagation of incremental block release and the ultimate full failure of the cliff, have each been shown to be fundamentally directly controlled by structure, which remains a vital control in understanding how cliffed coasts have changed in the past and will change in the future.     

Tree-ring dates on two pre-Little Ice Age advances in Glacier Bay National Park and Preserve, Alaska, USA

Two interstadial tree ring-width chronologies from Geikie Inlet, Glacier Bay Southeast, Alaska were built from 40 logs. One of these chronologies has been calendar dated to AD 224–999 (775 yr) crossdating with a living ring-width chronology from Prince William Sound, Alaska. Trees in this chronology were likely killed through inundation by sediments and meltwater from the advancing Geikie Glacier and its tributaries ca. AD 850. The earlier tree-ring chronology spans 545 yr and is a floating ring-width series tied to radiocarbon ages of about 3000 cal yr BP. This tree-ring work indicates two intervals of glacial expansion by the Geikie Glacier system toward the main trunk glacier in Glacier Bay between 3400 and 3000 cal yr BP and again about AD 850. The timing of both expansions is consistent with patterns of ice advance at tidewater glaciers in other parts of Alaska and British Columbia about the same time, and with a relative sea-level history from just outside Glacier Bay in Icy Strait. This emerging tree-ring dated history builds on previous radiocarbon-based glacial histories and is the first study to use tree-ring dating to assign calendar dates to glacial activity for Glacier Bay.     

A process-based approach to estimate lodgepole pine (Pinus contorta Dougl.) distribution in the Pacific Northwest under climate change

Lodgepole pine (Pinus contorta Dougl.) is a widely distributed species in the Pacific Northwest of North America. The extent that the current distribution of this species may be altered under a changing climate is an important question for managers of wood supply as well as those interested in conservation of subalpine ecosystems. In this paper, we address the question, how much might the current range of the species shift under a changing climate? We first assessed the extent that suboptimal temperature, frost, drought, and humidity deficits affect photosynthesis and growth of the species across the Pacific Northwest with a process-based model (3-PG). We then entered the same set of climatic variables into a decision-tree model, which creates a suite of rules that differentially rank the variables, to provide a basis for predicting presence or absence of the species under current climatic conditions. The derived decision-tree model successfully predicted weighted presence and absence recorded on 12,660 field survey plots with an accuracy of ~70%. The analysis indicated that sites with significant spring frost, summer temperatures averaging <15°C and soils that fully recharged from snowmelt were most likely to support lodgepole pine. Based on these criteria, we projected climatic conditions through the twenty-first century as they might develop without additional efforts to reduce carbon emissions using the Canadian Climate Centre model (CGCM2). In the 30-year period centered around 2020, the area suitable for lodgepole pine in the Pacific Northwest was projected to be reduced only slightly (8%). Thereafter, however, the projected climatic conditions appear to progressively favor other species, so that by the last 30 years of twenty-first century, lodgepole pine could be nearly absent from much of its current range. We conclude that process-based models, because they are highly sensitive to seasonal variation in solar radiation, are well adapted to identify the importance of different climatic variables on photosynthesis and growth. These same variables, once indentified, and run through a decision-tree model, provide a reasonable approach to predict current and future patterns in a species?? distribution.     

Microstructural evolution during experimental albitization of K-rich alkali feldspar

Crystals of K-feldspar (C2/m), in contact with highly concentrated aqueous NaCl solutions at 500°C and 200 MPa, are pseudomorphically replaced by high albite (C[`1]) (C\bar{1}) as a result of an interface-coupled dissolution/reprecipitation process. The reaction occurs at an extremely sharp reaction front (<10 nm) and involves the complete breakdown of the initial framework structure. This results in the release of tetrahedrally incorporated elements such as Fe³⁺ and Ti⁴⁺ and a significant increase in Si/Al disorder across the reaction interface. The evolving microstructure is controlled by crystallographic relations between the phases. This leads to highly anisotropic, sawtooth-shaped intergrowths of albite and initial K-feldspar, resulting in the least structural misfit between the two framework structures. As a result, the newly formed interfaces appear to be semicoherent, and cracks across the reaction fronts even indicate elastic strain. The reaction produces 2 distinctive albite types (albite-1 and albite-2). Both are polycrystalline, with albite-2 showing significantly larger subgrain sizes. This indicates a secondary coarsening step driven by the reduction in interfacial energy within the polycrystalline replacement product. The reaction also produces a highly porous rim. However, the porosity is not evenly distributed resulting in a porous albite-1 and a non-porous albite-2 that mostly surrounds large, euhedral pores. Despite the substantial volume fraction of porosity in albite-1, no significant 3D interconnectivity could be detected, making the presence of a pervasive porosity unlikely. However, the result of coarsening is the continuous modification of the 3D porosity distribution. This could potentially provide a mechanism for fluid transport through the replacement rim until textural and chemical equilibration is achieved.     

Vent fluid chemistry of the Rainbow hydrothermal system (36°N, MAR): Phase equilibria and in situ pH controls on subseafloor alteration processes

The Rainbow hydrothermal field is located at 36°13.8′N-33°54.15′W at 2300 m depth on the western flank of a non-volcanic ridge between the South AMAR and AMAR segments of the Mid-Atlantic Ridge. The hydrothermal field consists of 10-15 active chimneys that emit high-temperature (∼365 °C) fluid. In July 2008, vent fluids were sampled during cruise KNOX18RR, providing a rich dataset that extends in time information on subseafloor chemical and physical processes controlling vent fluid chemistry at Rainbow. Data suggest that the Mg concentration of the hydrothermal end-member is not zero, but rather 1.5-2 mmol/kg. This surprising result may be caused by a combination of factors including moderately low dissolved silica, low pH, and elevated chloride of the hydrothermal fluid. Combining end-member Mg data with analogous data for dissolved Fe, Si, Al, Ca, and H₂, permits calculation of mineral saturation states for minerals thought appropriate for ultramafic-hosted hydrothermal systems at temperatures and pressures in keeping with constraints imposed by field observations. These data indicate that chlorite solid solution, talc, and magnetite achieve saturation in Rainbow vent fluid at a similar pH_(T,P) (400 °C, 500 bar) of approximately 4.95, while higher pH values are indicated for serpentine, suggesting that serpentine may not coexist with the former assemblage at depth at Rainbow. The high Fe/Mg ratio of the Rainbow vent fluid notwithstanding, the mole fraction of clinochlore and chamosite components of chlorite solid solution at depth are predicted to be 0.78 and 0.22, respectively. In situ pH measurements made at Rainbow vents are in good agreement with pH_(T,P) values estimated from mineral solubility calculations, when the in situ pH data are adjusted for temperature and pressure. Calculations further indicate that pH_(T,P) and dissolved H₂ are extremely sensitive to changes in dissolved silica owing to constraints imposed by chlorite solid solution-fluid equilibria. Indeed, the predicted correlation between dissolved silica and H₂ defines a trend that is in good agreement with vent fluid data from Rainbow and other high-temperature ultramafic-hosted hydrothermal systems. We speculate that the moderate concentrations of dissolved silica in vent fluids from these systems result from hydrothermal alteration of plagioclase and olivine in the form of subsurface gabbroic intrusions, which, in turn are variably replaced by chlorite + magnetite + talc ± tremolite, with important implications for pH lowering, dissolved sulfide concentrations, and metal mobility.     

A multi-proxy approach to assessing isolation basin stratigraphy from the Lofoten Islands, Norway

This study takes a comprehensive approach to characterizing the isolation sequence of Heimerdalsvatnet, a coastal lake in the Lofoten Islands, northern Norway. We use established methods and explore new techniques to assess changes in marine influence. Bathymetric and sub-bottom profiles were acquired to examine basin-wide sedimentation and a 5.8 m sediment core spanning the last 7800 cal yr BP was analyzed. We measured magnetic susceptibility, bulk organic matter properties, molecular biomarkers, diatom assemblages, and elemental profiles acquired by scanning X-ray fluorescence. These characteristics of the sediment reflect detailed changes in salinity and water column conditions as the lake was progressively isolated. Three distinct litho/chemo-stratigraphic units represent a restricted marine phase (7800-6500 cal yr BP), a transitional phase characterized by intermittent marine influence (6500-4900 cal yr BP), and complete isolation and freshwater sedimentation (4900 cal yr BP to present). Although there are uncertainties in the estimate of the threshold elevation of the lake, the timing of these phases generally corresponds with previous interpretations of the local relative sea-level history. This record captures sea-level regression following the Tapes transgression and supports the interpretation of a subsequent sea-level stillstand, dated in Heimerdalsvatnet from 6500 to 4900 cal yr BP.     

Titanoceratops ouranos, a giant horned dinosaur from the late Campanian of New Mexico

At the end of the Cretaceous, 65.5 million years ago, the giant ceratopsids Triceratops and Torosaurus dominated North America’s dinosaur fauna. The origins of these giant ceratopsids, the Triceratopsini, are poorly understood. This paper describes Titanoceratops ouranos, a giant ceratopsid from the late Campanian (73-74 Ma) of New Mexico, and the earliest known triceratopsin. The holotype was previously interpreted as an aberrant and exceptionally large specimen of Pentaceratops sternbergi, but the animal does not show the diagnostic features of Pentaceratops. Instead, cladistic analysis shows that Titanoceratops is the sister taxon of a clade formed by Eotriceratops, Triceratops, and Torosaurus. With an estimated mass of 6.5 tons, Titanoceratops is among the largest dinosaurs known from the Campanian of North America, and rivaled Triceratops in size. The recognition of Titanoceratops suggests that giant chasmosaurines evolved once, among the Triceratopsini, and that the group evolved large size five million years earlier than previously thought. The giant horned dinosaurs probably originated in the southern part of the North American continent during the Campanian but only became widespread during the Maastrichtian.     

Insights into the Hadean Earth from experimental studies of zircon

Geologists investigate the evolution of the atmosphere, crust, and mantle through time by direct study of the rock record. However, the Hadean eon (>3.85 Ga) has been traditionally viewed as inaccessible due to the absence of preserved rocks. The discovery of >4.0 Ga detrital zircons from Western Australia in the 1980s — coupled with the development of new micro-analytical capabilities — made possible new avenues of early Earth research. The prevailing view that emerged is that the early Earth may have contained a stable hydrosphere, water-saturated or (near watersaturated) granitic magmas, and volcanic emanations dominated by neutral gas species (e.g., CO₂, H₂O, and SO₂). The Hadean Earth may have been capable of supporting life ～200 Ma after accretion and perhaps earlier. Many of these models are formulated — or have been subsequently supported — by laboratory experiments of zircon. Important petrological variables such as temperature, pressure, oxygen fugacity, and component activities (e.g., SiO₂/TiO₂-activities) can be controlled. These experiments are fundamental for extrapolation to ‘deep time’ because they provide a means to understand primary chemistry preserved in ancient zircons. This review paper specifically focuses on zircon experimental studies (oxygen isotope fractionations, Ti-thermometry, and redox sensitive element incorporation into zircon), which have influenced our view of the very early Earth.     

Tephrochronology,magnetostratigraphy and mammalian faunas of Middle and Early Pleistocene sediments at two sites on the Old Crow River,northern Yukon Territory,Canada

Alluvial and lacustrine sediments exposed beneath late Pleistocene glaciolacustrine silt and clay at two sites along the Old Crow River, northern Yukon Territory, are rich in fossils and contain tephra beds. Surprise Creek tephra (SZt) occurs in the lower part of the alluvial sequence at CRH47 and Little Timber tephra (LTt) is present near the base of the exposure at CRH94. Surprise Creek tephra has a glass fission-track age of 0.17 ± 0.07 Ma and Little Timber tephra is 1.37 ± 0.12 Ma. All sediments at CRH47 have a normal remanent magnetic polarity and those near LTt at CRH94 have a reversed polarity — in agreement with the geomagnetic time scale. Small mammal remains from sediments near LTt support an Early Pleistocene age but the chronology is not so clear at CRH47 because of the large error associated with the SZt age determination. Tephrochronological and paleomagnetic considerations point to an MIS 7 age for the interglacial beds just below SZt at CRH47 and at Chester Bluffs in east-central Alaska, but mammalian fossils recovered from sediments close to SZt suggest a late Irvingtonian age, therefore older than MIS 7. Further studies are needed to resolve this problem.     

Diffusion and solubility of hydrogen and water in periclase

Cylinders of synthetic periclase single crystals were annealed at 0.15–0.5 GPa and 900–1200 °C under water-saturated conditions for 45 min to 72 h. Infrared spectra measured on the quenched products show bands at 3,297 and 3,312 cm^?1 indicating V _OH ^? centers (OH-defect stretching vibrations in a half-compensated cation vacancy) in the MgO structure as a result of proton diffusion into the crystal. For completely equilibrated specimens, the OH-defect concentration, expressed as H₂O equivalent, was calculated to 3.5 wt ppm H₂O at 1,200 °C and 0.5 GPa based on the calibration method of Libowitzky and Rossmann (Am Min 82:1111–1115, 1997). This value was confirmed via Raman spectroscopy, which shows OH-defect-related bands at identical wavenumbers and yields an H₂O equivalent concentration of about 9 wt ppm using the quantification scheme of Thomas et al. (Am Min 93:1550–1557, 2008), revised by Mrosko et al. (Am Mineral 96:1748–1759, 2011). Results of both independent methods give an overall OH-defect concentration range of 3.5–9 (+4.5/?2.6) ppm H₂O. Proton diffusion follows an Arrhenius law with an activation energy E _a = 280 ± 64 kJ mol^?1 and the logarithm of the pre-exponential factor logD_o (m² s^?1) = ?2.4 ± 1.9. IR spectra taken close to the rims of MgO crystals that were exposed to water-saturated conditions at 1,200 °C and 0.5 GPa for 24 h show an additional band at 3,697 cm^?1, which is related to brucite precipitates. This may be explained by diffusion of molecular water into the periclase, and its reaction with the host crystal during quenching. Diffusion of molecular water may be described by logD_H2O (m² s^?1) = ?14.1 ± 0.4 (2σ) at 1,200 °C and 0.5 GPa, which is ~ 2 orders of magnitude slower than proton diffusion at identical P-T conditions.