Search USGSSearch

SEARCH · Search USGS

Results for “Mathematical and Computer Modelling”

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.

At least 55 records · Page 3Linked to original sources

Estimating the magnitude and frequency of low flows of streams in Massachusetts

The report presents techniques used for estimating 7-day 2-year (7Q2) and 7-day 10-year (7Q10) flows at continuous-record streamflow-gaging stations, partial-record stations, and ungaged sites on streams in Massachusetts. A two-parameter log-normal probability distribution was used to compute low flows at 31 continuous-record stream-flow-gaging stations. Graphical and mathematical techniques were used to estimate low flows at partial-record stations. Regional regression techniques commonly are used to estimate 7Q2 and 7Q10 flows at ungaged stream sites. The development and application of a low-flow regression model is described. The model contains two parameters, mean and standard deviation, computed from separate regression equations that use total drainage area and basin relief as independent variables. The coefficients of determination of the mean and standard deviation regression equations are 0.964 and 0.960, respec- tively; the percentage of standard error of the regression of each equation are 35 and 34, respec- tively. The model is applicable for basins in Massachusetts having a drainage area of 5 to 150 square miles, a slope of less than 4 percent, and an area of stratified drift greater than 4 percent of the total drainage area. A computer program is presented that can be used to compute 7Q2 and 7Q10 flows and the 95-percent intervals of confidence and prediction for a given basin.

Water-Resources Investigations Report

Flood-frequency analysis for small natural streams in Georgia

Flood information from 104 project basins (0.1 to 20 square miles) and data available from 170 larger basins (20 to 1,500 square miles) are analyzed to provide planners and designers with relations for estimating the magnitude and frequency of flood-peak discharges on small natural streams in Georgia. More than 90 percent of the project activity was directed toward the acquisition and processing of data for more than two thousand rainfall-runoff periods occurring during 1964-74 at the 104 project stations. Storm rainfall data were computed at 5-minute intervals from records of six long-term (about 75 years) National Weather Service Stations. The U.S. Geological Survey mathematical rainfall-runoff model is used to synthesize long-term flood records for about 80 percent of the project stations. Station flood-frequency data are developed, using the log-Pearson Type III distribution with regionalized skew coefficients. Multiple regression analysis is used to define relations between flood-frequency station data for small and large streams and 10 physical and climatological basin characteristics. The analysis indicates the drainage-basin size is the most significant variable. Five regions having distinct flood-peak characteristics are delineated. The developed relations, expressed as equations and nomographs, are considered usable for virtually any site in Georgia where the drainage area is between 0.1 and 20 square miles, and the flow is natural. (Woodard-USGS)

Open-File Report

Preliminary implications of viscoelastic ray theory for anelastic seismic tomography models

The recent developments in general viscoelastic ray theory provide a rigorous mathematical framework for anelastic seismic tomography. They provide closed‐form solutions of forward ray‐tracing and simple inverse problems for anelastic horizontal and spherical layered media with material gradients. They provide ray‐tracing computation algorithms valid for all angles of incidence that account for changes in wave speed, attenuation, and trajectory of anelastic P and S body waves induced by anelastic boundaries. They account for theoretical predictions that seismic waves refract as inhomogeneous waves across anelastic boundaries for all angles of incidence, which in turn accounts for energy carried by plane waves along seismic boundaries at head wave critical angles and wide‐angle refracted (WAR) ray paths that are not predicted by elastic models. Exact viscoelastic ray‐tracing numerical results for various models provide examples that illustrate the effects of anelastic boundaries on the travel times and amplitudes of seismic waves. They show the effects are strongly dependent on angle of incidence. For near‐critical and wide angles of incidence the anelastic effects on travel times and amplitudes can be large and are not explained by elastic ray theory, but the effects on travel times can be relatively small and difficult to distinguish from those for elastic media for pre‐near‐critical angles of incidence. The results for some models indicate that reflected anelastic WAR waves may be observable at the surface and possibly account for some prominent seismic arrivals not explained by elasticity. These preliminary results suggest that the application of exact viscoelastic ray‐tracing computation algorithms to exploration and teleseismic data sets can reveal new insights regarding the properties and distribution of anelastic materials in the Earth.

Bulletin of the Seismological Society of America

A computer program to calculate the resistivity and induced polarization response for a three-dimensional body in the presence of buried electrodes

