Coseismic Response Research Articles

Variations in Temperature and Pressure in the “Reservoir–Well” System Triggered by Blasting Recovery of Iron Ore at the Kursk Magnetic Anomaly

This paper presents the results of precise measurements of temperature and pore pressure in the “reservoir–well” system during the development of iron ore deposits of the Kursk Magnetic Anomaly (KMA) via blasting. For the observation period from October 2021 to June 2024, variations in compressibility, permeability and temperature in the upper Albian-Cenomanian confined aquifer, which is used for district water supply, were determined. The general trend in a decrease in water temperature was traced (from 12 °C to 11.4 °C). It was accompanied by an increase in the hydrostatic head (from 3.7 m to 7.4 m). Water temperature in the upper aquifer was measured for 9 industrial explosions in the mine and for 30 explosions in the quarry. For one explosion in the mine and five explosions in the quarry the coseismic changes in water temperature with amplitudes of 0.06–0.09 °C were established, while changes in pore pressure in the “reservoir–well” system were 0.4–2.2 kPa. Local changes in the permeability of the reservoir in the vicinity of the well (the skin effect) are considered to be the main factor that controls the coseismic response of temperature during industrial explosions. As the reservoir permeability increases, the water temperature in the “reservoir–well” system can decrease and vice versa. The same pattern was observed according to regime measurements performed in 2022–2023. The recorded coseismic responses of water temperature in the upper aquifer in the high-frequency range are similar to the effects observed during propagation of seismic waves originated from earthquakes in the low-frequency range at different sites all over the world for the seismic energy density of 0.05–0.45 J/m3. The observed variations in aquifer temperature in the “reservoir–well” system under episodic dynamic impacts are of particular interest from the point of view of activating hydrogeochemical processes that accompany the development of iron ore deposits.

Capacity curves for monitored existing buildings and within- and between-building variability of secant stiffness

In this study, accelerometric data from seven Japanese buildings under long-term monitoring were analysed to explore the variability of the buildings’ co-seismic response over time and its within- and between-building components, using co-seismic capacity curves developed in acceleration-displacement-response-spectrum format. The data include the 2011 Tohoku Mw9.1 earthquake, which caused building damage of different levels of severity, and the time-varying actual capacity curves were analysed considering earthquakes before and after 2011. Result showed that the initial slope of the capacity curves reflects the amount of damage. The between-building and within-building components of the variability are discussed by comparing a single building and several buildings in the same class for several earthquakes. Finally, the epistemic uncertainty of seismic risk assessment studies is discussed in relation to the selection of a generic capacity model for all buildings in a single class.

Quasi-coseismic variations and geosphere coupling associated with the strong 2023 Turkey earthquakes

On 6 February 2023, two massive shallow earthquakes of Mw 7.8 and Mw 7.5 occurred successively in Turkey. Unlike earlier studies on pre-seismic anomalous signals, here we focus on the co-seismic changes in ocean, land, and atmospheric parameters associated with these two earthquakes. The quasi-coseismic variations have been clearly observed in ocean, atmosphere, and snow parameters from satellite and reanalysis datasets. Our results show a decline in the seawater temperature and salinity, and enhancement in the chlorophyll-a concentration in the eastern Mediterranean Sea associated with the earthquakes. Over the epicentral region, the total ozone column and snowfall have been observed to increase with the occurrence of earthquakes. After a detailed analysis of the hourly parameters (atmospheric pressure, air temperature, relative humidity, and wind speed), we found the atmospheric pressure variation caused by the ground motion of seismic waves is the possible reason for the co-seismic changes. Based on the results discussed in this paper, a model of multi-geosphere co-seismic response associated with the 2023 Turkey events is proposed. Additionally, our results show that the crustal motion in the high mountain region within the earthquake preparation zone could be more sensitive prior to the earthquake due to the enhancement of electric field as suggested by the theory of positive holes.

Kinematic deformation and intensity assessment of the 2021 Maduo MS7.4 earthquake in Qinghai revealed by high-frequency GNSS

