Terrestrial hot springs have existed throughout Earth’s history and house some of the most ancient evidence of life on our planet. These settings are known for their high habitability and preservation potential, and are extensively studied as analog environments since hot spring deposits are thought to exist on the surface of Mars. Hot spring water commonly precipitates silica that coats microbial life dwelling in the hot spring outflow streams. This process can entomb microorganisms and preserve microbial remains over long timescales and with high morphological fidelity. Here we present research carried out on modern and sub-recent remains of microbial filaments from amorphous (unaltered) silica deposits in Yellowstone National Park. This work suggests that various elements sequestered by hot spring-dwelling organisms during life are preserved in microbial remains and persist over > 10,000 years. We also present findings from microfossils preserved in mid-Paleozoic terrestrial hot spring deposits which also show sequestrations of select elements in microfossil remains, suggesting that certain elements may persist even after several hundred million years and substantial host rock alteration. These elemental concentrations may be indicative of metabolic functioning during life and have application as biosignatures. Recent developments in analytical instrumentation now allow for even extremely low trace elemental abundances to be detected and mapped, regardless of sample complexity. This work is especially relevant to the search for life on Mars, as evidence of impact-induced hydrothermal activity may exist near the rim of Jezero Crater and may be sampled by the Perseverance rover. As a primary objective of the Mars 2020 mission is to search for evidence of past life on Mars, we suggest the application of this analytical technique to be valuable for potential samples returned to Earth by future Mars Sample Return missions. Distribution Statement A. Approved for public release: distribution unlimited.

Introduction The complex geomorphology of Triton reflects its geological history [1]. The morphological heterogeneity on small scale has led to discerning three main complexes: cantaloupe terrains, equatorial plains, and south polar cap terrains. We analyse an area located in the east equatorial zone of Triton called Monad Regio (centred at 37°N, 2°E), characterized by the presence of walled plains. We produced a new geological map and a DEM (Digital Elevation Model) to recognize the main terrains and features in the study area. We used Voyager 2 imagery named c1139533 (600 m/px) [2], properly calibrated, filtered, and georeferenced using the Integrated Software for Imagers and Spectrometers (ISIS4) [3]. Results and conclusions We mapped the different geological units and main features according to differences in surface morphology (Fig.1). Terraced terrain covers most of the studied area. It shows a chaotic pattern characterized by several terraces, some of which lay in a parallel arrangement around some of the large depressions. These basins have areas ranging from 1300 to 2050 km2, and their degree of alteration is variable, with the features inferred to be more recent showing an inner minor basin within the main one. The most altered basins appear smoother, featureless, and shallower. Sizes and excavation depths estimated using DEM data of the observed basin features appear to be relatively homogeneous, which leads us to exclude an impact related origin. We argue that the origin of these depressions is linked to processes analogue to those described in the formation of terrestrial maar craters and possible explosion craters discussed on Titan [4]. Alternatively, diapirism may also explain the origin of such features. Further analysis could help to understand the nature and related processes that originated these basins. Acknowledgements G.M., C.C. and D.S. acknowledge support from the Italian Space Agency (2020-13-HH.0). References [1] Basilevsky A.T. et al.,1992. Adv. Space Res.,12(11), 123-132. [2] Smith, B. A., et al. (1989), Science, 246 (4936), 1422-1449. [3] Houck J.C. and DeNicola L.A. (2000) Astronomical Data Analysis Software and Systems IX, ASP Conference Series,216. [4] Mitri G., et al. (2019), Nature Geoscience, 12, 791,796.

