SEARCH · Search USGS

Results for “Nonlinear Processes in Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

USGS2011

Scaling of peak flows with constant flow velocity in random self-similar networks

A methodology is presented to understand the role of the statistical self-similar topology of real river networks on scaling, or power law, in peak flows for rainfall-runoff events. We created Monte Carlo generated sets of ensembles of 1000 random self-similar networks (RSNs) with geometrically distributed interior and exterior generators having parameters p i and p e , respectively. The parameter values were chosen to replicate the observed topology of real river networks. We calculated flow hydrographs in each of these networks by numerically solving the link-based mass and momentum conservation equation under the assumption of constant flow velocity. From these simulated RSNs and hydrographs, the scaling exponents β and φ characterizing power laws with respect to drainage area, and corresponding to the width functions and flow hydrographs respectively, were estimated. We found that, in general, φ > β, which supports a similar finding first reported for simulations in the river network of the Walnut Gulch basin, Arizona. Theoretical estimation of β and φ in RSNs is a complex open problem. Therefore, using results for a simpler problem associated with the expected width function and expected hydrograph for an ensemble of RSNs, we give heuristic arguments for theoretical derivations of the scaling exponents β (E) and φ (E) that depend on the Horton ratios for stream lengths and areas. These ratios in turn have a known dependence on the parameters of the geometric distributions of RSN generators. Good agreement was found between the analytically conjectured values of β (E) and φ (E) and the values estimated by the simulated ensembles of RSNs and hydrographs. The independence of the scaling exponents φ (E) and φ with respect to the value of flow velocity and runoff intensity implies an interesting connection between unit hydrograph theory and flow dynamics. Our results provide a reference framework to study scaling exponents under more complex scenarios of flow dynamics and runoff generation processes using ensembles of RSNs.

Nonlinear Processes in Geophysics↗

USGS2018

Statistical approach to neural network imaging of karst systems in 3D seismic reflection data

The current lack of a robust, standardized technique for geophysical mapping of karst systems can be attributed to both the complexity of the environment and prior technological limitations. Abrupt lateral variations in physical properties that are inherent to karst systems generate significant geophysical noise, challenging conventional seismic signal processing and interpretation. Modern application of neural networks to multi-attribute seismic interpretation now provide a semiautomated method for identifying and leveraging the nonlinear relationships exhibited among seismic attributes. The ambiguity generally associated with designing neural networks for seismic object detection can be reduced via statistical analysis of the extracted attribute data. A data-driven approach to selecting the appropriate set of input seismic attributes, as well as the locations and minimum number of training examples, provides a more objective and computationally efficient method for identifying karst systems using reflection seismology. This statistically optimized neural network technique is thoroughly demonstrated using three-dimensional seismic reflection data collected from the southeastern portion of the Florida carbonate platform. Several dimensionality reduction methods are applied and the resulting karst probability models are evaluated relative to one another based on both quantitative and qualitative criteria. Comparing the preferred model, using quadratic discriminant analysis, to previously available seismic object detection workflows demonstrates the karst-specific nature of the tool. Results suggest that the karst multi-attribute workflow presented is capable of approximating the structural boundaries of karst systems with more accuracy and efficiency than a human counterpart or previously presented seismic interpretation schemes. This objective technique, using solely three-dimensional seismic reflection data, likely represents the most practical approach to mapping karst systems for subsequent hydrogeological modeling.

Interpretation↗

USGS2026

Thinking outside the rocks: Subsurface water storage, topography, and land cover are key modulators of large-scale riverine dissolved silicon dynamics

