Regional-Scale Excess Ar wave in a Barrovian type metamorphic belt, eastern Tibetan Plateau

Laser step heating ⁴⁰Ar/³⁹Ar analysis of biotite and muscovite single crystals from a Barrovian type metamorphic belt in the eastern Tibetan plateau yielded consistent cooling ages of ca. 40 Ma in the sillimanite zone with peak metamorphic temperatures higher than 600°C and discordant ages from 46 to 197 Ma in the zones with lower peak temperatures. Chemical Th‐U‐Total Pb Isochron Method (CHIME) monazite (65 Ma) and sensitive high mass‐resolution ion microprobe (SHRIMP) apatite (67 Ma) dating give the age of peak metamorphism in the sillimanite zone. Moderate amounts of excess Ar shown by biotite grains with ages of 46 to 94 Ma at metamorphic grades up to the high‐grade part of the kyanite zone probably represent incomplete degassing during metamorphism. In contrast, the high‐grade part of the kyanite zone yields biotite ages of 130 to 197 Ma. The spatial distribution of these older ages in the kyanite zone along the sillimanite zone boundary suggests they reflect trapped excess argon that migrated from higher‐grade regions. The most likely source is muscovite that decomposed to form sillimanite. The zone with extreme amounts of excess argon preserves trapped remnants of an ‘excess argon wave’. We suggest this corresponds to the area where biotite cooled below its closure temperature in the presence of an elevated Ar wave. Extreme excess Ar is not recognized in muscovite suggesting that the entrapment of the argon wave by biotite took place when the rocks had cooled down to temperatures lower than the closure temperature of muscovite. The breakdown of phengite during ultrahigh‐pressure (UHP) metamorphism may be a key factor in accounting for the very old apparent ages seen in many UHP metamorphic regions. This is the first documentation of a regional Ar‐wave spatially associated with regional metamorphism. This study also implies that resetting of the Ar isotopic systems in micas can require temperatures up to 600°C; much higher than generally thought.     

基于应变空间多机构CG模型的地震液化大变形数值分析

砂土液化导致的地基侧向大变形是地震中许多重要的工程设施和建筑物破坏的主要原因之一。简要介绍了可进行液化大变形分析的散粒体材料本构模型--应变空间多机构CG模型,基于FLIP ROSE程序平台,建立了预测和研究倾斜地基砂土液化导致侧向大变形的二维有限元数值分析方法。采用该模型对相同工况的土工动态离心模型试验进行了模拟,通过对比超孔隙水压力、剪切波水平加速度以及地基侧向位移发现,数值预测与试验结果吻合良好,从而验证了该有限元数值分析模型的可靠性。最后利用该数值分析模型预测了倾斜率不同的地基受到相同剪切波作用时,倾斜地基不同深度产生的侧向位移。预测结果显示,随着地基深度的减小,倾斜率对于地震液化导致倾斜地基侧向大变形的影响越来越显著。     

Behavior of zircon in the upper-amphibolite to granulite facies schist/migmatite transition, Ryoke metamorphic belt, SW Japan: constraints from the melt inclusions in zircon

Behavior of zircon at the schist/migmatite transition is investigated. Syn-metamorphic overgrowth is rare in zircon in schists, whereas zircon in migmatites has rims with low Th/U that give 90.3 ± 2.2 Ma U–Pb concordia age. Between inherited core and the metamorphic rim, a thin, dark-CL annulus containing melt inclusion is commonly developed, suggesting that it formed contemporaneous with the rim in the presence of melt. In diatexites, the annulus is further truncated by the brighter-CL overgrowth, suggesting the resorption and regrowth of the zircon after near-peak metamorphism. Part of the zircon rim crystallized during the solidification of the melt in migmatites. Preservation of angular-shaped inherited core of 5–10 μm in zircon included in garnet suggests that zircon of this size did not experience resorption but developed overgrowths during near-peak metamorphism. The Ostwald ripening process consuming zircon less than 5–10 μm is required to form new overgrowths. Curved crystal size distribution pattern for fine-grained zircons in a diatexite sample may indicate the contribution of this process. Zircon less than 20 μm is confirmed to be an important sink of Zr in metatexites, and ca. 35-μm zircon without detrital core are common in diatexites, supporting new nucleation of zircon in migmatites. In the Ryoke metamorphic belt at the Aoyama area, monazite from migmatites records the prograde growth age of 96.5 ± 1.9 Ma. Using the difference of growth timing of monazite and zircon, the duration of metamorphism higher than the amphibolite facies grade is estimated to be ca. 6 Myr.     

Electrical conduction and polarization of calcite single crystals

The electrical conductivity and polarization properties of calcite single crystals with three orientations, namely, a (00.1) plane perpendicular to the crystallographic c axis (10.0) plane parallel to the crystallographic c axis, and a (10.4) cleavage plane, were studied by both complex impedance and thermally stimulated depolarization current (TSDC) measurements. Conductivities for (00.1)-, (10.0)-, and (10.4)-oriented single calcite crystals at 400–600?°C were 1.16?×?10^?7?–?1.05?×?10^?5, 7.40?×?10^?8?–?4.27?×?10^?6, and 4.27?×?10^?7?–?2.86?×?10^?5 Ω^?1 m^?1, respectively, and the activation energies for conduction were 112, 103, and 101?kJ?mol^?1, respectively. The TSDC spectra verified the electrical polarizability of calcite crystals. The activation energy for depolarization, estimated from TSDC spectra, of the (00.1)-, (10.0)-, and (10.4)-oriented calcite substrates were 112, 119, and 114?kJ?mol^?1, respectively. Considering the correlation between the processes of conduction and electric polarization, we proposed the mechanisms of conduction and polarization in calcite on the assumption of oxide ion transport.     

