Subsurface Geometry Research Articles

Potential field modelling and U–Pb geochronology reveal the pluton emplacement dynamics of the Lower Devonian Tarnagulla Granodiorite, southeast Australia

The Tarnagulla Granodiorite in the western Lachlan Fold Belt is poorly exposed; however, magnetic data reveal its distinct U shape (in plan view) and patterns characterised by multiple nested cusp-shaped anomalies of contrasting magnetic intensity. We interpret this magnetic response to reflect incremental assembly of the pluton above a persistent magma point source that migrated laterally. Here, we use potential field data, geophysical modelling and high-precision zircon U–Pb dating to constrain the shape of the intrusion at depth and deduce the migration path of the magmatism epicentre. Potential field profile modelling across three traverses, constrained by petrophysics, resolves the subsurface geometry of the Tarnagulla Granodiorite. In cross-section, the intrusion is funnel-shaped, with a lower boundary that varies from ∼2 km depth in the north to ∼6 km depth in the southwest, reaching a maximum depth of ∼8 km, which may coincide with a magma feeder zone. New precise U–Pb geochronology provides ages that confirm the interpreted magnetic phases were emplaced sequentially along a clockwise intrusion path that progressed stepwise from the north to the southwest between 406.62 ± 0.27 Ma and 399.3 ± 1.1 Ma. We attribute migration of the epicentre of magmatism within the Tarnagulla Granodiorite to syn-intrusion (Early Devonian) reactivations of the underlying Tarnagulla Fault, with relative displacements between the fault hanging wall (which hosts the pluton) and footwall (which we interpret to have hosted a persistent point-source magma feeder zone during pluton emplacement). Our model explains the intrusion’s U-shaped magnetic response and reveals critical information about the paleo-stress and paleo-strain conditions during emplacement, characterised by three main displacement vectors and fault slip magnitudes. The geophysical characteristics, path of emplacement and mineralisation distribution all suggest that progressive fractional crystallisation contributed to the formation of an enriched final pulse of magma that resulted in Mo, Cu and W mineralisation hosted in the Moliagul Granodiorite. KEY POINTS High-precision U–Pb zircon geochronology reveal that the Tarnagulla Granodiorite was emplaced along a clockwise pathway from the northeast to the southwest between 406.62 ± 0.27 Ma and 399.3 ± 1.1 Ma. Active faulting during pluton emplacement is the suggested explanation for the Tarnagulla Granodiorites’ U-shaped magnetic response. Fractional crystallisation throughout pluton emplacement produced the varying internal geophysical response of the pluton as well as the enrichment of Mo and W at Mt Moliagul.

Clutter Reduction by Estimation of Echoes Direction of Arrival in Distributed Radar Sounders in Formation Flying

Spaceborne radar sounders are High Frequency (HF) / Very High Frequency (VHF) nadir-looking sensors devoted to subsurface investigations. Their data interpretation can be severely hindered by off-nadir surface clutter. Recent literature showed that the clutter suppression capabilities of this class of systems can be greatly enhanced by deploying an array of orbiting sensors in formation flight synthesizing a narrow radar antenna beam. In this paper, we assess the capability of distributed radar sounding to discriminate clutter from subsurface returns by exploiting Direction of Arrival (DOA) estimation techniques. This is achieved by first outlining an approach for designing and evaluating the distributed radar sounder DOA estimation performance as function of the radar system parameters (e.g, inter-sensor distance) and external noise factors such as ionospheric scintillations. Then, the theory is complemented by radar simulations of several acquisitions over Greenland assuming a variety of subsurface geometries. The simulations confirm that clutter discrimination through DOA estimation is a viable approach to further improve the array capability in disambiguation of subsurface echoes from surface ones.

Determining Grain Boundary Position and Geometry from EBSD Data: Limits of Accuracy.

As the feature size of crystalline materials gets smaller, the ability to correctly interpret geometrical sample information from electron backscatter diffraction (EBSD) data becomes more important. This paper uses the notion of transition curves, associated with line scans across grain boundaries (GBs), to correctly account for the finite size of the excitation volume (EV) in the determination of the geometry of the boundary. Various metrics arising from the EBSD data are compared to determine the best experimental proxy for actual numbers of backscattered electrons that are tracked in a Monte Carlo simulation. Consideration of the resultant curves provides an accurate method of determining GB position (at the sample surface) and indicates a significant potential for error in determining GB position using standard EBSD software. Subsequently, simple criteria for comparing experimental and simulated transition curves are derived. Finally, it is shown that the EV is too shallow for the curves to reveal subsurface geometry of the GB (i.e., GB inclination angle) for most values of GB inclination.

