Quantication of magnetic coercivity components by the analysis of acquisition curves of isothermal remanent magnetisation

A new method of analysing isothermal remanent magnetisation (IRM) acquisition curves based on cumulative log Gaussian analysis [Robertson and France, Phys. Earth Planet. Inter. 82 (1994) 223-234] is proposed. It is based on the curve fitting of the IRM acquisition curve versus the logarithm of the applied field with: (i) the acquisition curve on a linear scale, (ii) the acquisition curve expressed as a gradient, and (iii) the acquisition curve on a probability scale. Even when a sample is not saturated, its magnetic properties can be defined, although with less certainty. The number of magnetic components required for an optimal fit to a measured IRM acquisition curve is evaluated statistically. The method discriminates on the basis of different mineral coercivity. Therefore, additional rock-magnetic tests are still required to separate minerals with similar coercivities

Derivation of a near-surface damping model for the Groningen gas field

Seismic damping of near-surface deposits is an important input to site-response analysis for seismic hazard assessment. In Groningen, the Netherlands, gas production from a reservoir at 3 km depth causes seismicity. Above the gas field, an 800 m thick layer of unconsolidated sediments exist, which consists of a mixture of sand, gravel, clay and peat strata. Shear waves induced at 3 km depth experience most of their anelastic attenuation in these loose sediments. A good estimate of damping is therefore crucial for modelling realistic ground-motion levels. In Groningen, we take advantage of a large network of 200 m deep vertical arrays to estimate damping from recordings of the induced events. As a first step, we apply seismic interferometry by deconvolution to estimate local transfer functions over these vertical arrays. Subsequently, two different methods are employed. The first is the 'upgoing' method, where the amplitude decay of the retrieved upgoing wave is used. The second is the 'up-down' method, where the amplitude difference between retrieved up- and downgoing waves is utilized. For the upgoing method, the amplitude of the upgoing direct wave is affected by both elastic and anelastic effects. In order to estimate the anelastic attenuation, it is necessary to remove the elastic amplification first. Despite the fact that elastic compensation could be determined quite accurately, non-physical damping values were estimated for a number of boreholes. Likely, the underlying cause was small differences in effective response functions of geophones at different depths. It was found that the up-down method is more robust. With this method, elastic propagation corrections are not needed. In addition, small differences in in situ geophone response are irrelevant because the up- and downgoing waves retrieved at the same geophone are used. For the 1-D case, we showed that for estimating the local transfer function, the complex reverberations need to be included in the interferometric process. Only when this is done, the transfer function does not contain elastic transmission loss and Q estimation can be made without knowing the soil profile in detail. Uncertainty in the estimated damping was found from the signal-to-noise ratio of the estimated transfer function. The Q profiles estimated with the up-down method were used to derive a damping model for the top 200 m of the entire Groningen field. A scaling relation was derived by comparing estimated Q profiles with low-strain damping profiles that were constructed using published models for low-strain damping linked to soil properties. This scaling relation, together with the soil-property-based damping model, allowed up-scaling of the model to each grid-cell in the Groningen field. For depths below 200 m, damping was derived from the attenuation of the microseism over Groningen. The mean damping model, over a frequency band between 2 and 20 Hz, was estimated to be 2.0 per cent (0-50 m depth), 1.3 per cent (50-100 m), 0.66 per cent (100-150 m), 0.57 per cent (150-200 m) and 0.5 per cent (200-580 m)

Liquefaction hazard of the Groningen region of the Netherlands due to induced seismicity

The operator of the Groningen gas field is leading an effort to quantify the seismic hazard and riskof the region due to induced earthquakes, includingoverseeing one of the most comprehensive liquefaction hazard studies performedgloballyto date. Due tothe unique characteristics of the seismic hazard and the geologic deposits in Groningen, efforts first focused on developing relationships for a Groningen-specific liquefaction triggering model. The liquefaction hazard was then assessedusing a Monte Carlo method, wherein a range of credibleevent scenarios were considered in computingliquefaction damage-potentialhazard curves. Thiseffort entailed the use of a regional stochastic seismic source model,ground motion prediction equation,site response model,and geologic model that were developed as part of the broader regional seismic hazardassessment.“No-to-Minor Surficial Liquefaction Manifestations”arepredicted for mostsites across the study areafor a 75-year return period. The only sites where “Moderate Surficial Liquefaction Manifestations” are predicted are in the town of Zandeweer, with only some of the sites in the townbeing predicted to experience this severityof liquefactionfor thisreturn period

A database of ground-motion recordings, site profiles, and amplification factors from the Groningen gas field in the Netherlands