Induced fabric under cyclic and rotational loads in a strain space multiple mechanism model for granular materials

The strain space multiple mechanism model idealizes the behavior of granular materials on the basis of a multitude of virtual simple shear mechanisms oriented in arbitrary directions. Within this modeling framework, the virtual simple shear stress is defined as a quantity dependent on the contact distribution function as well as the normal and tangential components of interparticle contact forces, which evolve independently during the loading process. In other terms, the virtual simple shear stress is an intermediate quantity in the upscaling process from the microscopic level (characterized by contact distribution and interparticle contact forces) to the macroscopic stress. The stress space fabric produces macroscopic stress through the tensorial average. Thus, the stress space fabric characterizes the fundamental and higher modes of anisotropy induced in granular materials. Herein, the induced fabric is associated with monotonic and cyclic loadings, loading with the rotation of the principal stress, and general loading. Upon loading with the rotation of the principal stress axis, some of the virtual simple shear mechanisms undergo loading whereas others undergo unloading. This process of fabric evolution is the primary cause of noncoaxiality between the axes of principal stresses and strains. Although cyclic behavior and behavior under the rotation of the principal stress axis seem to originate from two distinct mechanisms, the strain space multiple mechanism model demonstrates that these behaviors are closely related through the hysteretic damping factor. Copyright © 2011 John Wiley & Sons, Ltd.     

基于能量平衡法的黑河中游灌溉渠道蒸发量估算

干旱区绿洲农田需要通过渠道引水进行灌溉,在此过程中蒸发损失不可避免,尤其是在渠系广泛分布的情况下,蒸发损失量大,如何估算这部分量对于干旱区水资源管理至关重要. 以甘肃张掖绿洲灌溉渠道为研究对象,采用Mihara所提供的基于能量平衡理论的经验方程,对渠道蒸发量进行了估算.结果表明：蒸发率在晴天时白天大于夜间,阴天时白天和夜间相差不大,气温和相对湿度是影响蒸发率的主要因子.累计蒸发量在晴天时白天是夜间的2倍左右,阴天时白天和夜间相差不大,蒸发率是直接原因.     

Boron in metapelites controlled by the breakdown of tourmaline and retrograde formation of borosilicates in the Yanai area,Ryoke metamorphic belt,SW Japan

Tourmaline-out isograd formed by the breakdown of tourmaline is defined in the upper amphibolite-facies metapelites in the Yanai area, Ryoke metamorphic belt, SW Japan. The rim composition of tourmaline progressively becomes aluminous with ascending metamorphic grade, and the chemical zoning of tourmaline is controlled by ^X□AlNa_–1Mg_–1 and MgTi^YAl_–2 vectors in low- to medium-grade zones where muscovite is stable, whereas it is controlled by Mg(OH)^YAl_–1O_–1, CaMgO^X□_–1 ^YAl_–1(OH)_–1 and MgTi^YAl_–2 vectors in further higher–grade, muscovite-unstable zones. The size of tourmaline increases drastically where breakdown of muscovite+quartz takes place, probably due to the growth of tourmaline during breakdown of muscovite. On the high-temperature side of the tourmaline-out isograd, depletion of whole-rock boron is observed. Escape of boron-bearing melt or the fluid evolved from the melt during its crystallization probably caused this depletion, although locally trapped, boron-bearing melt or fluid formed irregularly shaped tourmaline and dumortierite during retrograde metamorphism.     

A Numerical Simulation on the Mesoscale Dynamics of the Spring Bloom in the Sea of Japan

The northward migration of spring bloom was observed in the Sea of Japan from April to May 1997 by the Ocean Color and Temperature Scanner (OCTS) on board the Advanced Earth Observing Satellite (ADEOS). This phenomenon is well simulated with a numerical ecosystem model coupled with a hydrodynamic model. The hydrodynamic model is the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model (MOM). The ecosystem model consists of five components: dissolved inorganic nitrogen (DIN), dissolved inorganic phosphorus (DIP), phytoplankton, zooplankton and detritus. Results of the numerical ecosystem model suggest that the mesoscale development of the spring bloom in the Sea of Japan is related to that of sea water temperature, and that the bloom is limited by the depletion of DIN. This revised version was published online in August 2006 with corrections to the Cover Date.     

Sea level response to wind field fluctuation around the tip of the Izu Peninsula

The mechanism of a characteristic sea level response (barotropic coastal ocean response) to wind field fluctuation around the tip of the Izu Peninsula observed during the middle of December 2000 to the middle of January 2001 was investigated based on three types of numerical experiments using the Princeton Ocean Model with various parameters. The response was characterized by the relaxation of sea level falling (rising) during eastward upwelling (westward downwelling) favorable wind regime. Analyses of quasi-realistic numerical model results in terms of the vertically integrated momentum balances and vorticity balance for the barotropic mode revealed that: 1) development/abatement of two anomalous circulations generated around the tip of the Izu Peninsula controls the sea level response through the acceleration/deceleration of a quasi-geostrophic barotropic coastal current between the circulations; 2) nonlinear vorticity advection by the Kuroshio Current and by the coastal current, coupled with vorticity diffusion, decelerates the quasi-geostrophic coastal current in the latter half of the wind regimes, which induces the relaxation of sea level rise/fall. The results of the quasi-realistic numerical experiment suggest that an analysis of the vorticity balance for the barotropic mode contributes to a better understanding of sea level responses to wind in coastal regions with strong currents and complex topography. In addition, a numerical experiment with idealized spatially uniform density stratification and a quasi-realistic wind field shows that if the Kuroshio Current had been shifted far offshore from the Izu Peninsula during the observation period, westward propagating continental shelf waves would have controlled the coastal sea level response.