Geophysical Constraints on the Crustal Architecture of the Transtensional Warm Springs Valley Fault Zone, Northern Walker Lane, Western Nevada, USA

AbstractThe Walker Lane is a zone of distributed transtension where normal faults are overprinted by strike‐slip motion. We use two newly acquired, high‐resolution seismic reflection profiles and a reprocessed Consortium for Continental Reflection Profiling (COCORP) deep crustal reflection profile to assess the subsurface geometry of the Holocene‐active, transtensional Warm Springs Valley fault zone (WSVFZ) near Reno, Nevada, USA. Our multiscale observations extend to 12 km depth and suggest that the WSVFZ is more complex in the subsurface than implied by late Pleistocene surface fault traces. Two 4‐km‐long high‐resolution profiles image to a depth of ∼2 km and reveal moderately dipping reflections and truncations, some of which project to mapped scarps formed in late Pleistocene surfaces. The shallow lines are collocated with COCORP profile NV 08 along ∼40°N latitude. Reanalysis of the COCORP data reveals previously unidentified coherent reflections to a depth of ∼12 km and a previously mapped ∼30° west‐dipping fault at 8–12 km. From these seismic profiles, the WSVFZ is not a simple, subvertical fault zone extending through the entire seismogenic crust. Instead, the reflections are consistent with a zone of steeply and moderately dipping faults that simplify and steepen with depth before intersecting a mid‐crustal, low‐angle (∼25–30°) fault. The complex fault geometry of the WSVFZ implies that crustal shear is accommodated by a mix of dipping and subvertical faults in the transtensional northern Walker Lane. If so, transtensional fault zones may present challenges to paleoseismic and geodetic studies and require careful treatment when included in seismic hazard analyses.

Integrating pseudo-3D ERT and DEM to map the near-surface structures and morphology of the Piparo mud volcano, Trinidad

Mud volcanoes are often associated with violent eruptions which are hazardous and pose significant risks to lives and livelihoods. A pseudo-3D electrical resistivity tomography (ERT) method that is integrated with digital elevation model (DEM) is used to map and monitor the near-surface structures and morphology of a mud volcano in Piparo, Trinidad. The mudflow within the vent is characterized by very low resistivity of less than 2.5 Ωm. Most of the mudflow within the near-surface is located north of the main crater. A magnitude 6.9 earthquake, located near the northern coast of Venezuela, caused the vent to enlarge, the mud volcano to be displaced towards the north, the crater region to sink by 0.2 m, and the flanks to rise by 0.7 m. This study demonstrates that the integrated pseudo-3D ERT and DEM provide means to monitor the dynamic geological interplay between the surface and subsurface geometry and structures.

Linking the Surface and Subsurface in River Deltas—Part 1: Relating Surface and Subsurface Geometries

AbstractRiver deltas are densely populated regions of the world with vulnerable groundwater reserves. Contamination of these groundwater aquifers via saline water intrusion and pollutant transport is a growing threat due to both anthropogenic and climate changes. The arrangement and composition of subsurface sediment is known to have a significant impact on aquifer contamination; however, developing accurate depictions of the subsurface is challenging. In this work, we explore the relationship between surface and subsurface properties and identify the metrics most sensitive to different forcing conditions. To do so, we simulate river delta evolution with the rule‐based numerical model, DeltaRCM, and test the influence of input sand fraction and steady sea level rise (SLR) on delta evolution. From the model outputs, we measure a variety of surface and subsurface metrics chosen based on their applicability to imagery and modeling results. The Kullback‐Leibler (KL) divergence is then used to quantitatively gauge which metrics are most indicative of the imposed forcings. Both qualitative observations and the KL divergence analysis suggest that estimates of subsurface connectivity can be constrained using surface information. In particular, more variable shoreline roughness values and higher surface wetted fraction values correspond to increased subsurface connectivity. These findings complement traditional methods of estimating subsurface structure in river‐dominated delta systems and represent a step toward the identification of a direct link between surface observations and subsurface form.

Linking the Surface and Subsurface in River Deltas—Part 2: Relating Subsurface Geometry to Groundwater Flow Behavior

