Exhumation of eclogite and blueschist (Cyclades,Greece): Pressure–temperature evolution determined by thermobarometry and garnet equilibrium modelling

High‐P rocks such as eclogite and blueschist are metamorphic markers of palaeo‐subduction zones, and their formation at high‐P and low‐T (HP–LT) conditions is relatively well understood since it has been the focus of numerous petrological investigations in the past 40 years. The tectonic mechanisms controlling their exhumation back to the surface are, however, diverse, complex and still actively debated. Although the Cycladic Blueschist Unit (CBU, Greece) is among the best worldwide examples for the preservation of eclogite and blueschist, the proposed P–T evolution followed by this unit within the Hellenic subduction zone is quite different from one study to another, hindering the comprehension of exhumation processes. In this study, we present an extensive petrological data set that permits refinement of the shape of the P–T trajectory for different subunits of the CBU on Syros. High‐resolution quantitative compositional mapping has been applied to support the thermobarometric investigations, which involve semi‐empirical thermobarometry, garnet equilibrium modelling and P–T isochemical phase diagrams. The thermodynamic models highlight the powerful use of reactive bulk compositions approximated from local bulk compositions. The results are also combined with Raman spectrometry of carbonaceous material (RSCM) to retrieve the metamorphic peak temperature distribution at the scale of the island. A major result of this study is the good agreement between all the independent thermobarometric methods, permitting reconstruction of the prograde and retrograde P–T trajectories. Garnet compositional zoning was used to retrieve prograde, peak and retrograde growth stages in line with the results of the P–T isochemical phase diagrams, RSCM temperature and peak‐pressure crystallization of the garnet–omphacite–phengite assemblage. Our results are consistent with previous thermobarometric estimates from other occurrences of CBU rocks (Tinos, Andros), suggesting a multistage exhumation process with (1) early syn‐orogenic exhumation within the subduction channel, (2) isobaric heating at mid‐crustal depths (~10–12 kbar) following thermal re‐equilibration of the lithosphere from a cold syn‐orogenic regime in the subduction zone to a warmer post‐orogenic regime in the back‐arc domain and (3) exhumation and cooling related to a post‐orogenic phase of extension following slab retreat. Expanding to the general aspects of subduction zones, we suggest that such metamorphic evolution of HP–LT units should be regarded as a characteristic feature of exhumation driven by slab rollback.     

Numerical simulation of the heavy rainfall caused by a convection band over Korea: a case study on the comparison of WRF and CReSS

Natural Hazards - This study investigates the capability of two numerical models, namely the Weather Research and Forecasting (WRF) and Cloud Resolving Storm Simulator (CReSS), to simulate the...     

Geotechnical parameter estimation in tunnelling using relative convergence measurement

Accurate estimation of geotechnical parameters is an important and difficult task in tunnel design and construction. Optimum evaluation of the geotechnical parameters have been carried out by the back‐analysis method based on estimated absolute convergence data. In this study, a back‐analysis technique using measured relative convergence in tunnelling is proposed. The extended Bayesian method (EBM), which combines the prior information with the field measurement data, is adopted and combined with the 3‐dimensional finite element analysis to predict ground motion. By directly using the relative convergence as observation data in the EBM, we can exclude errors that arise in the estimation of absolute displacement from measured convergence, and can evaluate the geotechnical parameters with sufficient reliability. The proposed back‐analysis technique is applied and validated by using the measured data from two tunnel sites in Korea. Copyright © 2005 John Wiley & Sons, Ltd.     

A simulation using a hybrid method for predicting fault zones ahead of a tunnel face

By monitoring and analysing three‐dimensional absolute displacement during tunnel excavation, it is possible to predict discontinuity zones ahead of tunnel face. This paper presents results of the three‐dimensional displacement obtained from three‐dimensional simulations using hybrid method for assessing the effect of discontinuity zones on monitoring data. The influence of fault zones and in situ initial stresses on the three‐dimensional tunnel movement is evaluated. Copyright © 2002 John Wiley & Sons, Ltd.     

Propagation characteristics of historical tsunamis that attacked the east coast of Korea

In this study, a numerical modeling system based on the dispersion–correction finite difference scheme equipped with a grid-nesting scheme is constructed. The model is applied to simulate the propagation of three historical tsunami events that attacked the east coast of Korea. The calculated free-surface displacements for the cases of the 1983 Akita and the 1993 Okushiri tsunamis are compared with the observations at four tidal stations along the east coast of Korea. The comparison shows that the results agree well with the observations. The analyses of the simulated results show that underwater topography, such as submerged rises and ridges, plays an important role in the propagation of tsunamis in this region.     

Underestimation of roughness in rough rock joints

