Geology and Geophysics (G&G)
Permanent URI for this collection
Scientists in the G&G Department seek knowledge of the structure, composition, and dynamics of the earth’s interior, the origin and evolution of the earth’s crust, controls on ocean and climate change on time scales of decades to 100 million years, and processes of mass and energy transfer at the land-sea interface.
Browse
Recent Submissions
1 - 20 of 1848
-
ArticleNeurodiversity: An important axis of diversity in ocean sciences(Oceanography Society, 2023-12-05)Neurodiversity refers to variations in the human brain that affect information processing; it includes conditions, or “neurotypes,” such as autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), dyslexia, dyscalculia, and dyspraxia, among others. Neurodiversity can be conceptualized as significant differences in the ways that individuals process information; such differences may concern written or verbal language, sensory information, body language, or social interactions. These differences have been historically viewed within the medical model of disability, for example, as deficits in ability through a diagnosed condition, often associated with a goal of curing or managing the condition.
-
ArticleDeeper and stronger North Atlantic Gyre during the Last Glacial Maximum(Nature Research, 2024-07-10)Subtropical gyre (STG) depth and strength are controlled by wind stress curl and surface buoyancy forcing1,2. Modern hydrographic data reveal that the STG extends to a depth of about 1 km in the Northwest Atlantic, with its maximum depth defined by the base of the subtropical thermocline. Despite the likelihood of greater wind stress curl and surface buoyancy loss during the Last Glacial Maximum (LGM)3, previous work suggests minimal change in the depth of the glacial STG4. Here we show a sharp glacial water mass boundary between 33° N and 36° N extending down to between 2.0 and 2.5 km—approximately 1 km deeper than today. Our findings arise from benthic foraminiferal δ18O profiles from sediment cores in two depth transects at Cape Hatteras (36–39° N) and Blake Outer Ridge (29–34° N) in the Northwest Atlantic. This result suggests that the STG, including the Gulf Stream, was deeper and stronger during the LGM than at present, which we attribute to increased glacial wind stress curl, as supported by climate model simulations, as well as greater glacial production of denser subtropical mode waters (STMWs). Our data suggest (1) that subtropical waters probably contributed to the geochemical signature of what is conventionally identified as Glacial North Atlantic Intermediate Water (GNAIW)5,6,7 and (2) the STG helped sustain continued buoyancy loss, water mass conversion and northwards meridional heat transport (MHT) in the glacial North Atlantic.
-
ArticleRecent advances in vertical temperature profiler instrumentation and flux estimation methods facilitate groundwater – Surface water exchange studies in environments with strong discharge zones(Elsevier, 2024-06-23)Groundwater fluxes to many surface water systems are spatially heterogeneous with discharge focused into discrete, high-flux zones. Quantifying fluxes in these preferential discharge zones is critical to a range of surface water habitat and water quality processes, but characterization can be difficult due to short-scale spatial and temporal variability. Passive heat-as-a-tracer methods employing vertical temperature profiler (VTP) data can provide the necessary spatial and temporal resolution, but upward fluid flow strongly attenuates the thermal signals used for estimating fluxes. In preferential discharge zones it becomes difficult to measure the signals in the subsurface and the flux parameter can become insensitive in the analysis models, leading to large uncertainties. We use data from a high-flux site of contaminant-loaded groundwater discharge to the Quashnet River on Cape Cod, Massachusetts, USA, to demonstrate how recent advances in VTP instrumentation that allow for the acquisition of high-resolution (0.001 °C) temperature data at short (1 cm) offsets near the ground surface, combined with advances in flux estimation methods that exploit the information content of the high-resolution data, facilitate heat-as-a-tracer approaches for characterizing groundwater-surface water exchanges and make it possible to obtain accurate and statistically robust results in a preferential discharge zone with a specific discharge of ∼1 m/d.
-
ArticleNon-equilibrium scour evolution around an emerged structure exposed to a transient wave(MDPI, 2024-06-05)The present study evaluates the performance of two numerical approaches in estimating non-equilibrium scour patterns around a non-slender square structure subjected to a transient wave, by comparing numerical findings with experimental data. This study also investigates the impact of the structure’s positioning on bed evolution, analyzing configurations where the structure is either attached to the sidewall or positioned at the centerline of the wave flume. The first numerical method treats sediment particles as a distinct continuum phase, directly solving the continuity and momentum equations for both sediment and fluid phases. The second method estimates sediment transport using the quadratic law of bottom shear stress, yielding robust predictions of bed evolution through meticulous calibration and validation. The findings reveal that both methods underestimate vortex-induced near-bed vertical velocities. Deposits formed along vortex trajectories are overestimated by the first method, while the second method satisfactorily predicts the bed evolution beneath these paths. Scour holes caused by wave impingement tend to backfill as the flow intensity diminishes. The second method cannot sufficiently capture this backfilling, whereas the first method adequately reflects the phenomenon. Overall, this study highlights significant variations in the predictive capabilities of both methods in regard to the evolution of non-equilibrium scour at low Keulegan–Carpenter numbers.
-
ArticleMeltwater orientations modify seismic anisotropy in temperate ice(American Geophysical Union, 2024-07-04)Seismology is increasingly used to infer the magnitude and direction of glacial ice flow. However, the effects of interstitial meltwater on seismic properties remain poorly constrained. Here, we extend previous studies on seismic anisotropy in temperate ices to consider the role of melt preferred orientation (MPO). We used the ELLE numerical toolbox to simulate microstructural shear deformation of temperate ice with variable MPO strength and orientation, and calculated the effective seismic properties of these numerical ice-melt aggregates. Our models demonstrate that even 3.5% melt volume is sufficient to rotate fast directions by up to 90°, to increase Vp anisotropy by up to +110%, and to modify Vs anisotropy by −9 to +36%. These effects are especially prominent at strain rates ≥3.17 × 10−12 s−1. MPO may thus obscure the geophysical signatures of temperate ice flow in regions of rapid ice discharge, and is therefore pivotal for understanding ice mass loss.
-
ArticleInhabitation of bathyal hydrocarbon seeps by early-branching benthic foraminifera: Implications for Neo-Proterozoic ecosystem functioning(Cushman Foundation for Foraminiferal Research, 2024-07-26)Little is known about mechanisms allowing agglutinated and thecate (i.e., organic-walled) foraminifera associated with extreme environments to survive. Here, we present cytological observations of organic-walled and agglutinated monothalamid and milioline tubothalamids collected from sulfur-oxidizing microbial mats of Gulf of Mexico hydrocarbon seeps. This is the first ultrastructural analysis of allogromids, a basal group of organic-walled thecate foraminifera, from seeps. We show most allogromids contained numerous ingested bacteria of various forms; organelle abundance and distribution varied among their morphotypes. Saccaminids (agglutinated monothalamids) had abundant phagocytosed bacteria, including putative methanotrophic and sulfur-oxidizing bacteria, indicating foraminiferal activity in seeps. A porcelaneous tubothalamid morphotype contained phagocytosed bacteria of one morphology, suggesting food selectivity. Increasing our knowledge of the cytology and ecology of these modern representatives of early-evolving foraminifera could help elucidate their evolutionary history. Thus, we augment understanding of extremophile foraminifera, adding to our burgeoning understanding of microeukaryote protists and microfossil adaptations.
-
ArticleIntensified currents associated with benthic storms underneath an eddying jet(American Geophysical Union, 2024-07-08)Benthic storms are episodes of intensified near-bottom currents capable of sediment resuspension in the deep ocean. They typically occur under regions of high surface eddy kinetic energy (EKE), such as the Gulf Stream. Although they have long been observed, the mechanism(s) responsible for their formation and their relationships with salient features of the deep ocean, such as bottom mixed layers (BMLs) and benthic nepheloid layers (BNLs), remain poorly understood. Here we conduct idealized experiments with a primitive-equation model to explore the impacts of the unforced instability of a surface-intensified jet on near-bottom flows of a deep zonal channel. Vertical resolution is increased near the bottom to represent the bottom boundary layer. We find that the unstable near-surface jet develops meanders and evolves into alternating, deep-reaching cyclones and anticyclones. Simultaneously, EKE increases near the bottom due to the convergence of vertical eddy pressure fluxes, leading to near-bottom currents comparable to those observed during benthic storms. These currents in turn form BMLs with thickness of O(100 m) from enhanced velocity shears and turbulence production near the bottom. The deep cyclonic eddies transport fluid particles both laterally and vertically, from near the bottom through the entire BML and may contribute to the formation of the lower part of BNLs. A sloping bottom reduces both the intensity of the near-bottom currents and the extent of vertical transport. Overall, our study highlights a significant response of the abyssal environment to near-surface current instability, with potential implications for sediment transport in the deep ocean.
-
ArticleA reciprocity‐based efficient method for improved source parameter estimation of submarine earthquakes with hybrid 3‐D teleseismic Green’s functions(American Geophysical Union, 2024-05-08)Accurate source parameters of global submarine earthquakes are essential for understanding earthquake mechanics and tectonic dynamics. Previous studies have demonstrated that teleseismic P coda waveform complexities due to near-source 3-D structures are highly sensitive to source parameters of marine earthquakes. Leveraging these sensitivities, we can improve the accuracy of source parameter inversion compared to traditional 1-D methods. However, modeling these intricate 3-D effects poses significant computational challenges. To address this issue, we propose a novel reciprocity-based hybrid method for computing 3-D teleseismic Green's functions. Based on this method, we develop a grid-search inversion workflow for determining reliable source parameters of moderate-sized submarine earthquakes. The method is tested and proven on five Mw5+ earthquakes at the Blanco oceanic transform fault (OTF) with ground truth locations resolved by a local ocean bottom seismometer array, using ambient noise correlation and surface-wave relocation techniques. Our results show that fitting P coda waveforms through 3-D Green's functions can effectively improve the source location accuracy, especially for the centroid depth. Our improved centroid depths indicate that all the five Mw5+ earthquakes on the Blanco transform fault ruptured mainly above the depth of 600°C isotherm predicted by the half-space cooling model. This finding aligns with the hypothesis that the rupture zone of large earthquakes at OTFs is confined by the 600°C isotherm. However, it is noted that the Blanco transform fault serves as a case study. Our 3-D source inversion method offers a promising tool for systematically investigating global oceanic earthquakes using teleseismic waves.
-
ArticleNitrogen isotope homogenization of dissolved ammonium through depth and 15N enrichment of ammonium during the incorporation into expandable layer silicates occurred in organic-rich marine sediment from Guaymas Basin, Gulf of California(Elsevier, 2024-06-01)Sedimentary nitrogen isotopic ratios are used as a proxy for ancient biogeochemical cycles on Earth's surface. It is generally accepted that sediment hole tops record primary signatures because organic nitrogen (ON) is predominant in this part of the hole. In contrast to such early to middle diagenetic stages, it is well known that heavier nitrogen isotope 15N tends to enrich in sedimentary rocks during later diagenetic and metamorphic stages. However, there are some critical gaps in our understanding of nitrogen isotopic alteration associated with abiotic processes during early-middle diagenesis. In this study, we examined the isotope ratios of ammonium nitrogen in interstitial water (IW) and total nitrogen (TN), including exchangeable ammonium and mineral nitrogen, in the solid-phase of organic-rich-sediment recovered by International Ocean Discovery Program (IODP) Expedition 385 cores drilled in the Guaymas Basin, Gulf of California, that contained ammonium-rich IW. The isotopic ratios (δ15N value) of TN are the most variable with depth compared to any other type of nitrogen. This variation can be interpreted as reflecting changes in the water mass environment in the basin caused by glacial–interglacial climate changes, modifying the δ15N values of the marine primary producers. Thus, the δ15N value of TN is a proxy for environmental change in the basin, while each component of TN shows different trends. The δ15N values of IW and exchangeable ammonium did not exhibit significant changes with depth, but the latter values are about 3 ‰ enriched in 15N. This may be due to advective transport of solute into adjacent layers followed by the formation of an isotopic equilibrium between IW and exchangeable ammonium in the case of fast sediment accumulation rate. The δ15N value of exchangeable ammonium is the highest among the other types of nitrogen with one exception, where the δ15N value of TN is the highest. The calculated δ15N values of ON based on mass balance are almost the same as those of associated TN in the shallow sediment layers (< 150 m below seafloor), but the difference in the δ15N values of TN and ON are significant in the deeper layers, where proportions of ON contents are <50%. In particular, in the layer where the δ15N value of TN is the highest, that of ON shows an even higher value and the difference reaches 3.5 ‰. The δ15N values of mineral nitrogen are similar to that of IW ammonium except the surface layers. Under such conditions, when δ15N value of TN is intermediate between those of mineral nitrogen and exchangeable ammonium, calculated δ15N value of ON is close to that of TN. On the other hand, if δ15N value of TN is out of the range between mineral nitrogen and exchangeable ammonium, it causes further difference in δ15N value of ON. It means that the fluctuation of δ15N values of TN is reduced relative to those of ON through depth. It has been considered that δ15N value of TN in sediment is similar to that of ON, and changes in the δ15N value of TN due to diagenesis are limited, but in such environment ON fluctuations over depth may be slightly underestimated.
-
ArticleDevelopment of a quantum cascade laser absorption spectrometer for simultaneous measurement of 13C-18O and 18O-18O clumping in CO2(Wiley, 2024-06-18)Dual clumped isotope paleothermometry determines carbonate formation temperatures by measuring the frequency of 13C–18O (∆638) and 18O–18O (∆828) pairs in carbonates. It resolves isotopic kinetic biases and thus enables more accurate paleotemperature reconstructions. However, high-precision measurements of 18O–18O clumping using current techniques requires large sample sizes and long acquisition times. We developed a mid-infrared isotope ratio laser spectrometer (IRLS) for simultaneous measurement of the isotopologue ratios ∆638 and ∆828 in gas-phase carbon dioxide (CO2) at room temperature. Our IRLS uses a single laser scanning from 2290.7 to 2291.1 cm−1 and a 31 m pathlength optical cell, and it simultaneously measures the five isotopologues required for calculating ∆638 and ∆828: 16O12C16O, 16O13C16O, 16O12C18O, 16O13C18O, and 18O12C18O. In addition, our IRLS can measure 16O12C17O, enabling ∆17O analysis. At ~20°C and a CO2 pressure of ~2 Torr, our IRLS system achieved precisions of 0.128‰ and 0.140‰ within 20 s for abundances of the clumped isotopologues 16O13C18O and 18O12C18O, respectively, and precisions of 0.267‰, 0.245‰, and 0.128‰ for 16O12C16O, 16O13C16O, and 16O12C18O. This yielded precisions of 0.348‰ (∆638) and 0.302‰ (∆828) within 25 s. Simulated sample–reference switching highlights the potential of our system and the need for further development. We demonstrated simultaneous measurements of ∆638 and ∆828 in CO2 to precisions of <0.35‰ within 25 s using a room-temperature, single-laser IRLS. Future developments on better resolving 16O12C16O and 16O13C16O peaks and system temperature control could further improve the measurement precision.
-
ArticleElastic stress coupling between supraglacial lakes(American Geophysical Union, 2024-05-10)Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.
-
ArticleGlobal ocean cooling of 2.3°C during the last glacial maximum(American Geophysical Union, 2024-05-08)Quantitative constraints on past mean ocean temperature (MOT) critically inform our historical understanding of Earth's energy balance. A recently developed MOT proxy based on paleoatmospheric Xe, Kr, and N2 ratios in ice core air bubbles is a promising tool rooted in the temperature dependences of gas solubilities. However, these inert gases are systematically undersaturated in the modern ocean interior, and it remains unclear how air-sea disequilibrium may have changed in the past. Here, we carry out 30 tracer-enabled model simulations under varying circulation, sea ice cover, and wind stress regimes to evaluate air-sea disequilibrium in the Last Glacial Maximum (LGM) ocean. We find that undersaturation of all three gases was likely reduced, primarily due to strengthened high-latitude winds, biasing reconstructed MOT by −0.38 ± 0.37°C (1σ). Accounting for air-sea disequilibrium, paleoatmospheric inert gases indicate that LGM MOT was 2.27 ± 0.46°C (1σ) colder than the pre-industrial era.
-
ArticleSubducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia(American Association for the Advancement of Science, 2024-06-07)The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.
-
ArticlePatterns of causative faults of normal earthquakes in the fluid‐rich outer rise of Northeastern Japan, constrained with 3D teleseismic waveform modeling(American Geophysical Union, 2024-06-15)Accurate earthquake source parameters are crucial for understanding plate tectonics, yet, it is difficult to determine these parameters precisely for offshore events, especially for outer-rise earthquakes, as the limited availability of direct P or S wave data sets from land-based seismic networks and the unsuitability of simplified 1D methods for the complex 3D structures of subducting systems. To overcome these challenges, we employ an efficient hybrid numerical simulation method to model these 3D structural effects on teleseismic P/SH and P-coda waves and determine the reliable centroid locations and focal mechanisms of outer-rise normal-faulting earthquakes in northeastern Japan. Two M6+ events with reliable locations from ocean bottom seismic observations are utilized to calibrate the 3D velocity structure. Our findings indicate that 3D synthetic waveforms are sensitive to both event location, thanks to bathymetry and water reverberation effects, and the shallow portion of the lithospheric structure. With our preferred velocity model, which has Vs ∼16% lower than the global average, event locations are determined with uncertainties of <5 km for horizontal position and <1 km for depth. The refined event locations in a good match between one of the nodal strikes and the high-resolution bathymetry, enabling the determination of the causative fault plane. Our results reveal that trench-ward dipping normal faults are more active, with three parallel to the trench as expected, while five are associated with the abyssal hills. The significant velocity reduction in the uppermost lithosphere suggests abundant water migrating through active normal faults, enhancing both mineral alteration and pore density.
-
ArticleNEWTS1.0: numerical model of coastal erosion by waves and transgressive scarps(European Geosciences Union, 2024-04-30)Models of rocky-coast erosion help us understand the physical phenomena that control coastal morphology and evolution, infer the processes shaping coasts in remote environments, and evaluate risk from natural hazards and future climate change. Existing models, however, are highly complex, are computationally expensive, and depend on many input parameters; this limits our ability to explore planform erosion of rocky coasts over long timescales (thousands to millions of years) and over a range of conditions. In this paper, we present a simplified cellular model of coastline evolution in closed basins through uniform erosion and wave-driven erosion. Uniform erosion is modeled as a constant rate of retreat. Wave erosion is modeled as a function of fetch, the distance over which the wind blows to generate waves, and the angle between the incident wave and the shoreline. This reduced-complexity model can be used to evaluate how a detachment-limited coastal landscape reflects climate, sea-level history, material properties, and the relative influence of different erosional processes.
-
ArticleSignatures of wave erosion in Titan’s coasts(American Association for the Advancement of Science, 2024-05-13)The shorelines of Titan’s hydrocarbon seas trace flooded erosional landforms such as river valleys; however, it is unclear whether coastal erosion has subsequently altered these shorelines. Spacecraft observations and theoretical models suggest that wind may cause waves to form on Titan’s seas, potentially driving coastal erosion, but the observational evidence of waves is indirect, and the processes affecting shoreline evolution on Titan remain unknown. No widely accepted framework exists for using shoreline morphology to quantitatively discern coastal erosion mechanisms, even on Earth, where the dominant mechanisms are known. We combine landscape evolution models with measurements of shoreline shape on Earth to characterize how different coastal erosion mechanisms affect shoreline morphology. Applying this framework to Titan, we find that the shorelines of Titan’s seas are most consistent with flooded landscapes that subsequently have been eroded by waves, rather than a uniform erosional process or no coastal erosion, particularly if wave growth saturates at fetch lengths of tens of kilometers.
-
ArticlePlasmid-borne biosynthetic gene clusters within a permanently stratified marine water column(MDPI, 2024-05-02)Plasmids are mobile genetic elements known to carry secondary metabolic genes that affect the fitness and survival of microbes in the environment. Well-studied cases of plasmid-encoded secondary metabolic genes in marine habitats include toxin/antitoxin and antibiotic biosynthesis/resistance genes. Here, we examine metagenome-assembled genomes (MAGs) from the permanently-stratified water column of the Cariaco Basin for integrated plasmids that encode biosynthetic gene clusters of secondary metabolites (smBGCs). We identify 16 plasmid-borne smBGCs in MAGs associated primarily with Planctomycetota and Pseudomonadota that encode terpene-synthesizing genes, and genes for production of ribosomal and non-ribosomal peptides. These identified genes encode for secondary metabolites that are mainly antimicrobial agents, and hence, their uptake via plasmids may increase the competitive advantage of those host taxa that acquire them. The ecological and evolutionary significance of smBGCs carried by prokaryotes in oxygen-depleted water columns is yet to be fully elucidated.
-
ArticleSISALv3: a global speleothem stable isotope and trace element database(Copernicus Publications, 2024-04-26)Palaeoclimate information on multiple climate variables at different spatiotemporal scales is becoming increasingly important to understand environmental and societal responses to climate change. A lack of high-quality reconstructions of past hydroclimate has recently been identified as a critical research gap. Speleothems, with their precise chronologies, widespread distribution, and ability to record changes in local to regional hydroclimate variability, are an ideal source of such information. Here, we present a new version of the Speleothem Isotopes Synthesis and AnaLysis database (SISALv3), which has been expanded to include trace element ratios and Sr isotopes as additional, hydroclimate-sensitive geochemical proxies. The oxygen and carbon isotope data included in previous versions of the database have been substantially expanded. SISALv3 contains speleothem data from 365 sites from across the globe, including 95 Mg/Ca, 85 Sr/Ca, 52 Ba/Ca, 25 U-Ca, 29 P-Ca, and 14 Sr-isotope records. The database also has increased spatiotemporal coverage for stable oxygen (892) and carbon (620) isotope records compared with SISALv2 (which consists of 673 and 430 stable oxygen and carbon records, respectively). Additional meta information has been added to improve the machine-readability and filtering of data. Standardized chronologies are included for all new entities along with the originally published chronologies. Thus, the SISALv3 database constitutes a unique resource of speleothem palaeoclimate information that allows regional to global palaeoclimate analyses based on multiple geochemical proxies, permitting more robust interpretations of past hydroclimate and comparisons with isotope-enabled climate models and other Earth system and hydrological models. The database can be accessed at https://doi.org/10.5287/ora-2nanwp4rk (Kaushal et al., 2024).
-
ArticleRotoBOD─quantifying oxygen consumption by suspended particles and organisms(American Chemical Society, 2024-05-08)Sinking or floating is the natural state of planktonic organisms and particles in the ocean. Simulating these conditions is critical when making measurements, such as respirometry, because they allow the natural exchange of substrates and products between sinking particles and water flowing around them and prevent organisms that are accustomed to motion from changing their metabolism. We developed a rotating incubator, the RotoBOD (named after its capability to rotate and determine biological oxygen demand, BOD), that uniquely enables automated oxygen measurements in small volumes while keeping the samples in their natural state of suspension. This allows highly sensitive rate measurements of oxygen utilization and subsequent characterization of single particles or small planktonic organisms, such as copepods, jellyfish, or protists. As this approach is nondestructive, it can be combined with several further measurements during and after the incubation, such as stable isotope additions and molecular analyses. This makes the instrument useful for ecologists, biogeochemists, and potentially other user groups such as aquaculture facilities. Here, we present the technical background of our newly developed apparatus and provide examples of how it can be utilized to determine oxygen production and consumption in small organisms and particles.
-
ArticleThe influence of viscous slab rheology on numerical models of subduction(European Geosciences Union, 2024-05-07)Numerical models of subduction commonly use diffusion and dislocation creep laws from laboratory deformation experiments to determine the rheology of the lithosphere. The specific implementation of these laws varies from study to study, and the impacts of this variation on model behavior have not been thoroughly explored. We run simplified 2D numerical models of free subduction in SULEC, with viscoplastic slabs following (1) a diffusion creep law, (2) a dislocation creep law, and (3) both simultaneously, as well as several variations of model 3 with reduced resistance to bending. We compare the results of these models to a model with a constant-viscosity slab to determine the impact of the implementation of different lithospheric flow laws on subduction dynamics. In creep-governed models, higher subduction velocity causes a longer effective slab length, increasing slab pull and asthenospheric drag, which, in turn, affect subduction velocity. Numerical and analogue models implementing constant-viscosity slabs lack this feedback but still capture morphological patterns observed in more complex models. Dislocation creep is the primary deformation mechanism throughout the subducting lithosphere in our models. However, both diffusion creep and dislocation creep predict very high viscosities in the cold core of the slab. At the trench, the effective viscosity is lowered by plastic failure, rendering effective slab thickness the primary control on bending resistance and subduction velocity. However, at depth, plastic failure is not active, and the viscosity cap is reached in significant portions of the slab. The resulting high slab stiffness causes the subducting plate to curl under itself at the mantle transition zone, affecting patterns in subduction velocity, slab dip, and trench migration over time. Peierls creep and localized grain size reduction likely limit the stress and viscosity in the cores of real slabs. Numerical models implementing only power-law creep and neglecting Peierls creep are likely to overestimate the stiffness of subducting lithosphere, which may impact model results in a variety of respects.