Occurrence of Halogen-rich Phlogopite in Late Cenozoic Volcanic Rocks in the Japanese Arcs

Abstract. Halogen-rich phlogopite occurs in the groundmass of andesite and dacite lavas from Late Tertiary to Quaternary volcanoes associated with native sulfur and limonite deposits (Shiretoko-Iwozan, Hachimantai, Adatara, Omeshidake, Masaki) and hydrothermal ore deposits (Harukayama, Muineyama, Hishikari) in Japan. The F contents of the halogen-rich phlogopite range from 3.6 to 5.7 wt%, corresponding to atomic F/(F+C1+OH) ratios ranging from 0.45 to 0.69. On the other hand, the Cl contents of the halogen-rich phlogopite are around 0.2 wt%. The atomic Mg/(Mg+Fe) ratios range from 0.69 to 0.83.
The fluorine intercept value [IV(F)] defined by Munoz (1984) of the phlogopites ranges from 0.79 to 3.17, and the chlorine intercept value [IV(Cl)] ranges from -7.11 to -7.77. The observed IV(F) of the phlogopites broadly overlap the range of the IV(F) for biotites from porphyry copper deposits. On the other hand, the observed IV(Cl) are significantly lower than the IV(Cl) for biotites from porphyry copper deposits. Whereas the F contents of the phlogopite appear more prominent compared to the Cl contents, the calculation of halogen intercept values revealed that the phlogopites are enriched in Cl with respect to the element distribution effect of Mg-Fe substitution. Since the degree of Cl enrichment of the phlogopite is more significant compared to that of biotite in porphyry copper deposits, the phlogopites are considered to have formed under the condition of significantly high activity of halogens. Hydrothermal ore deposits may be formed in magmatic hydrothermal system associated with volcanoes where halogen-rich phlogopite is formed by hypersaline fluid.     

Primary Ore Mineral Assemblage and Fluid Inclusion Study of the Batu Hijau Porphyry Cu-Au Deposit, Sumbawa, Indonesia

Abstract. The Batu Hijau porphyry Cu‐Au deposit, Sumbawa Island, Indonesia, is associated with a tonalitic intrusive complex. The temperature‐pressure condition of mineralization at the Batu Hijau deposit is discussed on the basis of fluid inclusion microthermometry. Then, the initial Cu‐Fe sulfide mineral assemblage is discussed. Bornite and chalcopyrite are major copper ore minerals associated with quartz veinlets. The quartz veinlets have been classified into ‘A’ veinlets associated with bornite, digenite, chalcocite and chalcopyrite, ‘B’ veinlets having chalcopyrite bornite along vuggy center‐line, rare ‘C’ chalcopyrite‐quartz veinlets, and late ‘D’ veinlets consisting of massive pyrite and quartz (Clode et al., 1999). Copper and gold mineralization is associated with abundant ‘A’ quartz veinlets. Abundant fluid inclusions are found in veinlet quartz consisting mainly of gas‐rich inclusions and polyphase inclusions throughout the veinlet types. The hydrothermal activity occurred in temperature‐pressure conditions of aqueous fluid immiscibility into hypersaline brine and dilute vapor. The halite dissolution (Tm[halite]) and liquid‐vapor homogenization (Th) temperatures of the polyphase inclusions in veinlet quartz range from 270 to 472d?C and from 280 to 454d?C, respectively. The estimated salinity ranges from 36 to 47 wt% (NaCl equiv.). The apparent pressures lower than 300 bars are estimated to have been along the liquid‐vapor‐halite curve for the fluid inclusions having the Th lower than the Tm that trapped the brine saturated with halite, or at slightly higher pressure relative to liquid‐vapor‐halite curve for the fluid inclusions having the Th higher than the Tm that trapped the brine unsaturated with halite. The actual temperature and pressure during the hydrothermal activity at the Batu Hijau deposit are estimated to have been around 300d?C and 50 bars. At such temperature‐pressure conditions, the principal and initial Cu‐Fe sulfide mineral assemblages are thought to be chalcopyrite + bornite solid solution (bnss) for the chalcopyrite‐bearing assemblage, and chalcocite‐digenite solid solution and bnss for the chalcopyrite‐free assemblage.     

Evolution of Hydrothermal System at the Dizon Porphyry Cu-Au Deposit, Zambales, Philippines