Rapid acquisition of the kinematic deformation field and seismic intensity distribution of large earthquakes is crucial for postseismic emergency rescue, disaster assessment, and future seismic risk research. The advancement of GNSS observation and data processing makes it play an important role in this field, especially the high-frequency GNSS. We used the differential positioning method to calculate the 1 HZ GNSS data from 98 sites within 1000 km of the MS7.4 Maduo earthquake epicenter. The kinematic deformation field and the distribution of the seismic intensity by using the peak ground velocity derived from displacement waveforms were obtained. The results show that: 1) Horizontal coseismic response deformation levels ranging from 25 mm to 301 mm can be observed within a 1000 km radius from the epicenter. Coseismic response deformation on the east and west sides shows bilateral asymmetry, which markedly differs from the symmetry presented by surface rupture. 2) The seismic intensity obtained through high-frequency GNSS and field investigations exhibits good consistency of the scope and orientation in the high seismic intensity area, although the former is generally slightly smaller than the latter. 3) There may exist obstacles on the eastern side of the seismogenic fault. The Maduo earthquake induced a certain tectonic stress loading effect on the western Kunlun Pass-Jiangcuo fault (KPJF) and Maqin-Maqu segment, resulting in higher seismic risk in the future.

Reconciling the conflicting extent of overriding plate deformation before and during megathrust earthquakes in South America, Sunda and northeast Japan

SUMMARYWe aim to better understand the overriding plate deformation during the megathrust earthquake cycle. We estimate the spatial patterns of interseismic GNSS velocities in South America, Southeast Asia and northern Japan and the associated uncertainties due to variations in network density and observation uncertainties. Interseismic velocities with respect to the overriding plate generally decrease with distance from the trench with a steep gradient up to a ‘hurdle’, beyond which the gradient is distinctly lower and velocities are small. The hurdle is located 500–1000 km away from the trench for the trench-perpendicular velocity component, and either at the same distance or closer for the trench-parallel component. Significant coseismic displacements were observed beyond these hurdles during the 2010 Maule, 2004 Sumatra–Andaman, and 2011 Tohoku earthquakes. We hypothesize that both the interseismic hurdle and the coseismic response result from a mechanical contrast in the overriding plate. We test our hypothesis using physically consistent, generic, 3-D finite element models of the earthquake cycle. Our models show a response similar to the interseismic and coseismic observations for a compliant near-trench overriding plate and an at least five times stiffer overriding plate beyond the contrast. The model results suggest that hurdles are more prominently expressed in observations near strongly locked megathrusts. Previous studies inferred major tectonic or geological boundaries and seismological contrasts located close to the observed hurdles in the studied overriding plates. The compliance contrast probably results from thermal, compositional and thickness contrasts and might cause the observed focusing of smaller-scale deformation like backthrusting.

Responses of well-aquifer system to four earthquakes and implication for contaminant transport in Changan landfill, China

Evaluation of aquifer response to earthquakes is important to understanding the evolution of aquifer properties and contaminant transport. In Changan Landfill, a multi parameter monitoring system was installed to collect data every 12 h. Principal Component Analysis (PCA), energy density and wavelet analysis, and tidal analysis were used to reveal the response mechanisms to four earthquakes, sensitivity of various parameters to seismic stimuli, and permeability evolution, respectively. The results showed that the chemical characteristics of two aquifers were distinct. The shallow aquifer was polluted by leachate, and the deep aquifer was less polluted. Coseismic responses were dominated by the deep aquifer and the relative contribution of each aquifer remained stable. Coseismic chemical changes were dominated by the opening and closing of fractures, and the long-term evolution was controlled by seasonality. The deep aquifer had a higher sensitivity to seismic stimuli than the shallow aquifer. Tidal analysis showed that a general permeability decrease during this monitoring period reduced contaminant transport, but some contaminants crossed the aquitard between the two aquifers after the Zizhong earthquake, due to an increase in vertical permeability.

Probabilistic seismic source inversion of the 1886 Charleston, South Carolina, earthquake from macroseismic evidence: A major updating

