川中地区茅口组两期流体叠合控制下的白云石化模式

岩相学特征、地化分析揭示:川中地区下二叠统茅口组白云岩存在三种白云石:①平直晶面细晶白云石(δ13C_PDB=3.06‰,δ18O_PDB=-6.81‰;Fe:1×10-6,Mn:未检出,Sr:150×10-6);②非平直晶面粗晶鞍状白云石(δ13C_PDB=3.22‰,δ18O_PDB=-7.82‰;Fe:149×10-6,Mn:185×10-6,Sr:85×10-6);③非平直晶面细晶白云石(δ13C_PDB=3.49‰,δ18O_PDB=-9.45‰;Th=123℃,S=133~139‰NaCl;Fe:58×10-6,Mn:59×10-6,Sr:76×10-6)。研究结果表明存在两期白云石化流体:①早期压实作用形成的埋藏白云石化流体(T=37.4℃;S=29.8‰NaCl);②晚期受构造控制的热液白云石化流体(Th=114.8℃;S=153~226‰NaCl)。明确了三种白云石的成因:平直晶面细晶白云石是由早期埋藏流体交代泥晶基质形成的;非平直晶面细晶白云石是平直晶面细晶白云石在受到热液流体改造后所形成的;而非平直晶面粗晶鞍状白云石则是由热液流体直接沉淀出的。基于研究结果建立了相应的白云化模式。     

High‐pressure metamorphism in Taiwan: from oceanic subduction to arc‐continent collision?

Abstract The Taiwan orogen has been the focus of a number of models of mountain building processes, but little attention has been paid to high‐pressure (HP) metamorphic rocks that are found as exotic blocks intermingled within the deepest units of the mountain belt. In this study, we re‐appraise from updated petrological and thermodynamic databases the physical conditions of HP metamorphism in Taiwan, and we combine our findings with available geochronological data to estimate the thermal history of these rocks. Our results indicate that peak metamorphic conditions of ～550 °C and 10–12 kbar have been followed by a rapid isothermal decompression, with exhumation possibly as rapid as burial. These units have subsequently been stored at a pressure of ～3 kbar for ～4–5 Myr, before their final exhumation, probably facilitated by the accretion of passive margin sequences during the Late Cenozoic collision. Therefore, HP units in Taiwan maintain a record of processes at depth from the early stages of oceanic subduction to the present arc‐continent collision.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

Submarine groundwater discharge revealed by aerial thermal infrared imagery: a case study on Jeju Island,Korea

Submarine groundwater discharge (SGD) is a global phenomenon that carries large volumes of groundwater and dissolved chemical species such as nutrient, metals, and organic compounds to coastal zones. We report the influence of SGD on the coastal waters of Jeju Island, Korea, using high‐resolution aerial thermal infrared (TIR) mapping techniques and field investigations. An aircraft‐based system was implemented using a cost‐effective TIR camera for aerial TIR mapping. Ground‐based calibrations and system integration with GPS/IMU (global positioning system/inertial measurement unit) were performed for the aerial systems. The aerial surveys showed distinct low‐temperature signatures of SGD along the coasts of Jeju Island, revealing large groundwater inputs from the coastal aquifers to the ocean. Multiple aerial surveys over a range of seasons and tidal stages revealed that SGD rates dynamically affect the sea surface temperature (SST) of the coastal zone. The in‐situ measurements supported that SGD has a substantial influence on the coastal water chemistry as well as SST. Our observations highlight the extent to which aerial‐based TIR mapping can serve as a powerful tool for studying SGD and other coastal processes. Copyright © 2016 John Wiley & Sons, Ltd.     

A coupled discrete element model for the simulation of soil and water flow through an orifice

Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd.     

A western North Pacific tropical cyclone intensity prediction scheme

