Collision, subduction and accretion events in the Philippines: A synthesis

Abstract The Philippines preserves evidence of the superimposition of tectonic processes in ancient and present‐day collision and subduction zone complexes. The Baguio District in northern Luzon, the Palawan–Central Philippine region and the Mati–Pujada area in southeastern Mindanao resulted from events related to subduction polarity reversal leading to trench initiation, continent‐arc collision and autochthonous oceanic lithosphere emplacement, respectively. Geological data on the Baguio District in Northern Luzon reveal an Early Miocene trench initiation for the east‐dipping Manila Trench. This followed the Late Oligocene cessation of subduction along the west‐dipping proto‐East Luzon Trough. The Manila Trench initiation, which is modeled as a consequence of the counter‐clockwise rotation of Luzon, is attributed to the collision of the Palawan microcontinental block with the Philippine Mobile Belt. In the course of rotation, Luzon onramped the South China Sea crust, effectively converting the shear zone that bounded them into a subduction zone. Several collision‐related accretionary complexes (e.g. Romblon, Mindoro) are present in the Palawan–Central Philippine region. The easternmost collision zone boundary is located east of the Romblon group of islands. The Early Miocene southwestward shift of the collision boundary from Romblon to Mindoro started to end by the Pliocene. Continuous interaction between the Palawan microcontinental block and the Philippine Mobile Belt is presently taken up again along the collisional boundary east of the Romblon group of islands. The Mati–Pujada Peninsula area, on the other hand, is underlain by the Upper Cretaceous Pujada Ophiolite. This supra‐subduction zone ophiolite is capped by chert and pelagic limestones which suggests its derivation from a relatively deep marginal basin. The Pujada Ophiolite could be a part of a proto‐Molucca Sea plate. The re‐interpretation of the geology and tectonic settings of the three areas reaffirm the complex geodynamic evolution of the Philippine archipelago and addresses some of its perceived geological enigmas.     

Assessing traffic noise impact based on probabilistic and fuzzy approaches under uncertainty

Traffic noise can cause severe sound pollution for human communities. This paper proposes a hybrid approach to assess traffic noise impact under uncertainty. There are many factors influencing traffic noise level, but only three traffic parameters, namely, traffic flow, traffic speed and traffic component, are highly uncertain. These uncertain parameters are represented by probability distributions, and Monte Carlo simulations are performed to generate a noise distribution after considering about other certain influencing factors. Fuzzy set and binary fuzzy relations as well as probability analysis method are applied to identify the predicted traffic noise impacts in qualitatively and quantitatively. The applicability of this proposed technique is demonstrated using a case study.     

Slab melt as metasomatic agent in island arc magma mantle sources, Negros and Batan (Philippines)

Abstract Two new cases of association of adakites with ‘normal’ island arc lavas and transitional adakites are recognized in the islands of Batan and Negros in northern and central Philippines, respectively. The Batan lavas are related to the subduction of the middle Miocene portion of the South China Sea basin along the Manila trench; those of Negros come from the almost aseismic subduction of the middle Miocene Sulu Sea crust along the Negros trench. The occurrence of the Batan adakites is consistent with previous findings showing adakitic glass inclusions within minerals of mantle xenoliths associated with Batan arc lavas. The similarity of adakite ages (1.09 Ma) and that of the metasomatized xenoliths (1 Ma) suggests that both are linked to the same slab‐melting and metasomatic event. Earlier Sr, Pb and Nd‐isotopic studies, however, also reveal the presence of an important sediment contribution to the Batan lava geochemistry. Thus, the role played by slab melts, assumed to have mid‐ocean ridge basalts‐like (MORB) isotopic characteristics, in enriching the Batan subarc mantle is largely masked by the sediment input. The Negros adakites are present only in Mount Cuernos, the volcanic center nearest to the Negros trench. Batch partial melting calculations show that the Negros adakites could be derived from a garnet amphibolitic source with normal‐MORB (N‐MORB) geochemistry. This is supported by the MORB‐like isotopic characteristics of the Mount Cuernos lavas. The volcanic rocks from the other volcanoes consist of normal arc and transitional adakitic lavas that have slightly higher Sr‐ and Pb‐isotopic ratios, probably due to slight sediment input. Mixing of adakites and normal arc lavas to produce transitional adakites is only partly supported by trace element geochemistry and not by field evidence. The transitional adakites can be modeled as partial melts of an adakite‐enriched mantle. Trace element enrichment of non‐adakitic lavas could reflect the interaction of their mantle source with uprising slab melts, as metasomatic mantle minerals scavenge certain trace elements from the adakitic fluids. Therefore, in arcs beneath which thick (up to 2 km) continent‐derived detrital sediments are involved in subduction, like in Batan, the sediment signature can overwhelm the slab melt input. In arcs like Negros where slow subduction could cause a more efficient scraping of thinner (approximately 1 km) detrital sediments, the contribution of slab melts is easier to detect.     

Low resistivity zones at contacts of igneous intrusions emplaced in organic-rich formations and their implications on fluid flow and petroleum systems: A case study in the northern Neuquén Basin,Argentina

