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.     

Seasonal variations in infiltrability of moss‐dominated biocrusts on aeolian sand and loess soil in the Chinese Loess Plateau

Biocrust effects on soil infiltration have attracted increasing attention in dryland ecosystems, but their seasonal variations in infiltrability have not yet been well understood. On the Chinese Loess Plateau, soil infiltrability indicated by saturated hydraulic conductivity (K_s) of biocrusts and bare soil, both on aeolian sand and loess soil, was determined by disc infiltrometer in late spring (SPR), midsummer (SUM), and early fall (FAL). Then their correlations with soil biological and physiochemical properties and water repellency index (RI) were analysed. The results showed that the biocrusts significantly decreased K_s both on sand during SPR, SUM, and FAL (by 43%, 66%, and 35%, respectively; P < .05) and on loess (by 42%, 92%, and 10%, respectively; P <.05). As compared with the bare soil, the decreased K_s in the biocrusted surfaces was mostly attributed to the microorganism biomass and also to the increasing content of fine particles and organic matter. Most importantly, both the biocrusts and bare soil exhibited significant (F ≥ 11.89, P ≤ .003) seasonal variations in K_s, but their patterns were quite different. Specifically, the K_s of bare soil gradually decreased from SPR to SUM (32% and 42% for sand and loess, respectively) and FAL (29% and 39%); the K_s of biocrusts also decreased from SPR to SUM (59% and 92%) but then increased in FAL (36% and 588%). Whereas the seasonal variations in K_s of the biocrusts were closely correlated with the seasonal variations in RI, the RI values were not high enough to point at hydrophobicity. Instead of that, the seasonal variations of K_s were principally explained by the changes in the crust biomass and possibly by the microbial exopolysaccharides. We conclude that the biocrusts significantly decreased soil infiltrability and exhibited a different seasonal variation pattern, which should be carefully considered in future analyses of hydropedological processes.     

Pattern and chronology of glacial Lake Peace shorelines and implications for isostacy and ice‐sheet configuration in northeastern British Columbia,Canada

Recognition of positions of glacial lakes along the margin of continental ice sheets is critical in reconstructing ice configuration during deglaciation. Advances in remote sensing technology (e.g. LiDAR) have enabled the generation of accurate digital‐elevation models (DEMs) that reveal unprecedented geomorphic detail. Combined with geographical information systems, these tools have considerably advanced the mapping and correlation of geomorphic features such as relict shorelines. Shorelines of glacial Lake Peace (GLP) developed between the Laurentide and Cordilleran ice sheets in northeastern British Columbia and northwestern Alberta. Shoreline mapping from high resolution DEMs produced more than 55 500 elevation data points from 3231 shorelines, enabling the identification of four major phases of GLP: Phase I (altitude 960–990 m a.s.l.); Phase II (890–915 m a.s.l.); Phase III (810–865 m a.s.l.); and Phase IV (724–733 m a.s.l.). The timing of Phase II of GLP is estimated by two optical ages of <16.0±2.5 and 14.2±0.5 ka BP. Extensive mapping of the shorelines allows for measuring of glacial isostatic adjustment as ice retreated. Shorelines currently dip to the northeast at around 0.4–0.5 m km^?1. This slope reflects the asynchronous retreat of the Cordilleran (CIS) and Laurentide (LIS) ice sheets. The relative uplift in the southwest of the study area within the Rocky Mountains and foothills suggests that the Late Wisconsinan (MIS 2) CIS persisted in the foothill after the LIS lost mass and retreated, or that the Late Wisconsinan CIS was very thick and caused deep crustal loading, which resulted in more uplift in the southwest before reaching equilibrium during, or shortly after deglaciation.     

Assessing the effectiveness of imperviousness on stormwater runoff in micro urban catchments by model simulation

Imperviousness, considered as a critical indicator of the hydrologic impacts of urbanization, has gained increasing attention both in the research field and in practice. However, the effectiveness of imperviousness on rainfall–runoff dynamics has not been fully determined in a fine spatiotemporal scale. In this study, 69 drainage subareas <1 ha of a typical residential catchment in Beijing were selected to evaluate the hydrologic impacts of imperviousness, under a typical storm event with a 3‐year return period. Two metrics, total impervious area (TIA) and effective impervious area (EIA), were identified to represent the impervious characteristics of the selected subareas. Three runoff variables, total runoff depth (TR), peak runoff depth (PR), and lag time (LT), were simulated by using a validated hydrologic model. Regression analyses were developed to explore the quantitative associations between imperviousness and runoff variables. Then, three scenarios were established to test the applicability of the results in considering the different infiltration conditions. Our results showed that runoff variables are significantly related to imperviousness. However, the hydrologic performances of TIA and EIA were scale dependent. Specifically, with finer spatial scale and the condition heavy rainfall, TIA rather than EIA was found to contribute more to TR and PR. EIA tended to have a greater impact on LT and showed a negative relationship. Moreover, the relative significance of TIA and EIA was maintained under the different infiltration conditions. These findings may provide potential implications for landscape and drainage design in urban areas, which help to mitigate the runoff risk. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrologic and water isotope characterization of a regulated Canadian Shield river basin

Stable water isotope surveys have increasingly been integrated into river basins studies, but fewer have used them to evaluate impact of hydropower regulation. This study applies hydrologic and water isotope survey approaches to a Canadian Shield river basin with both regulated and natural flows. Historical streamflow records were used to evaluate the influence of three hydroelectric reservoirs and unregulated portions of the basin on downstream flows and changes in water level management implemented after an extreme flood year (1979). In 2013, water isotope surveys of surface and source waters (e.g., rainfall, groundwater, snowmelt) were conducted to examine spatial and temporal variation in contributions to river flow. Seasonal changes in relative groundwater contribution were assessed using a water‐isotope mass balance approach. Within the basin, two regulated reservoirs exhibited inverted hydrographs with augmented winter flows, whereas a third exhibited a hydrograph dominated by spring snowmelt. In 2013, spatial variation in rain‐on‐snow and air temperatures resulted in a critical lag in snowmelt initiation in the southern and northern portions of the basin resulting in a dispersed, double peak spring hydrograph, contrasting with 1979 when a combination of rain‐on‐snow and coincident snowmelt led to the highest flood on record. Although eastern basin reservoirs become seasonally enriched in δ¹⁸O and δ²H values, unregulated western basin flows remain less variable due to groundwater driven baseflow with increasing influence downstream. Combined analysis of historical streamflow (e.g., flood of 1979, drought of 2010) and the 2013 water isotope surveys illustrate extreme meteorological conditions that current management activities are unable to prevent. In this study, the influence of evaporative fractionation on large surface water reservoirs provides important evidence of streamflow partitioning, illustrating the value of stable water isotope tracers for study of larger catchments.