Electrical characteristics of winter lightning

To make clear the characteristics of winter lightning flashes, the long term observation of winter lightning on the west coast of Japan has been continued using magnetic links, digital recording systems for the current oscillograms, field mills, still cameras and video camera systems for the images of lightning channels. Of the 66 magnetic links records exceeding 2 kA, 73% were negative in polarity, and 27% were positive. Median peak values for winter flashes were 17 kA for negative flashes and 45 kA for positive flashes, respectively. 145 current waveforms for the winter lightning flashes, which have current amplitudes exceeding 1 kA, have been obtained by the shunt systems and/or the coils system. They show that winter lightning flashes often have a very long duration or continuing current, and sometimes have a very large amplitude exceeding 200 kA in positive flashes. As one example flash on 9 January 1987, the maximum current amplitude was +280 kA, the maximum current derivative 1.0 × 10¹⁰ As⁻¹, the total charge +400 C and the action integral 1.5 × 10⁷ A²s. The winter lightning current waveforms are classified into three types: single stroke flashes, monopolar multiple stroke flashes and bipolar flashes. Moreover, each flash is categorized as positive or negative, single peak or multiple peak, and with or without a continuing current.     

The equatorial electrojet and associated currents as seen in Magsat data

Magsat data are re-examined with regard to the presence and character of fields due to the equatorial electrojet and meridional currents at dawn and dusk local times. Dip-latitude organized field variations at dawn are:

• 1.(1) extremely weak,
• 2.(2) extremely variable with longitude,
• 3.(3) inconsistent with the pattern expected from a line or narrow sheet current.

It is shown that the use of Magsat dusk data can ‘contaminate’ a main field model, introducing apparent equatorial electrojet effects into the dawn data.Fields due to the equatorial electrojet and (presumably) associated meridional currents are clearly present in the dusk data. They show a variation with longitude which is apparently associated with the longitudinal variation of the strength, or square of the strength, of the main field in the E-region. Also evident is a variation with time of the year, although data are available for only a six month period. The meridional currents are generally minimum during January and February and maximum either during November and December or March and April, depending upon longitude. The E-region horizontal currents are minimum in November and December and maximum in March and April, except for − 30° to −90° longitude when the maximum occurs in January and February.Assuming that field gradients in local time are considerably smaller than field gradients in dip-latitude, current densities are estimated to be 1–3.6μA/m² for the horizontal current at 110km and about 10–20 × 10⁻⁹ A/m² for the vertical currents at 400km altitude. These results confirm and extend earlier results of Takeda and Maeda.Most models of the electrojet system in the literature disagree severely with these measurements either because their scope is inadequate or because of the wind system they assume. Those models which best describe the data invoke an eastward wind and/or an eastward electric field at dusk local time.     

An experimental investigation of atmospheric electricity and lightning activity to be performed during the descent of the Huygens Probe onto Titan

No terrestrial-like electrical activity was observed during the Voyager 1 flyby of Titan on 12 November 1980, in spite of a predicted global lightning energy dissipation rate of 4 × 10⁻⁶ Wm⁻². This lack of evidence does not, however, rule out the existence of electrical discharges with magnitudes, rates of occurrence and spectral characteristics drastically different from those known on Earth, owing to large dissimilarities between the temperatures, chemical compositions and, especially, electrical conductivities of the two atmospheres. Towards the end of the year 2004, the ESA probe Huygens will be jettisoned from the NASA Saturn orbiter, Cassini. This probe will descend onto Titan and perform in situ measurements during a period of 3 h, from an altitude of 170 km down to the satellite surface where the atmospheric pressure reaches 1.6 × 10⁵ Pa. The Huygens scientific payload will include a set of instruments entirely dedicated to the detection of lightning and to the characterization of the electrical properties of the atmosphere and surface. An electric antenna will search for natural emissions in the frequency range 0–10 kHz, at altitudes lower than those of ionized layers opaque to electromagnetic waves, and measure the magnitude of static electric fields due to charge separation. The conductivity of the atmosphere and the existence of free electrons will be checked during the whole descent with a combination of quadrupolar and relaxation probes; a microphone will also record acoustic phenomena associated with electrical discharges and atmospheric processes. The impedance of the surface will be evaluated from the measurements collected with a radar during the descent and a quadrupolar array after touch down.     

