Paleomagnetism of Middle Proterozoic mafic intrusions and Upper Proterozoic (Nankoweap) red beds from the Lower Grand Canyon Supergroup, Arizona

Paleomagnetic data from lavas and dikes of the Unkar igneous suite (16 sites) and sedimentary rocks of the Nankoweap Formation (7 sites), Grand Canyon Supergroup (GCSG), Arizona, provide two primary paleomagnetic poles for Laurentia for the latest Middle Proterozoic (ca. 1090 Ma) at 32°N, 185°E (dp=6.8°, DM=9.3°) and early Late Proterozoic (ca. 850–900 Ma) at 10°S, 163°E (dp=3.5°, DM=7.0°). A new ⁴⁰Ar/³⁹Ar age determination from an Unkar dike gives an interpreted intrusion age of about 1090 Ma, similar to previously reported geochronologic data for the Cardenas Basalts and associated intrusions. The paleomagnetic data show no evidence of any younger, middle Late Proterozoic tectonothermal event such as has been revealed in previous geochronologic studies of the Unkar igneous suite. The pole position for the Unkar Group Cardenas Basalts and related intrusions is in good agreement with other ca. 1100 Ma paleomagnetic poles from the Keweenawan midcontinent rift deposits and other SW Laurentia diabase intrusions. The close agreement in age and position of the Unkar intrusion (UI) pole with poles derived from rift related rocks from elsewhere in Laurentia indicates that mafic magmatism was essentially synchronous and widespread throughout Laurentia at ca. 1100 Ma, suggesting a large-scale continental magmatic event. The pole position for the Nankoweap Formation, which plots south of the Unkar mafic rocks, is consistent with a younger age of deposition, at about 900 to 850 Ma, than had previously been proposed. Consequently, the inferred 200 Ma difference in age between the Cardenas Basalts and overlying Nankoweap Formation provides evidence for a third major unconformity within the Grand Canyon sequence.     

Crustal thickness, discontinuity depth, and upper mantle structure beneath southern Africa: constraints from body wave conversions

The technique of receiver function analysis is applied to the study of crustal and upper mantle structures beneath the Kaapvaal craton in southern Africa and its surroundings. Seismic data were recorded by the seismic array of 82 sites deployed from April 1997 to April 1999 across southern Africa, as well as a dense array of 32 sites near Kimberley, in operation from December 1998 to June 1999. Arrival times for phases converted at the Moho are used to determine crustal thickness. The Moho depth in the south–western section of the craton was found to vary between 37 and 40 km, except for one station that recorded a depth of 43 km (SA23). Farther north along the western block of the craton (into Botswana) the depth increases up to 43 km. The depth increases even further in the north–eastern section of the craton, where results vary from 40 to 52 km. Just north of the Kaapvaal craton, in the neighbouring Zimbabwe craton, the crustal thickness drops significantly. The results obtained there varied from 36 to 40 km. For the Kimberley area, using the dense array, the Moho depth was found to be 37.3 km. Arrivals of the Ps and Ppps phases were used to determine the Poisson’s ratio in the region. This was found to be 0.26±0.01. Arrivals of phases from the 410 and 660 km mantle discontinuities are used to interpret the relative positions of these discontinuities, as well as for comparison of mantle temperatures and seismic velocities in the region with global averages. In the Kimberley area the 410 and 660 km discontinuities were found at their expected depth, implying that mantle temperatures in the region are close to the global average. The seismic velocities above the ‘410’ were found up to 5% faster than the averages from the global iasp91 model, which is fast even by Precambrian standards. In other sections of the Kaapvaal craton, the velocities are also faster than global averages, but not as fast as beneath Kimberley. In these sections, the ‘410’ is also slightly elevated, while the ‘660’ is depressed, which implies a slightly lower mantle temperature relative to the global average. Beneath the Kaapvaal craton we find evidence suggesting the presence of a zone with a reduced wavespeed gradient at an upper bound of approximately 300 km, which may mark the lower chemical boundary of the craton.     

Reproducibility of river water quality measurements: inter-agency comparisons for quality assurance

We assessed the reproducibility of river state-of-environment (SoE) water quality measurements in the Wellington Region, New Zealand (NZ). Field staff from GWRC and NIWA conducted 29 side-by-side water sampling and in-situ measurements at six river sites of diverse water quality for 12 variables measured routinely in river SoE monitoring across NZ. Field measurements of water temperature, dissolved oxygen, electrical conductivity and visual clarity agreed closely with strong numerical similarity (within 10%). Numerical similarity ranged widely for laboratory measurements, from strong for nitrate-nitrite-nitrogen to weak for turbidity, dissolved reactive phosphorus, and ammoniacal-nitrogen. Numerical agreement was very weak for laboratory pH (which is problematic) and E. coli–which is ‘tolerable’ for many applications given good correlation (R?=?0.94) over a 2000-fold concentration range. The findings of our inter-agency comparison have contributed to quality assurance recommendations in the NZ National Environmental Monitoring Standard (NEMS) for water quality.     

