Geology ReportsSearch

SEARCH · Geology Reports

Results for “Frontiers in Ecology and Evolution”

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

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

1,673 records · Page 83Linked to original sources

Evaluating mass flow meter measurements from chambers for greenhouse gas emissions from orphan wells and other point sources

This study evaluates the performance of a rigid gas flux chamber equipped with a mass flow meter (MFM) for measuring gas emissions from leaking orphan wells and similar pressure-driven gas point sources. We conducted a series of laboratory and field experiments to evaluate the sensitivity, stability, and dynamic range of an MFM chamber system and found an optimal method for sealing the chamber to the ground to isolate the emission source. From these results, we estimate the effects of different soil gas permeabilities on measurements and identify the uncertainty of environmental processes that can impact measurements. Simulations of an MFM chamber are compared to those of a dynamic flux chamber to contrast the data derived with both methodologies and illustrate the potential for measuring high variability leaks with the MFM chamber. Using a low flow resistance MFM and a chamber well-sealed to the ground, it is possible to measure leaks down to 1.08 x 10 -3 cubic meters per hour (m 3 h −1 ) (refenced to 25°/1 atm), corresponding to 0.77 grams per hour (g h −1 ) methane or 2.11 g h −1 carbon dioxide, with a mean uncertainty of 0.89 % relative standard deviation. Environmental processes such as heated gas inside the chamber from solar gain, wind blowing across the chamber vent, and changing humidity in the chamber, can cause variation in MFM measurements. Over 11 d of continuous monitoring under varying weather conditions, the standard deviation of the environmentally sourced signals was found to be 7.40 x 10 -3 m 3 h −1 (equivalent to or 5.27 g h −1 methane or 14.45 g h −1 carbon dioxide). Strategies to obtain the highest quality data from MFM chambers include burying the edges of the chamber below the surface sufficiently deep to seal the chamber edges against gas flow and soaking the dirt with water to lower the chances of escaping gases, while monitoring the gas flow and adjusting the chamber seal to achieve a maximum flow rate.

Atmospheric Measurement Techniques

Using a temporary emigration model to estimate abundance of stream fishes from hybrid removal surveys with and without block nets

Monitoring programs are often faced with a decision to allocate resources into either robust spatiotemporal coverage to estimate a population index (e.g., not true abundance) or confirming closed sampling conditions (e.g., with block nets) for an unbiased population estimate at the cost of spatiotemporal coverage. However, making accurate and precise abundance estimates at robust spatiotemporal scales is possible when combining open and closed sampling designs with integrated modeling techniques. We used simulations and a case study of backpack electrofishing surveys in the Santa Ana River, California to test the efficacy of an integrated abundance model (temporary emigration model, TE) to estimate abundance of fishes using removal sampling methods with a hybrid sampling design (sampling with and without block nets during removal sampling). We found that the TE model performed well under most modeling scenarios (sample size, amount of closure violation, number of samples collected during closure), although at least a few samples with block nets were necessary for all parameters to be estimable. When applied to fish surveys in the Santa Ana River, we found that catch of the fishes fit to the TE model (Santa Ana Sucker, Arroyo Chub, Channel Catfish, Largemouth Bass, Yellow Bullhead) showed little evidence that the closure assumption was violated when block nets were not used. Additionally, we found that the abundance of non-native fishes negatively affected the abundance of the native Santa Ana Sucker, which was also found to adversely affect the native fish’s access to critical habitat consisting of gravel and cobble substrate. Our results indicate that the TE model presents a viable solution to common sampling problems that impact many monitoring programs, where precise and accurate population estimates can be made at large spatiotemporal scales even when most samples violate the closure assumption.

California

The U.S. Geological Survey National Water Quality Network—Groundwater—2023

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, eight of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,114 wells in the network.

conterminous United States

2024 Crustal Deformation Modeling Workshop report

The 2024 Crustal Deformation Modeling Workshop was held June 10–14 at the Colorado School of Mines. The workshop included two days of tutorials on using PyLith for crustal deformation modeling, followed by three days of science talks and discussions. The workshop focused on four primary themes: ● Constraining long-term fault slip rates and their uncertainties using geodetic and geologic data; ● Earthquake cycle modeling with a focus on constraining models using seismic and geodetic data; ● Interaction of fluids and faulting; and ● Separating contributions of surface loading and tectonic loading in crustal deformation. The complete agenda is available on the CIG website.

Conference Paper

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California

Map of topographic lineaments interpreted as recent surface ruptures along the Bennett Valley and Southern Maacama Fault Zones, Sonoma County, California

This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.

California

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications

Cross-fade sampling: Extremely efficient Bayesian inversion for a variety of geophysical problems

This paper introduces cross-fade sampling, a computationally efficient Markov Chain Monte Carlo simulation method that uses a semi-analytical approach to quickly solve Bayesian inverse problems that do not themselves have an analytical solution. Cross-fading is efficient in two ways. First, it requires fewer samples to obtain the same quality simulation of the target probability density function (PDF). Secondly, it is much faster to evaluate the posterior probability of each sample than conventional sampling methods for simulating Bayesian posterior PDFs. Conventional methods require evaluating the prior probability (which describes your a priori constraints) and data likelihood (which describes the fit between the observations and the predictions of the model) for each sample model. However, cross-fading does not require evaluating the data likelihood, meaning that ‘big data’ can be fit with zero additional computational cost. Further, the cross-fading approach can be used to calculate the marginal likelihood associated with a model design, facilitating model comparison and Bayesian model averaging. Topics covered in this paper include derivation of the cross-fade approach and how it can be used to simulate Bayesian posterior PDFs and compute the marginal likelihood, discussion of the class of problems to which cross-fading can be applied (with examples from earthquake statistics, earthquake ground motion modelling, volcanic eruption forecasting, and finite fault slip modelling), demonstration of efficiency relative to existing sampling methods and discussion of how cross-fading can be used to account for prediction errors (i.e. epistemic errors) as part of the geophysical inverse problem.