Riverine dissolved silicon (DSi) dynamics reflect integrated geologic, hydrologic, climatic, and ecological controls. We compiled annual DSi data for 337 rivers across four continents and trained interpretable machine-learning models to predict concentrations and yields from 28 watershed variables. Both models reproduced testing data ( R 2 = 0.85 for concentration and 0.96 for yield) and withheld-site validation ( R 2 = 0.91 and 0.93). Lithology, especially volcanic rock fraction, strongly controlled DSi while subsurface storage, topography, and land cover further shaped DSi dynamics. DSi concentrations and yields exhibited nonlinear responses to basin slope, recession-curve slope, proportion of open-water cover, and nutrient availability. Concentrations showed sharper threshold responses to hydrologic and biotic variables, whereas yields varied more gradually with climate and lithology. These results provide a framework for forecasting DSi under land cover and climate change and for embedding realistic, nonlinear processes in mechanistic models.

Geophysical Research Letters↗

USGS2021

Characteristics of frequent dynamic triggering of microearthquakes in Southern California

Dynamic triggering of earthquakes has been reported at various fault systems. The triggered earthquakes are thought to be caused either directly by dynamic stress changes due to the passing seismic waves, or indirectly by other nonlinear processes that are initiated by the passing waves. Distinguishing these physical mechanisms is difficult because of the general lack of high‐resolution earthquake catalogs and robust means to quantitatively evaluate triggering responses, particularly, delayed responses. Here we use the high‐resolution Quake Template Matching catalog in Southern California to systematically evaluate teleseismic dynamic triggering patterns in the San Jacinto Fault Zone and the Salton Sea Geothermal Field from 2008 to 2017. We develop a new statistical approach to identify triggered cases, finding that approximately 1 out of every 5 global M w ≥ 6 earthquakes dynamically trigger microearthquakes in Southern California. The triggering responses include both instantaneous and delayed triggering, showing a highly heterogeneous pattern and indicating possible evolving triggering thresholds. We do not observe a clear peak ground velocity triggering threshold that can differentiate triggering earthquakes from non‐triggering events, but there are subtle differences in the frequency content that may possibly differentiate the earthquakes. In contrast to the depth distribution of background seismicity, the identified triggered earthquakes tend to concentrate at the edges of the seismogenic zones. Although instantaneously triggered earthquakes are likely a result of dynamic Coulomb stress changes, the cases of delayed dynamic triggering are best explained by nonlinear triggering processes, including cyclic material fatigue, accelerated transient creep, and stochastic frictional heterogeneities.

California↗

USGS2014

Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model

Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

Hydrology and Earth System Sciences↗

USGS2009

Inversion of multichannel geophysical data with projected kernels

Statistical de‐noising methods such as Principal Component Analysis modify data in a way not constrained by physics. In much the same way as frequency‐filtered data must incorporate altered frequency content into numerical interpretation, so must statistically rotated data include the rotation operator in inversion processes. We propose a method of accounting for statistical reduction of data in non‐linear inversions (such as for electromagnetic data) through an incorporation of the rotation operator into the kernels and sensitivity matrix. We show a generalized linear inversion to demonstrate the necessity of rotating the inversion kernels in this abstract and will present a nonlinear parametric example from an unexploded ordnance application.

Conference Paper↗

USGS1996

Factors affecting suspended-solids concentrations in South San Francisco Bay, California

Measurements of suspended-solids concentration (SSC) were made at two depths at three sites in South San Francisco Bay (South Bay) to determine the factors that affect SSC. Twenty-eight segments of reliable and continuous SSC time series data longer than 14 days were collected from late 1991 or 1992 through September 1993. Spectral analysis and singular spectrum analysis were used to relate these data segments to time series of several potential forcing factors, including diurnal and semidiurnal tides, the spring-neap tidal cycle, wind shear, freshwater runoff, and longitudinal density differences. SSC is greatest during summer, when a landward wind shear is applied to South Bay by the afternoon sea breeze. About one half the variance of SSC is caused by the spring-neap cycle, and SSC lags the spring-neap cycle by about 2 days. Relatively short duration of slack water limits the duration of deposition of suspended solids and consolidation of newly deposited bed sediment during the tidal cycle, so suspended solids accumulate in the water column as a spring tide is approached and slowly deposit as a neap tide is approached. Perturbations in SSC caused by wind and local runoff from winter storms during the study period were usually much smaller than SSC variations caused by the spring-neap cycle. Variations of SSC at the study sites at tidal timescales are tidally forced, and nonlinear physical processes are significant. Advective transport dominates during spring tides, when water with higher SSC due to wind wave resuspension is advected to the main channel from shallow water, but, during neap tides, advective transport is less significant. The findings of this and other studies indicate that the tidally averaged transport of suspended solids responds to seasonal variations of wind shear in South Bay.