Three-dimensional induced polarization and resistivity modeling for buried electrode configurations can be achieved by adapting surface integral techniques for surface electrode configurations to buried electrodes. Modification of. the surface technique is accomplished by considering the additional mathematical terms required to express-the changes in the electrical potential and geometry caused by placing the source and receiver electrodes below the surface. This report presents a listing of a computer program to calculate the resistivity and induced polarization response from a three-dimensional body for buried electrode configurations. The program is designed to calculate the response for the following electrode configurations: (1) hole-to-surface array with a buried bipole source and a surface bipole receiver, (2) hole-to-surface array with a buried pole source and a surface bipole receiver, (3) hole-to-hole array with a buried, fixed pole source and a moving bipole receiver, (4) surface-to-hole array with a fixed pole source on the surface and a moving bipole receiver in the borehole, (5) hole-to-hole array with a buried, fixed bipole source and a buried, moving bipole receiver, (6) hole-to-hole array with a buried, moving bipole source and a buried, moving bipole receiver, and (7) single-hole, buried bipole-bipole array. Input and output examples are given for each of the arrays.

Open-File Report

Uncertainty quantification for environmental models

Environmental models are used to evaluate the fate of fertilizers in agricultural settings (including soil denitrification), the degradation of hydrocarbons at spill sites, and water supply for people and ecosystems in small to large basins and cities—to mention but a few applications of these models. They also play a role in understanding and diagnosing potential environmental impacts of global climate change. The models are typically mildly to extremely nonlinear. The persistent demand for enhanced dynamics and resolution to improve model realism [17] means that lengthy individual model execution times will remain common, notwithstanding continued enhancements in computer power. In addition, high-dimensional parameter spaces are often defined, which increases the number of model runs required to quantify uncertainty [2]. Some environmental modeling projects have access to extensive funding and computational resources; many do not. The many recent studies of uncertainty quantification in environmental model predictions have focused on uncertainties related to data error and sparsity of data, expert judgment expressed mathematically through prior information, poorly known parameter values, and model structure (see, for example, [1,7,9,10,13,18]). Approaches for quantifying uncertainty include frequentist (potentially with prior information [7,9]), Bayesian [13,18,19], and likelihood-based. A few of the numerous methods, including some sensitivity and inverse methods with consequences for understanding and quantifying uncertainty, are as follows: Bayesian hierarchical modeling and Bayesian model averaging; single-objective optimization with error-based weighting [7] and multi-objective optimization [3]; methods based on local derivatives [2,7,10]; screening methods like OAT (one at a time) and the method of Morris [14]; FAST (Fourier amplitude sensitivity testing) [14]; the Sobol' method [14]; randomized maximum likelihood [10]; Markov chain Monte Carlo (MCMC) [10]. There are also bootstrapping and cross-validation approaches.Sometimes analyses are conducted using surrogate models [12]. The availability of so many options can be confusing. Categorizing methods based on fundamental questions assists in communicating the essential results of uncertainty analyses to stakeholders. Such questions can focus on model adequacy (e.g., How well does the model reproduce observed system characteristics and dynamics?) and sensitivity analysis (e.g., What parameters can be estimated with available data? What observations are important to parameters and predictions? What parameters are important to predictions?), as well as on the uncertainty quantification (e.g., How accurate and precise are the predictions?). The methods can also be classified by the number of model runs required: few (10s to 1000s) or many (10,000s to 1,000,000s). Of the methods listed above, the most computationally frugal are generally those based on local derivatives; MCMC methods tend to be among the most computationally demanding. Surrogate models (emulators)do not necessarily produce computational frugality because many runs of the full model are generally needed to create a meaningful surrogate model. With this categorization, we can, in general, address all the fundamental questions mentioned above using either computationally frugal or demanding methods. Model development and analysis can thus be conducted consistently using either computation-ally frugal or demanding methods; alternatively, different fundamental questions can be addressed using methods that require different levels of effort. Based on this perspective, we pose the question: Can computationally frugal methods be useful companions to computationally demanding meth-ods? The reliability of computationally frugal methods generally depends on the model being reasonably linear, which usually means smooth nonlin-earities and the assumption of Gaussian errors; both tend to be more valid with more linear

SIAM News

Flood frequency of Mississippi streams

