Search

Simple Summary: Earthquakes are a common occurrence at The Geysers geothermal field in California which, over the years, have led to general land sinking in the area. In our study, we explore the correlation of geothermal data and seismicity and find that the injection rate exhibits a high correlation with the number of earthquakes. Additionally, we rely on satellite imagery which measures the extent of the land subsidence and develop a machine learning model that predicts future land subsidence, finding that it has a relatively low error rate. Moreover, we incorporate geothermal data into this model and find that it performs even better. All in all, this encourages the use of machine learning models in hazard mitigation in order to minimize the potential impacts that land subsidence might bring. The Geysers geothermal field in California is experiencing land subsidence due to the seismic and geothermal activities taking place. This poses a risk not only to the underlying infrastructure but also to the groundwater level which would reduce the water availability for the local community. Because of this, it is crucial to monitor and assess the surface deformation occurring and adjust geothermal operations accordingly. In this study, we examine the correlation between the geothermal injection and production rates as well as the seismic activity in the area, and we show the high correlation between the injection rate and the number of earthquakes. This motivates the use of this data in a machine learning model that would predict future deformation maps. First, we build a model that uses interferometric synthetic aperture radar (InSAR) images that have been processed and turned into a deformation time series using LiCSBAS, an open-source InSAR time series package, and evaluate the performance against a linear baseline model. The model includes both convolutional neural network (CNN) layers as well as long short-term memory (LSTM) layers and is able to improve upon the baseline model based on a mean squared error metric. Then, after getting preprocessed, we incorporate the geothermal data by adding them as additional inputs to the model. This new model was able to outperform both the baseline and the previous version of the model that uses only InSAR data, motivating the use of machine learning models as well as geothermal data in assessing and predicting future deformation at The Geysers as part of hazard mitigation models which would then be used as fundamental tools for informed decision making when it comes to adjusting geothermal operations. [ABSTRACT FROM AUTHOR]

Simple Summary: Excessive use of groundwater resources in the Central Valley in California has led to major land sinking over the years. In this study, we rely on satellite imagery to monitor and assess the extent of this sinking. Specifically, we use images from a satellite that emits microwaves which allows us to directly obtain the deformation at a specific time. Then, we apply a machine learning algorithm to the resulting data in an attempt to effectively predict short-term future deformation. We find that the algorithm applied has a low error when compared to the actual data. This shows that machine learning algorithms could be incorporated into models that assess potential hazards associated with land sinking. Land subsidence caused by excessive groundwater pumping in Central Valley, California, is a major issue that has several negative impacts such as reduced aquifer storage and damaged infrastructures which, in turn, produce an economic loss due to the high reliance on crop production. This is why it is of utmost importance to routinely monitor and assess the surface deformation occurring. Two main goals that this paper attempts to accomplish are deformation characterization and deformation prediction. The first goal is realized through the use of Principal Component Analysis (PCA) applied to a series of Interferomtric Synthetic Aperture Radar (InSAR) images that produces eigenimages displaying the key characteristics of the subsidence. Water storage changes are also directly analyzed by the use of data from the Gravity Recovery and Climate Experiment (GRACE) twin satellites and the Global Land Data Assimilation System (GLDAS). The second goal is accomplished by building a Long Short-Term Memory (LSTM) model to predict short-term deformation after developing an InSAR time series using LiCSBAS, an open-source InSAR time series package. The model is applied to the city of Madera and produces better results than a baseline averaging model and a one dimensional convolutional neural network (CNN) based on a mean squared error metric showing the effectiveness of machine learning in deformation prediction as well as the potential for incorporation in hazard mitigation models. The model results can directly aid policy makers in determining the appropriate rate of groundwater withdrawal while maintaining the safety and well-being of the population as well as the aquifers' integrity. [ABSTRACT FROM AUTHOR]

Nowcasting is a term originating from economics, finance, and meteorology. It refers to the process of determining the uncertain state of the economy, markets or the weather at the current time by indirect means. In this paper, we describe a simple two‐parameter data analysis that reveals hidden order in otherwise seemingly chaotic earthquake seismicity. One of these parameters relates to a mechanism of seismic quiescence arising from the physics of strain‐hardening of the crust prior to major events. We observe an earthquake cycle associated with major earthquakes in California, similar to what has long been postulated. An estimate of the earthquake hazard revealed by this state variable time series can be optimized by the use of machine learning in the form of the Receiver Operating Characteristic skill score. The ROC skill is used here as a loss function in a supervised learning mode. Our analysis is conducted in the region of 5° × 5° in latitude‐longitude centered on Los Angeles, a region which we used in previous papers to build similar time series using more involved methods (Rundle & Donnellan, 2020, https://doi.org/10.1029/2020EA001097; Rundle, Donnellan et al., 2021, https://doi.org/10.1029/2021EA001757; Rundle, Stein et al., 2021, https://doi.org/10.1088/1361-6633/abf893). Here we show that not only does the state variable time series have forecast skill, the associated spatial probability densities have skill as well. In addition, use of the standard ROC and Precision (PPV) metrics allow probabilities of current earthquake hazard to be defined in a simple, straightforward, and rigorous way. Plain Language Summary: Earthquake nowcasting refers to the determination of hazard for major earthquakes at the present time, the recent past, and the near future. Nowcasting is an idea borrowed from economics, markets, and meteorology, where it has been frequently used. In this paper, we show that there is order hidden within chaotic earthquake seismicity using a very simple transformation of the data. Small earthquakes appear to transition from unstable stick‐slip events that produce seismic waves, to stable sliding where no seismic waves are produced. Our hypothesis is that this transition is due to a material phenomenon called strain‐hardening, that is frequently observed in laboratory rock mechanics experiments. The result is a state variable time series, computed over the last 51 years in California, that strongly resembles the long‐anticipated cycle of stress accumulation and release. Using supervised machine learning techniques, we can optimize the two‐parameter model. From that optimized model, we can rigorously calculate the probability of current hazard from major earthquakes. Extending these methods, we can also compute spatial hazard as well. The result is a new method for assessing earthquake hazard that may be useful for a variety of applications. Key Points: "Chaotic" seismicity contains hidden structure in the form of state variable time seriesStandard data science methods can be used to convert the time series to probabilitiesBoth temporal and spatial probabilities can be computed [ABSTRACT FROM AUTHOR]