The source of the 1886 Charleston, South Carolina earthquake influences the computed seismic hazard of the Southeastern U.S. and thus impacts public policy and engineering practice. However, because the 1886 rupture predated seismic instruments, its source is highly uncertain. This study presents probabilistic seismic-source inversions of the Charleston earthquake from liquefaction evidence and historical intensity reports. Using the latest predictive models and a novel inversion approach, we seek to constrain the magnitude, location, and orientation of the 1886 rupture. Probability distributions of rupture magnitude are conditioned on both the “Woodstock Fault” – a commonly inferred source of the 1886 event – and on an unknown source, wherein the uncertainties of fault location and orientation are considered. These distributions are compared to the Mw6.7-Mw7.5 distribution adopted by the U.S. National Seismic Hazard Model Project (NSHMP). Collectively, the results do not provide strong support for the hypothesized Woodstock Fault. This is not to say the Woodstock Fault does not exist or could not have caused the 1886 rupture, but rather, that the position of the 1886 source model cannot be constrained by the data and models studied herein, given the large uncertainties inherent to each. While this is at odds with the underlying assumption of prior studies, the results nonetheless generally uphold the magnitude distribution adopted by the NSHMP. The largest uncertainties inherent to this distribution are identified and could be diminished in the future. Finally, we note that the inversion methodology used here is not specific to any region, or to certain types of evidence, but can be applied to any seismic zone and to any co-seismic response. This methodology allows for uncertainty to be accounted for in a more complete and transparent manner when inverting seismic source parameters from macroseismic data. Of course, any limitations, biases, or unmodeled uncertainties inherent to these data must be understood, and their implications acknowledged, as further discussed herein.

Responses of Water Temperature and Level to Large Earthquakes in Tibet

Well water levels can reflect the stress placed on a confined subsurface aquifer system in a similar manner to a strain meter. Based on observations of the geophysical field in Lhasa combined with digital data recorded at an underground fluid well at the Lhasa geomagnetic station in recent years, we comprehensively analyzed the characteristics of co-seismic changes caused by 14 different-magnitude M ≥ 5 earthquakes recorded in the well. The results show that (1) the co-seismic changes in water temperature and water level are different; the water level exhibits oscillation-type changes, while the water temperature variations indicate first heating and subsequent recovery. (2) The co-seismic changes are related to the epicentral distance, magnitude and focal depth of the earthquake. The closer the epicenter is to the well, the earlier the co-seismic responses occur, but the time interval between the co-seismic changes in the water level and temperature differs. (3) The co-seismic ratio of the water temperature is higher than that of the water level; this may be related to faulty water level instrumentation or segmented records.

A New Mechanism for Earthquake‐Enhanced Permeability

AbstractEarthquakes may enhance the permeability of Earth's crust, but the mechanisms for such changes are not entirely clear. Recent studies have revealed a distinct coseismic response at large epicentral distances to great earthquakes and at close distances to major earthquakes , where both the horizontal and vertical permeabilities increased by orders of magnitude and the increases stayed at the co‐seismic level for many years, which cannot be explained by the existing models. To explain such prolong increases of the vertical permeability I propose a new mechanism, supported by laboratory and field observations of earthquake‐triggered hydrogeologic processes, that great earthquakes may pressurize and liquefy groundwater systems and inject the liquefied sands into fractures to keep them open. The proposed mechanism is testable and implies prolonged upward migration of hot fluids from deep basins, which may be an important process for subsurface transport of heat and solutes and may also be useful to understand the safety of groundwater and the security of underground repositories.

Multi Field Coupled Coseismic Changes of the Jiashi Ms = 6.4 Earthquake of 19 January 2020, Based on Ground Temperature Observation

The coseismic geothermal changes of ground temperature observed at observatories near the epicenter of the 2020 Jiashi Ms = 6.4 earthquake in China, provide a unique opportunity to study heat generation and conduction in rock. Here, evolutions of rock temperature at the Xikeer, Jiashizongchang, and Gedaliang observatories, which are located at epicentral distances of 1.4, 27.42, and 50 km respectively, were analyzed. Significant coseismic geothermal changes of 0.0432 °C were observed at the Xikeer observatory at the depth of 33.38 m, at which clear diurnal variations can be observed. Smaller changes of ~0.0001 °C were observed at the depths of 12.3 and 22.8 m at the Xikeer observatory and 22.3 m at the Jiashizongchang observatory. The stress transfer induced by the coseismic rupture induced a rise in local ground temperature, but the magnitude of the change was relatively small. The larger amplitude change at the Xikeer observatory was caused by fluid infiltration. We note that diurnal variation has been recorded at the Gedaliang observatory, but the coseismic response is no longer in existence. The temperature increases at the hypocentral area were higher than expected in the ground due to the coseismic stress transfer, but the change attenuated rapidly with distance.