Abstract. Evolution of hydrothermal system from initial porphyry Cu mineralization to overlapping epithermal system at the Dizon porphyry Cu‐Au deposit in western central Luzon, Zambales, Philippines, is documented in terms of mineral paragen‐esis, fluid inclusion petrography and microthermometry, and sulfur isotope systematics. The paragenetic stages throughout the deposit are summarized as follows; 1) stockwork amethystic quartz veinlets associated with chalcopyrite, bornite, magnetite and Au enveloped by chlorite alteration overprinting biotite alteration, 2) stockwork quartz veinlets with chalcopyrite and pyrite associated with Au and chalcopyrite and pyrite stringers in sericite alteration, 3) stringer quartz veinlets associated with molybdenite in sericite alteration, and 4) WNW‐trending quartz veins associated with sphalerite and galena at deeper part, while enargite and stibnite at shallower levels associated with advanced argillic alteration. Chalcopyrite and bornite associated with magnetite in quartz veinlet stockwork (stage 1) have precipitated initially as intermediate solid solution (iss) and bornite solid solution (bnss), respectively. Fluid inclusions in the stockwork veinlet quartz consist of gas‐rich inclusions and polyphase inclusions. Halite in polyphase inclusions dissolves at temperatures ranging from 360d?C to >500d?C but liquid (brine) and gas (vapor) do not homogenize at <500d?C. The maximum pressure and minimum temperature during the deposition of iss and bnss with stockwork quartz veinlets are estimated to be 460 bars and 500d?C. Fluid inclusions in veinlet stockwork quartz enveloped in sericite alteration (stage 2) consist mainly of gas‐rich inclusions and polyphase inclusions. In addition to the possible presence of saturated NaCl crystals at the time of entrapment of fluid inclusions that exhibit the liquid‐vapor homogenization temperatures lower than the halite dissolution temperatures in some samples, wide range of temperatures of halite dissolution and liquid‐vapor homogenization of polyphase inclusions from 230d?C to >500d?C and from 270d?C to >500d?C, respectively, suggests heterogeneous entrapment of gaseous vapor and hypersaline brine. The minimum pressure and temperature are estimated to be about 25 bars and 245d?C. Fluid inclusions in veinlet quartz associated with molybdenite (stage 3) are dominated by gas‐rich inclusions accompanied with minor liquid‐rich inclusions that homogenize at temperatures between 350d?C and 490d?C. Fluid inclusions in vuggy veinlet quartz associated with stibnite (stage 4) consist mainly of gas‐rich inclusions with subordinate polyphase inclusions that do not homogenize below 500d?C. Fluid inclusions in veinlet quartz associated with galena and sphalerite (stage 4) are composed of liquid‐rich two‐phase inclusions, and they homogenize into liquid phase at temperatures ranging widely from 190d?C to 300d?C (suggesting boiling) and the salinity ranges from 1.0 wt% to 3.4 wt% NaCl equivalent. A pressure of about 15 bars is estimated for the dilute aqueous solution of 190d?C from which veinlet quartz associated with galena and sphalerite precipitated. In addition to a change in temperature‐pressure regime from lithostatic pressure during the deposition of iss and bnss with stockwork quartz veinlets to hydrostatic pressure during fracture‐controlled quartz veinlet associated with galena and sphalerite, a decrease in pressure is supposed to have occurred due to unroofing or removal of the overlying piles during the temperature decrease in the evolution of hydrothermal system. The majority of the sulfur isotopic composition of sulfides ranges from ±0 % to +5 %. Sulfur originated from an iso‐topically uniform and homogeneous source, and the mineralization occurred in a single hydrothermal system.     

Deep structure of the northeastern Japan arc and its implications for crustal deformation and shallow seismic activity

Seismic tomography studies in the northeastern Japan arc have revealed the existence of an inclined sheet-like seismic low-velocity and high-attenuation zone in the mantle wedge at depths shallower than about 150 km. This sheet-like low-velocity, high-attenuation zone is oriented sub-parallel to the subducted slab, and is considered to correspond to the upwelling flow portion of the subduction-induced convection. The low-velocity, high-attenuation zone reaches the Moho immediately beneath the volcanic front (or the Ou Backbone Range) running through the middle of the arc nearly parallel to the trench axis, which suggests that the volcanic front is formed by this hot upwelling flow. Aqueous fluids supplied by the subducted slab are probably transported upward through this upwelling flow to reach shallow levels beneath the Backbone Range where they are expelled from solidified magma and migrate further upward. The existence of aqueous fluids may weaken the surrounding crustal rocks, resulting in local contractive deformation and uplift along the Backbone Range under the compressional stress field of the volcanic arc. A strain-rate distribution map generated from GPS data reveals a notable concentration of east–west contraction along the Backbone Range, consistent with this interpretation. Shallow inland earthquakes are also concentrated in the upper crust of this locally large contraction deformation zone. Based on these observations, a simple model is proposed to explain the deformation pattern of the crust and the characteristic shallow seismic activity beneath the northeastern Japan arc.     

