Geology ReportsSearch

USGS · 70015356

Geology of the Zambales ophiolite, Luzon, Philippines

Abstract

The Zambales ophiolite of western Luzon, Philippines, exposes a typical succession of basalt flows, diabasic dikes, gabbro and tectonized harzburgite. The age established by limiting strata is late Eocene. Lack of evidence of thrust faulting and the general domal disposition of the lithologie units indicate that the ophiolitic rocks are exposed by uplift. Highly complex internal layered structures within the complex are related to processes developed during formation of the ophiolite and the Zambales ophiolite may be one of the least disturbed (by emplacement) ophiolitic masses known. The exposed mass trends north and the upper surface plunges at low angles (a few degrees) to the north and south. The chemistry and composition of the rocks in the northwest part of the Zambales area (Acoje block) is distinct from that in the southeastern segment (Coto block). The Acoje block, according to Evans (1983) and Hawkins and Evans (1983), resembles (on a chemical basis) arc-tholeiite series rocks from intra-island arcs and the rocks in the Coto block are typical back-arc basin rock series. The present writer believes that the ophiolite composes a single genetic unit and that the changes in composition are the result of changes that took place during the initial formation. The gabbro probably formed below a spreading center in an elongate, in cross section, V-shaped, magma chamber. The gabbro is estimated by the writer to be less than 2 km thick and may be less than 1 km in places. Numerous erosional windows through the gabbro in the northern and eastern side of the Zambales area show that the gabbro remaining in those areas is likely to be only a few hundred meters thick. Harzburgite is exposed to a depth of about 800 m in the Bagsit River area and this may be the deepest part of the ophiolite accessible for study on which there is any control on depth. A transitional zone, about 200 m thick lying between the gabbro and harzburgite, is composed of serpentinized dunite. Commonly the dunite contains disseminated sulfide minerals and at the Acoje Mines, platinum-group elements. A compositional layering within the gabbro is in places cumulate in the lower part of the unit but may have formed by nucleation higher up on the relatively steep sides of the magma chamber. A widespread gneissic banding in the gabbro forms large mappable structures which are many times more complex than is the disposition of the major rock units. These structures are believed to be the result of extensive slumping in the magma chamber. The structure produced by the cumulate layering merges with the gneissic banding, commonly without discernible change in attitude. This tectonic layered structure crosses the gabbro-peridotite boundary at any angle without seeming to disturb the original rock distribution. At greater depths below the boundary (ca. 800 m), the harzburgite contains low dipping banding, which probably reflects the result of differential movement within the mantle. Chromite occurs almost exclusively in a zone that generally lies no more than 200–300 m below the gabbro-peridotite boundary. Refractory-grade chromite is found in this zone below the olivine gabbro in the Goto block and as low-grade metallurgical grade chromite below norite in the Acoje block. At Acoje Mines the chromite is present in layers in dunite, which the writer interprets as being distributed in a zone along the gently dipping (ca. 25°) gabbro-peridotite boundary. The steeply dipping (ca. 60–80 ° ) individual layers lie en echelon along the boundary at an angle (ca. 50 ° ) to the contact. At Coto the chromite forms large discontinuous masses in the lowest dunite and in the uppermost harzburgite. Except for the chromite present as layers at Acoje, the regional tectonic layering crosses the chromite deposits without structural deviation. The chromite deposits and associated peridotite may be cumulate in origin, but have been modified to such an extent that cumulate textures are generally obliterated. The angle of repose of cumulate layers in the Acoje area and in the Coto block dip towards each other raising the possibility that the Zambales area may contain the relics of a spreading center. Initial emplacement of the Zambales ophiolite took place by uplift and the ultramafic portion was exposed to erosion in the earliest Miocene or late Oligocene. Submergence of some of the ophiolite followed the previous uplift and on the west side of the Zambales Range submergence of several kilometers is indicated. Final emergence appears to have taken place in Pliocene or Pleistocene time by block uplift and areas of greatest uplift closely conform to the present topographic surface.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 14.71257942771895° to 16.379018641910505° latitude; 120.01436353828899° to 120.64963935647586° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Darwin L. Rossman, G.C. Castanada, G.C. Bacuta. 1989. Geology of the Zambales ophiolite, Luzon, Philippines. https://doi.org/10.1016/0040-1951(89)90366-1

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics

Fairweather transform boundary Oligocene to present orogenesis: Fairweather Range vertical extrusion and rotation of the Yakutat microplate at ca. 3 Ma

Oblique-slip along transform fault boundaries is often partitioned between a strike-slip system and thrust faults that accommodate contraction. However, topography along the Yakutat-North American transform (Fairweather fault), is asymmetric with low-terrain above active thrusts on the western, Yakutat side of the transform and high topography on the continental side with peaks >4500 m (Mount Fairweather: 4671 m) to the west of the Border Ranges fault, limited recorded earthquakes >M4, and no apparent reverse faults to generate the highest terrain. In this study we compile, for the first time, published U-Pb zircon, 40 Ar/ 39 Ar and K-Ar (hornblende, muscovite, and biotite) and U-Th/He and fission-track (zircon and apatite) bedrock ages (109) from 75 samples to investigate the exhumation history of the Fairweather Range region, complemented by a published detrital sample (ZFT and AFT) and 13 new 40 Ar/ 39 Ar (hornblende, biotite, and K-feldspar) ages on 9 bedrock samples from both sides of the Fairweather fault. Additionally, we examined published seismicity and geodetic data of the Fairweather region and assessed if plate paleo-vectors correlate with the cooling history of the Fairweather Range. Cooling age, seismic, and block-motion patterns indicate the Fairweather Range has been vertically extruded between the Fairweather and the Border Ranges faults as a coherent block since ca. 25 Ma. The pre-6 Ma Pacific plate motion (N30°W) aligns with the N33°W strike of the Fairweather Fault whereas a hypothetical pre-6 Ma Yakutat microplate paleo-vector of (N39°W) does not: indicating a post-6 Ma timing for Yakutat microplate counter-clockwise rotation (9°). We infer that rotation and impingement of the Yakutat microplate along the Fairweather fault at ca. 3 Ma led to the development of the Fairweather restraining bend and increased cooling rates. The resultant thickened Fairweather welt and the ∼30 km thick southeast end of the Yakutat microplate compounded double-indenter tectonics into Alaska's southeast convergent corner

Tectonophysics

A scaling relationship for the width of secondary deformation around strike-slip faults

Simple mechanical arguments suggest that slip along interlocked, rough faults, damages surrounding rocks. The same arguments require that the scale of secondary damage is proportional to the size of geometric irregularities along the main fault. This relationship could apply at all scales, but has, so far, been difficult to observe at the 10s to 100 s of km scales of large, natural faults, often because large-scale deformation is distributed across wide, complex plate-boundary fault systems, like the San Andreas Fault. The geometry and geology of another large-scale plate-boundary strike slip fault—the Queen Charlotte Fault (QCF)—is, in contrast, especially simple. Here, we show that observations of secondary deformation are well-aligned with predictions of stress variations caused by geometric irregularities along the QCF, suggesting a geometric relationship between primary fault geometry and secondary deformation. The analytic stress solution reveals that the highest stresses and highest likelihood of failure are confined to a zone of influence (ZOI) with a width quantified by ZOI = λ / 2 π "> ZOI=λ/2π , where λ is the wavelength of geometric variations along the main fault. This simple model is consistent with ∼100-km-scale observations along the QCF and can theoretically be used to predict the width of secondary deformation at all scales.

Queen Charlotte Fault