How does water near clay mineral surfaces influence the rock physics of shales?

Clays and clay‐bearing rocks like shale are extremely water sensitive. This is partly due to the interaction between water and mineral surfaces, strengthened by the presence of nanometer‐size pores and related large specific surface areas. Molecular‐scale numerical simulations, using a discrete‐element model, show that shear rigidity can be associated with structurally ordered (bound or adsorbed) water near charged surfaces. Building on these and other molecular dynamics simulations plus nanoscale experiments from the literature, the water monolayer adjacent to hydrophilic solid surfaces appears to be characterised by shear stiffness and/or enhanced viscosity. In both cases, elastic wave propagation will be affected by the bound or adsorbed water. Using a simple rock physics model, bound water properties were adjusted to match laboratory measured P‐ and S‐wave velocities on pure water‐saturated kaolinite and smectite. To fit the measured stress sensitivity, particularly for kaolinite, the contribution from solid‐grain contact stiffness needs to be added. The model predicts, particularly for S‐waves, that viscoelastic bound water could be a source of dispersion in clay and clay‐rich rocks. The bound‐water‐based rock physics model is found to represent a lower bound to laboratory‐measured velocities obtained with shales of different mineralogy and porosity levels.     

T-matrix approach to seismic forward modelling in the acoustic approximation

Forward seismic modelling in the acoustic approximation, for variable velocity but constant density, is dealt with. The wave equation and the boundary conditions are represented by a volume integral equation of the Lippmann-Schwinger (LS) or Fredholm type. A T-matrix (or transition operator) approach from quantum mechanical potential scattering theory is used to derive a family of linear and nonlinear approximations (cluster expansions), as well as an exact numerical solution of the LS equation. For models of 4D anomalies involving small or moderate contrasts, the Born approximation gives identical numerical results as the first-order t-matrix approximation, but the predictions of an exact T-matrix solution can be quite different (depending on spatial extention of the perturbations). For models of fluid-saturated cavities involving large or huge contrasts, the first-order t-matrix approximation is much more accurate than the Born approximation, although it does not lead to significantly more time-consuming computations. If the spatial extention of the perturbations is not too large, it is practical to use the exact T-matrix solution which allows for arbitrary contrasts and includes all the effects of multiple scattering.     

A Compilation of Radiocarbon Dates from Disko Bugt,Central West Greenland/Meteorological Observations in 1996 at the Arctic Station,Qeqertarsuaq (Godhavn), Central West Greenland/A Discussion on Pingos in Mellemfjord,Disko, Central West Greenland/Open System Pingos in Mellemfjord,Disko, Central West Greenland: A Reply to Gurney and Worsley

188 new and previously published radiocarbon dates on Holocene material from Disko Bugt, central West Greenland, are presented together with relevant informations (laboratory number, place name, dated material, geografic coordinates, altitude and δ¹³C-value) about the dates.

In October 1990 an automatic weather station was established at the Arctic station (65 °15′N,53 °31′W), Qeqertarsuaq (Godhavn), Central West Greenland. The Station register parameters each 20 min., and the parameters have been described in an earlier paper in this journal by Nielsen et al. (1995). The present paper summarises main points of the climate during 1996.

Open system, or more correctly, hydraulic pingos, are genetically poorly understood. A continuing problem concerns their need for a perennial groundwater supply (intra- or sub-permafrost). This has to be maintained despite the existence of continuous permafrost in many areas where they are located. Recent work on Disko Island has suggested a new type of hydraulic pingo developing only in a “marsh environment”. It is argued that the marsh setting is not relevant to the formation of these features and that they are simply hydraulic pingos.

Abstract

A group of marsh initiated open system pingo remnants from the Iterdlagssûp kûgssua valley mouth, in Mellemfjord, Disko Island, Central West Greenland was described in Christiansen (1995). Gurney and Worsley (1997) state that the location of this group of pingo remnants in the Iterdlagssûp kûgssua valley mouth is of no relevance to their genesis, and that they presumably were the result of an assumed late Holocene sea level regression, causing permafrost to be established in the valley bottom. In this reply the arguments by Gurney and Worsley (1997) on the Iterdlagssûp kûgssua valley pingos are commented, and it is argued by way of sea level information, frost penetration and water supply that the special setting must indeed have caused pingo initiation and growth. Furthermore, the area has experienced a relative sea level rise during the late Holocene.     