Riming electrification mechanism for charge generation within a thundercloud

A time dependent model for the rate of growth of the electric field within a thundercloud by the process of graupel formation (riming electrification) has been presented. The parameters of the model are: p₀, the precipitation intensity; q, the charge acquired by the graupel in each collision it makes with an ice crystal and 〈p〉, the charge transfer efficiency. Sets of values for the field growth and the maximum field have been obtained by varying these parameters.It is found that the estimated field for the initiation of a lightning flash (3.4 × 10⁵Vm⁻¹) within a time of about 1200 s can be achieved, taking a reasonable value of q = −1.67 × 10⁻¹⁴ C only for P₀ ≳ 5.56 × 10⁻⁶ms⁻¹(20 mmh⁻¹) and 〈p〉 ≳ 0.5. The maximum attainable electric field, E_max, itself is not sensitive to the value of q, within a reasonable range of variation in it (for a given p₀), but its rate of growth is: it grows faster, if g is larger.     

The EISCAT mesospheric measurements during the CAMP campaign

A brief outline is given of the experimental technique used during the Cold Arctic Mesopause Project to record the first D-region ion line spectra with the EISCAT incoherent scatter radar. The data analysis shows that echoes from mesospheric heights between about 70 km and 90 km can be detected during disturbed periods of enhanced electron density during particle precipitation events. Electron density profiles were determined which show a fairly high density, up to 5 × 10¹⁰ m⁻³ in the upper D-region. The measured meridional winds were lower than 10 m s⁻¹. A fit of the measured height profile of spectral width to temperature and neutral density models yielded a measured temperature profile in good agreement with simultaneous rocket data. The mesopause temperature was determined to be as low as 130 K. This detailed analysis of the spectral width profile indicates that below about 77–80 km the ratio of negative ions to electrons exceeded unity. Finally, some discussions are added on the limitations and significance of these first mesosphere observations.     

Constraining groundwater flow in the glacial drift and saginaw aquifers in the Michigan Basin through helium concentrations and isotopic ratios

³He and ⁴He concentrations in excess of those in water in solubility equilibrium with the atmosphere by up to two and three orders of magnitude are observed in the shallow Glacial Drift and Saginaw aquifers in the Michigan Basin. A simplified He transport model shows that in situ production is negligible and that most He excesses have a source external to the aquifer. Simulated results show that ³He and ⁴He fluxes entering the bottom of the Saginaw aquifer are 7.5 × 10^?14 and 6.1 × 10^?7 cm³STPcm^?2 yr^?1, both of which are lower than fluxes entering the underlying Marshall aquifer, 1.0 × 10^?13 and 1.6 × 10^?6 cm³STPcm^?2 yr^?1 for ³He and ⁴He, respectively. In contrast, He fluxes entering the Saginaw aquifer are higher than fluxes entering the overlying Glacial Drift aquifer of 5.2 × 10^?14 and 1.5 × 10^?7 cm³STPcm^?2 yr^?1 for ³He and ⁴He, respectively. The unusually high He fluxes and their decreasing values from the lower Marshall to the upper Glacial Drift aquifer strongly suggest the presence of an upward cross‐formational flow, with increasing He dilution toward the surface by recharge water. These fluxes are either comparable to or far greater than He fluxes in deeper aquifers around the world. Model simulations also suggest an exponential decrease in the horizontal groundwater velocity with recharge distance. Horizontal velocities vary from 13 to 2 myr^?1 for the Saginaw aquifer and from 18 to 6 myr^?1 for the Marshall aquifer. The highly permeable Glacial Drift aquifer displays a greater velocity range, from 250 to 5 myr^?1. While Saginaw ⁴He ages estimated based on the simulated velocity field display an overall agreement with ¹⁴C ages, ¹⁴C and ⁴He ages in the Glacial Drift and Marshall aquifers deviate significantly, possibly due to simplifications introduced in the He transport model leading to calculation of first‐order approximation He ages and high uncertainties in Glacial Drift ¹⁴C ages.     

E-region nitric oxide concentrations inferred from the 26 February 1979 eclipse expedition

