The variation and controls of mean transit times in Australian headwater catchments

Determining mean transit times in headwater catchments is critical for understanding catchment functioning and understanding their responses to changes in landuse or climate. Determining whether mean transit times (MTTs) correlate with drainage density, slope angle, area, or land cover permits a better understanding of the controls on water flow through catchments and allows first-order predictions of MTTs in other catchments to be made. This study assesses whether there are identifiable controls on MTTs determined using ³H in headwater catchments of southeast Australia. Despite MTTs at baseflow varying from a few years to >100 years, it was difficult to predict MTTs using single or groups of readily-measured catchment attributes. The lack of readily-identifiable correlations hampers the prediction of MTTs in adjacent catchments even where these have similar geology, land use, and topography. The long MTTs of the Australian headwater catchments are probably in part due to the catchments having high storage volumes in deeply-weathered regolith, combined with low recharge rates due to high evapotranspiration. However, the difficulty in estimating storage volumes at the catchment scale hampers the use of this parameter to estimate MTTs. The runoff coefficient (the fraction of rainfall exported via the stream) is probably also controlled by evapotranspiration and recharge rates. Correlations between the runoff coefficient and MTTs in individual catchments allow predictions of MTTs in nearby catchments to be made. MTTs are shorter in high rainfall periods as the catchments wet up and shallow water stores are mobilized. Despite the contribution of younger water, the major ion geochemistry in individual catchments commonly does not correlate with MTTs, probably reflecting heterogeneous reactions and varying degrees of evapotranspiration. Documenting MTTs in catchments with high storage volumes and/or low recharge rates elsewhere is important for understanding MTTs in diverse environments.     

Insights into agricultural influences and weathering processes from major ion patterns

Karst areas and their catchments pose a great challenge for protection because fast conduit flow results in low natural attenuation of anthropogenic contaminants. Studies of the hydrochemistry of karst sources and river solutes are an important tool for securing and managing water resources. A study of the geochemical downriver evolution of the Wiesent River and its tributaries, located in a typical karst terrain, revealed unexpected downstream decreases of nitrate with maximum mean values of 30 mg/L at the source to minimum values of 18 mg/L near the river mouth. This trend persisted over the length of the river even though increased agricultural activities are evident in the downstream section of the catchment. This pattern is caused by fertilizer inputs via diffusive and fast conduits flow from karst lithology in the upstream area that may have reached the river's source even from beyond the hydrological catchment boundaries. Further downstream, these influences became diluted by tributary inputs that drain subcatchments dominated by claystone and sandstone lithologies that increased potassium and sulphate concentrations. Our findings indicate that bedrock geology remains the dominant control on the major ion chemistry of the Wiesent River and that agricultural influences are strongest near the headwaters despite increased land use further downstream, due to long‐term storage and accumulation in karst aquifers. This feature may not be unique to the Wiesent River system, as carbonates cover significant portions of the Earth's surface and subsequent work in other river systems could establish whether such patterns are ubiquitous worldwide.     

Calibration of oxygen isotope fractionation and calcite-corundum thermometry in emery at Naxos,Greece

Corundum (Crn), including sapphire, occurs in emery pods surrounded by marble on the island of Naxos, Greece. The emery formed from bauxite deposited in karst that was metamorphosed to 400–700°C at 20–15 Ma. Many of these rocks initially appeared well suited for refractory accessory mineral (RAM) thermometry, which uses oxygen isotope fractionation between a RAM – corundum – and a modally dominant phase with faster diffusion of oxygen – calcite (Cc) – to determine peak metamorphic temperatures. However, previous attempts at oxygen isotope thermometry were confounded by highly variable fractionations (Δ¹⁸O) measured at mm-scale and the uncertain calibration of Δ¹⁸O(Cc-Crn) versus temperature. Secondary ion mass spectrometry (SIMS) permits in situ analysis of δ¹⁸O in corundum and calcite at the 10-μm scale in adjacent grains where textures suggest peak metamorphic equilibrium was attained. SIMS analyses of adjacent mineral pairs in eight rocks yield values of Δ(Cc-Crn) that systematically decrease from 7.2 to 2.9‰ at higher metamorphic grade. Pairing these data with independent temperature estimates from mineral isograds yields an empirical calibration of 1,000 lnα(Cc-Crn) = 2.72 ± 0.3 × 10⁶/T² (T in K). The new fractionations (2.7‰ at 1,000 K) are significantly smaller than those calculated from the modified increment method (6.5‰ at 1,000 K; Zheng, Geochimica et Cosmochimica Acta, 1991, 55:2299–2307; Zheng, Mineral Mag, 1994, 58A:1000–1001), which yield unreasonably high temperatures of 630 to 1,140°C when applied to the new Naxos data. The new calibration of Δ(Cc-Crn) can be combined with published fractionations to calculate A-factors for corundum versus a range of 14 other minerals. These new fractionation factors can be used for thermometry or to constrain the genesis of corundum. A compilation of gem corundum δ¹⁸O values shows that many igneous sapphires, including important deposits of basalt-associated sapphire, are mildly elevated in δ¹⁸O relative to the calculated range in equilibrium with mantle values (4.4–5.7‰) and formed from evolved magmas.     

Titanite dates crystallization: Slow Pb diffusion during super‐solidus re‐equilibration

Titanite can be found in rocks of wide compositional range, is reactive, growing or regrowing during metamorphic and hydrothermal events, and is generally amenable to U–Pb geochronology. Experimental evidence suggest that titanite has a closure temperature for Pb ranging from 550 to 650°C, and thus titanite dates are commonly interpreted as cooling ages. However, this view has been challenged in recent years by evidence from natural titanite which suggests the closure temperature may be significantly higher (up to 800°C). Here, we investigate titanite in an enclave of migmatitic gneiss included within a granite intrusion. The titanite crystals exhibit textural features characteristic of fluid‐mediated mass transfer processes on length scales of <100 µm. These textural features are associated with variation in both Pb concentrations and distinct U–Pb isotopic compositions. Zr‐in‐titanite thermometry indicates that modification of the titanite occurred at temperatures in excess of 840°C, in the presence of a high‐T silicate melt. The Pb concentration gradients preserved in these titanite crystals are used to determine the diffusivity of Pb in titanite under high‐T conditions. We estimate diffusivities ranging from 2 × 10^?22 to 5 × 10^?25 m²/s. These results are significantly lower than experimental data predict yet are consistent with other empirical data on natural titanites, suggesting that Pb diffusivity is similar to that of Sr. Thus our data challenge the wide‐held assumption that U–Pb titanite dates only reflect cooling ages.     

Influence of length effect on embankment slope reliability in 3D

Embankment slopes composed of spatially variable soils have a variety of different failure modes that are affected by the correlation distances of the material properties and the geometry and total length of the slope. This paper examines the reliability of soil slopes for embankments of different length and uses parallel computing to analyse very long embankments (up to 100 times the embankment height) for a clay soil characterised by a spatially varying undrained shear strength. Based on a series of analyses using the 3D random finite element method (RFEM), it is first shown that the reliability of slopes of various length can be efficiently computed by combining simple probability theory with a detailed 3D RFEM analysis of a representative shorter slope of length 10 times the slope height. RFEM predictions of reliability indices for longer slopes are then compared with results obtained using Vanmarcke's (1977a) simplified 3D method and Calle's (1985) extended 2D approach. It is shown that these methods can give significantly different results, depending on the horizontal scale of fluctuation relative to the slope length, with RFEM predicting a lower slope reliability than the Vanmarcke and Calle solutions in all cases. The differences in the solutions are evaluated and attributed to differences in the assumed and computed failure surface geometries.     

Neogene upper‐crustal cooling of the Olympus range (northern Aegean): major role of Hellenic back‐arc extension over propagation of the North Anatolia Fault Zone

The North Anatolian Fault Zone (NAFZ) is one of the most hazardous active faults on Earth, yet its Pliocene space‐time propagation across the north Aegean domain remains poorly constrained. We use low‐temperature multi‐thermochronology and inverse thermal modelling to quantify the cooling history of the upper crust across the Olympus range. This range is located in the footwall of a system of normal faults traditionally interpreted as resulting from superposed Middle–Late Miocene N–S stretching, related to the back‐arc extension of the Hellenic subduction zone, and a Pliocene‐Quaternary transtensional field, attributed to the south‐westward propagation of the NAFZ. We find that accelerated exhumational cooling occurred between 12 and 6 Ma at rates of 15–35 °C Ma^?1 and decreased to <3 °C Ma^?1 by 8–6 Ma. The absence of significant Plio‐Pleistocene cooling across Olympus suggests that crustal exhumation there is driven by late Miocene back‐arc extension, while the impact of the NAFZ remains limited.     