Journal of Geophysical Research C: Oceans↗

USGS2010

Generalizing a nonlinear geophysical flood theory to medium-sized river networks

The central hypothesis of a nonlinear geophysical flood theory postulates that, given space-time rainfall intensity for a rainfall-runoff event, solutions of coupled mass and momentum conservation differential equations governing runoff generation and transport in a self-similar river network produce spatial scaling, or a power law, relation between peak discharge and drainage area in the limit of large area. The excellent fit of a power law for the destructive flood event of June 2008 in the 32,400-km 2 Iowa River basin over four orders of magnitude variation in drainage areas supports the central hypothesis. The challenge of predicting observed scaling exponent and intercept from physical processes is explained. We show scaling in mean annual peak discharges, and briefly discuss that it is physically connected with scaling in multiple rainfall-runoff events. Scaling in peak discharges would hold in a non-stationary climate due to global warming but its slope and intercept would change.

Geophysical Research Letters↗

USGS1999

Multichannel analysis of surface waves

The frequency‐dependent properties of Rayleigh‐type surface waves can be utilized for imaging and characterizing the shallow subsurface. Most surface‐wave analysis relies on the accurate calculation of phase velocities for the horizontally traveling fundamental‐mode Rayleigh wave acquired by stepping out a pair of receivers at intervals based on calculated ground roll wavelengths. Interference by coherent source‐generated noise inhibits the reliability of shear‐wave velocities determined through inversion of the whole wave field. Among these nonplanar, nonfundamental‐mode Rayleigh waves (noise) are body waves, scattered and nonsource‐generated surface waves, and higher‐mode surface waves. The degree to which each of these types of noise contaminates the dispersion curve and, ultimately, the inverted shear‐wave velocity profile is dependent on frequency as well as distance from the source. Multichannel recording permits effective identification and isolation of noise according to distinctive trace‐to‐trace coherency in arrival time and amplitude. An added advantage is the speed and redundancy of the measurement process. Decomposition of a multichannel record into a time variable‐frequency format, similar to an uncorrelated Vibroseis record, permits analysis and display of each frequency component in a unique and continuous format. Coherent noise contamination can then be examined and its effects appraised in both frequency and offset space. Separation of frequency components permits real‐time maximization of the S/N ratio during acquisition and subsequent processing steps. Linear separation of each ground roll frequency component allows calculation of phase velocities by simply measuring the linear slope of each frequency component. Breaks in coherent surface‐wave arrivals, observable on the decomposed record, can be compensated for during acquisition and processing. Multichannel recording permits single‐measurement surveying of a broad depth range, high levels of redundancy with a single field configuration, and the ability to adjust the offset, effectively reducing random or nonlinear noise introduced during recording. A multichannel shot gather decomposed into a swept‐frequency record allows the fast generation of an accurate dispersion curve. The accuracy of dispersion curves determined using this method is proven through field comparisons of the inverted shear‐wave velocity ( υ s ">

Geophysics↗

USGS2006

Nonlinear inversion of potential-field data using a hybrid-encoding genetic algorithm

Using a genetic algorithm to solve an inverse problem of complex nonlinear geophysical equations is advantageous because it does not require computer gradients of models or "good" initial models. The multi-point search of a genetic algorithm makes it easier to find the globally optimal solution while avoiding falling into a local extremum. As is the case in other optimization approaches, the search efficiency for a genetic algorithm is vital in finding desired solutions successfully in a multi-dimensional model space. A binary-encoding genetic algorithm is hardly ever used to resolve an optimization problem such as a simple geophysical inversion with only three unknowns. The encoding mechanism, genetic operators, and population size of the genetic algorithm greatly affect search processes in the evolution. It is clear that improved operators and proper population size promote the convergence. Nevertheless, not all genetic operations perform perfectly while searching under either a uniform binary or a decimal encoding system. With the binary encoding mechanism, the crossover scheme may produce more new individuals than with the decimal encoding. On the other hand, the mutation scheme in a decimal encoding system will create new genes larger in scope than those in the binary encoding. This paper discusses approaches of exploiting the search potential of genetic operations in the two encoding systems and presents an approach with a hybrid-encoding mechanism, multi-point crossover, and dynamic population size for geophysical inversion. We present a method that is based on the routine in which the mutation operation is conducted in the decimal code and multi-point crossover operation in the binary code. The mix-encoding algorithm is called the hybrid-encoding genetic algorithm (HEGA). HEGA provides better genes with a higher probability by a mutation operator and improves genetic algorithms in resolving complicated geophysical inverse problems. Another significant result is that final solution is determined by the average model derived from multiple trials instead of one computation due to the randomness in a genetic algorithm procedure. These advantages were demonstrated by synthetic and real-world examples of inversion of potential-field data. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗

USGS2026

Constraining the permeability and outer-rise hydration at the Central America margin

Subduction zones are important drivers of Earth's geochemical evolution; the presence of water within the subduction system has substantial implications for geohazards, including the style of seismicity and the volume of magma generated at arc volcanoes. We utilize estimates of pore fluid distributions in the outer-rise region of the Central American margin imaged using controlled-source electro-magnetic surveys and novel two-dimensional forward models simulating poro-elastic deformation of slab bending to constrain the background permeability structure of the Cocos Plate and within outer-rise fault damage zones. The model results suggest that a relatively low background permeability structure is required to match the observed porosity at depths greater than 3 km (less than ~10 −21 m 2 . Similarly, increasing the permeability within the damage zones of outer-rise faults by more than a factor of 10 also over-predicts observed porosity distributions at depths greater than 3 km. However, the simulations consistently underpredict pore fluid distributions between depths of 1–3 km, which may reflect nonlinear elastic or anelastic deformation processes not included in the simulations or inherent uncertainty in the conversion of resistivity measurements to porosity distributions. The best-fitting depth-permeability relationship is relatively low, which has important implications for the timing of pore fluid release down-dip of the trench, as it restricts shallow fluid release, likely leading to higher pore fluid pressures, more widespread hydraulic fracturing, and deeper release of pore fluid along the plate interface.

Journal of Geophysical Research, Solid Earth↗

USGS2018

Degradation of 100‐m‐scale rocky ejecta craters at the InSight Landing Site on Mars and implications for surface processes and erosion rates in the hesperian and amazonian

Rocky ejecta craters (RECs) at the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) landing site on Elysium Planitia, Mars, provide constraints on crater modification and rates for the Hesperian and Amazonian. The RECs are between 10 m and 1.2 km in diameter and exhibit five classes of preservation. Class 1 represents pristine craters with sharp rims and abundant ejected rocks. From Classes 2 to 5, rims become more subdued, craters are infilled, and the ejecta become discontinuously distributed. High‐Resolution Imaging Science Experiment digital elevation models indicate a maximum depth to diameter ratio of ~0.15, which is lower than pristine models for craters of similar size. The low ratio is related to the presence of a loosely consolidated regolith and early‐stage eolian infill. Rim heights have an average height to diameter ratio of ~0.03 for the most pristine class. The size‐frequency distribution of RECs, plotted using cumulative and differential methods, indicates that crater classes within the diameter range of 200 m to 1.2 km are separated by ~100 to 200 Myr. Smaller craters degrade faster, with classes separated by <100 Myr. Rim erosion can be entirely modeled by nonlinear diffusional processes using the calculated timescales and a constant diffusivity of 8 × 10 −7 m 2 /year for craters 200 to 500 m in diameter. Diffusion models only partly capture depth‐related degradation, which requires eolian infill. Depth degradation and rim erosion rates are 10 −2 to 10 −3 m/Myr, respectively. The rates are consistent with relatively slow modification that is typical of the last two epochs of Martian history.