AbstractUnderstanding subsurface structure and groundwater flow in deltaic aquifers is essential for evaluating the vulnerability of groundwater resources in delta systems. Deltaic aquifers contain coarse‐grained paleochannels that preserve a record of former surface river channels as well as fine‐grained floodplain deposits. The distribution of these deposits and how they are interconnected control groundwater flow and contaminant transport. In this work, we link depositional environments of deltaic aquifers to stratigraphic (static) and flow and transport (dynamic) connectivity metrics. Numerical models of deltaic stratigraphy were generated using a reduced‐complexity numerical model (DeltaRCM) with different input sand fractions (ISF) and rates of sea‐level rise (SLR). The groundwater flow and advective transport behavior of these deltas were simulated using MODFLOW and MODPATH. By comparing the static and dynamic metrics calculated from these numerical models, we show that groundwater behavior can be predicted by particular aspects of the subsurface architecture, and that horizontal and vertical connectivity display different characteristics. We also evaluate relationships between connectivity metrics and two environmental controls on delta evolution: ISF and SLR rate. The results show that geologic setting strongly influences both static and dynamic connectivity in different directions. These results provide insights into quantitatively differentiated subsurface hydraulic behavior between deltas formed under different external forcing (ISF and SLR rate) and they are a potential link in using information from delta surface networks and depositional history to predict vulnerability to aquifer contamination.

Geophysical investigation of the Orle River fracture system in Igarra Township, southwestern Nigeria

The subsurface structure and propagation geometry of the fracture system controlling the Orle River Channel in the Igarra Township, Southwestern Nigeria, have been investigated via a multi-method geophysical survey. The goal was to delineate the nature, distribution, and spatial propagation geometry of the fracture system and evaluate its potential to serve as storage and distribution features for groundwater within adjoining areas. 2D electrical resistivity, total field magnetic, and co-planar loop conductivity measurements were collected along four traverses using an ABEM1000 Terrameter unit, a GEMS Magnetometer, and an EM34 Co-planar loop electromagnetic system. Traverses were established to run across and parallel to the river channel. The observed electrical resistivity field data were inverted for subsurface 2D resistivity structure using a commercially available 2.5D finite element modelling inversion software. Magnetic field intensity data and ground conductivity data were presented against station positions. Three subsurface layers were delineated at the survey site; (1) surficial humus-rich and wet top-soil, (2) a thin poorly developed weathering layer, and (3) the fresh bedrock which occurs as relatively shallow levels and often outcrops. The River channel is controlled by multiple fractures usually located at or near the contracts between contrasting rock types. Fracture dip is usually in the northerly direction but conjugates, dipping southwards also occur in the most northerly extremes of the imaged fracture system. Upstream the fracture path is wider and along with the imaged overburden, isjuxtaposed northwards of the channel axis suggesting a much broader river channel in the geologic past. Low aperture fractures imaged tangential to channel axis likely serve to funnel surface and groundwater from the channels to the surrounding areas. Where such low aperture fractures can be delineated, they offer the best chances of groundwater abstraction within adjoining areas, particularly during the dry season.

Space–time geological model of the Quaternary syntectonic fill of a foreland basin (Po basin, Northern Italy)

Quantitative models of large and widely populated Quaternary foreland basins, which host relevant resources like groundwater, heat, and fertile soils, may be computed if space–time models of the basin fills are available. The paper presents the combination of geometrical (space) and evolutionary (time) constraints to model Quaternary stratigraphy in the southern Po basin in Lombardy, disentangling the tectonic and climatic controls. Five Quaternary high-rank unconformities and seven intermediate-rank unconformities were recognised, down-traced, and correlated from surface to subsurface. Within this framework, the surface facies associations of high-, intermediate-, and low-rank stratigraphic units were compared to the litho-textural associations of their subsurface equivalents obtained from borehole logs and published geophysical images, to draw the best fitting model, which was constrained to the Quaternary geological evolution.Discrete, long-term thrusting, wrenching, folding, and extensional collapsing stages are inferred to have shaped the Gelasian, intra-Calabrian, Early–Middle Pleistocene, Middle–Late Pleistocene and Latest Pleistocene–Holocene high-rank composite unconformities. These surfaces are cut into synsedimentary, polyphasic, buried anticlines and correspond to conformable boundaries in the adjacent depocentres. The intermediate-rank, low-angle unconformities composing the high-rank ones, bound the stratigraphic units forming the overall regressive fill of the basin. These surfaces might correspond to steady tectonic uplift/subsidence periods when deposition was mostly controlled by Pleistocene climate changes since the Gelasian.The hierarchic space/time structure, translated into hard geometrical constraints (a fence diagram of cross-sections) and relative chronology of the stratigraphic relationships (intersect/erode, onlap and offset) allows computing the geological quantitative model using a potential field interpolation method. The model building phase and subsequent inspection permit critically analysing the current interpretations on the evolution of the Po basin during the Quaternary compared with other foreland settings. The 4-D model integrates the surface–subsurface geometry with the tectono-sedimentary evolution in a multi-scale structure that may be exported for hydrogeological, geothermal and environmental modelling.