Igneous sills and laccoliths emplaced in sedimentary basins may significantly impact petroleum systems, both positively and negatively. Igneous intrusions provide heat to maturate regionally immature organic-rich host rocks, act as fractured reservoirs hosting commercial accumulations of hydrocarbons, and form structures affecting fluid flow and trapping at different scales. Nevertheless, the petrophysical implications of igneous intrusions on their host rock are poorly known. In this study, we analyse 200 wells in the Río Grande Valley oil field, Neuquén basin, Argentina, where the main reservoirs are in fractured igneous sills. This dataset represents a globally unique possibility to characterize the igneous–host rock interaction using both wireline logs and core material. We identify a systematic Contact Low Resistivity Zone (CLRsZ) at both the upper and lower contacts of the sills emplaced in the organic-rich Vaca Muerta and Agrio Formations. We characterize the nature of these CLRsZ and their petrophysical properties by integrating resistivity and gamma ray well logs, petrographic analyses, petrophysical tests and geochemical analyses. The low resistivity signal of the CLRsZ is dominantly carried by massive-sulphide deposits, mainly pyrite, observed both in the host rock and the chilled margin of the sills. Well log images and porosity-permeability analysis on core plugs show that both the sills and their associated CLRsZ can act as carrier for fluid flow and reservoir for hydrocarbons storage. The thickness of the upper and lower CLRsZ correlates linearly with the thickness of the sill, and the volume of both the upper and lower CLRsZ represents ca. 40% with respect to the volume of their associated sill. The thickness of the CLRsZ represents ca. 13% of the thickness of contact aureole induced by the sills. In the CLRsZ, a great proportion of kerogen was transformed to hydrocarbon, so that CLRsZ were restricted to the innermost contact aureole of the sills. Our results show that the CLRsZ can have major implications on fluid flow and should be considered in reservoir models in volcanic basins hosting sills emplaced in organic-rich formations.     

Greenland and Antarctica Ice Sheet Mass Changes and Effects on Global Sea Level

Thirteen years of GRACE data provide an excellent picture of the current mass changes of Greenland and Antarctica, with mass loss in the GRACE period 2002–2015 amounting to 265 ± 25 GT/year for Greenland (including peripheral ice caps), and 95 ± 50 GT/year for Antarctica, corresponding to 0.72 and 0.26 mm/year average global sea level change. A significant acceleration in mass loss rate is found, especially for Antarctica, while Greenland mass loss, after a corresponding acceleration period, and a record mass loss in the summer of 2012, has seen a slight decrease in short-term mass loss trend. The yearly mass balance estimates, based on point mass inversion methods, have relatively large errors, both due to uncertainties in the glacial isostatic adjustment processes, especially for Antarctica, leakage from unmodelled ocean mass changes, and (for Greenland) difficulties in separating mass signals from the Greenland ice sheet and the adjacent Canadian ice caps. The limited resolution of GRACE affects the uncertainty of total mass loss to a smaller degree; we illustrate the “real” sources of mass changes by including satellite altimetry elevation change results in a joint inversion with GRACE, showing that mass change occurs primarily associated with major outlet glaciers, as well as a narrow coastal band. For Antarctica, the primary changes are associated with the major outlet glaciers in West Antarctica (Pine Island and Thwaites Glacier systems), as well as on the Antarctic Peninsula, where major glacier accelerations have been observed after the 2002 collapse of the Larsen B Ice Shelf.     

Investigations on Aral Sea Regressions from Mirabilite Deposits and Remote Sensing

Remote sensing techniques including radar (Topex/Poseidon, Jason-1 and Envisat) and laser altimetry (Icesat), and moderate resolution spectro-radiometer (MODIS) images, are used to estimated current level and surface extent time variations of the Aral Sea. During the Holocene several phases of regression occurred, leading to desiccation of the Aral Sea. During the last 50 years, Aral Sea has drastically shrunk due to intense use of river’s water for irrigation purposes. It is currently separated into four distinct water bodies, namely, the Small Aral in the North, the Tchebas Bay in the North West, and the South West and the South East basins. The Kulandy strait connected the SW and SE basins until very recent times. These basins are now almost separated and salinity becomes very high (140–180 g/l) in the Eastern part. Rubanov discovered past deposits of mirabilite in the years 1970–1980. We investigate the significance of these deposits in the light of current evolution of the four water bodies that constitute the heritage of Aral Sea contemporary desiccation. Using remote sensing techniques, we have attempted to calculate the water balance of south Aral Sea during the last 3 years. We conclude in strong probability that the Kulandy strait carries water most of the time from the Eastern Basin to the Western Basin. We have demonstrated that it should have been the same process in the past to explain the Mirabilite deposit, but unfortunately, due to recent artificial water monitoring of the Aral Sea (dam in the Berg’s strait, new reservoirs in the Amu Darya’s delta), it is impossible to make definitive conclusion from actual Aral Sea water balance.     

Evaluation of the controls affecting the quality of spatial data derived from historical aerial photographs

