Was the M<Subscript>w</Subscript> 7.5 1952 Kern County,California, earthquake induced (or triggered)?

Several recent studies have presented evidence that significant induced earthquakes occurred in a number of oil-producing regions during the early and mid-twentieth century related to either production or wastewater injection. We consider whether the 21 July 1952 M_w 7.5 Kern County earthquake might have been induced by production in the Wheeler Ridge oil field. The mainshock, which was not preceded by any significant foreshocks, occurred 98 days after the initial production of oil in Eocene strata at depths reaching 3 km, within ~1 km of the White Wolf fault (WWF). Based on this spatial and temporal proximity, we explore a potential causal relationship between the earthquake and oil production. While production would have normally be expected to have reduced pore pressure, inhibiting failure on the WWF, we present an analytical model based on industry stratigraphic data and best estimates of parameters whereby an impermeable splay fault adjacent to the main WWF could plausibly have blocked direct pore pressure effects, allowing the poroelastic stress change associated with production to destabilize the WWF, promoting initial failure. This proof-of-concept model can also account for the 98-day delay between the onset of production and the earthquake. While the earthquake clearly released stored tectonic stress, any initial perturbation on or near a major fault system can trigger a larger rupture. Our proposed mechanism provides an explanation for why significant earthquakes are not commonly induced by production in proximity to major faults.     

Comparison of spectral and spatial methods for a Moho recovery from gravity and vertical gravity-gradient data

In global studies investigating the Earth’s lithospheric structure, the spectral expressions for the gravimetric forward and inverse modeling of the global gravitational and crustal structure models are preferably used, because of their numerical efficiency. In regional studies, the applied numerical schemes typically utilize the expressions in spatial form. Since the gravity-gradient observations have a more localized support than the gravity measurements, the gravity-gradient data (such as products from the Gravity field and steady-state Ocean Circulation Explorer - GOCE - gravity-gradiometry satellite mission) could preferably be used in regional studies, because of reducing significantly the spatial data-coverage required for a regional inversion or interpretation. In this study, we investigate this aspect in context of a regional Moho recovery. In particular, we compare the numerical performance of solving the Vening Meinesz-Moritz’s (VMM) inverse problem of isostasy in spectral and spatial domains from the gravity and (vertical) gravity-gradient data. We demonstrate that the VMM spectral solutions from the gravity and gravity-gradient data are (almost) the same, while the VMM spatial solutions differ from the corresponding spectral solutions, especially when using the gravity-gradient data. The validation of the VMM solutions, however, reveals that the VMM spatial solution from the gravity-gradient data has a slightly better agreement with seismic models. A more detailed numerical analysis shows that the VMM spatial solution formulated for the gravity gradient is very sensitive to horizontal spatial variations of the vertical gravity gradient, especially in vicinity of the computation point. Consequently, this solution provides better results in regions with a relatively well-known crustal structure, while suppressing errors caused by crustal model uncertainties from distant zones. Based on these findings we argue that the gravity-gradient data are more suitable than the gravity data for a regional Moho recovery.     

Experimental study of deformable connection consisting of buckling‐restrained brace and rubber bearings to connect floor system to lateral force resisting system

This paper presents experimental and numerical studies of a full‐scale deformable connection used to connect the floor system of the flexible gravity load resisting system to the stiff lateral force resisting system (LFRS) of an earthquake‐resistant building. The purpose of the deformable connection is to limit the earthquake‐induced horizontal inertia force transferred from the floor system to the LFRS and, thereby, to reduce the horizontal floor accelerations and the forces in the LFRS. The deformable connection that was studied consists of a buckling‐restrained brace (BRB) and steel‐reinforced laminated low‐damping rubber bearings (RB). The test results show that the force–deformation responses of the connection are stable, and the dynamic force responses are larger than the quasi‐static force responses. The BRB+RB force–deformation response depends mainly on the BRB response. A detailed discussion of the BRB experimental force–deformation response is presented. The experimental results show that the maximum plastic deformation range controls the isotropic hardening of the BRB. The hardened BRB force–deformation responses are used to calculate the overstrength adjustment factors. Details and limitations of a validated, accurate model for the connection force–deformation response are presented. Numerical simulation results for a 12‐story reinforced concrete wall building with deformable connections show the effects of including the RB in the deformable connection and the effect of modeling the BRB isotropic hardening on the building seismic response. Copyright © 2016 John Wiley & Sons, Ltd.     

Spillage sedimentation on large river floodplains

Active deposition across the floodplains of large rivers arises through a variety of processes; collectively these are here termed ‘spillage sedimentation’. Three groups of 11 spillage sedimentation styles are identified and their formative processes described. Form presences on large river floodplains show different combinations of active spillage styles. Only some large floodplains have prominent levees; some have coarse splays; many have accessory channel dispersion and reworking, while still‐water sedimentation in lacustrine environments dominates some lower reaches. Infills are also commonly funnelled into prior, and often linear, negative relief forms relating to former migration within the mainstream channel belt. Shuttle Radar Topography Mission (SRTM) and Landsat 8 data are used to map spillage form types and coverage along a 1700 km reach of the Amazon that has an active floodplain width of up to 110 km with a systematic character transformation down‐valley. Spillage forms associated directly with mainstream processes rarely account for more than 5% of the floodplain deposits. There is a marked decrease in floodplain point bar complexes (PBC) over 1700 km downstream (from 34% to 5%), and an increase in the prevalence of large water bodies (2% to 37%) and accompanying internal crevasses and deltas (0% to 5%). Spillage sedimentation is likely within the negative relief associated with these forms, depending on mainstream sediment‐laden floodwater inputs. Spillage style dominance depends on the balance between sediment loadings, hydrological sequencing, and morphological opportunity. Down‐river form sequences are likely to follow gradient change, prior up‐river sediment sequestration and the altered nature of spilled loads, but also crucially, local floodplain relief and incident water levels and velocities at spillage times. Considering style distribution quantitatively, as a spatially distributed set of identifiable forms, emphasizes the global variety to spillage phenomena along and between large rivers. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Predicting gravel bed river response to environmental change: the strengths and limitations of a regime‐based approach