Geoelectrical signatures for detecting water-bearing zones in a micro-watershed of granitic terrain from Southern India

This communication represents the application of geoelectrical signatures such as Dar-Zarrouk parameters, namely longitudinal conductance (S) and transverse resistance (T), which play a vital role in deciphering subsurface geometry and detect water-bearing zones in hard rock area. A micro-watershed is situated in granite terrain of Southern India, where the variation of bedrock depth and location of water-bearing zones are identified using these signatures, and achieved from the results of 47 Vertical Electrical Soundings (VESs) data. Pumping wells are recommended for the drilling, and the high yielded wells are confirmed in the water-bearing zones. The observed borehole lithology is very closely correlated with the interpreted geoelectric layer parameters. It indicates that the water-bearing zones exit in highly weathered zones having a resistivity range of 50–100 Ω-m, and 100–300 Ω-m in fractured granites at 12 to 22 m depths. The estimated T and S-values reveal the conductive mediums within the host granitic rock, whereas low S and comparatively high T might be indicative possible water-bearing/shear zone. The aquifer properties at the selective wells are compared to the aquifer potentiality. The high permeability values have been observed to confirm the presence of a shear zone, which is more productive for providing a good yielding bore well. Therefore, this approach could be replicated for groundwater exploration in any hard rock terrain.

Test on the Reliability of the Subsurface Fault Geometry Estimated by Deformed River Terraces Along the Bailang River, North Front of the Qilian Shan (North West China)

The subsurface fault geometry is the base for understanding a process of crust deformation and mountain building. Based on kinematic models for fault-related folds, a geomorphic method is recently applied to estimate the subsurface fault geometry, while the validation on its reliability is lacking. In this study, we surveyed a suit of river terrace surfaces across an active fold at the north front of the Qilian Shan. According to the deformation geometry of the terraces, the fold deformation is interpreted by a listric fault fold model, and based on this kinematic model, the fault geometry underlying the fold is estimated. In comparison between the estimated fault geometry and a seismic reflection profile, we found that the decollement depth and the back thrust are highly consistent with each other. Although some small fault bends or internal shearing cannot be estimated solely by the terrace deformation, the overall fault geometry is successfully revealed by the terrace deformation. Using this fault geometry and the terrace dating results, the region deformation kinematics are re-evaluated, which suggest that the dip slip (in a rate of 1.8 ± 0.4 mm/a) along the decollement is mainly accommodated by two structures, one is the blind-back-thrust fault within the piggy basin in a dip-slip rate of 0.9 ± 0.3 mm/a and another is the thrust and fold at the west portion of the Yumu Shan range.

Toward an Integrated and Sustainable Water Resources Management in Structurally-Controlled Watersheds in Desert Environments Using Geophysical and Remote Sensing Methods

Sustainable water resources management in desert environment has yet to be reached due to the limited hydrological datasets under such extreme arid conditions. In the Eastern Sahara, the tectonic activity associated with the opening of the Red Sea adds more complexity to developing sustainable water management by creating multiple aquifers within subsided half-grabens along the Red Sea extension. To overcome these difficulties, a two-fold approach is adopted including integrated remote sensing and geoelectrical methods using Wadi Al-Ambagi watershed in the Eastern Desert of Egypt as a test site. First, the total discharge is estimated as 15.7 × 106 m3 following the application of a uniform storm of 10 mm effective precipitation, which exceeds the storage capacity of existing mitigation measures (5.5 × 106 m3), and thus additional dams are required. Second, the subsurface geometry of alluvium and sedimentary aquifers, within subsided blocks in the Arabian–Nubian shield (ANS), is delineated using 1D direct current and 2D electrical-resistivity tomography (ERT). Findings indicate that significant thicknesses of more than 80 m of permeable sedimentary units occur within the subsided blocks. Therefore, the scarce water resources can be managed by controlling the flash floods and suggesting proper dam sites at the location of thick alluvium and sedimentary rocks, where aquifers can be recharged representing a sustainable source for freshwater. The proposed approach is transferable and can be applied in similar arid rift-related watersheds in Saudi Arabia and worldwide.