Assessment of Erythemal UV Level in Nepal Based on Solar UV Estimates from Total Ozone Mapping Spectrometer

Nepal lies on the southern slope of Himalaya in Asia. In a width ranging between 150 and 250 km, the altitude varies greatly from about 100 m at its southern border to a maximum of 8848 min the northern part. Like the variation in altitude, climatic condition varies quite a lot. Long-term monthly mean erythemal UV daily dose values for Nepal are evaluated using Total Ozone Mapping Spectrometer (TOMS) estimation from the time of its overpass between 1996 and 2003. The results are presented as summer and winter maps of mean UV levels in each satellite grid. The mean winter erythemal UV daily dose ranges between 2.1 and 3.6 kJ m-2 whereas summer values are found to lie between 4.6 and 9.7 kJ m-2. The altitude variation increases the UV levels by about 0.2 kJ km-1 in winter months, and 0.9 kJ km-1 in summer. A multiyear monthly average erythemal daily dose in most of the areas shows that the summer value is about three times higher than that in winter. Although year-to-year variation is not pronounced in high- and mid-elevation regions, UV levels seemed to decrease from 1997 to 2002 in the southern part of the country in the low elevation region by about 5.35%. Due to the combined effects of the altitude, low ozone concentration in the troposphere, and thin air, surface UV radiation at higher altitudes is found to be higher than in the surrounding regions.     

Seismic characterization of reservoirs with multiple fracture sets using velocity and attenuation anisotropy data

Knowledge about the spatial distribution of the fracture density and the azimuthal fracture orientation can greatly help in optimizing production from fractured reservoirs. Frequency-dependent seismic velocity and attenuation anisotropy data contain information about the fractures present in the reservoir. In this study, we use the measurements of velocity and attenuation anisotropy data corresponding to different seismic frequencies and azimuths to infer information about the multiple fracture sets present in the reservoir. We consider a reservoir model with two sets of vertical fractures characterized by unknown azimuthal fracture orientations and fracture densities. Frequency-dependent seismic velocity and attenuation anisotropy data is computed using the effective viscoelastic stiffness tensor and solving the Christoffel equation. A Bayesian inversion method is then applied to measurements of velocity and attenuation anisotropy data corresponding to different seismic frequencies and azimuth to estimate the azimuthal fracture orientations and the fracture densities, as well as their uncertainties. Our numerical examples suggest that velocity anisotropy data alone cannot recover the unknown fracture parameters. However, an improved estimation of the unknown fracture parameters can be obtained by joint inversion of velocity and attenuation anisotropy data.     

The Zn abundance and isotopic composition of diatom frustules,a proxy for Zn availability in ocean surface seawater

We have developed cleaning methods for extracting diatomopal from bulk marine sediment samples, for measurement of both zinc (Zn) abundance and isotope composition. This cleaning technique was then applied to a set of Holocene core-top samples from the Southern Ocean. The measured δ⁶⁶Zn (reported relative to the JMC_Lyon standard) and Zn/Si ratios from the Southern Ocean diatom_opal samples range from 0.7 to 1.5‰, and from 14 to 0.9 μmol/mol, respectively. The Zn abundance and isotope composition data show a clear correlation with opal burial rates and other oceanographic parameters. In common with previous work, we interpret the systematic changes in the Zn/Si ratio to be linked to the variability in the concentrations of bioavailable Zn in the ambient surface seawater where the diatom opal is formed. This variability is likely to be primarily controlled by the degree to which Zn is taken up into phytoplankton biomass. The observed systematic pattern in the δ⁶⁶Zn compositions of the diatom_opal core-top samples is, similarly, likely to reflect changes in the δ⁶⁶Zn composition of the ambient Zn in the surface waters above the core-top sites, which is progressively driven towards isotopically heavier values by preferential incorporation of the lighter isotopes into phytoplankton organic material. Thus, the systematic relationship between Zn isotopes and abundance observed in the core-top diatom_opal samples suggests a potential tool for investigating the biogeochemical cycling of Zn in the past surface ocean for down-core diatom_opal material. In this respect, it may be possible to test hypotheses that attribute variations in atmospheric CO₂ on glacial–interglacial timescales to the degree to which trace metals limited primary productivity in HNLC zones.     

