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.     

Partial melting kinetics of plagioclase-diopside pairs

Partial melting experiments on plagioclase (An₆₀) and diopside have been carried out using pairs of large crystals to investigate textures and kinetics of melting. The experiments were done at one atmosphere pressure as a function of temperature (1,190–1,307° C) and time (1.5–192 h). Melting took place mainly at the plagioclase-diopside contact planes. Reaction zones composed of fine mixtures of calcic plagioclase and melt were developed from the surface of the plagioclase crystal inward. There exists a critical temperature, below which only a few % melting can occur over the duration of the experiments. This sluggish melting is caused by slow NaSi-CaAl diffusion in plagioclase, because the plagioclase crystal must change its composition to produce albite-rich cotectic melts. Diffusion in the solid also affects the chemical composition of the melts. During initial melting, potassium is preferentially extracted from plagioclase because K-Na diffusion in plagioclase is faster than that of NaSi-CaAl. This also causes a shift in the cotectic compositions. Above the critical temperature, on the other hand, melting is promoted by a metastable reaction in which the plagioclase composition does not change, and which produces melts with compositional gradients along the original An₆₀-diopside tie line. The critical temperature is determined by the intersection of the cotectic and the An₆₀-diopside tie line. Interdiffusion coefficients of plagioclase-diopside components in the melt are estimated from melting rates above the critical temperature by using a simplified steady-state diffusion model (e.g., 10^–8 cm²/sec at 1,300° C).Many examples of reaction zones due to partial melting have been described as spongy or fingerprint-like textures in xenoliths. Metastable melting above the critical temperature is considered to take place in natural melting where there is a high degree of melting. However, we cannot exclude the possibility of disequilibrium created by sluggish melting controlled by diffusion in the minerals. If melting occurs close to the solidus, this process can be important even for partial melting in the upper mantle.     

Seasonal variations of the surface currents in the Tsushima Strait

Seasonal variations of the surface currents in the Tsushima Strait were investigated by analyzing the monthly mean surface currents measured with HF radar. Several new features of the surface currents have been found. One notable feature is the large, complicated seasonal variation in the current structure in the eastern channel of the strait. For example, in the southeastern and northwestern regions of the channel, southwestward countercurrents are found in summer while southeastward acrossshore currents are found in autumn and winter. The wind-driven flow (Ekman flow) as well as surface geostrophic currents are responsible for these complicated variations of the surface currents. To quantify each variation of the flow and current, the wind-driven flow was calculated from the monthly wind (more precisely, the friction velocity) using the monthly speed factor and deflection angle estimated in our previous study, and the surface geostrophic currents were then estimated by subtracting the wind-driven flow from the measured surface currents. It was found that the acrossshore currents are the wind-driven flow, and that the surface geostrophic currents flow almost in the along-shore direction, indicating the validity of the decomposition of the surface velocity into the wind-driven flow and the geostrophic currents using the speed factor and deflection angle. A real-vector empirical orthogonal function (EOF) analysis of the surface geostrophic currents shows a pair of eddies in the lee of Tsushima and Iki Islands as the first mode, which indicates that the southwestward countercurrents in the eastern channel are formed primarily by the incoming Tsushima Warm Current.     

Metallogenesis of Porphyry Cu Deposits of the Western Luzon Arc, Philippines: K-Ar ages, SO3 Contents of Microphenocrystic Apatite and Significance of Intrusive Rocks

Abstract: K–Ar ages of the following porphyry Cu deposits in the western Luzon arc are determined: Lobo-Boneng (10.50.4 Ma), Santo Niño (9.50.3 Ma), Black Mountain (2.10.1 Ma), Dizon (2.50.2 Ma) and Taysan (7.30.2 Ma). Microphenocrys-tic apatite in the late Cenozoic intermediate to silicic intrusions associated with porphyry Cu deposits in the western Luzon arc contains sulfur as SO₃ detectable by electron probe microanalyzer. Sulfur is supposed to have been accommodated dominantly as oxidized species in oxidizing hydrous magmas that generated porphyry Cu deposits. Likewise, such high SO₃ contents in microphenocrystic apatite are common characteristics of the intermediate to silicic magmatism of the western Luzon arc, from tonalitic rocks of the Luzon Central Cordillera of about 15 Ma to an active magmatism at Mount Pinatubo. Thus, the western Luzon arc has been generating porphyry Cu mineralization associated with oxidizing hydrous intermediate to silicic magmatism related to eastward subduction, since Miocene to the present day. Intermediate to silicic rocks since 15 Ma to present-day western Luzon arc generally show high whole-rock Sr/Y ratio ranging from 20 to 184. However, porphyry Cu deposit is not necessarily related to the rocks that show higher Sr/Y ratios compared to the other barren rocks in the western Luzon arc. The characteristics of the intermediate to silicic magma associated with porphyry Cu deposit are not attributed to the composition of the source material of the magma, but to the properties defined by the high activity of oxidized species of sulfur in the fluid phase that is encountered during the generation of intermediate to silicic magmas.     

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.     

Dissolution kinetics of plagioclase in the melt of the system diopside-albite-anorthite,and origin of dusty plagioclase in andesites