Analysis of gravity anomaly decomposition and depth to basement, case study: Cenozoic Bogor Basin, Indonesia.

The Cenozoic Bogor Basin situated on the SE edge part of Sunda shelf is presumed to have hydrocarbon potential on its turbidite deep-water play. The deep subsurface geometry of this basin may play an important role, yet unexplored, to the hydrocarbon exploration. The gravity method is advantageous to illuminate the subsurface structure on the arbitrary depth and various sources. Eight-hundred-thirty-eight points of ground-based gravity survey were collected on roughly one kilometer spacing in the North of Bandung to Pamanukan region covering the Cenozoic Bogor Basin on the Subang-Purwakarta segment to generate complete Bouguer anomaly (CBA) map. This study examines the two robust methods of gravity anomaly decomposition (i.e., polynomial trend surface and upward continuation) by using multiple parameters to match the geological background. Radially averaged power spectrum was used to estimate the depth of anomalous source corresponds to the top of basement layer and resulting 4 km basement depth in North West Java basin and below 1 km on Tangkuban Parahu volcanic zone. The resulting estimated depth of Cenozoic Bogor Basin was evaluated by the matched a priori published data on those two areas and revealed the depth of depocenter (deepest sediment) on the Cenozoic Bogor Basin is up to 9 km, quite a deeper extent than previously assumed.

Unveiling buried aeolian landscapes: reconstructing a late Holocene dune environment using 3D ground‐penetrating radar

ABSTRACTAcross the UK, sandy beaches and dunes protect coastal infrastructure from waves and extreme water levels during large‐scale storms, while providing important habitats and recreational opportunities. Understanding their long‐term evolution is vital in managing their condition in a changing climate. Recently, ground‐penetrating radar (GPR) methods have grown in popularity in geomorphological applications, yielding centimetre‐scale resolution images of near‐surface stratigraphy and structure, thus allowing landscape evolution to be reconstructed. Additionally, abrupt changes in palaeo‐environments can be visualized in three dimensions. Although often complemented by core data, GPR allows interpretations to be extended into areas with minimal ground‐truth control. Nonetheless, GPR data interpretation can be non‐intuitive and ambiguous, and radargrams may not initially resemble the expected subsurface geometry. Interpretation can be made yet more onerous when handling the large 3D data volumes that are facilitated with modern GPR technology. Here we describe the development of novel semi‐automated GPR feature‐extraction tools, based on ‘edge detection’ and ‘thresholding’ methods, which detect regions of increased GPR reflectivity which can be applied to aid in the reconstruction of a range Quaternary landscapes. Since reflectivity can be related to lithological and/or pore fluid changes, the 3D architecture of the palaeo‐landscape can be reconstructed from the features extracted from a geophysical dataset. We present 500 MHz GPR data collected over a buried Holocene coastal dune system in North Wales, UK, now reclaimed for use as an airfield. Core data from the site, reaching a maximum depth 2 m, suggest rapid vertical changes from sand to silty‐organic units, and GPR profiles suggest similar lateral complexity. By applying thresholding methods to GPR depth slices, these lateral complexities are effectively and automatically mapped. Furthermore, automatic extraction of the local reflection power yields a strong correlation with the depth variation of organic content, suggesting it is a cause of reflectivity contrast. GPR‐interpolated analyses away from core control thus offer a powerful proxy for parameters derived from invasive core logging. The GPR data collected at Llanbedr airfield highlight a complex dune system to a depth of 2.8 m, probably deposited in several phases over ~700 years, similar to elsewhere in North Wales.

Slip Rate of the Danghe Nan Shan Thrust Fault from 10Be Exposure Dating of Folded River Terraces: Implications for the Strain Distribution in Northern Tibet