A western North Pacific tropical cyclone (TC) intensity prediction scheme (WIPS) is developed based on TC samples from 1996 to 2002 using the stepwise regression technique, with the western North Pacific divided into three sub-regions: the region near the coast of East China (ECR), the South China Sea region (SCR), and the far oceanic region (FOR). Only the TCs with maximum sustained surface wind speed greater than 17.2 m s⁻¹ are used in the scheme. Potential predictors include the climatology and persistence factors, synoptic environmental conditions, potential intensity of a TC and proximity of a TC to land. Variances explained by the selected predictors suggest that the potential intensity of a TC and the proximity of a TC to land are significant in almost all the forecast equations. Other important predictors include vertical wind shear in ECR, 500-hPa geopotential height anomaly at the TC center, zonal component of TC translation speed in SCR, intensity change of TC 12 or 24 h prior to initial time, and the longitude of TC center in FOR.     

Lithium isotope systematics of lavas from the Cook–Austral Islands: Constraints on the origin of HIMU mantle

Radiogenic isotope variations in lavas from the Cook–Austral volcanic chain have delineated three distinct mantle sources: a HIMU component, a depleted component (DM), and an enriched component (EM). To better constrain the mantle sources for South Pacific hot spot volcanism, we determined lithium isotopic compositions of lavas from Raivavae, Rapa, Mangaia and Tubuai of the volcanic chain. The study includes whole rock and mineral analyses. In general, δ⁷Li of most olivines resemble bulk rock composition whereas clinopyroxenes are variably lighter. This implies that clinopyroxene is more susceptible to diffusion-induced fractionation, in agreement with previous studies. Olivine δ⁷Li values span a narrower range than whole rock values, and do not extend to the very heavy compositions previously reported in HIMU bulk lavas. This discrepancy likely reflects alteration of bulk lavas, and suggests that Li-isotope analyses of bulk lavas should be interpreted with caution. Olivines from HIMU lavas have heavy δ⁷Li values (up to 6.2‰), and extend beyond the range reported for fresh MORB. Correlations between Li-isotopes and radiogenic isotopes suggest that the heavy δ⁷Li values in HIMU olivines are a source characteristic and not the result of post-magmatic alteration. Although the Li-isotope composition of recycled, dehydrated oceanic crust is currently under debate, our results suggest that HIMU lavas derive from a source containing recycled dehydrated oceanic crust, and that the “heavy” Li-isotope signature of altered oceanic crust is partially preserved during passage through the subduction factory.     

Stable axisymmetric SPH formulation with no axis singularity

The axisymmetric formulation of the governing equations for geomechanics in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two forms of SPH discretization for the motion equations, which are labeled as form I and form II, are proposed, and the methods to compute the hoop stress and strain terms including hoop strain rate and the acceleration introduced by the hoop stress are compared. To avoid possible singularity problem near the axis of symmetry, a perfectly smooth contact along with ghost particles are applied to prevent the real particles from overly approaching the axis of symmetry to remove this potential singularity. In addition, the Mohr–Coulomb constitutive model is implemented into the SPH formulation in describing soil behavior. Four numerical tests are carried out to validate and compare the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions and results from FEM analysis. The performance in these comparisons suggests that SPH II with hoop terms computed through direct hoop method is more stable than the others, and the adoption of contact for the symmetric axis is efficient in eliminating the singularity problem. Copyright © 2015 John Wiley & Sons, Ltd.     

青藏高原片麻岩穹窿与找矿前景

片麻岩穹窿是研究俯冲-折返和碰撞-折返造山过程的重要窗口。已查明的大量青藏高原片麻岩穹窿(群)分布在古特提斯和新特提斯大洋俯冲-折返以及地体碰撞-折返过程中。松潘-甘孜造山带中雅江甲基卡片(麻)岩穹窿的三叠纪变质片岩的含矿伟晶岩脉中发现了超大型锂矿床,揭示片(麻)岩穹窿构造与同构造花岗岩、含矿伟晶岩脉以及大型印支滑脱带在时空和成因上有天然联系,为片麻岩穹窿的找矿前景提供了范例。