Automated near-field deformation detection from mobile laser scanning for the 2014 Mw 6.0 South Napa earthquake

Abstract Quantifying off-fault deformation in the near field remains a challenge for earthquake monitoring using geodetic observations. We propose an automated change detection strategy using geometric primitives generated using a deep neural network, random sample consensus and least squares adjustment. Using mobile laser scanning point clouds of vineyards acquired after the magnitude 6.0 2014 South Napa earthquake, our results reveal centimeter-level horizontal ground deformation over three kilometers along a segment of the West Napa Fault. A fault trace is detected from rows of vineyards modeled as planar primitives from the accumulated coseismic response, and the postseismic surface displacement field is revealed by tracking displacements of vineyard posts modeled as cylindrical primitives. Interpreted from the detected changes, we summarized distributions of deformation versus off-fault distances and found evidence of off-fault deformation. The proposed framework using geometric primitives is shown to be accurate and practical for detection of near-field off-fault deformation.

Various far-field hydrological responses during 2015 Gorkha earthquake at two distant wells

Aquifer hydraulic parameter can change during earthquakes. Continuous monitoring of the response of water level to seismic waves or solid Earth tides provides an opportunity to document how earthquakes influence hydrological properties. Here, we use data of two groundwater wells, Dian-22 (D22) and Lijiang (LJ) well, in southeast Tibet Plateau in response to the 2015 Mw 7.8 Gorkha earthquake to illustrate hydrological implications. The coherences of water level and seismic wave before and after the far-field earthquake show systematic variations, which may confirm the coseismic dynamic shaking influence at high frequencies (f > 8 cpd). The tidal response of water levels in these wells shows abrupt coseismic changes of both phase shift and amplitude ratio after the earthquake, which may be interpreted as an occurrence in the vertical permeability of a switched semiconfined aquifer in the D22 well, or an enhancement unconfined aquifer in the LJ well. Using the continuous short-term transmissivity monitoring, we show that the possible coseismic response for about 10 days and instant healing after 10 days to the causal earthquake impact. Thus, the dynamic shaking during the Gorkha earthquake may have caused the short-term aquifer responses by reopening of preexisting vertical fractures and later healing at epicentral distances about 1500 km.

Relative Sea Level Changes and Morphotectonic Implications Triggered by the Samos Earthquake of 30th October 2020

On 30th October 2020, the eastern Aegean Sea was shaken by a Mw = 7.0 earthquake. The epicenter was located near the northern coasts of Samos island. This tectonic event produced an uplift of the whole island as well as several cases of infrastructure damage, while a small tsunami followed the mainshock. Underwater and coastal geological, geomorphological, biological observations and measurements were performed at the entire coast revealing a complex character for the uplift. At the northwestern part of the island, maximum vertical displacements of +35 ± 5 cm were recorded at the northwestern tip, at Agios Isidoros. Conversely, the southeastern part was known for its subsidence through submerged archaeological remains and former sea level standstills. The 2020 underwater survey unveiled uplifted but still drowned sea level indicators. The vertical displacement at the south and southeastern part ranges between +23 ± 5 and +8 ± 5 cm suggesting a gradual fading of the uplift towards the east. The crucial value of tidal notches, as markers of co-seismic events, was validated from the outcome of this study. The co-seismic response of Samos coastal zone to the 30th October earthquake provides a basis for understanding the complex tectonics of this area.

Modelling coseismic displacements of fracture systems in crystalline rock during large earthquakes: Implications for the safety of nuclear waste repositories

Crystalline rock has been tested/selected by many countries for developing underground nuclear waste repositories at ~500 m depth to achieve geological isolation of high-level, long-lived radioactive waste. During the assessment period of up to one million years, large earthquakes may occur around the repository and trigger coseismic displacements along secondary fractures, jeopardising the integrity of the buffer-waste canister system. It is, therefore, of great importance to understand the coseismic behaviour of the repository site during large earthquakes. In this paper, we develop a finite element method-based seismo-mechanical model to simulate the response of fractured rocks subject to both in-situ stresses and seismic activities. The model can well capture the fracture displacement behaviour under dynamic loadings involving alternating regimes of contact loss, partial slip and total sliding. We model the earthquake-induced displacement field according to the seismic source theory in combination with a generic source time function. We apply the model to the nuclear waste repository site at Forsmark, Sweden, and analyse the coseismic responses of both single fractures and fracture networks during a potential post-glacial earthquake with a moment magnitude of Mw = 5.6. The shear dislocation of a single fracture is strongly dependent on the fracture length and dip angle, while the displacement pattern of a fracture network is dominated by its “backbone” structures. We observe that significant coseismic shear displacement occurs if the fractured rock is close to the hypocentre and located in the dilational quadrant of the primary fault due to reduced shear resistance. However, the earthquake-induced shear displacement decreases drastically with the increased distance to the hypocentre and a distance of ~700 m may be needed for a fracture up to 100 m long to not displace beyond 5 cm. If two repeated earthquakes would occur, a distance of ~1200 m may be necessary.