Geophysical Journal International

Using geologic mapping to understand temporal and spatial relations of closely clustered to concurrent latest Holocene surface ruptures on two intersecting faults, south-central Mojave Desert, California

The Pinto Mountain Fault Zone (PMFZ) marks a major structural boundary between east-oriented sinistral faults of the eastern Transverse Ranges (to the south) and northwest-oriented dextral faults of the south-central Mojave Desert (to the north). These structural fault systems comprise sinistral and dextral deformational domains of the Eastern California Shear Zone (ECSZ) that intersect one another in the Copper Mountain and Twentynine Palms areas. The U.S. Geological Survey (USGS) is conducting detailed geologic mapping and geochronologic investigations designed to clarify geometric, kinematic, and temporal relations among the two domains, that are focused on the central portion of the left-lateral PMFZ near its intersection with the right-lateral Copper Mountain Fault (CMF) and Mesquite Lake Fault Zone (MLFZ).

California

Digitizer Suite: The Albuquerque Seismological Laboratory Digitizer Testing Suite

Laboratory testing of digitizers and seismometers helps ensure that prior to deployment the instrumentation can produce high quality data and is operating within specifications. In this work we detail the software package called: the Albuquerque Seismological Laboratory (ASL) Digitizer Test Suite. This Java software package provides several algorithms to verify various performance parameters of digitizers commonly used for recording analog seismic instruments. The goal of these tests is not to be exhaustive, but to identify common failures that could compromise the integrity of seismic data being recorded on the digitizer. For example, Sandia National Laboratories (e.g., Slad and Merchant, 2018) routinely do comprehensive testing of digitizers for various monitoring missions. While these tests reports are valuable for comprehensively characterizing a recording system, it would be resource intensive to conduct such tests on every seismic recorder used in a network. We focus on tests that include ways to estimate the sensitivity, timing, self-noise, and clip-level of the digitizer, as well as the fidelity of the signal being recorded. The software is publicly available and provides a way for the community to verify the integrity of a digitizer using a minimum amount of outside equipment.

Seismological Research Letters

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone

Field observations and logs from the Rose Hip trench exposure across a north-facing scarp within the Seattle Fault Zone, southern Bainbridge Island, Washington

The Seattle Fault Zone is an approximately 70-km-long, east-west-trending zone of south-dipping blind reverse faults within the Puget lowland region in Washington. Because of the proximity, the Seattle Fault Zone poses a significant earthquake hazard to the Puget sound and Seattle metropolitan regions. We present preliminary mapping and trench-site information from a paleoseismic investigation across a newly identified active fault scarp located within the hanging wall of the Seattle Fault Zone on southern Bainbridge Island, Washington. The trench exposed monoclinally folded Miocene bedrock, fractured and faulted glacial-related deposits, and laminated lacustrine deposits capped by slope-derived colluvium. The observations from this investigation record late Pleistocene to Holocene north-vergent folding and faulting along this new fault scarp.

Washington

Conventional hydrothermal power-producing systems of the Great Basin, USA

As part of the update to the electric-grade conventional hydrothermal assessment of the Great Basin, USA, Monte Carlo analyses of identified resources within explored regions will be performed to make estimates of discovered resources and associated uncertainty. Analyses use conditional statistics where estimates are conditioned upon a hydrothermal favorability map, allowing for the likelihood that more resources exist in regions of higher hydrothermal favorability. For these analyses, a dataset of identified hydrothermal systems is compiled, and the new compilation is described herein. Recognizing that a single hydrothermal system may be developed with multiple power plants, and that the hydrothermal upflow zone may be several kilometers across with many measurements characterizing a single hydrothermal system, a procedure was developed and employed to create clusters of points (power plants, measurements, etc.) that are associated with a single system, and a new central point was defined as the best estimator of the center of the hydrothermal system. Hydrothermal systems were uniquely identified by grouping electric-grade hydrothermal measurements and operating power plants within a distance of 10 km. Groups that are >10 km apart are assumed to be different electric-grade hydrothermal systems. While 10 km was used as the threshold, most systems were significantly further apart, and most points within groups were typically within 5 km of each other. A well measurement was considered an electric-grade measurement of a hydrothermal system if it had two properties: a measured temperature of >85 °C and evidence of hydrothermal convection. Other points that were added to the dataset are locations of operating powerplants or locations that have been classified as an electric-grade hydrothermal resource by either the U.S Geological Survey (USGS) or the Great Basin Center for Geothermal Energy. After all points are assigned to systems, new points were computed with the goal of identifying the center of the throat of the hydrothermal upflow zone. If operating powerplants exist for a system, then the arithmetic average of all power plant locations is used. Otherwise, if USGS made an estimate, that location is used. In the absence of both powerplants or USGS estimates, the arithmetic average of all electric-grade measurement locations is used. An example is shown of how these newly compiled locations might be ranked for uncertainty analyses, where higher confidence is assumed if measured temperature is higher and there are many supporting measurements indicating an electric-grade resource. In summary, 28 systems have operating power plants, an additional 78 systems are known identified electric-grade hydrothermal resources, and 99 new systems were identified as probable electric-grade systems with varying levels of confidence. These 205 locations are shown as a function of a recent hydrothermal favorability map, conceptually illustrating the conditional statistics that can be used to make estimates of the undiscovered resources of the Great Basin. An accompanying data release provides summaries of developed capacity by system and USGS estimates of likely total capacity and associated uncertainty.

Arizona, California, Idaho, Nevada, Oregon, Utah

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş