Discretizing the fracture-matrix interface to simulate solute transport

This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases.     

Airborne laser altimetry survey of Glaciar Tyndall, Patagonia

The first airborne laser altimetry measurements of a glacier in South America are presented. Data were collected in November of 2001 over Glaciar Tyndall, Torres del Paine National Park, Chilean Patagonia, onboard a Twin Otter airplane of the Chilean Air Force. A laser scanner with a rotating polygon-mirror system together with an Inertial Navigation System (INS) were fixed to the floor of the aircraft, and used in combination with two dual-frequency GPS receivers. Together, the laser–INS–GPS system had a nominal accuracy of 30 cm after data processing. On November 23rd, a total of 235 km were flown over the ablation area of Glaciar Tyndall, with 5 longitudinal tracks with a mean swath width of 300 m, which results in a point spacing of approximately 2 m both along and across track. A digital elevation model (DEM) generated using the laser altimetry data was compared with a DEM produced from a 1975 map (1:50,000 scale — Instituto Geográfico Militar (IGM), Chile). A mean thinning of − 3.1 ± 1.0 m a^− 1 was calculated for the ablation area of Glaciar Tyndall, with a maximum value of − 7.7 ± 1.0 m a^− 1 at the calving front at 50 m a.s.l. and minimum values of between − 1.0 and − 2.0 ± 1.0 m a^− 1 at altitudes close to the equilibrium line altitude (900 m a.s.l.). The thinning rates derived from the airborne survey were similar to the results obtained by means of ground survey carried out at  600 m of altitude on Glaciar Tyndall between 1975 and 2002, yielding a mean thinning of − 3.2 m a^− 1 [Raymond, C., Neumann, T.A., Rignot, E., Echelmeyer, K.A., Rivera, A., Casassa, G., 2005. Retreat of Tyndall Glacier, Patagonia, over the last half century. Journal of Glaciology 173 (51), 239–247.]. A good agreement was also found between ice elevation changes measured with laser data and previous results obtained with Shuttle Radar Topography Mission (SRTM) data. We conclude that airborne laser altimetry is an effective means for accurately detecting glacier elevation changes in Patagonia, where an ice thinning acceleration trend has been observed during recent years, presumably in response to warming and possibly also drier conditions.     

Mineralogy of Erromango lavas (New Hebrides): Evidence of an early stage of fractionation in island arc basalts

The orogenic volcanic rocks of Erromango island are divided into a Quaternary theoleiitic group and an older Pliocene one showing calc-alkalic affinity. The microprobe mineralogy of these lavas agrees with this geochemical distinction and indicates a marked iron enrichment trend in the recent tholeiitic lavas, whereas titanomagnetite occurs as early formed phenocrysts (Al, Cr, Mg-rich) only in Pliocene basalts. Some Ni and Mg-rich olivine phenocrysts (Fo88) and unusual chemically zoned clinopyroxene megacrysts are described from Pliocene basic lavas and interpreted as deep-fractionated minerals. The large compositional range of these latter (diopside core to augite rim, through salite intermediate zone) leads us to discuss the physical conditions and the corresponding cation substitutions prevailing during their formation. Diopside cores crystallized at depths of about 20 km from a relatively ‘primitive’ melt. The formation of salite is considered to be the result of concomitant fractionation of Mg-rich phases and increase of f_O2. An early stage of fractionation of Cr-bearing diopside + Ni-bearing olivine would account for the typical Cr and Ni depletion of Erromango lavas.     

Magnetic phase diagrams of Mg-, Fe-, (Cu + Ti)- and Cu-substituted braunite Mn<Superscript>2+</Superscript>Mn<Subscript>6</Subscript><Superscript>3+</Superscript>SiO<Subscript>12</Subscript>

 The magnetic behavior of the Jahn-Teller structure braunite, (Mn²⁺ _1−yM _y )(Mn³⁺ _{6− x} M_x)SiO₁₂, is strongly influenced by the incorporation of elements substituting manganese. Magnetic properties of well-defined synthetic samples were investigated in dependence on the composition. The final results are presented in magnetic phase diagrams. To derive the necessary data, ac susceptibility and magnetization of braunites with the substitutional elements M = Mg, Fe, (Cu+Ti) and Cu were measured. Whereas the antiferromagnetic ordering temperature, T _N , of pure braunite is hardly affected by the substitution of nonmagnetic Mg, it is rapidly suppressed by the substitution of magnetic atoms at the Mn positions. Typically for a concentration (x, y) ≥ 0.7 of the substituted elements, a spin glass phase occurs in the magnetic phase diagrams. Additionally, for the braunite system with Fe³⁺ substitutions, we observe in the concentration range 0.2 < x< 0.7 a double transition from the paramagnetic state, first to the antiferromagnetic state, followed by a transition to a spin glass state at lower temperatures. The unusual change of the magnetic properties with magnetic substitution at the Mn positions is attributed to the peculiar antiferromagnetic structure of braunite, which has been resolved recently. Received: 19 April 2001 / Accepted: 6 September 2001     

A place for place: Modelling and analysing platial representations