Passage zones are stratigraphic surfaces found in littoral settings separating deposits diagnostic of subaqueous environments from overlying sequences of subaerial deposits. In glaciovolcanic settings, passage zone surfaces are unequivocal records of the heights and depths of paleo-englacial lakes at a specific point in time and space thereby informing on the presence and nature of the enclosing ice sheet. Kima Kho, a Pleistocene glaciovolcano (i.e. tuya) features multiple and diverse passage zones. The basaltic volcano comprises 4 main stratigraphic packages: i) subaqueously and subaerially deposited lapilli tuffs forming a central tephra cone and representing an explosive onset to the eruption, ii) subaqueously deposited, steeply-inclined beds of tuff breccia dominated by pillow lava fragments, iii) stacked sheets of subaerial lavas , and iv) dykes and sills (I) intruding all units. Stratigraphic and geochemical relationships suggest that Kima Kho volcanism was continuous and 40Ar/39Ar geochronometry on 3 samples yields a mean age of 1949 ±63 ka. Three temporally distinct passage zones record the interplay between growth of the volcanic edifice, syn-volcanic melting of the enclosing ice sheet, and fluctuations in the depth of the englacial lake. The earliest passage zone is expressed in two different ways indicating a transition to effusive eruption: (i) within pyroclastic deposits of the tephra cone (< 1800 masl), and (ii) by pillow lava tuff breccia deposits overlain by subaerial lavas. Together they record a peak, sustained lake depth of 320-340 m that constrains the enclosing ice sheet to a minimum thickness of ~400 m and a minimum radial extent >7 km relative to Kima Kho. Two subsequent passage zones, also defined by sequences of subaerial lavas resting on dipping beds of pillow lava tuff breccias, occur at lower elevations: 1690-1640 masl and 1740-1720 masl, respectively. The latter two passage zones indicate a major draining of the englacial lake followed by refilling to depths of 230-180 m and 260-280 m, respectively. The substantial decline in lake level between passage zones suggests a massive, catastrophic deluge (i.e. jökulhlaup) of 1-2 km3. Lastly, the reconstructed evolution of Kima Kho demands the presence of a regionally extensive, cold-based ice sheet on the Kawdy plateau at ~1.9 Ma.

Oceanic Transform Faults (TF) comprise first order discontinuities bounded between mid-ocean ridge spreading centres. TF mainly accommodate strike slip motion, separating lithospheric plates of different age and thermal structure. Oceanic TF are intriguing in that they do not produce earthquakes as large as might be expected given their long length, with seismic slip corresponding only to a small fraction of the total tectonic slip. The relative geologic simplicity of oceanic TF means that they are an important analogue for more hazardous continental TF, with high potential for improving insights into the earthquake cycle. We investigate the earthquake properties along Chain, a ~300 km long TF in the equatorial MAR by combining both microseismic and teleseismic data. We use the ~1-year microseismicity data (total of 812 events) gathered during the PI-LAB (Passive Imaging of the Lithosphere-Asthenosphere Boundary) experiment and EURO-LAB (Experiment to Unearth the Rheological Lithosphere-Asthenosphere Boundary). We perform cluster analysis in multi-dimensional phase space, consisting of various seismic (epicentral coordinates, magnitude) and geophysical (gravity anomalies, bathymetry, tidal height) parameters. We investigate potential triggering mechanisms, including tidal, static and dynamic stresses. We extend our analysis back in time by considering stronger earthquakes (MW>~5.0) from Global Centroid Moment Tensor (GCMT) since 1976. We find three unique, 50-100 km long clusters or segments from our analysis going from east to west, separated by seismic gaps. Microseismic activity is highest at the eastern segment of Chain where there is the largest positive flower structure, negative rMBA gravity anomaly but very few M>5.5 events. The western segment has reduced seismicity rates relative to the eastern, and is associated with a positive rMBA and a few small flower structures. The central segment is bounded between two seismic gaps and demonstrates relatively high activity rates in the middle. Our result suggests that trans-pression of highly altered mantle/crust and/or high pore pressure due to hydrothermal fluid circulation in the eastern flower structure enhances seismic activity. Overall, we find the existence of consecutive locking and creeping segments, with some of the patches exhibiting hybrid behaviour, potentially causing their sporadic activation/reactivation.