AbstractThe northeastward motion of the Tibetan Plateau along the Altyn Tagh strike‐slip fault causes thrust faulting in three parallel mountain ranges (Qilian Shan, Daxue Shan, and Danghe Nan Shan) in the plateau interior, and leads to NNE‐directed crustal shortening and plateau growth. While slip rates at the plateau margin (i.e., along the Qilian Shan and the Altyn Tagh fault) are well constrained, rates of thrust faulting and the strain distribution in the plateau interior remain poorly resolved. Here, we use field investigations, a high‐resolution DEM, and 10Be exposure dating to quantify the shortening rate across the Danghe Nan Shan thrust fault from fluvial terraces, which are deformed by a growing NNE‐vergent anticline. 10Be exposure ages from two terrace levels range from 70 ± 5 to 92 ± 7 ka. When combined with uplift values of 37–68 m along the fold hinge, the 10Be ages yield a mean uplift rate of 0.6 ± 0.2 mm/year. Using the cross‐sectional area of the fold and the subsurface geometry of the listric thrust fault, we obtain a shortening rate of 0.8 ± 0.2 mm/year, which is consistent with the rate of elastic strain accumulation recorded by GPS data. Together with published fault slip rates and GPS data, our results indicate that northern Tibet experiences NNE‐directed shortening at a rate of ∼5 mm/year between the Qaidam Basin and the Hexi Corridor. In the plateau interior, this shortening is accommodated by several range‐bounding thrust faults and closely coupled with the eastward decrease in the slip rate of the Altyn Tagh fault.

The interplay of stratal and vertical migration pathways in shallow hydrocarbon plumbing systems

AbstractHydrocarbon plumbing systems have been extensively documented in the past two decades using high‐resolution 3D seismic data, exploiting the ability of seismic imaging techniques to reveal the subsurface geometry of gas charged sediments. In this paper, we present a detailed study of a hydrocarbon plumbing system from the South China Sea, that involves both vertical and lateral (stratal) hydrocarbon migration in Miocene to Recent clastic sediments that comprise multilayer stacking of thinly layered clays, silts and sands. We show that a transtensive fault system that provides lateral seal for fault‐dip traps of deep Miocene reservoirs, and offers a vertical pathway for migration to shallower silty units. These silty units in turn form a ‘spillway’ in a regional, northward migration path. This path involves filling each shallow fault‐dip trap to spill point towards the fault tips, with stratal migration forced around the outer flanks of the fault‐related folds. Successive fill‐to‐spill leads to a continuous trail of amplitude anomalies that merge into a continuous, larger, gas‐charged anomaly pattern. The migrating gas finally accumulates within a zone bounded by a large boundary fault with full juxtaposition seal. The pattern of anomaly distribution suggests that the hydrocarbon migration has been active in the Late Pleistocene and is probably continuing at the present day. Hence this plumbing system may be one of very few examples described to date in which dynamic hydrocarbon migration pathways have been directly imaged by seismic data.

Inferring the Subsurface Geometry and Strength of Slow‐Moving Landslides Using 3‐D Velocity Measurements From the NASA/JPL UAVSAR

AbstractThe hazardous impact and erosive potential of slow‐moving landslides depends on landslide properties including velocity, size, and frequency of occurrence. However, constraints on size, in particular, subsurface geometry, are lacking because these types of landslides rarely fully evacuate material to create measurable hillslope scars. Here, we use pixel offset tracking with data from the NASA/JPL Uninhabited Aerial Vehicle Synthetic Aperture Radar to measure the three‐dimensional surface deformation of 134 slow‐moving landslides in the northern California Coast Ranges. We apply volume conservation to infer the actively deforming thickness, volume, geometric scaling, and frictional strength of each landslide. These landslides move at average rates between ∼0.1–2 m/yr and have active areas of ∼6.10 × 103–2.35 × 106 m2, inferred mean thicknesses of ∼1.1–25 m, and volumes of ∼7.01 × 103–9.75 × 106 m3. The best fit volume‐area geometric scaling exponent is γ ∼ 1.2–1.5, indicating that these landslides fall between typical soil and bedrock landslide scaling. A rollover in the scaling relationship suggests that the largest landslide complexes in our data set become large primarily by increasing in area rather than thickness. In addition, the slow‐moving landslides display scale‐dependent frictional strength, such that large landslide tend to be weaker than small landslides. This decrease in frictional strength with landslide size is likely because larger landslides are composed of higher proportions of weak material. Our work shows how state of the art remote sensing techniques can be used to better understand landslide processes and quantify their contribution to landscape evolution and hazards to human safety.

