Structure and evolution of the volcanic rift zone at Ponta de São Lourenço,eastern Madeira

Ponta de São Lourenço is the deeply eroded eastern end of Madeira’s east–west trending rift zone, located near the geometric intersection of the Madeira rift axis with that of the Desertas Islands to the southeast. It dominantly consists of basaltic pyroclastic deposits from Strombolian and phreatomagmatic eruptions, lava flows, and a dike swarm. Main differences compared to highly productive rift zones such as in Hawai’i are a lower dike intensity (50–60 dikes/km) and the lack of a shallow magma reservoir or summit caldera. ⁴⁰Ar/³⁹Ar age determinations show that volcanic activity at Ponta de São Lourenço lasted from >5.2 to 4 Ma (early Madeira rift phase) and from 2.4 to 0.9 Ma (late Madeira rift phase), with a hiatus dividing the stratigraphy into lower and upper units. Toward the east, the distribution of eruptive centers becomes diffuse, and the rift axis bends to parallel the Desertas ridge. The bending may have resulted from mutual gravitational influence of the Madeira and Desertas volcanic edifices. We propose that Ponta de São Lourenço represents a type example for the interior of a fading rift arm on oceanic volcanoes, with modern analogues being the terminations of the rift zones at La Palma and El Hierro (Canary Islands). There is no evidence for Ponta de São Lourenço representing a former central volcano that interconnected and fed the Madeira and Desertas rifts. Our results suggest a subdivision of volcanic rift zones into (1) a highly productive endmember characterized by a central volcano with a shallow magma chamber feeding one or more rift arms, and (2) a less productive endmember characterized by rifts fed from deep-seated magma reservoirs rather than from a central volcano, as is the case for Ponta de São Lourenço.     

Effects of pore fluid pressure on the seismic response of a fractured carbonate reservoir

An effective medium model for the stress-dependent seismic properties of fractured reservoirs is developed here on the basis of a combination of a general theory of viscoelastic waves in rock-like composites with recently published formulae for deformation of communicating and interacting cavities (interconnected pores/cracks and fractures at finite concentration) under drained loading. The inclusion-based model operates with spheroidal cavities at two different length scales; namely, the microscopic scale of the pores and (grain-boundary) cracks, and the mesoscopic scale of the fractures (controlling the flow of fluid). The different cavity types can in principle have any orientation and aspect ratio, but the microscopic pores/cracks and mesoscopic fractures were here assumed to be randomly and vertically oriented, respectively. By using three different aspect ratios for the relatively round pores (representing the stiff part of the pore space) and a distribution of aspect ratios for the relatively flat cracks (representing the compliant part of the pore space), we obtained a good fit between theoretical predictions and ultrasonic laboratory measurements on an unfractured rock sample under dry conditions. By using a single aspect ratio for the mesoscopic fractures, we arrived at a higher-order microstructural model of fractured porous media which represents a generalization of the first-order model developed by Chapman et al. (2002,2003). The effect of cavity size was here modelled under the assumption that the characteristic time for wave-induced (squirt) flow at the scale of a particular cavity (pore/crack vs. fracture) is proportional with the relevant scale-size. In the modelling, we investigate the effect of a decreasing pore pressure with constant confining pressure (fixed depth), and hence, increasing effective pressure. The analysis shows that the attenuation-peak due to the mesoscopic fractures in the reservoir will move downward in frequency as the effective pressure increases. In the range of seismic frequencies, our modelling indicates that the P-wave velocities may change by more than 20% perpendicular to the fractures and close to 10% parallel to the fractures. In comparison, the vertical S-wave velocities change by only about 5% for both polarization directions (perpendicular and parallel to the fractures) when the effective pressure increases from 0 to 15 MPa. This change is mainly due to the overall change in porosity with pressure. The weak pressure dependence is a consequence of the fact that the S waves will only sense if the fractures are open or not, and since all the fractures have the same aspect ratio, they will close at the same effective pressure (which is outside the analysed interval). Approximate reflection coefficients were computed for a model consisting of the fractured reservoir embedded as a layer in an isotropic shale and analysed with respect to variations in Amplitude Versus Offset and aZimuth (AVOZ) properties at seismic frequencies for increasing effective pressure. For the P-P reflections at the top of the reservoir, it is found that there is a significant dependence on effective pressure, but that the variations with azimuth and offset are small. The lack of azimuthal dependence may be explained from the approximate reflection coefficient formula as a result of cancellation of terms related to the S-wave velocity and the Thomson’s anisotropy parameter δ. For the P-S reflection, the azimuthal dependence is larger, but the pressure dependence is weaker (due to a single aspect ratio for the fractures). Finally, using the effective stiffness tensor for the fractured reservoir model with a visco-elastic finite-difference code, synthetic seismograms and hodograms were computed. From the seismograms, attenuation changes in the P wave reflected at the bottom of the reservoir can be observed as the effective pressure increases. S waves are not much affected by the fractures with respect to attenuation, but azimuthal dependence is stronger than for P waves, and S-wave splitting in the bottom reservoir P-S reflection is clearly seen both in the seismograms and hodograms. From the hodograms, some variation in the P-S reflection with effective pressure can also be observed.     