Nitric oxide (NO) concentrations have been determined from the analysis of positive ion composition data obtained by AFGL for eclipse and post-eclipse conditions near Red Lake, Canada. Values of about 3 × 10⁸ cm⁻³ for 105–110 km and about 3 × 107 cm⁻³ for 90 km have been established. A residual ion-pair production rate of about 50 cm⁻³ s⁻¹ is estimated for eclipse totality at 110 km. A factor of two uncertainty is thought to be appropriate for all deduced values. The calculated NO concentrations appear to be within the range of typical variations for this season (late winter) and latitude (51°N).     

Positive and negative ion composition measurements in the D- and E-regions during the 26 February 1979 solar eclipse

Two rockets bearing quadrupole mass spectrometers capable of measuring both positive and negative ion composition were launched from Red Lake, Canada, during the solar eclipse. Both instruments had liquid helium cryopumps and shock-attaching conical samplers. The payloads also contained two Gerdien condensers to measure total positive and negative ion concentrations and ion mobilities. Attitude control systems aligned the payloads with the velocity vector throughout ascent and descent. The first rocket was launched so that the D-region was in darkness 35 ± 8 s on the upleg and about 150 ± 15 s on the downleg for the study of ionospheric decay processes. The second rocket was fired after totality into 75% solar illumination for the study of ionospheric recovery. The positive ion composition above 105 km exhibited a strongly increasing NO⁺/O₂⁺ ratio with time after second contact due to O₂⁺ charge transfer with NO and a sharply diminished ionization rate. However, in both nights, the ionization below 105 km was created mainly by energetic particle deposition as exemplified by the increased ion concentrations and the composition signatures of a particle event: asignificant enhancement of O₂⁺ below 105 km and large amounts of H₅O₂⁺ ions in the D-region which result from the O₂⁺ clustering scheme. H₅O₂ was the major ion in the upper D-region while H₇O⁺₃, H₉O₄⁺ and H₅O₂⁺ were dominant ions at lower altitudes. Numerous minor species were also detected. The negative ion distributions in both flights exhibited a distinct shelf at 83 ± 2 km, decreasing by more than an order of magnitude by 90 km and with minima near 75 km. In the 75–90 km range, a significant percentage of the negative ions had masses exceeding 160 a.m.u. Comparisons are made with prior negative ion measurements during similar daytime auroral zone absorption (AZA) events. Two striking characteristics of the precipitating particles were apparent from these and past observations in daytime AZA events: there is a near absence of low energy electrons capable of ionizing above about 105 km and there is'a significant spatial and/or temporal variability in the electron flux. This paper is devoted principally to a presentation of the ion composition measurements and associated uncertainties.     

Magnetopause plasma transients: signatures in the auroral ionosphere

The magnetopause and adjacent boundary layers of the Earth's magnetosphere play important roles in transferring momentum and energy from the solar wind to the magnetosphere-ionosphere system. The details of the different boundary processes, their ionospheric signatures and relative importance are not well known at present. Particle precipitation, field-aligned current, auroral emission, ionospheric ion drift and ground magnetic perturbations are among the low-altitude parameters that show signatures of various plasma processes in the LLBL and the magnetopause current layer. Magnetic merging events, Kelvin-Helmholtz waves, and pressure pulses excited by the variable solar wind/magnetosheath plasma are examples of boundary phenomena that may be coupled to the ionosphere via field-aligned currents. In this paper, attention is focussed on a specific category of auroral activity occurring in the cusp/cleft region predominantly during the southward directed interplanetary magnetic field (IMF). Co-ordinated observations from the ground and satellites in polar orbit have been used to study the temporal/spatial development of the events in relation to the background patterns of particle precipitation and ionospheric convection as well as the field-aligned current and ion drift characteristics of the individual events. The auroral phenomenon is characterized by a sequence of elongated forms moving laterally into the polar cap. Spatial scales of major events repeating every 5–10 min are ∼200 km (N-S) times 300–1000 km (E-W). Smaller scale auroral structures with more irregular occurrence rates are observed at times. The preliminary evidence suggests that the motion pattern is regulated by the IMF orientation, that is, the direction of longitudinal motion along the polar cap boundary is determined by the IMF B_Y polarity. The examples reported here occurred within 1000–1400 MLT, near the zero point potential line separating the morning and post-noon convection cells. During nonzero IMF B_Y the auroral structures are associated with channels of enhanced zonal ionospheric ion flow and Birkeland current sheets of opposite polarity, imbedded within the larger scale IMF B_Y-related cusp-mantle current system. These characteristics are discussed in relation to model predictions of ionospheric signatures of magnetopause plasma transients, with particular emphasis placed on impulsive magnetic merging events.     