Modeling of Empirical Transfer Functions with 3D Velocity Structure

ABSTRACTEmpirical transfer functions (ETFs) between seismic records observed at the surface and depth represent a powerful tool to estimate site effects for earthquake hazard analysis. However, conventional modeling of site amplification, with assumptions of horizontally polarized shear waves propagating vertically through 1D layered homogeneous media, often poorly predicts the ETFs, particularly, in which large lateral variations of velocity are present. Here, we test whether more accurate site effects can be obtained from theoretical transfer functions (TTFs) extracted from physics-based simulations that naturally incorporate the complex material properties. We select two well-documented downhole sites (the KiK-net site TKCH05 in Japan and the Garner Valley site, Garner Valley Downhole Array, in southern California) for our study. The 3D subsurface geometry at the two sites is estimated by means of the surface topography near the sites and information from the shear-wave profiles obtained from borehole logs. By comparing the TTFs to ETFs at the selected sites, we show how simulations using the calibrated 3D models can significantly improve site amplification estimates as compared to 1D model predictions. The primary reason for this improvement in 3D models is redirection of scattering from vertically propagating to more realistic obliquely propagating waves, which alleviates artificial amplification at nodes in the vertical-incidence response of corresponding 1D approximations, resulting in improvement of site effect estimation. The results demonstrate the importance of reliable calibration of subsurface structure and material properties in site response studies.

Measurements of Iceberg Melt Rates Using High‐Resolution GPS and Iceberg Surface Scans

AbstractIncreasing freshwater input to the subpolar North Atlantic through iceberg melting can influence fjord‐scale to basin‐scale ocean circulation. However, the magnitude, timing, and distribution of this freshwater have been challenging to quantify due to minimal direct observations of subsurface iceberg geometry and melt rates. Here we present novel in situ methods capturing iceberg change at high‐temporal and ‐spatial resolution using four high‐precision GPS units deployed on two large icebergs (>500 m length). In combination with measurements of surface and subsurface geometry, we calculate iceberg melt rates between 0.10 and 0.27 m/d over the 9‐day survey. These melt rates are lower than those proposed in previous studies, likely due to using individual subsurface iceberg geometries in calculations. In combining these new measurements of iceberg geometry and melt rate with the broad spatial coverage of remote sensing, we can better predict the impact of increasing freshwater injection from the Greenland Ice Sheet.

The 18 March 2020 M 5.7 Magna, Utah, Earthquake: Strong-Motion Data and Implications for Seismic Hazard in the Salt Lake Valley

Abstract The 2020 oblique normal-faulting M 5.7 Magna mainshock has provided the best dataset of recorded strong ground motions for an earthquake within the Wasatch Front region, Utah, and the larger Basin and Range Province. We performed a preliminary evaluation of the strong motion and broadband data from this earthquake and compared the data with the Next Generation Attenuation - West2 Project (NGA-West2) ground-motion models (GMMs). The highest horizontal peak ground acceleration (PGA) recorded was 0.43g (geometric mean of the two horizontal components) at a station located above the rupture plane at a rupture distance of 8 km. Eleven stations recorded PGAs >0.20g. Most of these stations are located on the deep sedimentary deposits within the Salt Lake Valley, and all are at rupture distances <20 km. The data compare favorably with the NGA-West2 GMMs, although the expected variability was observed. PGAs exceed the GMM predictions at the closest distances for the source model that we used. The area of the strongest ground shaking encompassed the town of Magna, where some of the heaviest damage occurred. A significant implication of the 2020 Magna earthquake for seismic hazards in the Salt Lake Valley arises from the possibility that this earthquake occurred on the Salt Lake City segment of the Wasatch fault. If so, then the dip of this fault segment must decrease with depth to ≤30°–35°, as proposed by Pang et al. (2020)—at least along the northern part of the segment where the earthquake occurred. Because of the lack of information about the subsurface geometry of the Wasatch fault zone, modeling of this fault zone in seismic hazard analyses has assumed a moderate dip of 50°±15°. Assuming a more shallowly dipping fault results in higher estimates of ground shaking in future large earthquakes on this fault. Alternative interpretations of the Magna earthquake are that it occurred (1) on an auxiliary fault within the Wasatch fault zone or (2) on a listric section of the northern Salt Lake City segment that is not representative of the geometry of the whole fault segment.