Sedimentary recycling in arc magmas: geochemical and U–Pb–Hf–O constraints on the Mesoproterozoic Suldal Arc, SW Norway

The Hardangervidda-Rogaland Block within southwest Norway is host to ~1.52 to 1.48 Ga continental building and variable reworking during the ~1.1 to 0.9 Ga Sveconorwegian orogeny. Due to the lack of geochronological and geochemical data, the timing and tectonic setting of early Mesoproterozoic magmatism has long been ambiguous. This paper presents zircon U–Pb–Hf–O isotope data combined with whole-rock geochemistry to address the age and petrogenesis of basement units within the Suldal region, located in the centre of the Hardangervidda-Rogaland Block. The basement comprises variably deformed grey gneisses and granitoids that petrologically and geochemically resemble mature volcanic arc lithologies. U–Pb ages confirm that magmatism occurred from ~1,521 to 1,485 Ma, and conspicuously lack any xenocrystic inheritance of distinctly older crust. Hafnium isotope data range from εHf_(initial) +1 to +11, suggesting a rather juvenile magmatic source, but with possible involvement of late Palaeoproterozoic crust. Oxygen isotope data range from mantle-like (δ¹⁸O ~5 ‰) to elevated (~10 ‰) suggesting involvement of low-temperature altered material (e.g., supracrustal rocks) in the magma source. The Hf–O isotope array is compatible with mixing between mantle-derived material with young low-temperature altered material (oceanic crust/sediments) and older low-temperature altered material (continent-derived sediments). This, combined with a lack of xenoliths and xenocrysts, exposed older crust, AFC trends and S-type geochemistry, all point to mixing within a deep-crustal magma-generation zone. A proposed model comprises accretion of altered oceanic crust and the overlying sediments to a pre-existing continental margin, underthrusting to the magma-generation zone and remobilisation during arc magmatism. The geodynamic setting for this arc magmatism is comparable with that seen in the Phanerozoic (e.g., the Sierra Nevada and Coast Range batholiths), with compositions in the Suldal Sector reaching those of average upper continental crust. As within these younger examples, factors that drive magmatism towards the composition of the average continental crust include the addition of sedimentary material to magma source regions, and delamination of cumulate material. Underthrusting of sedimentary materials and their subsequent involvement in arc magmatism is perhaps a more widespread mechanism involved in continental growth than is currently recognised. Finally, the Suldal Arc magmatism represents a significant juvenile crustal addition to SW Fennoscandia.     

Reassembling the Paleogene–Eocene North Atlantic igneous province: New paleomagnetic constraints from the Isle of Mull, Scotland

The paleomagnetic data sets from the British Tertiary Igneous Province (BTIP) have recently been criticized as being unreliable and discordant with data from elsewhere in the North Atlantic Igneous Province (NAIP) [Riisager et al. Earth Planet. Sci. Lett. 201 (2002) 261–276; Riisager et al. Earth Planet. Sci. Lett. 214 (2003) 409–425]. We offer new paleomagnetic data for the extensive lava flow sequence on the Isle of Mull, Scotland, and can confirm the paleomagnetic pole positions emanating from important earlier studies. Our new north paleomagnetic pole position for Eurasia at 59 ± 0.2 Ma has latitude 73.3°N, longitude 166.2°E (dp/dm = 5.2/7.0).A re-evaluation and an inter-comparison of the paleomagnetic database emanating from the NAIP were carried out to test for sub-province consistency. We find a general agreement between the Eurasian part of NAIP (BTIP and Faeroes) and East Greenland data. However a compilation of West Greenland data displays a large and unexplained dispersion. We speculate on if this is related to different sense of block rotation of the Tertiary West Greenland constituents. Combining all data from the NAIP constituents, give a pole position at 75.0°N, 169.9°E (N = 25, K = 84.3, A₉₅ = 3.2) in Eurasian reference frame.