Measurement of O2(a1Δg) emission in a total solar eclipse

The altitude distribution of the oxygen infrared atmospheric bands at 1.27 μm was measured during the total solar eclipse of 26 February 1979. The ozone concentration profile has been derived from these airglow measurements and indicates that at 85 km the concentration at totality was 7 × 1.7 cm⁻³, with no well defined upper layer. This reduced concentration, which is typical of summertime conditions, was probably due to perturbations in the mesospheric chemistry and transport induced by a winter warming event that was in progress at the time of the eclipse. At 60 km the ozone concentration, 2.7 × 10¹⁰ cm⁻³, was enhanced above that normally measured. This increase may also have been caused by the stratospheric warming event but the effects of a particle precipitation event, which was also in progress during the eclipse, may be important.     

Observations of auroral E-region plasma waves and electron heating with EISCAT and a VHF radar interferometer

Two radars were used simultaneously to study naturally occurring electron heating events in the auroral E-region ionosphere. During a joint campaign in March 1986 the Cornell University Portable Radar Interferometer (CUPRI) was positioned to look perpendicular to the magnetic field to observe unstable plasma waves over Tromsø, Norway, while EISCAT measured the ambient conditions in the unstable region. On two nights EISCAT detected intense but short lived (< 1 min) electron heating events during which the temperature suddenly increased by a factor of 2–4 at altitudes near 108 km and the electron densities were less than 7 × 10⁴ cm⁻³. On the second of these nights CUPRI was operating and detected strong plasma waves with very large phase velocities at precisely the altitudes and times at which the heating was observed. The altitudes, as well as one component of the irregularity drift velocity, were determined by interferometric techniques. From the observations and our analysis, we conclude that the electron temperature increases were caused by plasma wave heating and not by either Joule heating or particle precipitation.     

Analytical and numerical models of hydrothermal fluid flow at fault intersections

Fault intersections are the locus of hot spring activity and Carlin‐type gold mineralization within the Basin and Range, USA. Analytical and numerical solutions to Stokes equation suggest that peak fluid velocities at fault intersections increase between 20% and 47% when fracture apertures have identical widths but increase by only about 1% and 8% when aperture widths vary by a factor of 2. This suggests that fault zone intersections must have enlarged apertures. Three‐dimensional finite element models that consider intersecting 10‐ to 20‐m wide fault planes resulted in hot spring activity being preferentially located at fault zone intersections when fault zones were assigned identical permeabilities. We found that the onset of convection at the intersections of the fault zones occurred in our hydrothermal model over a narrow permeability range between 5 × 10^?13 and 7 × 10^?13 m². Relatively high vertical fluid velocities (0.3–3 m year^?1) extended away from the fault intersections for about 0.5–1.5 km. For the boundary conditions and fault plane dimensions used, peak discharge temperatures of 112°C at the water table occurred with an intermediate fault zone permeability of 5 × 10^?13 m². When fault plane permeability differed by a factor of 2 or more, the locus of hot spring activity shifted away from the intersections. However, increasing the permeability at the core of the fault plane intersection by 40% shifted the discharge back to the intersections. When aquifer units were assigned a permeability value equal to those of the fault planes, convective rolls developed that extend about 3 km laterally along the fault plane and into the adjacent aquifer.     

Concentrations of H2O and NO in the mesosphere and the lower thermosphere at high latitudes

Water vapour and nitric oxide concentrations in the mesosphere and lower thermosphere were derived from infrared emission and positive ion composition measurements above northern Europe during the Energy Budget Campaign 1980. The experiments were performed at different levels of geomagnetic disturbance. Both water vapour and nitric oxide are highly variable. Water vapour mixing ratios between 0.2 ppm and 10 ppm were observed. The nitric oxide peak densities varied by more than a factor of ten. Maximum values of 2 × 10⁹cm⁻³ were obtained.     

On the role of dust in the summer mesopause