U-Pb and Pb-Pb ages of zircons from basaltic eucrites: Implications for early basaltic volcanism on the eucrite parent body

We have undertaken petrologic and SHRIMP U-Th-Pb isotopic studies on zircons from basaltic eucrites (Yamato [Y]-75011, Y-792510, Asuka [A]-881388, A-881467 and Padvarninkai) with different thermal and shock histories. Eucritic zircons are associated with ilmenite in most cases and have subhedral shapes in unmetamorphosed and metamorphosed eucrites. Some zircons in highly metamorphosed eucrites with granulitic texture occur alone in pyroxene, and typically have rounded to subrounded shapes due to recrystallization. Superchondritic Zr/Hf ratios of eucritic zircons indicate that they crystallized from incompatible element-rich melts after crystallization of ilmenite. Concentrations of uranium and thorium in zircons in the unmetamorphosed eucrite Y-75011 are higher than those in metamorphosed eucrites.The U-Pb systems of eucritic zircons are almost concordant but some zircon grains show reverse discordance. Radiogenic lead-loss up to 48% from zircons is observed in the shock-melted eucrite Padvarninkai. The ²⁰⁷Pb-²⁰⁶Pb ages of zircon in Y-75011 (4550 ± 9 Ma, n = 5) are nearly identical, within analytical uncertainty, to the ages of zircons from the metamorphosed eucrite Y-792510 (4545 ± 15 Ma, n = 13), the highly metamorphosed eucrites A-881388 (4555 ± 54 Ma, n = 5) and A-881467 (4558 ± 13 Ma, n = 8), and the shock-melted eucrite Padvarninkai (4555 ± 13 Ma, n = 18). The averaged ²⁰⁷Pb-²⁰⁶Pb age of zircon from five eucrites analyzed in this study is 4554 ± 7 Ma (95% confidence limits, n = 49), indistinguishable from the averaged U-Pb age (4552 ± 9 Ma) of the same samples. Because of the high closure temperature of lead in zircon (T_closure = ∼1050°C with a cooling rate of 0.2°C/yr), the ²⁰⁷Pb-²⁰⁶Pb ages of eucritic zircon do not represent metamorphic ages but crystallization ages of extrusive lavas.This fact strongly suggests that volcanism of the eucrite parent body occurred at a very early stage of the Solar System history, 7-20 Ma after CAI formation (4567.2 ± 0.6 Ma), thus basaltic eucrites crystallized from parental magmas within a short interval following the differentiation of their parent body. The U-Pb ages of eucritic zircons are older than the U-Pb, Sm-Nd and Rb-Sr ages of some basaltic eucrites, which is consistent with differences in closure temperatures of each isotopic system, and suggests that thermal and shock metamorphism affected the isotopic systems of pyroxene, plagioclase and phosphates.     

Continuous measurement of whole‐tree water balance for studying urban tree transpiration

Street and garden trees in urban areas are often exposed to advection of strong vapour pressure deficit (VPD) air that can raise the whole‐tree transpiration rate (E_T), known as the oasis effect. However, urban trees tend to have small soil volume compared with natural conditions, and so they are believed to strongly regulate stomata. E_T characteristics of such urban trees have not been well understood because of a lack of reliable measurement methods. Therefore, we propose a novel weighing lysimeter method and investigate the whole‐tree water balance of an isolated container‐grown Zelkova serrata to examine (a) which biotic and abiotic factors determine E_T and (b) which spatial and temporal information is needed to predict E_T under urban conditions. Whole‐tree water balance and environmental conditions were measured from 2010 to 2012. Although leaf area substantially increased in the study period, daily E_T did not vary much. E_T increased with VPD almost linearly in 2010 but showed saturation in 2011 and 2012. Root water uptake lagged E_T by 40 min in 2012. These results suggest that the small planter box interfered with root growth and that hydraulic supply capacities did not increase sufficiently to support leaf area increase. From analysis of water balance, we believe that neglecting soil drought effects on street trees without irrigation in Japan will overestimate E_T over 4–5 sunny days at the longest. This is unlike previous studies of forest.     

Tropical synoptic-scale wave disturbances over the western Pacific simulated by a global cloud-system resolving model

