Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

A model for the Am (Km) planetary geomagnetic activity index and application to prediction

We studied the existence of dynamical stochastic relations in the evolution of the am index. A first analysis of the autocorrelation functions showed evidence of several seasonalities. We first used linear (ARMA) models, and it was found that these do not account for the whole internal dynamics of the data series. We then used various non-linear models to provide a better fit to reality. The forecast performances of the non-linear models are not significantly different from those of the linear model. We give a tentative explanation for the failure of the non-linear predictions. Finally, ARCH models were used in order to take into account the fact that the confidence interval for the predicted value depends on past observations.     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.     

Geothermal Characterization of the St. Lawrence Lowlands Sedimentary Basin,Québec,Canada

Defining temperature at depth to identify geothermal resources relies on the evaluation of the Earth heat flow based on equilibrium temperature measurements as well as thermal conductivity and heat generation rate assessment. Such high-quality geothermal data can be sparse over the region of interest. This is the case of the St. Lawrence Lowlands sedimentary basin covering 20,000 km² to the south of Québec, Canada, and enclosing only three wells up to a depth of 500 m with equilibrium heat flow measurements. However, more than 250 oil and gas exploration wells have been drilled in this area, providing for this study (parce que c'est 93 sinon) 81 locations with bottom-hole temperature up to a depth of 4300 m, however, not at equilibrium. Analyzing these data with respect to the deep geothermal resource potential of this sedimentary basin requires evaluating the thermal conductivity and heat generation rate of its geological units to properly extrapolate temperature downward. This was done by compiling literature and recent thermal conductivity measurements in outcrop and core samples as well as new heat generation rate estimates from spectral gamma ray logs to establish a first thermal assessment of geological units deep down into the basin. The mean thermal conductivity of the thermal units varies from 2.5 to 6.3 W/m·K, with peak values in the basal sandstones, while the heat generation rate varies from 1.6 to 0.3 µW/m³, decreasing from the upper caprocks toward the base of the sequence. After correcting the bottom-hole temperatures for drilling disturbance with the Harrison correction and subsequently for paleoclimate variations, results indicate a mean geothermal gradient of 23.1 °C/km, varying from 14 to 40 °C/km. Evaluating the basin thermal state from oil and gas data is a significant challenge facilitated by an understanding of its thermal properties.     

Nouvelles données pour l’histoire des fluctuations holocènes du niveau du lac de Neuchatel (Suisse): la séquence sédimentaire de Marin-Les Piécettes

Resume. Les fouilles archéologiques réalisées de 1998 à 2002 sur le site néolithique de Marin-Les Piécettes sur la rive nord-ouest du lac de Neuchatel ont offert l’opportunité d’observer deux séquences sédimentaires intéressantes pour la reconstitution des variations holocènes du niveau du lac, en particulier au cours de la période correspondant au changement de cours de l’Aar en direction du nord-est, vers 5650–5500 cal BP. La géométrie des corps sédimentaires et les changements lithologiques permettent de reconna?tre plusieurs phases majeures de bas niveaux relatifs du plan d’eau vers 6870, 5660, 5060–3990, 3580, 2280, 1640 et 1170 cal BP. Le site néolithique de Marin-Les Piécettes s’est développé vers 5454–5433 cal BP à la faveur d’un court épisode de bas niveau du lac qui a interrompu une importante phase de haut niveau au début du Subboréal. Les événements paléohydrologiques identifiés à Marin-Les Piécettes apparaissent cohérents avec ceux reconstitués sur le site de Montilier au bord du lac de Morat, ainsi qu’avec les variations paléohydrologiques reconstituées à partir des lacs du Jura, du Plateau suisse et des Préalpes du Nord fran?aises. Ceci suggère qu’ils témoignent d’un signal climatique. La séquence sédimentaire de Marin-Les Piécettes met enfin en évidence la formation d’un cordon littoral sur la rive nord-ouest du lac juste après 5660 cal BP. Manuscrit re?u le 22 septembre 2004 Révision acceptée le 11 mars 2005     

The extreme diversity of uranium deposits

Most available classifications of uranium deposits are based on the characteristics of the host rocks or on the morphology of the ore deposits. The aim of the present paper is to propose the basis for a genetic classification of these deposits. After a short introduction on the geochemical behavior of uranium in fluids and silicate melts and on the main uranium fractionation mechanisms operating in uranium-rich peraluminous, metaluminous, and peralkaline melts, the most recent metallogenic models of the main types of uranium deposits are shortly reviewed.     

Petrological and experimental evidence for differentiation of water-rich magmas beneath St. Kitts,Lesser Antilles

St. Kitts lies in the northern Lesser Antilles, a subduction-related intraoceanic volcanic arc known for its magmatic diversity and unusually abundant cognate xenoliths. We combine the geochemistry of xenoliths, melt inclusions and lavas with high pressure–temperature experiments to explore magma differentiation processes beneath St. Kitts. Lavas range from basalt to rhyolite, with predominant andesites and basaltic andesites. Xenoliths, dominated by calcic plagioclase and amphibole, typically in reaction relationship with pyroxenes and olivine, can be divided into plutonic and cumulate varieties based on mineral textures and compositions. Cumulate varieties, formed primarily by the accumulation of liquidus phases, comprise ensembles that represent instantaneous solid compositions from one or more magma batches; plutonic varieties have mineralogy and textures consistent with protracted solidification of magmatic mush. Mineral chemistry in lavas and xenoliths is subtly different. For example, plagioclase with unusually high anorthite content (An_≤100) occurs in some plutonic xenoliths, whereas the most calcic plagioclase in cumulate xenoliths and lavas are An₉₇ and An₉₅, respectively. Fluid-saturated, equilibrium crystallisation experiments were performed on a St. Kitts basaltic andesite, with three different fluid compositions (XH₂O = 1.0, 0.66 and 0.33) at 2.4 kbar, 950–1025 °C, and fO₂ = NNO ? 0.6 to NNO + 1.2 log units. Experiments reproduce lava liquid lines of descent and many xenolith assemblages, but fail to match xenolith and lava phenocryst mineral compositions, notably the very An-rich plagioclase. The strong positive correlation between experimentally determined plagioclase-melt Kd_Ca–Na and dissolved H₂O in the melt, together with the occurrence of Al-rich mafic lavas, suggests that parental magmas were water-rich (> 9 wt% H₂O) basaltic andesites that crystallised over a wide pressure range (1.5–6 kbar). Comparison of experimental and natural (lava, xenolith) mafic mineral composition reveals that whereas olivine in lavas is predominantly primocrysts precipitated at low-pressure, pyroxenes and spinel are predominantly xenocrysts formed by disaggregation of plutonic mushes. Overall, St. Kitts xenoliths and lavas testify to mid-crustal differentiation of low-MgO basalt and basaltic andesite magmas within a trans-crustal, magmatic mush system. Lower crustal ultramafic cumulates that relate parental low-MgO basalts to primary, mantle -derived melts are absent on St. Kitts.     

Embryo dune development drivers: beach morphology,growing season precipitation,and storms

For development of embryo dunes on the highly dynamic land–sea boundary, summer growth and the absence of winter erosion are essential. Other than that, however, we know little about the specific conditions that favour embryo dune development. This study explores the boundary conditions for early dune development to enable better predictions of natural dune expansion. Using a 30 year time series of aerial photographs of 33 sites along the Dutch coast, we assessed the influence of beach morphology (beach width and tidal range), meteorological conditions (storm characteristics, wind speed, growing season precipitation, and temperature), and sand nourishment on early dune development. We examined the presence and area of embryo dunes in relation to beach width and tidal range, and compared changes in embryo dune area to meteorological conditions and whether sand nourishment had been applied. We found that the presence and area of embryo dunes increased with increasing beach width. Over time, embryo dune area was negatively correlated with storm intensity and frequency. Embryo dune area was positively correlated with precipitation in the growing season and sand nourishment. Embryo dune area increased in periods of low storm frequency and in wet summers, and decreased in periods of high storm frequency or intensity. We conclude that beach morphology is highly influential in determining the potential for new dune development, and wide beaches enable development of larger embryo dune fields. Sand nourishment stimulates dune development by increasing beach width. Finally, weather conditions and non‐interrupted sequences of years without high‐intensity storms determine whether progressive dune development will take place. Copyright © 2017 John Wiley & Sons, Ltd.