Permafrost underlies approximately one fifth of the global land area and affects ground stability, freshwater runoff, soil chemistry, and surface‑atmosphere gas exchange. The depth of thawed ground overlying permafrost (active layer thickness, ALT) has broadly increased across the Arctic in recent decades, coincident with a period of increased streamflow, especially the lowest flows (baseflow). Mechanistic links between ALT and baseflow have recently been explored using linear reservoir theory, but most watersheds behave as nonlinear reservoirs. We derive theoretical nonlinear relationships between long‑term average saturated soil thickness η (proxy for ALT) and long-term average baseflow. The theory is applied to 38 years of daily streamflow data for the Kuparuk River basin on the North Slope of Alaska. Between 1983–2020, the theory predicts that η increased 0.11±0.17 [2σ] cm a-1, or 4.4±6.6 cm total. The rate of change nearly doubled to 0.20±0.24 cm a-1 between 1990–2020, during which time field measurements from CALM (Circumpolar Active Layer Monitoring) sites in the Kuparuk indicate η increased 0.31±0.22 cm a-1. The predicted rate of change more than doubled again between 2002–2020, mirroring a near doubling of observed ALT rate of change. The inferred increase in η is corroborated by GRACE (Gravity Recovery and Climate Experiment) satellite gravimetry, which indicates that terrestrial water storage increased ~0.80±3.40 cm a-1, ~56% higher than the predicted increase in η. Overall, hydrologic change is accelerating in the Kuparuk River basin, and we provide a theoretical framework for estimating changes in active layer water storage from streamflow measurements alone.

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H2, CH4) or permanent (CO2, radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2nm-5μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

We propose a thermal-fluid system model for hollow formation. A subsurface heat source (typically impact-related) produces volatiles from LRM and drives them to the surface. Volatiles generated through heating of LRM are likely S and S-bearing gases produced by thermal decomposition of sulfides heated by the impact process. C-bearing volatiles, such as CH4 and other simple organics, and potentially fullerenes within LRM, may also be involved in proposed thermal-fluid systems responsible for hollow formation.

Lunar mare basalts are the products of their corresponding parent magma compositions, sourced from the lunar upper mantle. The lunar mantle has been repeatedly modeled through numeric simulations to reflect lunar magma ocean (LMO) crystallization, resulting in an early-stage anorthositic crust and immediately underlying, late-stage, KREEP-rich and ilmenite-rich layer. This negatively buoyant layer is expected to have induced mixing with the underlying mantle, potentially to the core-mantle boundary. The lunar mare basalts, in this context, reflect mantle sources that are variably mixed between pristine mantle compositions and the dense ilmenite-rich layer. In order to constrain the geometry of lunar mantle heterogeneity, we simultaneously examined multiple mare basalt characteristics to extract significant multivariate patterns that might lend insight into the nature of this mixing-induced heterogeneity. Using two fundamental machine learning approaches and a newly compiled database of Apollo basalt characteristics (ApolloBasaltDB), we conducted a preliminary investigation, holding the assumptions that 1) mare basalts are assumed to retain the majority of their original characteristics at the time of extrusion: texture, isotopic age, major element composition, mineral mode, and geographic occurrence; 2) negligible basalt alteration occurred due to the lack of an atmosphere; and 3) impact gardening did not have significant bearing on final geographic location of basalt samples based on our nearside spatial partitions. The results of cluster and principal component analyses over changing spatial basalt groupings suggest that lunar nearside changes in major element concentrations and mineral modes vary spatially. Al2O3 concentrations increase in diversity within the Procellarum KREEP Terrane (PKT) compared to older regions immediately exterior to the eastern PKT, while a general nearside trend appears to suggest that ilmenite (TiO2) diversity comes at the expense of plagioclase (Al2O3) diversity. Cluster analysis suggests PKT perimeter rifting may have tapped into increasingly Ti-rich sources as rifting proceeded SE to NW. By establishing such trends over varying spatial scales through multivariate processing, the changing strengths of these surface correlative patterns may indicate the changing conditions (including temporal) of the immediately underlying mantle at the time of extrusion.