Techniques have been developed for estimating future flood magnitudes having recurrence intervals as great as 100 years. Estimates for gaged sites were obtained by application of the log-Pearson Type III frequency distribution. Estimates for ungaged sites are defined by mathematical relations which may be solved using data from topographic maps. Nomographs were designed for making rapid estimates for floodfrequency values. Snythetic records generated by a digital computer Model of the rainfall-runoff process for drainage areas ranging from 0.04 to 4.35 square miles (0.1 to 11.3 square kilometres) were combined with long-term observed records to define the relation of basic characteristics to flood frequency on drainage areas as great as 6,630 square miles (17,200 square kilometres). Multiple-regression techniques indicated that drainage area, slope, and length were the most significant variables affecting flood frequency.

Mississippi

Approaches to highly parameterized inversion: Pilot-point theory, guidelines, and research directions

Pilot points have been used in geophysics and hydrogeology for at least 30 years as a means to bridge the gap between estimating a parameter value in every cell of a model and subdividing models into a small number of homogeneous zones. Pilot points serve as surrogate parameters at which values are estimated in the inverse-modeling process, and their values are interpolated onto the modeling domain in such a way that heterogeneity can be represented at a much lower computational cost than trying to estimate parameters in every cell of a model. Although the use of pilot points is increasingly common, there are few works documenting the mathematical implications of their use and even fewer sources of guidelines for their implementation in hydrogeologic modeling studies. This report describes the mathematics of pilot-point use, provides guidelines for their use in the parameter-estimation software suite (PEST), and outlines several research directions. Two key attributes for pilot-point definitions are highlighted. First, the difference between the information contained in the every-cell parameter field and the surrogate parameter field created using pilot points should be in the realm of parameters which are not informed by the observed data (the null space). Second, the interpolation scheme for projecting pilot-point values onto model cells ideally should be orthogonal. These attributes are informed by the mathematics and have important ramifications for both the guidelines and suggestions for future research.

Scientific Investigations Report

Subsurface storage of freshwater in South Florida; a digital model analysis of recoverability

As part of a study of the feasibility of recovering freshwater injected and stored underground in south Florida, a digital solute-transport model was used to investigate the relation of recovery efficiency to the variety of hydrogeologic conditions that could prevail in brackish artesian aquifers and to a variety of management alternatives. The analyses employed a modeling approach in which the control for sensitivity testing was a hypothetical aquifer considered representative of permeable zones in south Florida that might be used for storage of freshwater. Parameter variations in the tests represented possible variations in aquifer conditions in the area. The applicability of the analyses to south Florida limestone aquifers required the assumption that flow nonuniformities in those aquifers are small on the scale of volumes of water likely to be injected, and that their effect could be represented as hydrodynamic dispersion. Generally, it was shown that a loss of recovery efficiency is caused by (1) processes causing mixing of injected freshwater with native saline water (hydrodynamic dispersion), (2) processes causing the more or less irreversible displacement of the injected freshwater with respect to the well (buoyancy stratification, background hydraulic gradients, and interlayer dispersion), or (3) processes causing injection and withdrawal flow patterns to be dissimilar (directionally biased well-bore plugging, and dissimilar injection and withdrawal schedules in multiple-well systems). Other results indicated that recovery efficiency improves considerably with successive cycles, providing that each recovery phase ends when the chloride concentration of withdrawn water exceeds established criteria for potability (usually 250 milligrams per liter), and that freshwater injected into highly permeable or highly saline aquifers (such as the 'boulder zone') would buoy rapidly. Many hydrologic conditions were posed for model analysis. To have obtained comparable results with operational testing would have been more costly by orders of magnitude. The tradeoff is that the validity of results obtained from computer modeling is somewhat less certain. In particular, results must be qualified with observations that (1) the complex set of processes lumped as hydrodynamic dispersion is represented with a somewhat simplified mathematical approximation, and (2) other flow processes in limestone injection zones are as yet incompletely understood. Despite such reservations, the study is considered a practical example of the use of transport models in ground-water investigations.

Water Supply Paper

Carrying capacity of a population diffusing in a heterogeneous environment