Journal of Geophysical Research - Planets↗

USGS2019

Multidecadal geomorphic evolution of a profoundly disturbed gravel-bed river system—a complex, nonlinear response and its impact on sediment delivery

A 2.5-km3 debris avalanche during the 1980 eruption of Mount St. Helens reset the fluvial landscape of upper North Fork Toutle River valley. Since then, a new drainage network has formed and evolved. Cross-section surveys repeated over nearly 40 years at 16 locations along a 20-km reach of river valley document channel evolution, geomorphic processes, and their impacts on sediment delivery. We analyzed spatial and temporal changes in channel morphology using two new metrics: 1) a shape index that defines the degree of U-shaped or V-shaped valley geometry; and 2) an alluvial phase-space diagram that relates bed degradation or aggradation between consecutive surveys to increases or decreases in cross-section area. Metric relations reveal more diverse channel evolution than originally described by a simple, linear-response model of sequential channel initiation and incision; aggradation and widening; and subsequent episodic scour and fill with little change in bed elevation. Instead, vertical and lateral adjustments have been crucial processes intertwined throughout channel evolution. Channel evolution has followed a distinctly nonlinear and non-sequential trajectory, migrating through several phase spaces and involving varied combinations of (1) degradation and aggradation with widening and narrowing, (2) bed-level fluctuations with little change in cross-section area, and (3) changes in cross-section area with little change of bed elevation. Persistent channel widening and reworking of the channel bed presently drive elevated sediment delivery from this basin. Elevated sediment delivery is likely to persist until valley-floor widths greatly exceed that of the channel-migration corridor, and/or channel banks and valley walls stabilize.