Accurate descriptions of natural fault surfaces and associated fault rocks are important for understanding fault zone processes and properties. Slickensides–grooved polished surfaces that record displacement and wear along faults– develop measurable roughness and characteristic microstructures during fault slip. We quantify the roughness of natural slickensides from three different fault surfaces by calculating the surfaces power spectra and height distributions and analyze the microstructures formed above and below the slickensides. Slickenside surfaces exhibit anisotropic self-affine roughness with corresponding mean Hurst exponents in directions parallel– 0.53±0.07– and perpendicular –0.6±0.1– to slip, consistent with reports from other fault surfaces. Additionally, surfaces exhibit non-Gaussian height distributions, with their skewness and kurtosis roughness parameters having noticeable dependence on the scale of observation. Below the surface, microstructural analyses reveal that S-C-C’ fabrics develop adjacent to a C-plane-parallel principal slip zone characterized by a sharp decrease in clast size and a thin (≤100 µm) nanoparticulate-rich principal slip surface (PSS). These microstructures are present in most analyzed samples suggesting they commonly form during slickenside development regardless of lithology or tectonic setting. Our results suggest that 1) PSS likely arise by progressive localization along weaker oriented fabrics 2) deformation along PSS’s is energetic enough to comminute the rocks into nanometric grains, and 3) fault geometry can be further characterized by studying the height distributions of fault surfaces, which are likely to impact stress distributions and frictional responses along faults.

The new dual-frequency radar satellite ASAR-ISRO provide simultaneously and for the first time a wide wavelength spectrum which is critical to discriminate surface roughness based on different backscattering characteristics. Here we use such dual L-band and S-band airborne SAR system to characterize and map various volcanic areas in the Northern Cascades through their backscattering properties. Mapping volcanic flows (lava flows, pyroclastic currents, lahars) is vastly improved by using backscattering as a metric of surface roughness. Various types of volcanic flow deposits and surface textures are distinguished by their roughness measured with radar systems. The ability of radar systems to distinguish volcanic flow textures, represented by roughness, is a key factor in understanding the processes and timescales of flow emplacement. For instance, transition from pahoehoe to aa lava flows is associated with change in flux and steepness of topographic gradient. Therefore, lava textures are essential data for calibrating and improving lava flow emplacement codes. Similarly, the textures of high velocity, and more deadly, pyroclastic currents and lahars change along their flow paths, also revealing critical data about the mechanisms of flow emplacement. As with lava flows, we use the ASAR-ISRO L+S SAR system to characterize the run-out distances of pyroclastic currents and lahars where large numbers of blocks accumulate in such deposits (places where flow momentum was lost, whether due to friction or break-in-slope and vast quantities of blocks accumulate in a relatively small area compared to the total area inundated by the flow). Mapping volcanic flow textures in a variety of volcanic terrains will provide clues about modes and rates of emplacement, and change in these through time, in a way that is simply unavailable by traditional geologic mapping. The next generation of volcanic flow maps, used for hazard assessment, will rely on radar data to delineate these textures.

My earth science focused PhD research included the analysis of annual fluorescent laminae in cave stalagmites, under the supervision of karst hydrologist Peter Smart (Bristol, UK). At the time, the source of this fluorescence was uncertain, opening new research opportunities characterizing what is now understood to be soil-derived, water-soluble fluorescent dissolved and colloidal organic matter that is transported to the cave by vadose zone percolation waters. Aquatic organic matter fluorescence research benefitted at this time from significant laboratory analytical advances that were commercially driven, including Sony’s Blu-ray technology and the use of fluorescent labelling in the biomedical sciences. Faster analyses at increasingly higher energy excitation energies opened opportunities for novel fingerprinting of organic matter in hydrological science, including landfill leachates and sewage contamination. Today, hand-held fluorescence sensors can instantaneously determine microbial water quality. And back in the field of speleothem (cave deposit) science, annual geochemical laminae are now recognized to be widely preserved in regions where there is a seasonality in recharge, providing a precise chronology for the stalagmite paleoenvironmental archive. Currently, we are utilizing this precise chronology to generate high-resolution fire history records, where the fire proxy is water-soluble ash-derived elements transported from soil to cave by vadose zone percolation waters.