Numerous studies have been made to improve Barton's shear strength model for the quantification of rock joints. However, in these previous studies, the roughness and shear strength of the rock joint have been underestimated especially for relatively high undulated profiles (joint roughness coefficient (JRC) >14). The main factors of roughness underestimation in rough rock joints are investigated for the proper quantification of rock joint roughness. The aliasing effect and the roughness characteristics are analyzed by using artificial joint profiles and natural rock joint profiles. A 3D camera scanner is adopted to verify the main source of underestimation when using conventional measurement methods. Shear strength tests are carried out by using two types of shear apparatus to study the roughness mobilization characteristics, which may also affect the roughness underestimation. The results of joint roughness assessment, such as aliasing and undulation of waviness, show that the roughness can be underestimated in relatively rough joint profiles (JRC>14). At least two components of roughness parameters are needed to properly represent the joint roughness, for example, the amplitude and the inclination angle of joint asperity. Roughness mobilization is affected by both the normal stress and the asperity scale. Copyright © 2008 John Wiley & Sons, Ltd.     

Harpacticoid copepod response to epiphyte load variations in Posidonia oceanica (L.) Delile meadows

We conducted a field experiment to assess the response of phytal harpacticoids to nutrient‐driven increases of epiphyte load in Posidonia oceanica meadows. First, we evaluated differences in species richness, diversity and assemblage structure of phytal harpacticoids in P. oceanica meadows with differing epiphyte loads. Secondly, we conducted a field experiment where epiphyte load was increased through an in situ addition of nutrients to the water column and evaluated the responses of the harpacticoid assemblages. We predicted that there would be changes in the harpacticoid assemblages as a result of nutrient‐driven increases of epiphyte load, and that these changes would be of a larger magnitude in meadows of low epiphyte load. Our results show that the harpacticoid fauna (>500 μm) present in P. oceanica meadows in the Bay of Palma comprised taxa which are considered phytal and other less abundant ones previously described as sediment dwellers or commensal on other invertebrate species. Nutrient addition had an overall significant effect on epiphyte biomass and on harpacticoid abundance, diversity and assemblage structure, possibly as a response to the increased resources and habitat complexity provided by epiphytes. The abundance of dominant species at each location was favoured by nutrient addition and in some cases correlated with epiphytic biomass, although never strongly. This may indicate that structural complexity or diversity of the epiphytic cover might be more important than the actual epiphytic biomass for the harpacticoid species investigated. More species‐specific studies are necessary to ascertain this and clarify the relationships between harpacticoids and epiphytes in seagrass meadows. To our knowledge, this is the first account of harpacticoid species associated with P. oceanica leaves and the epiphytic community they harbour in the Mediterranean Sea.     

Reinforcement and Arching Effect of Geogrid-Reinforced and Pile-Supported Embankment on Marine Soft Ground

The effectiveness of constructing a geogrid-reinforced and pile supported embankment on soft ground to reduce differential settlement has been studied by pilot scale field tests and numerical analysis. Three-by-three pile groups with varying pile spacing were driven into a layer of soft ground, and a layer of geogrid was used as reinforcement over each pile group. Further, a 2-D numerical analysis has been conducted using the computer program FLAC 2D. The mechanisms of load transfer can be considered as a combination of embankment soil arching, geogrid tension, and stress transfer due to the difference in stiffness between pile and soft ground. Based on the pilot scale field tests and results of numerical analysis, we find that the geosynthetic reinforcement slightly interferes with soil arching, and helps reduce differential settlement of the soft ground. Also, the most effective load transfer and vertical stress reduction at the midspan between piles occurs when the pile cap spacing index D/b (D: pile cap spacing, b: diameter of pile) is 3.0.     

A condensation model for the formation of chondrules in enstatite chondrites

Abstract— It is proposed that the chondrules in enstatite chondrites formed near the Sun from rain‐like supercooled liquid silicate droplets and condensed Fe‐Ni alloys in thermodynamic equilibrium with a slowly cooling nebula. FeO formed and dissolved in the droplets in an initial stage when the nucleation of iron was blocked, and was later mostly reduced to unalloyed Fe. At high temperatures, the silicate droplets contained high concentrations of the less volatile components CaO and Al₂O₃. At somewhat lower temperatures the equilibrium MgO content of the droplets was relatively high. As cooling progressed, some droplets gravitated toward the Sun, and moved in other directions, depleting the region in CaO, Al₂O₃, and MgO and accounting for the relatively low observed CaO/SiO₂, Al₂O₃/SiO₂, and MgO/SiO₂ ratios in enstatite chondrites. At approximately 1400 K, the remaining supercooled silicate droplets crystallized to form MgSiO₃ (enstatite) with small amounts of olivine and a high‐SiO₂ liquid phase which became the mesostases. The high enstatite content is the result of the supercooled chondrules crystallizing at a relatively low temperature and relatively high total pressure. Finally, FeS formed at temperatures below 680 K by reaction of the condensed Fe with H₂S. All calculations were performed with the evaluated optimized thermodynamic databases of the FactSage thermodynamic computer system. The thermodynamic properties of compounds and solutions in these databases were optimized completely independently of any meteoritic data. Agreement of the model with observed bulk and phase compositions of enstatite chondrules is very good and is generally within experimental error limits for all components and phases.