Modeling Earthquake‐Induced Spring Discharge and Temperature Changes in a Fault Zone Hydrothermal System

AbstractAlthough different methods have been employed to study the earthquake‐induced changes in aquifer permeability (e.g., tidal/barometric pressure response and analytical models), identification of the mechanism is usually hindered by the limitation of application conditions or multiplicity of solutions. A numerical simulation that couples flow and heat transport is developed in the present study to explain the coseismic responses that followed two large earthquakes (Mw 7.0 Lijiang earthquake on 3 February 1996 and Mw 9.1 Sumatra earthquake on 26 December 2004) at Banglazhang #1 hot spring, China. The results show that decreased spring discharge and temperature are attributable to earthquake‐induced permeability decrease, which can be related to the blocking of narrow fracture apertures. Analysis of the model revealed that seismic waves may have caused heterogeneous changes in permeability within the fault zone hydrothermal system, which led to nonsynchronous changes in spring flow and temperature. Our study contributes to an enhanced understanding of the impact of the evolution of earthquake‐induced permeability in hydrothermal systems.

An improved method of Newmark analysis for mapping hazards of coseismic landslides

Abstract. Coseismic landslides can destroy buildings, dislocate roads, sever pipelines, and cause heavy casualties. It is thus important but challenging to accurately map the hazards posed by coseismic landslides. Newmark's method is widely applied to assess the permanent displacement along a potential slide surface and model the coseismic response of slopes. This paper proposes an improved Newmark analysis for mapping the hazards of coseismic landslides by considering the roughness and effect of the size of the potential slide surfaces. This method is verified by data from a case study on the 2014 Mw 6.1 (the United States Geological Survey) Ludian earthquake in Yunnan Province, China. Permanent displacements due to the earthquake ranged from 0 to 122 cm. The predicted displacements were compared with a comprehensive inventory of landslides triggered by the Ludian earthquake to map the spatial variation in the hazards of coseismic landslides using the certainty factor model. The confidence levels of coseismic landslides indicated by the certainty factors ranged from −1 to 0.95. A hazard map of the coseismic landslide was generated based on the spatial distribution of values of the certainty factor. A regression curve relating the predicted displacement and the certainty factor was drawn, and can be applied to predict the hazards of coseismic landslides for any seismic scenario of interest. The area under the curve was used to compare the improved and the conventional Newmark analyses, and revealed the improved performance of the former. This mapping procedure can be used to predict the hazards posed by coseismic landslides, and provide guidelines for decisions regarding the development of infrastructure and post-earthquake reconstruction.

Groundwater trace elements change induced by M5.0 earthquake in Yunnan

Geochemical compositions changes in groundwater in response to earthquakes have been documented for decades, and most such studies have documented the changes in major ions or geochemical gases (e.g., Radon). Changes in trace elements, however, have received less attention. In this paper, we collected groundwater samples from a wellbore in Yunnan Province weekly and analyzed both the major and trace chemical compositions from February 2018 to September 2018. Combined with well water level records, we found that the August 2018 Tonghai M5.0 earthquake caused changes in both water levels and trace elements but caused no significant changes in the stable isotope and major ions. Such a response indicates that trace elements are more sensitive to crustal strain than are the major elements and isotopes. And We considered that earthquake-induced shaking caused a breach in the hydraulic barrier of the isolated reservoir that led to the mixing of water with low concentrations of trace elements as the mechanism to explain the coseismic response. This study demonstrate that some trace elements are more sensitive to earthquake strain than the commonly used chemical constitutents and that trace element monitoring should considered to be a priority in the future to provide a complete assessment of the sensitivity of different chemical constitutents to earthquake stresses.