Surface heterogeneities below the spatial resolution of thermal infrared (TIR) instruments result in anisothermality and produce emissivity spectra with negative slopes at longer wavelengths. Sloped spectra arise from an incorrect assumption of either a uniform surface temperature or a maximum emissivity during the temperature-emissivity separation of radiance data. Surface roughness and lateral mixing of differing sub-pixel surface units result in spectral slopes that are distinct, with magnitudes proportional to the degree of temperature mixing. Routine Off-nadir Targeted Observations (ROTO) of the Thermal Emission Imaging Spectrometer (THEMIS) are used here for the first time to investigate anisothermality below the spatial resolution of THEMIS. The southern flank of Apollinaris Mons and regions within the Medusae Fossae Formation are studied using THEMIS ROTO data acquired just after local sunset. At higher emission angles, differing relative proportions of rocky and unconsolidated surface units are observed. This produces a range of sloped TIR emission spectra dependent on the magnitude of temperature differences within a THEMIS pixel. Spectral slopes and wavelength-dependent brightness temperature differences are forward-modeled for a series of two-component surfaces of varying thermal inertia values. This creates a thermophysical model suggesting a local rock abundance 6 times greater than currently published results and four orders of magnitude more sensitive than those relying on nadir data High-resolution visible images of these regions indicate a mixture of surface units from boulders to dunes, providing credence to the model.

High-conductivity anomalies of 0.1–1 S/m are widely distributed in the mid-lower crust of the Tibetan Plateau. Dehydration of amphibole-bearing rocks may play an important role in explaining these anomalies. To survey the anomalies’ origin, therefore, the electrical conductivities of amphibole-bearing samples, containing varying amphibole content, are measured at 1.5 GPa and 600–1300 K. Our experiments show that dehydration melting occurs at about 1100 K. Proton conduction and ionic conduction dominate the conduction mechanisms before and after dehydration melting, respectively. The dehydration melting of felsic rocks, containing 25 vol% of amphibole, is unable to account for the high-conductivity anomalies of 0.1–1 S/m. In contrast, the dehydration melting of garnet amphibolite, with an amphibole content higher than 60 wt%, can enhance the bulk conductivity to higher than 0.1 S/m under the lower-crust conditions beneath the Tibetan Plateau. The melt fraction of the garnet-amphibolite is estimated to be 3.8–36 vol% in the partial molten region based on a cube-model simulation.

The fractal dimension and multifractal spectrum can characterize the complexity of fracture sets. However, studies of impacts of fracture geometries on their fractal and multifractal characteristics are largely insufficient, especially for three-dimensional (3-D) fracture networks (natural fractures are always 3-D instead of 2-D). In this work, we construct 3-D stochastic discrete fracture networks with an open-source DFN software, HatchFrac. Systematical investigations are then conducted to study the impact of geometrical fracture properties and system sizes on the fractal and multifractal characteristics. The box-counting method is adopted to calculate the fractal dimension and multi-fractal descriptors. The fractal dimension, D, and the difference of the singularity exponent, ∆α, represent the fractal and multifractal patterns, respectively. Two critical (percolative and over-percolative) stages of fracture networks are considered. 3-D fracture networks share similar characteristics with 2-D fracture networks at percolation. However, results at an over-percolative stage are systematically different. At the first stage, fracture orientations (κ), lengths (a) and system sizes (L) have positive correlations with D and ∆α. D is weakly correlated with fracture positions (FD), meaning that the fractal dimension is insensitive to clustering effects. However, ∆α is strongly correlated with FD, implying that ∆α can characterize the heterogeneity caused by clustering effects. a and L are positively correlated with ∆α, and κ and FD have negative correlations. At stage two, the sensitivity results on D are similar to stage one, but a and L become negatively correlated with ∆α. Impacts of κ and FD become more significant.

We surveyed the subsurface structure in eastern Coprates and Capri Chasmata in the equatorial region using high-resolution visible images, digital terrain models, and radar sounding data. We identified subsurface reflectors in four areas of the chasmata. At the stratigraphic exposure on the chasmata walls, the corresponding depth of the reflector is ~60 m. The bulk dielectric constants of the layers above the reflectors are calculated as 3.4-4.0, suggesting a rock-air mixture with ~46.1% porosity, or a rock-air-ice mixture with ~21.2% water-ice fraction. Recent climate models suggest that water-ice is unstable on the surface around the equatorial regions. However, considering the recent high obliquity that occurred ~0.4 Ma and a slow diffusivity of water-ice, the existence of subsurface water-ice deeper than a few meters cannot be ruled out. If water-ice is actually contained in the layer, our results show the maximum volume of putative water-ice in the chasmata is 16.6 km.