The carrying capacity of the environment for a population is one of the key concepts in ecology and it is incorporated in the growth term of reaction-diffusion equations describing populations in space. Analysis of reaction-diffusion models of populations in heterogeneous space have shown that, when the maximum growth rate and carrying capacity in a logistic growth function vary in space, conditions exist for which the total population size at equilibrium (i) exceeds the total population that which would occur in the absence of diffusion and (ii) exceeds that which would occur if the system were homogeneous and the total carrying capacity, computed as the integral over the local carrying capacities, was the same in the heterogeneous and homogeneous cases. We review here work over the past few years that has explained these apparently counter-intuitive results in terms of the way input of energy or another limiting resource (e.g., a nutrient) varies across the system. We report on both mathematical analysis and laboratory experiments confirming that total population size in a heterogeneous system with diffusion can exceed that in the system without diffusion. We further report, however, that when the resource of the population in question is explicitly modeled as a coupled variable, as in a reaction-diffusion chemostat model rather than a model with logistic growth, the total population in the heterogeneous system with diffusion cannot exceed the total population size in the corresponding homogeneous system in which the total carrying capacities are the same.

Mathematics and Computers in Modern Science - Acou

Computational methodology to analyze the effect of mass transfer rate on attenuation of leaked carbon dioxide in shallow aquifers

Exsolution and re-dissolution of CO 2 gas within heterogeneous porous media are investigated using experimental data and mathematical modeling. In a set of bench-scale experiments, water saturated with CO 2 under a given pressure is injected into a 2-D water-saturated porous media system, causing CO 2 gas to exsolve and migrate upwards. A layer of fine sand mimicking a heterogeneity within a shallow aquifer is present in the tank to study accumulation and trapping of exsolved CO 2 . Then, clean water is injected into the system and the accumulated CO 2 dissolves back into the flowing water. Simulated exsolution and dissolution mass transfer processes are studied using both nearequilibrium and kinetic approaches and compared to experimental data under conditions that do and do not include lateral background water flow. The mathematical model is based on the mixed hybrid finite element method that allows for accurate simulation of both advection- and diffusion- dominated processes.

Acta Polytechnica

Fitting a structured juvenile-adult model for green tree frogs to population estimates from capture-mark-recapture field data

We derive point and interval estimates for an urban population of green tree frogs ( Hyla cinerea ) from capture–mark–recapture field data obtained during the years 2006–2009. We present an infinite-dimensional least-squares approach which compares a mathematical population model to the statistical population estimates obtained from the field data. The model is composed of nonlinear first-order hyperbolic equations describing the dynamics of the amphibian population where individuals are divided into juveniles (tadpoles) and adults (frogs). To solve the least-squares problem, an explicit finite difference approximation is developed. Convergence results for the computed parameters are presented. Parameter estimates for the vital rates of juveniles and adults are obtained, and standard deviations for these estimates are computed. Numerical results for the model sensitivity with respect to these parameters are given. Finally, the above-mentioned parameter estimates are used to illustrate the long-time behavior of the population under investigation.

Bulletin of Mathematical Biology

Documentation of a Regional Aquifer Simulation Model, RAQSIM, and a description of support programs applied in the Twin Platte-Middle Republican Study Area, Nebraska