Lower-tropospheric tropical synoptic-scale disturbances (TSDs) are associated with severe weather systems in the Asian Monsoon region. Therefore, exact prediction of the development and behavior of TSDs using atmospheric general circulation models is expected to improve weather forecasting for this region. Recent state-of-the art global cloud-system resolving modeling approaches using a Nonhydrostatic Icosahedral Atmospheric Model (NICAM) may improve representation of TSDs. This study evaluates TSDs over the western Pacific in output from an Atmospheric Model Intercomparison Project (AMIP)-like control experiment using NICAM. Data analysis compared the simulated and observed fields. NICAM successfully simulates the average activity, three-dimensional structures, and characteristics of the TSDs during the Northern summer. The variance statistics and spectral analysis showed that the average activity of the simulated TSDs over the western Pacific during Northern summer broadly captures that of observations. The composite analysis revealed that the structures of the simulated TSDs resemble the observed TSDs to a large degree. The simulated TSDs exhibited a typical southeast- to northwest-oriented wave-train pattern that propagates northwestward from near the equator around 150 ^° E toward the southern coast of China. However, the location of the simulated wave train and wave activity center was displaced northward by approximately a few degrees of latitude from that in the observation. This displacement can be attributed to the structure and strength of the background basic flow in the simulated fields. Better representation of the background basic states is required for more successful simulation of TSDs.     

Using a regional numerical weather prediction model for GNSS positioning over Brazil

The global navigation satellite system (GNSS) can provide centimeter positioning accuracy at low costs. However, in order to obtain the desired high accuracy, it is necessary to use high-quality atmospheric models. We focus on the troposphere, which is an important topic of research in Brazil where the tropospheric characteristics are unique, both spatially and temporally. There are dry regions, which lie mainly in the central part of the country. However, the most interesting area for the investigation of tropospheric models is the wet region which is located in the Amazon forest. This region substantially affects the variability of humidity over other regions of Brazil. It provides a large quantity of water vapor through the humidity convergence zone, especially for the southeast region. The interconnection and large fluxes of water vapor can generate serious deficiencies in tropospheric modeling. The CPTEC/INPE (Center for Weather Forecasting and Climate Studies/Brazilian Institute for Space Research) has been providing since July 2012 a numerical weather prediction (NWP) model for South America, known as Eta. It has yield excellent results in weather prediction but has not been used in GNSS positioning. This NWP model was evaluated in precise point positioning (PPP) and network-based positioning. Concerning PPP, the best positioning results were obtained for the station SAGA, located in Amazon region. Using the NWP model, the 3D RMS are less than 10 cm for all 24 h of data, whereas the values reach approximately 60 cm for the Hopfield model. For network-based positioning, the best results were obtained mainly when the tropospheric characteristics are critical, in which case an improvement of up to 7.2 % was obtained in 3D RMS using NWP models.     

Temporal variations in the current structure and volume transport of the Coastal Oyashio revealed by direct current measurement

Direct current measurements by a shipboard and bottom-mounted acoustic Doppler current profiler and concurrent hydrographic observations with a CTD were conducted off southeastern Hokkaido, Japan, between January and May 2005 to reveal temporal variations in the current structure and volume transport of the Coastal Oyashio (CO). The CO, which has a baroclinic jet structure with southwestward speeds exceeding 90 cm s^?1 and a width of 7–8 km, was associated with a surface-to-bottom density front and was formed on the offshore side of the shelf break. The volume transport of CO (T _CO) was estimated by integrating the fluxes of lower-density water that was trapped against the coast along the density front represented by the 26.2 σ _θ isopycnal line. This transport decreased monotonously from 0.79 Sv (1 Sv = 10⁶ m³ s^?1) in January to 0.21 Sv in March and subsequently to 0.12 Sv in May, possibly due to the decay of the East Sakhalin Current Water in the Okhotsk Sea. Accompanied by a decrease in T _CO, the location of the jet structure associated with the density front moved toward the coast while the maximum speed of the jet decreased and the tilt of the front became more horizontal. Consequently, more saline offshore Oyashio water flowed into the deep part of the shelf area, and the current structure altered from relatively barotropic in winter to baroclinic in spring. This study is the first to estimate the observed volume transport of the CO from direct current measurements.     

Limit states and failure mechanisms of viscous dampers and the implications for large earthquakes

Fluid viscous dampers are used to control story drifts and member forces in structures during earthquake events. These elements provide satisfactory performance at the design‐level or maximum considered earthquake. However, buildings using fluid viscous dampers have not been subjected to very large earthquakes with intensities greater than the design and maximum considered events. Furthermore, an extensive database of viscous damper performance during large seismic events does not exist. To address these issues, a comprehensive analytical and experimental investigation was conducted to determine the performance of damped structures subjected to large earthquakes. A critical component of this research was the development and verification of a detailed viscous damper mathematical model that incorporates limit states. The development of this model and the laboratory and simulation results conclude good correlation with the new model and the damper limit states and provide superior results compared with the typical damper model when considering near collapse evaluation of structures. Copyright © 2010 John Wiley & Sons, Ltd.