Goundwater investigations are increasingly becoming important for Indonesia. It is natural resources for sustainable development of a region. In this study, PT Sejahtera Alam Energy conducted an investigation of subsurface conditions, especially related to groundwater potential, carried out by a resistivity investigation method by means of Vertical Electrical Sounding (VES) or Geoelectricity in Embung Area, Baturaden Geothermal Project, Pandansari Village, Paguyangan Subdistrict, Brebes Regency. This method is one of the subsurface estimation methods that is considered suitable for water investigations as well as in terms of accuracy as well as in terms of low cost and faster implementation time compared to other geophysical methods.

Velocity is one of the fundamental data obtained from seismic and it is the direct behavior of the solids and fluids in lithosphere. Here we present an analysis of the derivative gradient of seismic velocity which can help identify the featured boundary of favorable petroleum accmulations such as sedimentary facies boundary, faults, and flow unit edges. The derivative gradient can be calculated both horizontally and vertically, thus it can help discriminate favorable targets in three-dimensional. We find the application of derivative gradient to detect or enhance edge is relatively mature in gravity and geomagnetic analysis but rarely mentioned in seismic nor in petroleum exploration. We believe we can make better use of the seismic velocity data by this means as it is quite efficient in pinpointing the favorable petroleum targets in subsurface with precisions of tens-of-meter scale, depending on the horizontal resolution of seismic survey.

Hearty, P. J., Rovere, A., Sandstrom, M. R., O’Leary, M. J., Roberts, D., & Raymo, M. E. (2020). Pliocene‐Pleistocene stratigraphy and sea‐level estimates, Republic of South Africa with implications for a 400 ppmv CO2 world. Paleoceanography and Paleoclimatology, 35, e2019PA003835. https://doi.org/10.1029/2019PA003835 The Mid-Pliocene Warm Period (MPWP, 2.9 to 3.3 Ma), along with older Pliocene (3.2 to 5.3 Ma) records, offers potential past analogues for our 400-ppmv world. The coastal geology of western and southern coasts of the Republic of South Africa expose an abundance of marine deposits of Pliocene and Pleistocene age. In this study, we report differential GPS elevations, detailed stratigraphic descriptions, standardized interpretations, and dating of relative sea-level indicators measured across ~700 km from the western and southern coasts of the Cape Provinces. Wave abrasion surfaces on bedrock, intertidal sedimentary structures, and in situ marine invertebrates including oysters and barnacles provide precise indicators of past sea levels. Multiple sea-level highstands imprinted at different elevations along South African coastlines were identified. Zone I sites average +32 ± 5 m (6 sites). A lower topographic Zone II of sea stands were measured at several sites around +17 ± 5 m. Middle and late Pleistocene sites are included in Zone III. Shoreline chronologies using 87Sr/86Sr ages on shells from these zones yield ages from Zone I at 4.6 and 3.0 Ma, and Zone II at 1.04 Ma. Our results show that polar ice sheets during the Plio-Pleistocene were dynamic and subject to significant melting under modestly warmer global temperatures. These processes occurred during a period when CO2 concentrations were comparable to our current and rapidly rising values above 400 ppmv.