Washington↗

USGS2011

Preface: Multiscale feedbacks in ecogeomorphology

Geomorphic systems are known to exhibit nonlinear responses to physical–biological feedbacks (Thornes, 1985; Baas, 2002; Reinhardt et al., 2010). These responses make understanding and/or predicting system response to change highly challenging. With growing concerns over ecosystem health, a pressing need exists for research that tries to elucidate these feedbacks (Jerolmack, 2008; Darby, 2010; National Research Council, 2010). A session was convened at the Fall 2008 meeting of the American Geophysical Union (AGU) to provide an outlet for some of this truly interdisciplinary and original research, which is central to understanding geomorphic and ecological dynamics. The session attracted over 39 contributions, which were divided into two well-attended oral sessions and a very busy poster session. This special issue presents new research from the AGU session, which highlights clear physical–biological feedbacks. The aim is to bring together contrasting perspectives on biological and geomorphic feedbacks in a diversity of physiographic settings, ranging from wetlands and estuaries, through rivers, to uplands. These papers highlight biological and physical feedbacks which involve the modulation or amplification of geomorphic processes. These papers will be of interest to a core geomorphology audience, and should also draw attention from the fields of ecohydraulics, hydroecology, ecohydrology, ecomorphology, biogeochemistry and biogeography, and biogeomorphology as well as the more traditional fields of hydrology, ecology and biology. In this preface to the special issue, we a) review past contributions to the emerging field of ecogeomorphology and related disciplines, b) provide some context for how this topical special issue came to fruition, and c) summarize the contributions to this special issue.

Geomorphology↗

USGS2007

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts↗

USGS2010

Numerical simulation of magmatic hydrothermal systems

