Changes in soil conditions before and after earthquakes at a repetitive soil liquefaction site in Taiwan

Beishih Village of Hsinhua Township in southern Taiwan is a unique location for studying soil liquefaction. Soil liquefaction was observed at the same site after earthquakes in 1946, 2010, and 2016, each of which had a Richter magnitude greater than six. This recurrence provides an opportunity for analyzing soil condition variations resulting from soil liquefaction. Seismic data sets were collected in 2011, 2014, 2016, and 2017. We used seismic refraction tomography and the multichannel analysis of surface waves to estimate P- and S-wave velocities. In S-wave velocity profiles, low shear velocity zones were located beneath sand volcanoes shortly after two earthquakes and disappeared 4 years after a 2010 earthquake. However, the P-wave velocity is less sensitive to soil condition changes, possibly because groundwater obscures the effect of soil liquefaction on velocity profiles. In addition, we used seismic wave velocities to determine the importance of soil properties such as Poisson’s ratio, shear modulus, and porosity to identify the cause of the low shear velocity zone. Notably, although porosity decreased after soil grain rearrangement, sand and clay mixing increased the Poisson’s ratio, reducing the shear modulus of the soil. In addition, a soil layer between 2 and 7 m and a deeper layer below 10 m that resulted in sand volcanoes were both liquefied. We also considered how the evaluation of soil liquefaction potential could be affected by long-term variations in soil conditions and changes resulting from liquefaction. The factor of safety was used to evaluate the liquefaction potential of the site. The results revealed that the assessment conducted long after the earthquake underestimated risk because the soil developed shear strength after the earthquake.

Generation of Gravity-Capillary Wind Waves by Instability of a Coupled Shear-Flow

The theory of surface wave generation, in viscous flows, is modified by replacing the linear-logarithmic shear velocity profile, in the air, with a model which links smoothly the linear and logarithmic layers through the buffer layer. This profile includes the effects of air flow turbulence using a damped mixing-length model. In the water, an exponential shear velocity profile is used. It is shown that this modified and coupled shear-velocity profile gives a better agreement with experimental data than the coupled linear-logarithmic, non smooth profile, (in the air)–exponential profile (in the water), widely used in the literature. We also give new insights on retrograde modes that are Doppler shifted by the surface velocity at the air-sea interface, namely on the threshold value of the surface current for the occurrence of a second unstable mode.

Experimental analysis of the effect of bed-load movement on flow hydraulic characteristics

In this study, the effect of bed-load movement on mean flow characteristics was evaluated in two rigid rectangular flumes. The experiments consisted of creating flow conditions carrying sediments with mean diameters of D50 = 0.5, 0.6, and 2.84 mm over both smooth and rough beds. Various sediment concentrations were injected at the upstream end of the flume at non-deposit injection rates to study the effect of various concentrations on flow resistance. The effect of sediment movement on flow resistance was examined by comparing the results with those of clear water flows (without sediment injection on both smooth and rough beds). The results showed that the sediment transport in maximum injection rate may increase the friction factor up to 50 and 58 percent for smooth bed, and up to about 75 and 80 percent in rough bed with mean diameter of 0.5 and 0.6 mm. Besides, for D50 = 2.84 mm, the friction factor decreased in smooth bed and increased up to 50 percent in rough bed. In general, it can be concluded that bed-load transport can be increased by the flow friction factor. The results also showed that bed-loads may decrease the average velocity and increase shear velocity with extraction of momentum from the flow, which both of mentioned factors may increase the flow friction factor. Raising the bed-load concentration in the flow may increase the elevation of the friction factor, approaching a constant value after reaching to the aggregation threshold and generation of bed forms.

A dimensionless framework for predicting submarine fan morphology

Limited observations of active turbidity currents at field scales challenges the development of theory that links flow dynamics to the morphology of submarine fans. Here we offer a framework for predicting submarine fan morphologies by simplifying critical environmental forcings such as regional slopes and properties of sediments, through densimetric Froude (ratio of inertial to gravitational forces) and Rouse numbers (ratio of settling velocity of sediments to shear velocity) of turbidity currents. We leverage a depth-average process-based numerical model to simulate an array of submarine fans and measure rugosity as a proxy for their morphological complexity. We show a systematic increase in rugosity by either increasing the densimetric Froude number or decreasing the Rouse number of turbidity currents. These trends reflect gradients in the dynamics of channel migration on the fan surface and help discriminate submarine fans that effectively sequester organic carbon rich mud in deep ocean strata.

Lifting of Tribocharged Grains by Martian Winds