This editorial presents a special collection of papers addressing the concept of place and its use in geographical information science (GIScience). The concept of place is a topic of increasing interest among GIScience scholars. First attempts to formalise platial information have been made and it is increasingly held that user‐generated data sets in particular are often more platial than spatial in nature. At the same time, and especially when compared to geometric spatial concepts, the concept of place is ambiguous, complex and difficult to capture in formal and analytical terms, suggesting the need for interdisciplinary approaches. This collection presents articles covering a wide range of place‐related aspects, including both conceptual and more applied contributions. In the present editorial we summarise these and comment on their individual contributions, and hope that the readership of Transactions in GIS will find the special collection inspiring and informative.     

Temporal Geochemical Evolution of Neogene Magmatism in the Baguio Gold–Copper Mining District (Northern Luzon, Philippines)

Baguio, in the Central Cordillera of Northern Luzon, is a district that displays porphyry copper and epithermal gold mineralization, associated with Early Miocene–Pliocene–Quaternary calc‐alkaline and adakitic intrusions. Systematic sampling, K‐Ar dating, major and trace elements, and Sr, Nd, Pb isotopic analyses of fresh magmatic rocks indicate three magmatic pulses: an Early Miocene phase (21.2–18.7 Ma), a Middle–Late Miocene phase (15.3–8 Ma) and finally a Pliocene–Quaternary event (3–1 Ma). The first phase emplaced evolved calc‐alkaline magmas, essentially within the Agno Batholith complex, and is thought to be related to the westward‐dipping subduction of the West Philippine Basin. After a quiescence period during which the Kennon limestone was deposited, magmatic activity resumed at 15.3 Ma, in connection with the start of the subduction of the South China Sea along the Manila Trench. It emplaced first petrogenetically related and relatively unradiogenic low‐K calc‐alkaline lavas and intermediate adakites. Temporal geochemical patterns observed from 15.3 to 1 Ma include progressive enrichment in K and other large ion lithophile elements, increase in radiogenic Sr and Pb and corresponding decrease in radiogenic Nd. These features are thought to reflect the progressive addition to the Luzon arc mantle wedge of incompatible elements largely inherited from South China Sea sediments. The origin of the long quiescence period, from 8 to 3 Ma, remains problematic. It might represent a local consequence of the docking of the Zambales ophiolitic terrane to Northern Luzon. Then, magmatic activity resumed at 3 Ma, emplacing chemically diversified rocks ranging from low K to high K and including a large proportion of adakites, especially during the Quaternary (dacitic plugs). The authors tentatively relate this diversity to the development of a slab tear linked with the subduction of the fossil South China Sea ridge beneath the Baguio area.     

Exploring aberrant bivalve shell ultrastructure and geochemistry as proxies for past sea water acidification

Throughout much of Earth's history, marine carbonates have represented one of the most important geological archives of environmental change. Several pivotal events during the Phanerozoic, such as mass extinctions or hyperthermal events have recently been associated with ocean acidification. Nevertheless, well‐defined geological proxies for past ocean acidification events are, at best, scarce. Here, experimental work explores the response of bivalve shell ultrastructure and isotope geochemistry (δ¹³C, δ¹⁸O and δ²⁶Mg) to stressful environments, in particular to sea water acidification. In this study, the common blue mussel, Mytilus edulis, was cultured (from early juvenile stages to one year of age) at four pH regimes (pH_NBS 7·2 to pH 8·0). Shell growth rate and ultrastructure of mainly the calcitic portion of the shells were compared between experimental treatments. Specimens exposed to low‐pH environments show patches of disordered calcitic fibre orientation in otherwise well‐structured shells. Furthermore, the electron backscattered diffraction analyses reveal that, under acidified conditions, the c‐axis of the calcite prisms exhibits a bimodal or multi‐modal distribution pattern. Similar shell disorder patterns have been reported from mytilids kept under naturally acidified sea water conditions. In contrast, this study found no evidence that different pH regimes affect shell carbon, oxygen or magnesium isotope ratios. Based on these observations, it is proposed that: (i) stressful environments, in this case low sea water pH, predictably affect bivalve biomineralization patterns; and (ii) these findings bear potential as a novel (petrographic) proxy for ancient sea water acidification. An assessment of the applicability of these data to well‐preserved fossil shell material from selected time intervals requires additional work.     

Novel Algorithms for Astronomical Plate Analyses

Powerful computers and dedicated software allow effective data mining and scientific analyses in astronomical plate archives. We give and discuss examples of newly developed algorithms for astronomical plate analyses, e.g., searches for optical transients, as well as for major spectral and brightness changes.     

Statistical Properties of the Immediate Aftershocks of the 15 October 2013 Magnitude 7.1 Earthquake in Bohol,Philippines

The aftershock records of the magnitude 7.1 earthquake that hit the island of Bohol in central Philippines on 15 October 2013 is investigated in the light of previous results for the Philippines using historical earthquakes. Statistics of interevent distances and interevent times between successive aftershocks recorded for the whole month of October 2013 show marked differences from those of historical earthquakes from two Philippine catalogues of varying periods and completeness levels. In particular, the distributions closely follow only the regimes of the historical distributions that were previously attributed to the strong spatio-temporal correlations. The results therefore suggest that these correlated regimes which emerged naturally from the analyses are strongly dominated by the clustering of aftershock events.