Stratigraphic architecture of Solander Basin records Southern Ocean currents and subduction initiation beneath southwest New Zealand

Solander Basin is characterized by subduction initiation at the Pacific‐Australia plate boundary, where high biological productivity is found at the northern edge of the Antarctic Circumpolar Current. Sedimentary architecture results from tectonic influences on accommodation space, sediment supply and ocean currents (via physiography); and climate influence on ocean currents and biological productivity. We present the first seismic‐stratigraphic analysis of Solander Basin based on high‐fold seismic‐reflection data (voyage MGL1803, SISIE). Solander Trough physiography formed by Eocene rifting, but basinal strata are mostly younger than ca. 17 Ma, when we infer Puysegur Ridge formed and sheltered Solander Basin from bottom currents, and mountain growth onshore increased sediment supply. Initial inversion on the Tauru Fault started at ca. 15 Ma, but reverse faulting from 12 to ca. 8 Ma on both the Tauru and Parara Faults was likely associated with reorganization and formation of the subduction thrust. The new seabed topography forced sediment pathways to become channelized at low points or antecedent gorges. Since 5 Ma, southern Puysegur Ridge and Fiordland mountains spread out towards the east and Solander Anticline grew in response to ongoing subduction and growth of a slab. Solander Basin had high sedimentation rates because (1) it is sheltered from bottom currents by Puysegur Ridge; and (2) it has a mountainous land area that supplies sediment to its northern end. Sedimentary architecture is asymmetric due to the Subtropical Front, which moves pelagic and hemi‐pelagic sediment, including dilute parts of gravity flows, eastward and accretes contourites to the shelf south of Stewart Island. Levees, scours, drifts and ridges of folded sediment characterize western Solander Basin, whereas hemi‐pelagic drape and secondary gravity flows are found east of the meandering axial Solander Channel. The high‐resolution record of climate and tectonics that Solander Basin contains may yield excellent sites for future scientific ocean drilling.     

Reconstruction of the Cenozoic deformation of the Bohai Bay Basin,North China

A well‐constrained plate deformation model may lead to an improved understanding of sedimentary basin formation and the connection between subduction history and over‐riding plate deformation. Building quantitative models of basin kinematics and deformation remains challenging often due to the lack of comprehensive constraints. The Bohai Bay Basin (BBB) is an important manifestation of the destruction of the North China Craton, and records the plate kinematic history of East Asia during the Cenozoic. Although a number of interpretations of the formation of the BBB have been proposed, few quantitative basin reconstruction models have been built to test and refine previous ideas. Here, we developed a quantitative deformation reconstruction of the BBB constrained with balanced cross‐sections and structural, stratigraphic and depositional age data. Our reconstruction suggests that the basin formation process was composed of three main stages: Paleocene‐early Eocene (65–42 Ma) extension initiation, middle Eocene‐early Oligocene (42–32.8 Ma) extension climax and post‐Oligocene (32.8–0 Ma) post‐extensional subsidence. The deformation of the BBB is spatially heterogeneous, and its velocity directions rotated clockwise during the basin formation process. The reconstruction supports the interpretation that the BBB formed via strike‐slip faulting and orthogonal extension and that the basin is classified as a composite extensional‐transtensional basin. We argue that the clockwise rotation of the basin velocity field was driven by the counter‐clockwise rotation in the direction of Pacific Plate subduction. The kinematics of the BBB imply that the Pacific Plate may have been sufficiently coupled to the over‐riding East Asian Plate during the critical period of Pacific Plate reorganization. The new reconstruction provides a quantitative basis for studies of deformation processes not only in the vicinity of the BBB, but also more broadly throughout East Asia.     

Another Year of Record Heat for the Oceans

Changes in ocean heat content(OHC), salinity, and stratification provide critical indicators for changes in Earth’s energy and water cycles. These cycles have been profoundly altered due to the emission of greenhouse gasses and other anthropogenic substances by human activities, driving pervasive changes in Earth’s climate system. In 2022, the world’s oceans, as given by OHC, were again the hottest in the historical record and exceeded the previous 2021 record maximum.According to IAP/CAS data, ...