Rivers respond to environmental changes such as climate shifts, land use changes and the construction of hydro‐power dams in a variety of ways. Often there are multiple potential responses to any given change. Traditionally, potential stream channel response has been assessed using simple, qualitative frameworks based largely on professional judgement and field experience, or using some form of regime theory. Regime theory represents an attempt to use a physically based approach to predict the configuration of stable channels that can transport the imposed sediment supply with the available discharge. We review the development of regime theory, and then present a specific regime model that we have created as a stand‐alone computer program, called the UBC Regime Model (UBCRM). UBCRM differs from other regime models in that it constrains its predictions using a bank stability criterion, as well as a pattern stability criterion; it predicts both the stable channel cross‐sectional dimensions as well as the number of anabranches that the stream must have in order to establish a stable channel pattern. UBCRM also differs from other models in that it can be used in a stochastic modelling mode that translates uncertainty in the input variables into uncertainty in the predicted channel characteristics. However, since regime models are fundamentally based on the concept of grade, there are circumstances in which the model does not perform well. We explore the strengths and weaknesses of the UBCRM in this paper, and we attempt to illustrate how the UBCRM can be used to augment the existing qualitative frameworks, and to help guide professionals in their assessments. Copyright © 2016 John Wiley & Sons, Ltd.     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.     

Observations of preburst heating and magnetic field changes in a coronal loop at 20 cm wavelength

The Very Large Array (VLA) has been used at 20 cm wavelength to study the evolution of a burst loop with 4 resolution on timescales as short as 10 s. The VLA observations show that the coronal loop began to heat up and change its structure about 15 min before the eruption of two impulsive bursts. The first of these bursts occurred near the top of the loop that underwent preburst heating, while the second burst probably occurred along the legs of an adjacent loop. These observations evoke flare models in which coronal loops twist, develop magnetic instabilities and then erupt. We also combine the VLA observations with GOES X-ray data to derive a peak electron temperature of T _e = 2.5 × 10⁷ K and an average electron density of N _e 1 × 10¹⁰ cm^–3 in the coronal loop during the preburst heating phase.     

A method for uncertainty constraint of catchment discharge and phosphorus load estimates

River discharge and nutrient measurements are subject to aleatory and epistemic uncertainties. In this study, we present a novel method for estimating these uncertainties in colocated discharge and phosphorus (P) measurements. The “voting point”‐based method constrains the derived stage‐discharge rating curve both on the fit to available gaugings and to the catchment water balance. This helps reduce the uncertainty beyond the range of available gaugings and during out of bank situations. In the example presented here, for the top 5% of flows, uncertainties are shown to be 139% using a traditional power law fit, compared with 40% when using our updated “voting point” method. Furthermore, the method is extended to in situ and lab analysed nutrient concentration data pairings, with lower uncertainties (81%) shown for high concentrations (top 5%) than when a traditional regression is applied (102%). Overall, for both discharge and nutrient data, the method presented goes some way to accounting for epistemic uncertainties associated with nonstationary physical characteristics of the monitoring site.     

Strong variation in weathering of layered rock maintains hillslope‐scale strength under high precipitation

The evolution of volcanic landscapes and their landslide potential are both dependent upon the weathering of layered volcanic rock sequences. We characterize critical zone structure using shallow seismic V_p and V_s profiles and vertical exposures of rock across a basaltic climosequence on Kohala peninsula, Hawai’i, and exploit the dramatic gradient in mean annual precipitation (MAP) across the peninsula as a proxy for weathering intensity. Seismic velocity increases rapidly with depth and the velocity–depth gradient is uniform across three sites with 500–600 mm/yr MAP, where the transition to unaltered bedrock occurs at a depth of 4 to 10 m. In contrast, velocity increases with depth less rapidly at wetter sites, but this gradient remains constant across increasing MAP from 1000 to 3000 mm/yr and the transition to unaltered bedrock is near the maximum depth of investigation (15–25 m). In detail, the profiles of seismic velocity and of weathering at wet sites are nowhere monotonic functions of depth. The uniform average velocity gradient and the greater depths of low velocities may be explained by the averaging of velocities over intercalated highly weathered sites with less weathered layers at sites where MAP > 1000 mm/yr. Hence, the main effect of climate is not the progressive deepening of a near‐surface altered layer, but rather the rapid weathering of high permeability zones within rock subjected to precipitation greater than ~1000 mm/yr. Although weathering suggests mechanical weakening, the nearly horizontal orientation of alternating weathered and unweathered horizons with respect to topography also plays a role in the slope stability of these heterogeneous rock masses. We speculate that where steep, rapidly evolving hillslopes exist, the sub‐horizontal orientation of weak/strong horizons allows such sites to remain nearly as strong as their less weathered counterparts at drier sites, as is exemplified by the 50°–60° slopes maintained in the amphitheater canyons on the northwest flank of the island. Copyright © 2017 John Wiley & Sons, Ltd.