Dob’s Linn (Scotland) is a location that has significantly influenced our understanding of how life evolved over the Ordovician to early Silurian. The current chronostratigraphic boundary between the Ordovician and Silurian periods is a Global Boundary Stratotype Section and Point (GSSP) at Dob’s Linn calibrated to 443.8±1.5 Ma, partly based on biostratigraphic markers, radiometric ages, and statistical modeling. Graptolites are used here as relative dating markers. We dated hundreds of zircon grains extracted from defined metabentonites from six horizons exposed at Dob’s Linn using Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS). Each zircon was imaged using cathodoluminescence, and most show igneous zoning with minimal alteration. Sample locations range from 42 meters above to 5 meters below the recognized GSSP for the Ordovician-Silurian. Samples were responsibly collected and analyzed for paleontology and geochemistry in other work. Overall, many 238U-206Pb zircon ages from the section are significantly younger than expected. The youngest zircon in sample DL7, located 5 meters below the GSSP, yielded a 238U-206Pb age of 402±12 Ma (±2s, 5% disc). Nineteen spots on zircons from this sample are younger than the presently assigned GSSP age, including more concordant results of 426±8 Ma (0.8% disc) and 435±5 Ma (0.2% disc). The youngest zircon in sample 19DL12, < 1 m below the GSSP, is 377±8 Ma (2% disc) with a more concordant age of 443±7 Ma (0.6% disc). A sample located directly on the GSSP (19DL09) yields 327±5 Ma (0.8% disc). Eight spots on zircons from this sample are also younger than the presently assigned GSSP age. We also dated two samples (DL24 and BRS23) 8 meters above the GSSP, and the youngest, most concordant zircon ages in these samples are 400±11 Ma (5% disc) and 421±9 Ma (0.4% disc), respectively. Overall, the U-Pb ages would re-assign the Dob’s Linn chronostratigraphic section to Silurian-Devonian. The young age results could be attributed to Pb loss due to hydrothermal alteration during the Acadian and Alleghenian orogenies. Future work will implement Chemical Abrasion Isotope Dilution Thermal Ionization Mass Spectrometry (CA-ID-TIMS) to obtain accurate U-Pb dating and evaluate the potential effects of Pb loss.

Despite many studies on seafloor hydrothermal systems conducted to date, the generation mechanism of seafloor massive sulfide (SMS) deposits is not yet fully understood. To elucidate this mechanism, this study clarifies the three-dimensional regional temperature distribution and fluid flow of a seafloor hydrothermal system of the Iheya North, middle Okinawa Trough. Lateral flow and boiling of hydrothermal fluids below the seafloor were the main features found by the simulation, leading to an interpretation of two-layered SMS deposit generation as follows. Hydrothermal fluids discharging from black smokers first formed the upper SMS deposits on the seafloor. Caprocks formed below the seafloor, and the above-mentioned occurrences were then induced under the caprocks. In the present system, vapor-rich hydrothermal fluids poor in metals are discharged from the vents as white smokers, whereas liquid-dominated hydrothermal fluids rich in metals flow laterally below the caprocks, forming lower SMS deposits tens of meters below the seafloor.

Little is known of the structure and seismic nature of icy satellites, planetary bodies that have an outer shell of ice instead of a rocky crust. One such body, Saturn's moon Titan, is unique because of its thick atmosphere, hydrologic cycle similar to Earth's, and salty subsurface ocean. The Dragonfly mission, set to launch in 2026, will use a robotic rotorcraft to transport various instruments, including a seismometer, to the Shangri-La dune field on Titan. Until then, Titan's seismic regime can be estimated by simulating wave propagation through assumed layering using source models. Ridge belts with low slope angles suggest fold-and-thrust belts possibly caused by fluid overpressures in Titan's icy shell. These have been observed in synthetic aperture radar images from Cassini and may be a likely site for Titan-quakes. Synthetic seismograms are calculated for Titan-quakes using Instaseis, a Python-based software sourced by Green's function databases computed by the axisymmetric spectral element method. Back-azimuth angles are calculated for different source-receiver configurations. Very high amplitude waves are interpreted as surface waves. The PP wave is the first arrival at distances greater than 35 degrees. At 45 degrees, S waves are not distinguishable.