Sensitivity of the response of well-aquifer systems to different periodic loadings: A comparison of two wells in Huize, China

Groundwater well can be used as a seismograph or as strainmeter, since the well water level is sensitive to seismic waves and other crustal strains through mechanical coupling between rock and groundwater. The use of the groundwater level in monitoring tectonic strain, however, is hindered by the enigmatic response discrepancy among wells (different responses may be recorded by even closely clustered wells). Therefore, understanding the factors that control the sensitivity of well-aquifer systems in response to different periodic loadings is important. In this paper, we analyzed two closely located wells that show no response to seismic waves and respond differently to earth tides and barometric pressure. We discuss these wells’ sensitivity of response to periodic loadings by using the spectrum, barometric response function, tidal response, and slug test methods. Our results show that the different response characteristics are due to the differences in permeability and confinement of the well-aquifer systems. Well-aquifer systems with a low permeability cannot record coseismic responses, even if they can record clear tidal and barometric signals. In addition, periodic loading signals in shallow aquifers will be easily masked by the interference from the surface, for example rainfall or pumping, which diminishes the signal-to-noise ratio of the tectonic activity. Thus, the hydrogeological and tectonic settings are the two key factors that should be considered in choosing a suitable monitoring well location for tectonic strain monitoring.

Coseismic Groundwater Temperature Response Associated with the Wenchuan Earthquake

Various physical, geophysical, geochemical, electrical and hydrological parameters are measured on the surface and in shallow/deep boreholes throughout mainland China to obtain early warning signals of impending earthquakes. Numerous wells are equipped with water level and temperature sensors for continuous observations of the water level and groundwater temperature. An analysis of water temperature data from boreholes equipped with water temperature sensors reveals that nearly half of the boreholes show coseismic response associated with the Wenchuan earthquake (Mw 7.9, 12 May, 2008). The coseismic response of the groundwater temperature cannot be differentiated from the groundwater flow or movement when the earthquake occurred, but there is no fixed relationship between the temperature variation and the water flow. At the same time, we observed that the rock temperature in dry wells can record the seismic events and even the pre-seismic abnormal information. The spatial distribution of the coseismic groundwater temperature response is random and irregular, which does not support the dislocation model of seismic faults at the regional or larger scale. Changes in the groundwater temperature are closely related to the borehole temperature gradient, lithology profile and geological environment of the borehole and depths of the aquifers. The mechanism of the coseismic groundwater temperature response can be explained by an enhanced permeability induced by an earthquake. The groundwater temperature increases if the temperature sensor in the borehole is located near the deep-circulating aquifer and decreases if the sensor is near the shallow-circulating aquifer when seismic waves arrive. The groundwater temperature may be slightly affected or even unchanged if the temperature sensor is far from the aquifer during the propagation of seismic waves. However, it was hard to conclude the changes of rock temperature observed in six dry wells.

Groundwater level response to the Wenchuan earthquake of May 2008

We have comprehensively analysed the co-seismic response of the groundwater levels of 280 wells in mainland China that were associated with the Wenchuan earthquake (Mw 7.9) that occurred on 12 May 2008. The observed co-seismic responses can be classified as step-like changes in 138 wells, variations in 69 wells and non-responses in 73 wells. After a quantitative analysis of spatial distribution, there was no spatially coherent signal found in the step-like changes (positive values indicate a step-like rise, and negative values indicate a step-like fall), even within 300 km of the epicenter. The amplitude and the phase shift of the M2-wave were compared between the pre- and post-earthquake conditions. The phase was forward and was concentrated at , regardless of the proximity of wells to the epicenter; however, the change in amplitude randomly increased or decreased. By computing the post-seismic groundwater recession, the characteristic times of the aquifers were . We concluded that by assuming a first-order approximation, i.e. a one-dimensional aquifer, the causal mechanism of the co-seismic response of the groundwater level was that the seismic waves enhanced rock permeability by clearing the facture-filling materials. The shapes of the co-seismic responses were determined by the well-aquifer system and the proximity of the wells to recharge or discharge areas (). The water level rose if the well was closer to the recharge area, the water level decline if the well was closer to the discharge area, and the water level oscillated if the well was farther from the recharge or discharge areas when the seismic wave was transmitted. The water level remained unchanged if the well did not penetrate any confined aquifer.