We propose that dust formed at the cool summer mesopause may have optical properties very different from that measured for bulk material of ice. The smallness of the dust and possible surface impurities may lead to high photoelectric yields and low workfunctions. For such reasons the dust in the summer mesopause may, at least occasionally, be charged to substantial positive surface potentials while pure ice, with its high photoelectric workfunction, would be charged to low and negative potentials by collisions with plasma particles. The presence of ‘dressed’ dust particles, with surface potentials of some volts, can lead to enhanced radar backscatter. We also suggest that the apparent reductions in electron density (‘bite-out’), which have been observed in the radar backscatter region, can be caused by the inability of an electrostatic probe to deflect the massive dust particles.The dust density which is required by our model to explain radar backscatter and electron bite-outs is of the order of 10 cm⁻³ for dust of radius above 5 × 10⁻⁶ cm.     

Effective gas permeability of Tight Gas Sandstones as a function of capillary pressure – a non‐steady‐state approach

Single‐ and two‐phase (gas/water) fluid transport in tight sandstones has been studied in a series of permeability tests on core plugs of nine tight sandstones of the southern North Sea. Absolute (Klinkenberg‐corrected) gas permeability coefficients (k_gas__inf) ranged between 3.8 × 10^?16 and 6.2 × 10^?19 m² and decreased with increasing confining pressure (10–30 MPa) by a factor 3–5. Klinkenberg‐corrected (intrinsic) gas permeability coefficients were consistently higher by factors from 1.4 to 10 than permeability coefficients determined with water. Non‐steady‐state two‐phase (He/water) flow experiments conducted up to differential pressures of 10 MPa document the dynamically changing conductivity for the gas phase, which is primarily capillary‐controlled (drainage and imbibition). Effective gas permeability coefficients in the two‐phase flow tests ranged between 1.1 × 10^?17 and 2.5 × 10^?22 m², corresponding to relative gas permeabilities of 0.03% and 10%. In the early phase of the nonstationary flow regime (before establishment of steady‐state conditions), they may be substantially (>50%) lower. Effective gas permeability measurements are affected by the following factors: (i) Capillary‐controlled drainage/imbibition, (ii) viscous–dynamic effects (iii) and slip flow.     

Diurnal variations of equatorial electrojet parameters derived from POGO solstitial data

The POGO electrojet data have been analysed for the winter and summer solstitial seasons of the two years, 1968 and 1969, respectively. Our analysis yielded a very large number of values (about 432), each of the electrojet half width, w, its peak current intensity, J₀, and its total eastward current, I₊, at 0900–1400 LT in December, and at 1000–1500 LT in June solstitial seasons, respectively. The all-longitude daytime values of the parameters are 246 ± 48 km for w,216 ± 60 A km⁻¹ for J₀, and (58 ± 8) × 10³A for I₊, in December of 1968 and 218 ± 19 km for w, 187 ± 20 A km⁻¹ for J₀, and (45 ± 3) × 10³ A for I₊, in June of 1969, respectively. We therefore present a diurnal study covering the entire Earth, in which for the first time, morning data earlier than 1000 LT are incorporated in the analysis. This has enabled us to chart a clearer picture of the temporal variations of electrojet parameters at two different solstices. This shows that all of the three parameters vary substantially with local time, in such a manner that J₀ and I₊ attain maximal values around local noon, while w is a minimum then, and therefore confirms the finding of Agu and Onwumechili.     

Gas breakthrough experiments on fine-grained sedimentary rocks

The capillary sealing efficiency of fine‐grained sedimentary rocks has been investigated by gas breakthrough experiments on fully water saturated claystones and siltstones (Boom Clay from Belgium, Opalinus Clay from Switzerland and Tertiary mudstone from offshore Norway) of different lithological compositions. Sand contents of the samples were consistently below 12%, major clay minerals were illite and smectite. Porosities determined by mercury injection lay between 10 and 30% while specific surface areas determined by nitrogen adsorption (BET method) ranged from 20 to 48 m² g ^{? 1}. Total organic carbon contents were below 2%. Prior to the gas breakthrough experiments the absolute (single phase) permeability (k_abs) of the samples was determined by steady state flow tests with water or NaCl brine. The k_abs values ranged between 3 and 550 nDarcy (3 × 10^?21 and 5.5 × 10^?19 m²). The maximum effective permeability to the gas‐phase (k_eff) measured after gas breakthrough on initially water‐saturated samples extended from 0.01 nDarcy (1 × 10^?23 m²) up to 1100 nDarcy (1.1 × 10^?18 m²). The residual differential pressures after re‐imbibition of the water phase, referred to as the ‘minimum capillary displacement pressures’ (P_d), ranged from 0.06 to 6.7 MPa. During the re‐imbibition process the effective permeability to the gas phase decreases with decreasing differential pressure. The recorded permeability/pressure data were used to derive the pore size distribution (mostly between 8 and 60 nm) and the transport porosity of the conducting pore system (10^‐5–10^‐2%). Correlations could be established between (i) absolute permeability coefficients and the maximum effective permeability coefficients and (ii) effective or absolute permeability coefficients and capillary sealing efficiency. No correlation was found between the capillary displacement pressures determined from gas breakthrough experiments and those derived theoretically by mercury injection.     