We quantitatively evaluate transducer-transducer one-source one-station active seismic waveform data, in order to monitor time-lapse changes of elastic and anelastic structure during deformation experiments in laboratory. The experiment data of dry and water-saturated sample are provided by Zaima and Katayama (2018, https://doi.org/10.1029/2018JB016377). A transducer receiver, at the mid-point of cylindrical rock sample, is located on the antipodal position of the transducer source, emitting compressional and shear waves. Due to the extremely underdetermined nature of inverse problem, we limit the number of unknowns to be four: global P- and S- wave velocities and their corresponding anelastic attenuation factors, which can represent the micro-cracks nucleation during the loading and before the appearance of the largest crack that causes the fracture. We first performed a trial-and-error search for a realistic boundary condition in three-dimensional seismic waveform modeling using spectral-element method, in order to fit the synthetic data with the observed waveforms. We then generated synthetic data for 6000 combinations of elastic and anelastic parameters, in order to conduct Monte-Carlo waveform inversion based on the cost functions using waveform misfit and zero-lag cross-correlation. We obtained the time-lapse changes in velocity and attenuation during the deformation, which are then linked to crack development. Compared with the wet experiment, the dry experiment has a larger change in both the velocity and attenuation. However, regardless of the configuration, global seismic wave speeds rise first and then decrease during the experiments. The quality factor shows roughly the same trend.

Carbon isotope (δ13C) records from marine sediments have been extensively used in Cenozoic chemostratigraphy. The early Paleogene interval in particular has received exceptional attention because negative carbon isotope excursions (CIEs) documented in the sedimentary record, e.g. at the Paleocene Eocene Thermal Maximum (PETM), ca ~56 Ma, are believed to reflect significant global carbon cycle perturbations during the warmest interval of the Cenozoic era. However, while bulk-carbonate δ13C values exhibit robust correlations across widely separated marine sedimentary sections, their absolute values and magnitude of CIEs vary spatially. Moreover, bulk-carbonates in open-marine environments are an ensemble of different components, each with a distinct isotopic composition. Consequently, a complete interpretation of the bulk δ13C record requires an understanding of co-evolution of these components. In this study, we dissect sediments, from early Paleogene interval, at ODP Site 1209, Shatsky Rise, Pacific Ocean to investigate how an evolving bulk-carbonate ensemble influences the overall carbon isotope record. A set of 45 samples were examined for δ13C and δ18O compositions, as bulk and individual size fractions. We find a significant increase in coarse-fraction abundance across PETM, driven by a changing community structure of calcifiers, modulating the size of CIE at Site 1209 and thus making it distinct from those recorded at other open-marine sites. These results highlight the importance of biogeography in marine stable-isotope record, especially when isotopic excursions are driven by climate- and/or carbon-cycle changes. In addition, community composition changes will alter the interpretation of weight percent coarse fraction as conventional proxy for carbonate dissolution.

For youth with limited role models in the STEM fields, and restricted summer research opportunities resulting from a lack of financial resources and academic connections, the opportunity to participate in academically connected, community based science research programs can be incredibly empowering. Providing these opportunities is critically important but it takes purposeful work, persistent outreach and strong community networks. We note that while providing these opportunities is incredibly rewarding, but there is a lot of work both up front and ongoing. This work is itself rewarding when networks are humming and enthusiasm for involvement is high, but it can be challenging when ceilings are hit and walls seem to arise unexpectedly. “Early Engagement in Research: Broadening participation through engagement in authentic science research” builds a regional network of summer research experiences for high school students underrepresented in STEM, starting from a successful model that has provided high school summer field research opportunities for New York City youth for over a decade (Secondary School Field Research Program). The program is developed around regional partnerships between various combinations of academic institutions and research centers, community environmental and education centers, state cooperative extensions, high schools and school networks, state and local park systems and land management groups. Each location has a unique approach, but all include some similar attributes. Each tackles an authentic science research issue that affects the local community such as microbiology in the local streams and microplastics in the local bays and biology, and each includes peer and near peer mentoring for the students along with a scientist mentor. Encouraging professional development of each student is central to the program. Technical instruction includes the use of scientific instruments and equipment, data recording and interpretation. Professional discussions include how to successfully read and dissect a science journal article, how to create and present a science poster and most importantly how build a network for themselves in STEM, and how to help us work with them to support the diversity that is needed for all of science to be inclusive and ultimately meet the needs of our future.