It is a long-standing open question whether electrification of wind-blown sand due to tribocharging—the generation of electric charges on the surface of sand grains by particle–particle collisions—could affect rates of sand transport occurrence on Mars substantially. While previous wind tunnel experiments and numerical simulations addressed how particle trajectories may be affected by external electric fields, the effect of sand electrification remains uncertain. Here we show, by means of wind tunnel simulations under air pressure of 20 mbar, that the presence of electric charges on the particle surface can reduce the minimal threshold wind shear velocity for the initiation of sand transport, u *ft, significantly. In our experiments, we considered different samples, a model system of glass beads as well as a Martian soil analog, and different scenarios of triboelectrification. Furthermore, we present a model to explain the values of u *ft obtained in the wind tunnel that is based on inhomogeneously distributed surface charges. Our results imply that particle transport that subsides, once the wind shear velocity has fallen below the threshold for sustained transport, can more easily be restarted on Mars than previously thought.

The evolution of basal mantle structure in response to supercontinent aggregation and dispersal

Seismic studies have revealed two Large Low-Shear Velocity Provinces (LLSVPs) in the lowermost mantle. Whether these structures remain stable over time or evolve through supercontinent cycles is debated. Here we analyze a recently published mantle flow model constrained by a synthetic plate motion model extending back to one billion years ago, to investigate how the mantle evolves in response to changing plate configurations. Our model predicts that sinking slabs segment the basal thermochemical structure below an assembling supercontinent, and that this structure eventually becomes unified due to slab push from circum-supercontinental subduction. In contrast, the basal thermochemical structure below the superocean is generally coherent due to the persistence of a superocean in our imposed plate reconstruction. The two antipodal basal thermochemical structures exchange material several times when part of one of the structures is carved out and merged with the other one, similarly to “exotic” tectonic terranes. Plumes mostly rise from thick basal thermochemical structures and in some instances migrate from the edges towards the interior of basal thermochemical structures due to slab push. Our results suggest that the topography of basal structures and distribution of plumes change over time due to the changing subduction network over supercontinent cycles.

Determining seismic shear-velocity from ambient noise sources at regional and local scales

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction. Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions. We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells. Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>

Using ambient seismic noise to study temporal and spatial surface wave velocity structures and ambient noise field characteristics of central South Island, New Zealand

<p>Ambient seismic noise is used to examine the spatial and temporal surface wave velocity structures and ambient seismic noise fields in the vicinity of different fault zone environments. This study focuses on two distinct regions of central South Island, New Zealand. The Canterbury Plains is a sedimentary basin with many minor faults, which was considered to have low seismic hazard prior to the 2010 – 2011 Canterbury earthquake sequence. We focus on the time period immediately following the 2010 Darfield earthquake, which ruptured the previously unmapped Greendale Fault. The second region of interest is the central Southern Alps. The locked portion of the Alpine Fault currently poses one of the largest seismic hazards for New Zealand. The wealth of data from both permanent and temporary seismic deployments in these regions make them ideal areas in which to assess the effectiveness of ambient noise for velocity modelling in regions surrounding faults at different stages of their seismic cycles. Temporal velocity changes are measured following the Mw 7.1 Darfield earthquake of 4 September 2010 in the Canterbury Plains. Nine-component cross-correlations are computed from temporary and permanent seismic stations lying on and surrounding the Greendale Fault. Using the Moving-Window Cross-Spectral method, surface wave velocity changes are calculated for the four months immediately following the earthquake until 10 January 2011, for 0.1 — 1.0 Hz. An average increase in seismic velocity of 0.14 ± 0.04 % is determined throughout the region, providing the first such estimate of postseismic relaxation rates in Canterbury. Depth analyses further showed that velocity changes are confined to the uppermost 5 km of the subsurface and we attribute this to postseismic relaxation via crack-healing of the Greendale Fault and throughout the surrounding region. Rayleigh and Love wave dispersion is examined throughout the Canterbury region. Multi-component cross-correlation functions are analysed for group and phase dispersion curves. These are inverted using frequency-time analysis for 2-D phase and group velocity maps of Rayleigh and Love waves. A high-velocity zone to the southeast of the region coincides with volcanic rocks of Banks Peninsula. Dispersion curves generated from the surface wave tomography are further inverted for one-dimensional shear velocity profiles. These models show a thin, low-velocity near surface layer consistent with the basin sediments, which thins towards the foothills of the Southern Alps. A near-surface damage zone is identified along the length of the Greendale Fault, with consistent reduced Vs velocities to depth of up to 5 km. Surface and shear wave velocity maps are computed for the central Southern Alps to image the seismic structure of the region. Tomographic surface maps at periods of 5 – 12 s are produced from dispersion measurements of three-component cross-correlation functions. At periods of 5 – 8 s a strong NE-SW trending velocity contrast highlights the Alpine Fault. One-dimensional shear velocity models, computed from the surface wave maps, are in agreement with previous models produced by other conventional methods. An analysis of surface wave amplitudes through signal-to-noise ratios of cross-correlations reveals strong directional effects. Calculated signal-to-noise ratios are up to eight times higher for surface waves travelling north-west than for waves travelling to the south or east. We attribute this to a combination of more energetic ocean wave signals from the Southern Ocean compared to the Tasman Sea.</p>