Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

Numerical modeling of three‐phase dissolution of underground cavities using a diffuse interface model

Natural evaporite dissolution in the subsurface can lead to cavities having critical dimensions in the sense of mechanical stability. Geomechanical effects may be significant for people and infrastructures because the underground dissolution may lead to subsidence or collapse (sinkholes). The knowledge of the cavity evolution in space and time is thus crucial in many cases. In this paper, we describe the use of a local nonequilibrium diffuse interface model for solving dissolution problems involving multimoving interfaces within three phases, that is, solid–liquid–gas as found in superficial aquifers and karsts. This paper generalizes developments achieved in the fluid–solid case, that is, the saturated case [1]. On one hand, a local nonequilibrium dissolution porous medium theory allows to describe the solid–liquid interface as a diffuse layer characterized by the evolution of a phase indicator (e.g., porosity). On the other hand, the liquid–gas interface evolution is computed using a classical porous medium two‐phase flow model involving a phase saturation, that is, generalized Darcy's laws. Such a diffuse interface model formulation is suitable for the implementation of a finite element or finite volume numerical model on a fixed grid without an explicit treatment of the interface movement. A numerical model has been implemented using a finite volume formulation with adaptive meshing (e.g., adaptive mesh refinement), which improves significantly the computational efficiency and accuracy because fine gridding may be attached to the dissolution front. Finally, some examples of three‐phase dissolution problems including density effects are also provided to illustrate the interest of the proposed theoretical and numerical framework. Copyright © 2014 John Wiley & Sons, Ltd.     

A major lava tube system from Undara Volcano,North Queensland

The Undara Volcano erupted 0.19 m.y. ago and formed lava fields covering 1,500 km² with a volume of approximately 23 km³. One of the flows extended 160 km on a gradient that averaged only 0.3°. This great length was a result of very high effusion rates, favourable topography and lava tube efficiency. The Undara lavas are rather uniform hawaiites. Lava temperatures are estimated to have been somewhat less than 1175–1220°C and viscosities greater than 10 to 30 Pa s. Long, apparently single lava tubes are well preserved in many places. They are marked by depressions, caves and long level ridges. A system of lava tubes extends for perhaps more than 100 km. The features of the lava tubes are comparable with those described elsewhere. Aligned depressions associated with caves appear to have formed contemporaneously. Most are much wider than the caves and probably represent collapsed lava ponds. The lava tubes appear to have formed by roofing over of lava channels. Close to lava tubes, the rocks developed strongly oxidised characteristics, such as oxidised olivine phenocrysts, ferric clinopyroxene and extensively developed hematite. Differentiated lava forms drips in some caves and is also oxidised.     

Attenuation relations for strong seismic ground motion in the Himalayan region

Strong motion data from various regions of India have been used to study attenuation characteristics of horizontal peak acceleration and velocity. The strong ground motion data base considered in the present work consists of various earthquakes recorded in the northern part of India since 1986 with magnitudes 5.7 to 7.2. Using these data, relations for horizontal peak acceleration and velocity, which are $$\begin{gathered} log_{10} a = 1.14 + 0.31M + 0.65log_{10} R \hfill \\ log_{10} v = 0.571 + 0.41M + 0.768log_{10} R \hfill \\ \end{gathered} $$ have been proposed wherea is the peak horizontal acceleration in cm/sec²,v is the peak horizontal velocity in mm/sec,M is body wave magnitude, andR is the hypocentral distance in km. The proposed relations are in reasonable agreement with the small amount of strong ground motion data available for the northern part of India. The present results will be useful in estimating strong ground motion parameters and in the earthquake resistant design in the Himalayan region.     

Gravel‐bed rivers: from particles to patterns

The papers in this special issue reflect several of the major themes and topics from the 7^th International Workshop on Gravel‐Bed Rivers. The papers focus primarily on aspects of bed material transport in gravel‐bed rivers and larger scale morpho‐dynamics. Research in gravel‐bed rivers is increasingly integrating processes over a wide range of temporal and spatial scales by combining field observation, lab experimentation, numerical modeling and theory testing in a range of river types, aided by new technological developments in particle tracking, computational modeling and high resolution spatial data. This is leading to greater understanding of the processes leading to distinctive morpho‐dynamics of river types and a more reliable basis for river management. Copyright © 2012 John Wiley & Sons, Ltd.     

The New Multichannel Radiospectrograph ARTEMIS-IV/HECATE, of the University of Athens