This paper is concerned with the fundamental controls affecting the quality of data derived from historical aerial photographs typically used in geomorphological studies. A short review is provided of error sources introduced into the photogrammetric workflow. Data‐sets from two case‐studies provided a variety of source data and hence a good opportunity to evaluate the influence of the quality of archival material on the accuracy of coordinated points. Based on the statistical weights assigned to the measurements, precision of the data was estimated a priori, while residuals of independent checkpoints provided an a posteriori measure of data accuracy. Systematic discrepancies between the two values indicated that the routinely used stochastic model was incorrect and overoptimistic. Optimized weighting factors appeared significantly larger than previously used (and accepted) values. A test of repeat measurements explained the large uncertainties associated with the use of natural objects for ground control. This showed that the random errors not only appeared to be much larger than values accepted for appropriately controlled and targeted photogrammetric networks, but also small undetected gross errors were induced through the ‘misidentification’ of points. It is suggested that the effects of such ‘misidentifications’ should be reflected in the stochastic model through selection of more realistic weighting factors of both image and ground measurements. Using the optimized weighting factors, the accuracy of derived data can now be more truly estimated, allowing the suitability of the imagery to be judged before purchase and processing. Copyright © 2010 John Wiley & Sons, Ltd.     

Large-scale numerical simulation of groundwater flow and solute transport in discretely-fractured crystalline bedrock

A large-scale fluid flow and solute transport model was developed for the crystalline bedrock at Olkiluoto Island, Finland, which is considered as potential deep geological repository for spent nuclear fuel. Site characterization showed that the main flow pathways in the low-permeability crystalline bedrock on the island are 13 subhorizontal fracture zones. Compared to other sites investigated in the context of deep disposal of spent nuclear fuel, most deep boreholes drilled at Olkiluoto are not packed-off but are instead left open. These open boreholes intersect the main fracture zones and create hydraulic connections between them, thus modifying groundwater flow. The combined impact of fracture zones and open boreholes on groundwater flow is simulated at the scale of the island. The modeling approach couples a geomodel that represents the fracture zones and boreholes with a numerical model that simulates fluid flow and solute transport. The geometry of the fracture zones that are intersected by boreholes is complex, and the 3D geomodel was therefore constructed with a tetrahedral mesh. The geomodel was imported into the numerical model to simulate a pumping test conducted on Olkiluoto Island. The pumping test simulation demonstrates that fracture-borehole intersections must be accurately discretized, because they strongly control groundwater flow. The tetrahedral mesh provides an accurate representation of these intersections. The calibrated flow model was then used for illustrative scenarios of radionuclide migration to show the impact of fracture zones on solute transport once the boreholes were backfilled. These mass transport simulations constitute base cases for future predictive analyses and sensitivity studies, since they represent key processes to take into consideration for repository performance assessment.     

Sedimentary carbon forms in relation to climate and phytoplankton biomass in a large,shallow, hard-water boreal lake

Carbon storage in lakes can have huge implications for the global carbon cycle, as lakes annually accumulate up to one half the amount of organic carbon buried in marine sediments. Yet little is known of the effect of recent climate change on carbon storage in lakes. We analyzed century-scale time series of climate variables (precipitation, temperature, NAO winter index) and profiles of sediment characteristics in a dated sediment core from shallow, eutrophic Lake Võrtsjärv, south Estonia. We used path analysis to evaluate the effect of climate conditions on phytoplankton biomass in the lake and accumulation of organic and inorganic carbon in the sediment. Changes in winter and spring climate influenced the lake’s phytoplankton growth significantly. Carbon pathways in hard-water Lake Võrtsjärv were influenced by both hydrological (most significant in colder periods) and biogeochemical processes. Increased nutrient and water input to Lake Võrtsjärv, anticipated with projected climate warming, favours greater in-lake productivity, larger accumulation of inorganic carbon in sediments, and an increase in organic carbon mineralisation, which fuels atmospheric greenhouse gas emissions from the lake.     

Copper Biosorption by Spent Coffee Ground: Equilibrium,Kinetics, and Mechanism

The potential use of spent coffee ground (SCG) for the removal of copper has been investigated as a low‐cost adsorbent for the biosorption of heavy metals. Adsorption batch experiments were conducted to determine isotherms and kinetics. The biosorption equilibrium data were found to fit well the Freundlich model and an experimental maximum biosorption capacity of copper ions 0.214 mmol/g was achieved. The biosorption kinetics of SCG was studied at different adsorbate concentrations (0.1–1.0 mM) and stirring speeds (100–400/min). The results showed an increase in the copper ion uptake with raising the initial metal concentration and the kinetic data followed the pseudo‐second order rate expression. The effect of stirring speed was a significant factor for the external mass transfer resistance at 100/min and coefficients were estimated by the Mathews and Weber model. Biosorption of copper ions onto SCG was observed to be related mainly with the release of calcium and hydrogen ions suggesting that biosorption performance by SCG can be attributed to ion‐exchange mechanism with calcium and hydrogen ions neutralizing the carboxyl and hydroxyl groups of the biomass.