The dynamic behavior of magmatic hydrothermal systems entails coupled and nonlinear multiphase flow, heat and solute transport, and deformation in highly heterogeneous media. Thus, quantitative analysis of these systems depends mainly on numerical solution of coupled partial differential equations and complementary equations of state (EOS). The past 2 decades have seen steady growth of computational power and the development of numerical models that have eliminated or minimized the need for various simplifying assumptions. Considerable heuristic insight has been gained from process-oriented numerical modeling. Recent modeling efforts employing relatively complete EOS and accurate transport calculations have revealed dynamic behavior that was damped by linearized, less accurate models, including fluid property control of hydrothermal plume temperatures and three-dimensional geometries. Other recent modeling results have further elucidated the controlling role of permeability structure and revealed the potential for significant hydrothermally driven deformation. Key areas for future reSearch include incorporation of accurate EOS for the complete H2O-NaCl-CO2 system, more realistic treatment of material heterogeneity in space and time, realistic description of large-scale relative permeability behavior, and intercode benchmarking comparisons. Copyright 2010 by the American Geophysical Union.

Reviews of Geophysics↗

USGS2022

Fast rupture of the 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California dynamically triggers seismicity in California

In the Gulf of California, Mexico, the relative motion across the North America-Pacific boundary is accommodated by a series of marine transform faults and spreading centers. About 40 M>6 earthquakes have occurred in the region since 1960. On 3 August 2009, an Mw 6.9 earthquake occurred near Canal de Ballenas in the region. The earthquake was a strike-slip event with a shallow hypocenter that is likely close to the seafloor. In contrast to an adjacent M7 earthquake, this earthquake triggered a ground-motion-based earthquake early warning algorithm being tested in southern California (∼600 km away). This observation suggests that the abnormally large ground motions and dynamic strains observed for this earthquake relate to its rupture properties. To investigate this possibility, we image the rupture process and resolve the slip distribution of the event using a P-wave back-projection approach and a teleseismic, finite-fault inversion method. Results from these two independent analyses indicate a relatively simple, unilateral rupture propagation directed along-strike in the northward direction. However, the average rupture speed is estimated around 4 km/s, suggesting a possible supershear rupture. The supershear speed is also supported by a Rayleigh wave Mach cone analysis, although uncertainties in local velocity structure preclude a definitive conclusion. The Canal de Ballenas earthquake dynamically triggered seismicity at multiple sites in California, with triggering response characteristics varying from location-to-location. For instance, some of the triggered earthquakes in California occurred up to 24 hours later, suggesting that nonlinear triggering mechanisms likely have modulated their occurrence.

Baja California↗

USGS2017

Interactions of estuarine shoreline infrastructure with multiscale sea level variability

Sea level rise increases the risk of storms and other short‐term water‐rise events, because it sets a higher water level such that coastal surges become more likely to overtop protections and cause floods. To protect coastal communities, it is necessary to understand the interaction among multiday and tidal sea level variabilities, coastal infrastructure, and sea level rise. We performed a series of numerical simulations for San Francisco Bay to examine two shoreline scenarios and a series of short‐term and long‐term sea level variations. The two shoreline configurations include the existing topography and a coherent full‐bay containment that follows the existing land boundary with an impermeable wall. The sea level variability consists of a half‐meter perturbation, with duration ranging from 2 days to permanent (i.e., sea level rise). The extent of coastal flooding was found to increase with the duration of the high‐water‐level event. The nonlinear interaction between these intermediate scale events and astronomical tidal forcing only contributes ∼1% of the tidal heights; at the same time, the tides are found to be a dominant factor in establishing the evolution and diffusion of multiday high water events. Establishing containment at existing shorelines can change the tidal height spectrum up to 5%, and the impact of this shoreline structure appears stronger in the low‐frequency range. To interpret the spatial and temporal variability at a wide range of frequencies, Optimal Dynamic Mode Decomposition is introduced to analyze the coastal processes and an inverse method is applied to determine the coefficients of a 1‐D diffusion wave model that quantify the impact of bottom roughness, tidal basin geometry, and shoreline configuration on the high water events

Journal of Geophysical Research: Oceans↗