We present the new solar radiospectrograph of the University of Athensoperating at the Thermopylae Station since 1996. Observations cover thefrequency range from 110 to 688 MHz. The radiospectrograph has a 7-meterparabolic antenna and two receivers operating in parallel. One is a sweepfrequency receiver and the other a multichannel acousto-optical receiver.The data acquisition system consists of a front-end VME based subsystem anda Sun Sparc-5 workstation connected through Ethernet. The two subsystems areoperated using the VxWorks real-time package. The daily operation is fullyautomated: pointing of the antenna to the sun, starting and stopping theobservations at pre-set times, data acquisition, data compression by`silence suppression', and archiving on DAT tapes. The instrument can beused either by itself to study the onset and evolution of solar radio bursts or in conjunction with other instruments including theNançay Decametric Array and the WIND/WAVES RAD1 and RAD2 low frequencyreceivers to study associated interplanetary phenomena.     

Evaluating the effect of land development on sediment transport using a probability density function

An important problem in sedimentation analysis is the development of a channel section that preserves, as best as possible, the current sedimentation regime even though the flood frequency tendencies have been altered due to land development within the catchment. In order to accomplish this task, a methodology is needed that estimates sediment transport capacity for various channel configurations. Such a procedure is described which allows the computation of the total sediment transport capacity for each of several T-year return frequency runoff hydrographs. This information is used to obtain an approximate probability distribution for the total sediment transport capacity, and the mean and standard deviation of this distribution are computed.Comparing the results for the catchment in its present state with a future developed state, using a selection of new channel parameters, indicates how to improve the channel to control changes in sedimentation due to development. The analysis procedure provides a basis for estimating a new channel configuration such that the new flow conditions retain, as best as possible, the existing condition sedimentation effects, and hence retain the natural sediment supply and transport trends even though runoff flow rates have changed due to land development within the catchment.The results of Wilson Creek are typical of the several sites examined, see Table 3 below. The T=2, T=5, T=25, and T=100 year values for total sediment transport capacity, in kilotons, are 6.9, 39.4, 61.3, and 96.7 with a mean of 17.1 and standard deviation of 19.3. After development with no change in the channel the respective values increase to: 17.9, 84.6, 128.1, and 258.0 with a mean of 39.1 and standard deviation of 44.3. A new channel can be constructed which will reduce these sediment transport capacity values, after development, to 5.2, 41.0, 62.0, and 124.8 with a mean of 17.4 and standard deviation of 22.0.     

Characterizing spatial and temporal variation in 18O and 2H content of New Zealand river water for better understanding of hydrologic processes

Time series of hydrogen and oxygen stable isotope ratios (δ²H and δ¹⁸O) in rivers can be used to quantify groundwater contributions to streamflow, and timescales of catchment storage. However, these isotope hydrology techniques rely on distinct spatial or temporal patterns of δ²H and δ¹⁸O within the hydrologic cycle. In New Zealand, lack of understanding of spatial and temporal patterns of δ²H and δ¹⁸O of river water hinders development of regional and national-scale hydrological models. We measured δ²H and δ¹⁸O monthly, together with river flow rates at 58 locations across New Zealand over a two-year period. Results show: (a) general patterns of decreasing δ²H and δ¹⁸O with increasing latitude were altered by New Zealand's major mountain ranges; δ²H and δ¹⁸O were distinctly lower in rivers fed from higher elevation catchments, and in eastern rain-shadow areas of both islands; (b) river water δ²H and δ¹⁸O values were partly controlled by local catchment characteristics (catchment slope, PET, catchment elevation, and upstream lake area) that influence evaporation processes; (c) regional differences in evaporation caused the slope of the river water line (i.e., the relationship between δ²H and δ¹⁸O in river water) for the (warmer) North Island to be lower than that of the (cooler, mountain-dominated) South Island; (d) δ²H seasonal offsets (i.e., the difference between seasonal peak and mean values) for individual sites ranged from 0.50‰ to 5.07‰. Peak values of δ¹⁸O and δ²H were in late summer, but values peaked 1 month later at the South Island sites, likely due to greater snow-melt contributions to streamflow. Strong spatial differences in river water δ²H and δ¹⁸O caused by orographic rainfall effects and evaporation may inform studies of water mixing across landscapes. Generally distinct seasonal isotope cycles, despite the large catchment sizes of rivers studied, are encouraging for transit time analysis applications.