A comprehensive database that has been used to develop ground motion models for induced earthquakes in the Groningen gas field is provided in a freely accessible online repository. The database includes more than 8500 processed ground motion recordings from 87 earthquakes of local magnitude ML between 1.8 and 3.6, obtained from a large network of surface accelerographs and borehole geophones placed at 50 m depth intervals to a depth of 200 m. The 5%-damped pseudo-acceleration spectra and Fourier amplitude spectra of the records are also provided. Measured shear-wave velocity (VS) profiles, obtained primarily from seismic Cone Penetration Tests (CPTs), are provided for 80 of the ∼100 recording stations. A model representing the regional dynamic soil properties is presented for the entire gas field plus a 5 km onshore buffer zone, specifying lithology, VS, and damping for all layers above the reference baserock horizon located at about 800 m depth. Transfer functions and frequency-dependent amplification factors from the reference rock horizon to the surface for the locations of the recording stations are also included. The database provides a valuable resource for further refinement of induced seismic hazard and risk modeling in Groningen as well as for generic research in site response of thick, soft soil deposits and the characteristics of ground motions from small-magnitude, shallow-focus induced earthquakes

Characterisation of ground motion recording stations in the Groningen gas field

The seismic hazard and risk analysis for the onshore Groningen gas field requires information about local soil properties, in particular shear-wave velocity (VS). A fieldwork campaign was conducted at 18 surface accelerograph stations of the monitoring network. The subsurface in the region consists of unconsolidated sediments and is heterogeneous in composition and properties. A range of different methods was applied to acquire in situ VS values to a target depth of at least 30 m. The techniques include seismic cone penetration tests (SCPT) with varying source offsets, multichannel analysis of surface waves (MASW) on Rayleigh waves with different processing approaches, microtremor array, cross-hole tomography and suspension P-S logging. The offset SCPT, cross-hole tomography and common midpoint cross-correlation (CMPcc) processing of MASW data all revealed lateral variations on length scales of several to tens of metres in this geological setting. SCPTs resulted in very detailed VS profiles with depth, but represent point measurements in a heterogeneous environment. The MASW results represent VS information on a larger spatial scale and smooth some of the heterogeneity encountered at the sites. The combination of MASW and SCPT proved to be a powerful and cost-effective approach in determining representative VS profiles at the accelerograph station sites. The measured VS profiles correspond well with the modelled profiles and they significantly enhance the ground motion model derivation. The similarity between the theoretical transfer function from the VS profile and the observed amplification from vertical array stations is also excellent

Developing a model for the prediction of ground motions due to earthquakes in the Groningen gas field

AbstractMajor efforts are being undertaken to quantify seismic hazard and risk due to production-induced earthquakes in the Groningen gas field as the basis for rational decision-making about mitigation measures. An essential element is a model to estimate surface ground motions expected at any location for each earthquake originating within the gas reservoir. Taking advantage of the excellent geological and geophysical characterisation of the field and a growing database of ground-motion recordings, models have been developed for predicting response spectral accelerations, peak ground velocity and ground-motion durations for a wide range of magnitudes. The models reflect the unique source and travel path characteristics of the Groningen earthquakes, and account for the inevitable uncertainty in extrapolating from the small observed magnitudes to potential larger events. The predictions of ground-motion amplitudes include the effects of nonlinear site response of the relatively soft near-surface deposits throughout the field.</jats:p

High‐resolution rock magnetic cyclostratigraphy in an Eocene flysch, Spanish Pyrenees

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95346/1/ggge1746.pd

Ground-motion prediction models for induced earthquakes in the Groningen gas field, the Netherlands

Small-magnitude earthquakes induced by gas production in the Groningen field in the Netherlands have prompted the development of seismic risk models that serve both to estimate the impact of these events and to explore the efficacy of different risk mitigation strategies. A core element of the risk modelling is ground-motion prediction models (GMPM) derived from an extensive database of recordings obtained from a dense network of accelerographs installed in the field. For the verification of damage claims, an empirical GMPM for peak ground velocity (PGV) has been developed, which predicts horizontal PGV as a function of local magnitude, ML; hypocentral distance, Rhyp; and the time-averaged shear-wave velocity over the upper 30 m, VS30. For modelling the risk due to potential induced and triggered earthquakes of larger magnitude, a GMPM for response spectral accelerations has been developed from regressions on the outputs from finite-rupture simulations of motions at a deeply buried rock horizon. The GMPM for rock motions is coupled with a zonation map defining frequency-dependent non-linear amplification factors to obtain estimates of surface motions in the region of thick deposits of soft soils. The GMPM for spectral accelerations is formulated within a logic-tree framework to capture the epistemic uncertainty associated with extrapolation from recordings of events of ML ≤ 3.6 to much larger magnitudes