The reactivity of iron oxides towards reductive dissolution with ascorbic acid in a shallow sandy aquifer (Rømø, Denmark)

The pool of iron oxides, available in sediments for reductive dissolution, is usually estimated by wet chemical extraction methods. Such methods are basically empirically defined and calibrated against various synthetic iron oxides. However, in natural sediments, iron oxides are present as part of a complex mixture of iron oxides with variable crystallinity, clays and organics etc. Such a mixture is more accurately described by a reactive continuum covering a range from highly reactive iron oxides to non-reactive iron oxide. The reactivity of the pool of iron oxides in sediment can be determined by reductive dissolution in 10 mM ascorbic acid at pH 3. Parallel dissolution experiments in HCl at pH 3 reveal the release of Fe(II) by proton assisted dissolution. The difference in Fe(II)-release between the two experiments is attributed to reductive dissolution of iron oxides and can be quantified using the rate equation J/m₀ = k′(m/m₀)^γ, where J is the overall rate of dissolution (mol s⁻¹), m₀ the initial amount of iron oxide, k′ a rate constant (s⁻¹), m/m₀ the proportion of undissolved mineral and γ a parameter describing the change in reaction rate over time. In the Rømø aquifer, Denmark, the reduction of iron oxides is an important electron accepting process for organic matter degradation and is reflected by the steep increase in aqueous Fe²⁺ over depth. Sediment from the Rømø aquifer was used for reductive dissolution experiments with ascorbic acid. The rate parameters describing the reactivity of iron oxides in the sediment are in the range k′ = 7·10⁻⁶ to 1·10⁻³ s⁻¹ and γ = 1 to 2.4. These values are intermediate between a synthetic 2-line ferrihydrite and a goethite. The rate constant increases by two orders of magnitude over depth suggesting an increase in iron oxide reactivity with depth. This increase was not captured by traditional oxalate and dithionite extractions.     

Geochemical processes and solute transport at the seawater/freshwater interface of a sandy aquifer

Geochemical processes occurring at a seawater/freshwater interface were studied in a shallow coastal siliclastic aquifer containing minor amounts of calcite. Data were collected from 106 piezometers in a 120-m transect from the coastline and landward. In the first 40 m from the coastline, a wedge of saltwater is intruding below the freshwater aquifer. The aquifer is strongly reduced with mineralization of organic matter by methanogenesis in the freshwater aquifer, and sulfate reduction dominating in the most seaward part of the saline aquifer. The spatial separation of cations in the aquifer indicated a slow freshening process where Ca²⁺ from freshwater displaced the marine cations Na⁺ and Mg²⁺ from the exchanger complex. The resulting loss of Ca²⁺ from solution decreases the saturation state for calcite and possibly causes calcite dissolution. A storm-flooding event was recorded where pulses of dense seawater sank through the fresh aquifer. As a result, the terminal electron accepting process switched from methanogenesis to sulfate reduction. The pulses of sinking seawater also triggered cation exchange reactions where Ca²⁺ was expelled from the exchanger by seawater Na⁺ and Mg²⁺. The released Ca²⁺ is being flushed from the aquifer by groundwater flow, and this export of Ca²⁺ will, in the long term, cause decalcification of the sediment. The water composition in the aquifer is in a transient state as the result of various processes that operate on different timescales. Oxidation of organic matter occurs continuously but at a rate decreasing on a geological time scale. The freshening of the aquifer operates on the timescale of a few years. The episodic flooding and sinking of seawater through the aquifer proceeds in the course of days to weeks, but occurs irregularly with years in between.     