The textures and kinetics of reaction between plagioclase and melts have been investigated experimentally, and origin of dusty plagioclase in andesites has been discussed. In the experiments plagioclase of different compositions (An₉₆, An₆₁, An₅₄, An₂₃, and An₂₂) surrounded by glasses of six different compositions in the system diopside-albite-anorthite was heated at temperatures ranging from 1,200 to 1,410° C for 30 min to 88 h. Textures were closely related to temperature and chemical compositions. A crystal became smaller and rounded above the plagioclase liquidus temperature of the starting melt (glass) and remained its original euhedral shape below the liquidus. Whatever the temperature, the crystal-melt interface became rough and often more complicated (sieve-like texture composed of plagioclase-melt mixture in the scale of a few m was developed from the surface of the crystal inward; formation of mantled plagioclase) if the crystal is less calcic than the plagioclase in equilibrium with the surrounding melt, and the interface remained smooth if the crystal is more calcic than the equilibrium plagioclase. From these results the following two types of dissolution have been recognized; (1) a crystal simply dissolves in the melt which is undersaturated with respect to the phase (simple dissolution), and a crystal is partially dissolved to form mantled plagioclase by reaction between sodic plagioclase and calcic melt (partial dissolution). The amount of a crystal dissolved and reacted increased proportional to the square root of time. This suggests that these processes are controlled by diffusion, probably in the crystal.Mantled plagioclase produced in the experiments were very similar both texturally and chemically to some of the so-called resorbed plagioclase in igneous rocks. Chemical compositions and textures of plagioclase phenocrysts in island-arc andesites of magma mixing origin have been examined. Cores of clear and dusty plagioclase were clacic (about An₉₀) and sodic (about An₅₀), respectively. This result indicates that dusty plagioclases were formed by the partial melting due to reaction between sodic plagioclase already precipitated in a dacitic magma and a melt of intermediate composition in a mixed magma during the magma mixing.     

Sulfur isotope study of Quaternary volcanic rocks from the Japanese Islands Arc

The sulfide and sulfate contents and their δ³⁴S values were determined in Quaternary volcanic rocks from the Japanese Islands Arc. The total sulfur contents are much lower (less than 40 ppm) and the δ³⁴S values are higher (+4.4 ± 2.1) than those of ocean-floor basalts (800 ± 100 ppm and +0.8 ± 0.5, respectively; Moore and Fabbi, 1971; Sakaiet al., 1982). Lateral variations of both sulfur content and δ³⁴S values were observed in the four volcanic belts in Japan. In the Northeast Japan belt, the sulfur content (30 ± 10 ppm) of the rocks in the inner zone (the Japan Sea side) is 3 to 5 times that in the outer zone (the Pacific side), although the δ³⁴S values of the two zones are almost the same (+4.3 ± 1.0). The δ³⁴S values for the two belts in West Japan are on the average 2%. higher than those of East Japan.This study suggests that the primary magmas that formed the island arc volcanic rocks are initially depleted in sulfur (<120 ppm) and enriched in ³⁴S ($\text{δ}{}^{34}\text{S: +5 ～ +7}$) compared to ocean-floor tholeiitic basalts which formed at mantle under oceanic region. This indicates that the upper-mantle is heterogeneous in sulfur content and isotope composition.     

Double-planed deep seismic zone and anomalous structure in the upper mantle beneath northeastern Honshu (Japan)

Seismicity located by using the most recent data obtained from the high-gain seismograph network of Tohoku University shows that the deep seismic zone beneath northeastern Honshu, Japan, is composed of two thin planes which are parallel to each other and are 30–40 km apart. Focal mechanisms derived from the earthquakes in the upper plane are reverse-faulting, or, some of them, down-dip compression. As a contrast, those in the lower plane are down-dip extension. The location of the upper boundary of the descending lithospheric slab, inferred from the arrival-time difference between ScS and ScSp waves and from the travel-time anomaly of intermediate-depth earthquakes observed at the small-scale seismic array, coincides exactly with the upper plane of the double-planed deep seismic zone. Anelasticity (1/Q) structure of the upper mantle consists of three distinct zones: a high-Q (Qs− 1500) inclined lithospheric slab, an intennediate-Q (Qs−350) land-side mantle between the Pacific coast and the volcanic front, and a low-Q (Qs − 100) land-side mantle between the volcanic front and the coast of the Japan Sea.The evidence obtained here provides valuable information as to the definition of the type of mechanism producing the plate motion beneath island arcs.     

Strong motion uncertainty determined from observed records by dense network in Japan

The variation of ground motions at specific stations from events in six narrow areas was inspected by using K-NET and KiK-net records. A source-area factor for individual observation stations was calculated by averaging ratios between observed values for horizontal peak acceleration and velocity, as well as acceleration response spectra for 5% damping, and predicted values using a ground-motion model (usually known as an attenuation relation) by Kanno et al. (Bull Seismol Soc Am, 96:879–897, 2006). Standard deviations between observed and predicted amplitudes after the correction factor are less than 0.2 on the logarithmic scale and decrease down to around 0.15 in the short-period range. Intra-event standard deviation clearly increases with decreasing distance due to differing paths around near source area. Standard deviations may increase with amplitude or decrease with magnitude; however, both amplitude and magnitude of the data are strongly correlated with distance. The standard deviation calculated in this study is obviously much smaller than that of the original ground-motion model, as epistemic uncertainties are minimized by grouping ground motions at specific stations. This result indicates that the accuracy of strong ground motion prediction could be improved if ground-motion models for specified region are determined individually. For this to be possible, it is necessary to have dense strong-motion networks in high-seismicity regions, such as K-NET and KiK-net.