Molecular transport of methane, ethane and nitrogen and the influence of diffusion on the chemical and isotopic composition of natural gas accumulations

Laboratory experiments have been performed to determine diffusion coefficients of natural gas components (methane, ethane and nitrogen) and isotope fractionation effects under simulated in situ pressure (up to 45 MPa effective stress) and temperature conditions (50–200°C) in water‐saturated pelitic and coarse‐grained rocks. Effective diffusion coefficients of molecular nitrogen (0.39 × 10^?11 to 21.6 × 10^?11 m² sec^?1 at 90°C) are higher than those for methane (0.18 × 10^?11 to 18.2 × 10^?11 m² sec^?1 at 90°C). Diffusive flux rates expressed in mass units are generally higher for N₂ than for CH₄. Both methane and (to a lesser extent) nitrogen diffusion coefficients decrease with increasing total organic carbon (TOC) content of the rock samples because of sorption processes on the organic matter. This effect decreases with increasing temperature. Effective diffusion coefficients increase upon a temperature increase from 50 to 200°C by a factor of four. Effective diffusion coefficients and steady‐state diffusive flux decrease with effective stress. Stationary diffusive fluxes drop by 50–70% for methane and 45–62% for nitrogen while effective diffusion coefficients are reduced by 38% (CH₄) and 32–48% (N₂), respectively. Isotope fractionation coefficients of diffusive transport are higher for methane (?1.56 and ?2.77‰) than for ethane (?0.84 and ?1.62‰). Application of the experimental results to geological systems show that diffusive transport has only a low transport efficiency. Significant depletion of natural gas reservoirs by molecular diffusion is only expected in cases of very poor caprock qualities (in terms of thickness and/or porosity) and over extended periods of geological time. Under these circumstances, the chemical and isotopic composition of a gas reservoir will change and maturity estimates based on these parameters may be deceptive. To account for these potential effects, nomograms have been developed to estimate diffusive losses and apply maturity corrections.     

Numerical model of pore‐pressure diffusion associated with the initiation of the 2010–2011 Guy–Greenbrier,Arkansas earthquakes

We model pore‐pressure diffusion caused by pressurized waste‐fluid injection at two nearby wells and then compare the buildup of pressure with the observed initiation and migration of earthquakes during the early part of the 2010–2011 Guy–Greenbrier earthquake swarm. Pore‐pressure diffusion is calculated using MODFLOW 2005 that allows the actual injection histories (volume/day) at the two wells to diffuse through a fractured and faulted 3D aquifer system representing the eastern Arkoma basin. The aquifer system is calibrated using the observed water‐level recovery following well shut‐in at three wells. We estimate that the hydraulic conductivities of the Boone Formation and Arbuckle Group are 2.2 × 10^?2 and 2.03 × 10^?3 m day^?1, respectively, with a hydraulic conductivity of 1.92 × 10^?2 m day^?1 in the Hunton Group when considering 1.72 × 10^?3 m day^?1 in the Chattanooga Shale. Based on the simulated pressure field, injection near the relatively conductive Enders and Guy–Greenbrier faults (that hydraulically connect the Arbuckle Group with the underlying basement) permits pressure diffusion into the crystalline basement, but the effective radius of influence is limited in depth by the vertical anisotropy of the hydraulic diffusivity. Comparing spatial/temporal changes in the simulated pore‐pressure field to the observed seismicity suggests that minimum pore‐pressure changes of approximately 0.009 and 0.035 MPa are sufficient to initiate seismic activity within the basement and sedimentary sections of the Guy–Greenbrier fault, respectively. Further, the migration of a second front of seismicity appears to follow the approximately 0.012 MPa and 0.055 MPa pore‐pressure fronts within the basement and sedimentary sections, respectively.