Lake metabolism scales with lake morphometry and catchment conditions

We used a comparative data set for 25 lakes in Denmark sampled during summer to explore the influence of lake morphometry, catchment conditions, light availability and nutrient input on lake metabolism. We found that (1) gross primary production (GPP) and community respiration (R) decline with lake area, water depth and drainage ratio, and increase with algal biomass (Chl), dissolved organic carbon (DOC) and total phosphorus (TP); (2) all lakes, especially small with less incident light, and forest lakes with high DOC, have negative net ecosystem production (NEP < 0); (3) daily variability of GPP decreases with lake area and water depth as a consequence of lower input of nutrients and organic matter per unit water volume; (4) the influence of benthic processes on free water metabolic measures declines with increasing lake size; and (5) with increasing lake size, lake metabolism decreases significantly per unit water volume, while depth integrated areal rates remain more constant due to a combination of increased light and nutrient limitation. Overall, these meta-parameters have as many significant but usually weaker relationships to whole-lake and benthic metabolism as have TP, Chl and DOC that are directly linked to photosynthesis and respiration. Combining water depth and Chl to predict GPP, and water depth and DOC to predict R, lead to stronger multiple regression models accounting for 57–63% of the variability of metabolism among the 25 lakes. It is therefore important to consider differences in lake morphometry and catchment conditions when comparing metabolic responses of lakes to human impacts.     

六元体系Li-Na-K-Mg-Cl-SO4-H2O多温相平衡预测

采用基于扩展UNIQUAC模型的AQSOL027软件针对水盐体系Li-Na-K-Mg-Cl-SO_4-H_2O的六个五元子体系进行相平衡理论计算,首次报道了多个五元子体系在多个温度条件下的相图,将计算相图与已有实验相图以及基于其它热力学模型的计算相图进行了比较。结果表明:AQSOL027在相图计算中具有较好的准确性和可靠性。AQSOL027是一款基于Excel、用于计算复杂多温水—盐体系Li-Na-K-Mg-Ca-H-Cl-(HSO_4/SO_4)-(CO_2/HCO_3/CO_3)-(H_3BO_3/B_4O_7/BO_2)-OH-H_2O的计算机软件,集成了扩展UNIQUAC模型和相平衡算法,该软件可结合SysCAD和Aspen过程模拟器,实现盐田工艺过程的计算机辅助设计。     

Meteorological observations 2005 and in a longer perspective at the Sermilik Station,Ammassalik Island,Southeast Greenland

Geografisk Tidsskrift, Danish Journal of Geography 107(2):59–68, 2007

Meteorological stations have been in operation since 1993 at Sermilik (65°40'N, 38°10'W), Ammassalik Island, southeast Greenland. This note presents meteorological observations for the year 2005 from the two meteorological stations: Station Nunatak (515 m a.s.l.) and Station Coast (25 m a.s.l.). The year 2005 was dominated by a temporal break down of Station Nunatak and defect instrument recordings at both stations. Average yearly measurements at both meteorological stations: Station Nunatak (1994–2004) and Station Coast (1998–2005) are presented to illustrate the climatic variability within the Mittivakkat Glacier catchment, and further compared with simultaneous standard synoptic meteorological DMI station climate observations (1994–2004) at the town Tasiilaq (Ammassalik) in order to establish transfer functions (linear regressions) and describe the spatial variability of climate variables.