Fault seal prediction of seismic-scale normal faults in porous sandstone: A case study from the eastern Gulf of Suez rift,Egypt

A study of normal faults in the Nubian Sandstone Sequence, from the eastern Gulf of Suez rift, has been conducted to investigate the relationship between the microstructure and petrophysical properties of cataclasites developed along seismic-scale faults (slip-surface cataclasites) and smaller displacement faults (deformation bands) found in their damage zones. The results help to quantify the uncertainty associated with predicting the fluid flow behaviour of seismic-scale faults by analysing small faults recovered from core, a common procedure in the petroleum industry. The microstructure of the cataclasites was analysed as well as their single-phase permeability and threshold pressure. Faulting occurred at a maximum burial depth of ∼1.2 km. The permeability of deformation band and slip-surface cataclasites varies over ∼1.5 orders of magnitude for a given fault. Our results suggest that the lowest measured deformation band permeabilities provide a good estimate for the arithmetic-mean permeability of the major slip-surface cataclasites. This is because the cataclastic permeability reduction is mostly established early in the deformation history. Stress at the time of faulting rather than final strain appears to be the critical factor determining fault rock permeability. For viable predictions it is important that the slip-surface cataclasites and deformation bands originate from the same host. On the other hand, a higher uncertainty is associated with threshold pressure predictions, as the arithmetic-mean slip-surface cataclasite threshold pressure exceeds the highest measured deformation band threshold pressure by at least a factor of 4.     

Intertidal spring dynamics and coastal nutrient loading revealed through geophysics,drone surveys,and in-situ monitoring

Coastal eutrophication poses an increasing risk to ecosystem health due to enhanced nutrient loading to the global coastline. Submarine groundwater discharge (SGD) represents a significant pathway for nitrate-nitrogen (NO₃-N) transport to the coast, but diffusive SGD transport is difficult to monitor directly, given the low flux rates and expansive discharge areas. In contrast, focused SGD from intertidal springs can potentially be sampled and directly gauged, providing unique insight into SGD and associated contaminant transport. Basin Head is a coastal lagoon in Prince Edward Island, Canada that is a federally protected ecosystem. Nitrate-nitrogen is conveyed from agricultural fields in the contributing watershed to the eutrophic lagoon via intertidal groundwater springs and groundwater-dominated tributaries. We used several field methods to characterize groundwater discharge, nutrient loading, and in-channel mixing associated with intertidal springs. The tributaries and intertidal springs were gauged and sampled to estimate a representative summer nitrate load to the lagoon. Our analysis revealed that NO₃-N export to the lagoon through tributaries and springs throughout summer 2023 was on average 401 kg N/month, with the combined spring loading comparable in magnitude to the combined tributary loading. We collected thermal infrared and visual imagery using drone surveys and found spatial overlap between cold-water plumes from the spring discharge and macroalgae blooms, indicating the local thermal and ecosystem impacts of the focused SGD. We also mapped the electrical resistivity (salinity) distribution in the water column around one large spring with electromagnetic geophysics at different tidal stages to reveal the three-dimensional spring plume dynamics. Results showed that the fresher spring water floated above the saline lagoon water with the brackish plume oriented in the direction of the tidal current. Collectively, our multi-pronged field investigations help elucidate the hydrologic, thermal, and nutrient dynamics of intertidal springs and the cascading ecosystem impacts.     

Homo floresiensis and the late Pleistocene environments of eastern Indonesia: defining the nature of the relationship

Evidence from Liang Bua, a limestone cave on the island of Flores in East Indonesia, provides a unique opportunity to explore the long term relationship between hominins and their environment. Occupation deposits at the site span ～95 ka and contain abundant stone artefacts, well preserved faunal remains and evidence for an endemic species of hominin: Homo floresiensis. Work at the site included detailed geomorphological and environmental analysis, which has enabled comparisons to be drawn between changes in the occupational intensity in the cave, using stone tool and faunal counts, and changes in the environmental conditions, using the characteristics of the sedimentary layers in the cave and speleothem records. These comparisons demonstrate that H. floresiensis endured rapidly fluctuating environmental conditions over the last ～100 ka, which influenced the geomorphological processes in the cave and their occupational conditions. The intensity of occupation in the cave changed significantly between 95 and 17 ka, with peaks in occupation occurring at 100–95, 74–61 and 18–17 ka. These correlate with episodes of channel formation and erosion in the cave, which in turn correspond with high rainfall, thick soils and high bio-productivity outside. In contrast, periods of low occupational intensity correlate with reduced channel activity and pooling associated with drier periods from 94 to 75 and 36 to 19 ka. This apparent link between intensity of hominin use of the cave and the general conditions outside relates to the expansion and contraction of the rainforest and the ability of H. floresiensis to adapt to habitat changes. This interpretation implies that these diminutive hominins were able to survive abrupt and prolonged environmental changes by changing their favoured occupation sites. These data provide the basis for a model of human–environment interactions on the island of Flores. With the addition of extra data from other sites on Flores, this model will provide a greater understanding of H. floresiensis as a unique human species.     

Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.