RAQSIM, a generalized flow model of a groundwater system using finite-element methods, is documented to explain how it works and to demonstrate that it gives valid results. Three support programs that are used to compute recharge and discharge data required as input to RAQSIM are described. RAQSIM was developed to solve transient, two-dimensional, regional groundwater flow problems with isotropic or anisotropic conductance. The model can also simulate radially-symmetric flow to a well and steady-state flow. The mathematical basis, program structure, data input and output procedures, organization of data sets, and program features and options of RAQSIM are discussed. An example , containing listings of data and results and illustrating RAQSIM 's capabilities, is discussed in detail. Two test problems also are discussed comparing RAQSIM 's results with analytical procedures. The first support program described, the PET Program, uses solar radiation and other climatic data in the Jensen-Haise method to compute potential evapotranspiration. The second support program, the Soil-Water Program, uses output from the PET Program, soil characteristics, and the ratio of potential to actual evapotranspiration for each crop to compute infiltration, storage, and removal of water from the soil zone. The third program, the Recharge-Discharge Program, uses output from the Soil-Water Program together with other data to compute recharge and discharge from the groundwater flow system. For each support program, a program listing and examples of the data and results for the Twin Platte-Middle Republican study are provided. In addition, a brief discussion on how each program operates and on procedures for running and modifying these programs are presented. (Author 's abstract)

Nebraska

Mathematical models of catchment behavior

After an examination of trends in the modeling of hydrologic systems, a review of some recent studies is given. The authors' preliminary studies on the feasibility and efficiency of the automatic evaluation of catchment model parameters by use of a digital computer are described and some results presented.

Journal of the Hydraulics Division

HyWaves: Hybrid downscaling of multimodal wave spectra to nearshore areas

Long-term and accurate wave hindcast databases are often required in different coastal engineering projects. The assessment of the nearshore wave climate is often accomplished by using downscaling techniques to translate offshore waves to coastal areas. However, dynamical downscaling approaches may incur huge computational cost. Additionally, the common use of bulk parameterizations are often not accurate for multidimensional waves. To overcome these limitations, we present a hybrid downscaling approach that combines mathematical algorithms (statistical downscaling) and numerical modeling (dynamical downscaling) over the individual spectral partitions. Every wave partition is downscaled and aggregated afterward by using principles of wave linear theory. By assuming linearity in the propagation of the wave celerity, the application of the method is limited from offshore to intermediate water depths. In addition, the method proposed uses a technique to simplify the spectral boundary conditions in complex domains. The methodology has been applied and validated in the island states of Samoa, American Samoa, Majuro, and Kwajalein, showing good skill at reproducing the spectral hourly time series of significant wave height, peak period, and peak direction. Moreover, an accurate representation of the observed energy spectrum was achieved. This study provides insight into the numerical approximation of the combined sea-swell states while improving the quality of fast spectral forecasting and early warning systems.

Ocean Modeling

Assessing local population vulnerability to wind energy development with branching process models: an application to wind energy development

Quantifying the impact of anthropogenic development on local populations is important for conservation biology and wildlife management. However, these local populations are often subject to demographic stochasticity because of their small population size. Traditional modeling efforts such as population projection matrices do not consider this source of variation whereas individual-based models, which include demographic stochasticity, are computationally intense and lack analytical tractability. One compromise between approaches is branching process models because they accommodate demographic stochasticity and are easily calculated. These models are known within some sub-fields of probability and mathematical ecology but are not often applied in conservation biology and applied ecology. We applied branching process models to quantitatively compare and prioritize species locally vulnerable to the development of wind energy facilities. Specifically, we examined species vulnerability using branching process models for four representative species: A cave bat (a long-lived, low fecundity species), a tree bat (short-lived, moderate fecundity species), a grassland songbird (a short-lived, high fecundity species), and an eagle (a long-lived, slow maturation species). Wind turbine-induced mortality has been observed for all of these species types, raising conservation concerns. We simulated different mortality rates from wind farms while calculating local extinction probabilities. The longer-lived species types (e.g., cave bats and eagles) had much more pronounced transitions from low extinction risk to high extinction risk than short-lived species types (e.g., tree bats and grassland songbirds). High-offspring-producing species types had a much greater variability in baseline risk of extinction than the lower-offspring-producing species types. Long-lived species types may appear stable until a critical level of incidental mortality occurs. After this threshold, the risk of extirpation for a local population may rapidly increase with only minimal increases in wind mortality. Conservation biologists and wildlife managers may need to consider this mortality pattern when issuing take permits and developing monitoring protocols for wind facilities. We also describe how our branching process models may be generalized across a wider range of species for a larger assessment project and then describe how our methods may be applied to other stressors in addition to wind.

Ecosphere

Stream temperature investigations: field and analytic methods

This document provides guidance to the user of the U.S. Fish and Wildlife Service’s Stream Network Temperature Model (SNTEMP). Planning a temperature study is discussed in terms of understanding the management objectives and ensuring that the questions will be accurately answered with the modeling approach being used. A sensitivity analysis of SNTEMP is presented to illustrate which input variables are most important in predicting stream temperatures. This information helps prioritize data collection activities, highlights the need for quality control, focuses on which parameters can be estimated rather than measured, and offers a broader perspective on management options in terms of knowing where the biggest temperature response will be felt. All of the major input variables for stream geometry, meteorology, and hydrology are discussed in detail. Each variable is defined, with guidance given on how to measure it, what kind of equipment to use, where to obtain it from another agency, and how to calculate it if the data are in a form other than that required by SNTEMP. Examples are presented for the various forms in which water temperature, discharge, and meteorological data are commonly found. Ranges of values for certain input variables that are difficult to measure of estimate are given. Particular attention is given to those variables not commonly understood by field biologists likely to be involved in a stream temperature study. Pertinent literature is cited for each variable, with emphasis on how other people have treated particular problems and on results they have found. Model calibration, verification, and validation steps are defined and outlines, with measures of "goodness-of-fit" given for comparing simulated stream temperatures with observed values. The question of how good is good enough is explored, and attention is given to the kinds of simulation and data reduction errors that one should be alert for. Some special cases dealing with ice and reservoir temperature are mentioned. Special attention is given to understanding micro-thermal habitats that act as important thermal refugia under low flow conditions; their causes, extent, and management implications are discussed. Alternative public domain stream and reservoir temperature models are contrasted with SNTEMP. A distinction is made between steady-flow and dynamic-flow models and their respective capabilities. Regression models are offered as an alternative approach for some situations, with appropriate mathematical formulas suggested. Appendices provide information on State and Federal agencies that are good data sources, vendors for field instrumentation, and small computer programs useful in data reduction.

Report

Stochastic simulation of spatially correlated geo-processes

In this study, developments in the theory of stochastic simulation are discussed. The unifying element is the notion of Radon projection in Euclidean spaces. This notion provides a natural way of reconstructing the real process from a corresponding process observable on a reduced dimensionality space, where analysis is theoretically easier and computationally tractable. Within this framework, the concept of space transformation is defined and several of its properties, which are of significant importance within the context of spatially correlated processes, are explored. The turning bands operator is shown to follow from this. This strengthens considerably the theoretical background of the geostatistical method of simulation, and some new results are obtained in both the space and frequency domains. The inverse problem is solved generally and the applicability of the method is extended to anisotropic as well as integrated processes. Some ill-posed problems of the inverse operator are discussed. Effects of the measurement error and impulses at origin are examined. Important features of the simulated process as described by geomechanical laws, the morphology of the deposit, etc., may be incorporated in the analysis. The simulation may become a model-dependent procedure and this, in turn, may provide numerical solutions to spatial-temporal geologic models. Because the spatial simu??lation may be technically reduced to unidimensional simulations, various techniques of generating one-dimensional realizations are reviewed. To link theory and practice, an example is computed in detail. ?? 1987 International Association for Mathematical Geology.

Mathematical Geology

The community code verification exercise for simulating sequences of earthquakes and aseismic slip (SEAS)

Numerical simulations of sequences of earthquakes and aseismic slip (SEAS) have made great progress over past decades to address important questions in earthquake physics. However, significant challenges in SEAS modeling remain in resolving multiscale interactions between earthquake nucleation, dynamic rupture, and aseismic slip, and understanding physical factors controlling observables such as seismicity and ground deformation. The increasing complexity of SEAS modeling calls for extensive efforts to verify codes and advance these simulations with rigor, reproducibility, and broadened impact. In 2018, we initiated a community code-verification exercise for SEAS simulations, supported by the Southern California Earthquake Center. Here, we report the findings from our first two benchmark problems (BP1 and BP2), designed to verify different computational methods in solving a mathematically well-defined, basic faulting problem. We consider a 2D antiplane problem, with a 1D planar vertical strike-slip fault obeying rate-and-state friction, embedded in a 2D homogeneous, linear elastic halfspace. Sequences of quasi-dynamic earthquakes with periodic occurrences (BP1) or bimodal sizes (BP2) and their interactions with aseismic slip are simulated. The comparison of results from 11 groups using different numerical methods show excellent agreements in long-term and coseismic fault behavior. In BP1, we found that truncated domain boundaries influence interseismic stressing, earthquake recurrence, and coseismic rupture, and that model agreement is only achieved with sufficiently large domain sizes. In BP2, we found that complexity of fault behavior depends on how well physical length scales related to spontaneous nucleation and rupture propagation are resolved. Poor numerical resolution can result in artificial complexity, impacting simulation results that are of potential interest for characterizing seismic hazard such as earthquake size distributions, moment release, and recurrence times. These results inform the development of more advanced SEAS models, contributing to our further understanding of earthquake system dynamics.

Seismological Research Letters