Geology ReportsSearch

SEARCH · Geology Reports

Results for “International Journal of Fracture”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

47 records · Page 3Linked to original sources

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International

Combining ash analyses with remote sensing to identify juvenile magma involvement and fragmentation mechanisms during the 2018/19 small eruption of Peteroa volcano (Southern Andes)

The Planchón Peteroa Volcanic Complex (PPVC) is located on the border of Chile and Argentina, and is one of the most active volcanic systems in the Andes. Holocene activity has included magma-water interaction with an evolving series of crater lakes, mainly sourced from Peteroa volcano. This study examines data from the 2018/19 eruption, together with the volcanic history of the PPVC, to elucidate the complex interplay between magmatic activity and summit water and ice. From February 2016 to mid-2019, three seismic swarms occurred in the PPVC, preceding the explosive eruption from September 2018 to April 2019. The activity originated from a small vent nested within the easternmost crater, the most active portion of the complex (Peteroa). The explosions interacted with a crater lake, producing ash plumes up to 2 km above the crater and building a small tephra cone. To investigate the eruption mechanisms, we performed remote sensing analysis of plume dispersal, thermal anomalies and ground deformation, and characterized the volcanic products, including grain size, componentry, morphology, internal textures, composition and mineralogy. Our results suggest that the precursory seismicity beginning in 2016 was related to the intrusion of a new magma batch that reached the surface during the 2018/19 eruption. The eruption was also preceded by thermal anomalies, geomorphic changes and increased hydrothermal activity at the surface, though without any ground deformation recognized through radar interferometry (InSAR). The eruption initially produced predominantly recycled ash (phreatic activity), then evolved to increasing proportions of juvenile magma (phreatomagmatic) by April 2019. The juvenile clasts had a trachyandesite composition (~59 wt% SiO 2 ), with vesicular and dense scoria containing plagioclase and pyroxene. The ash surfaces show external quenching cracks and step fractures consistent with phreatomagmatic fragmentation within the active crater lake. Textural characteristics also point to a slowly ascending batch of magma that was relatively viscous by the time it interacted with water in the crater lake. Notably, these juvenile particles are distinctive from the pre-2018 products. Ash erupted from 2010/11 did not contain recognizable juvenile material, and is inferred to have been a mainly phreatic eruption. Our findings suggest that the interplay between phreatic and phreatomagmatic eruptions fed by small magma batches intruding at shallow levels characterize much of the eruptive behavior of the PPVC during the last three decades. Multi-parametric assessment is a powerful tool to discriminate between phreatic and phreatomagmatic eruptions.

Southern Andes

Probabilistic 3-D time-lapse inversion of magnetotelluric data: Application to an enhanced geothermal system

Surface-based monitoring of mass transfer caused by injections and extractions in deep boreholes is crucial to maximize oil, gas and geothermal production. Inductive electromagnetic methods, such as magnetotellurics, are appealing for these applications due to their large penetration depths and sensitivity to changes in fluid conductivity and fracture connectivity. In this work, we propose a 3-D Markov chain Monte Carlo inversion of time-lapse magnetotelluric data to image mass transfer following a saline fluid injection. The inversion estimates the posterior probability density function of the resulting plume, and thereby quantifies model uncertainty. To decrease computation times, we base the parametrization on a reduced Legendre moment decomposition of the plume. A synthetic test shows that our methodology is effective when the electrical resistivity structure prior to the injection is well known. The centre of mass and spread of the plume are well retrieved.We then apply our inversion strategy to an injection experiment in an enhanced geothermal system at Paralana, South Australia, and compare it to a 3-D deterministic time-lapse inversion. The latter retrieves resistivity changes that are more shallow than the actual injection interval, whereas the probabilistic inversion retrieves plumes that are located at the correct depths and oriented in a preferential north-south direction. To explain the time-lapse data, the inversion requires unrealistically large resistivity changes with respect to the base model. We suggest that this is partly explained by unaccounted subsurface heterogeneities in the base model from which time-lapse changes are inferred.

Geophysical Journal International

Geometric controls on megathrust earthquakes

The role of subduction zone geometry in the nucleation and propagation of great-sized earthquake ruptures is an important topic for earthquake hazard, since knowing how big an earthquake can be on a given fault is fundamentally important. Past studies have shown subducting bathymetric features (e.g. ridges, fracture zones, seamount chains) may arrest a propagating rupture. Other studies have correlated the occurrence of great-sized earthquakes with flat megathrusts and homogenous stresses over large distances. It remains unclear, however, how subduction zone geometry and the potential for great-sized earthquakes ( M 8+) are quantifiably linked—or indeed whether they can be. Here, we examine the potential role of subduction zone geometry in limiting earthquake rupture by mapping the planarity of seismogenic zones in the Slab2 subduction zone geometry database. We build from the observation that historical great-sized earthquakes have preferentially occurred where the surrounding megathrust is broadly planar, and we use this relationship to search for geometrically similar features elsewhere in subduction zones worldwide. Assuming geometry exerts a primary control on earthquake propagation and termination, we estimate the potential size distribution of large ( M 7+) earthquakes and the maximum earthquake magnitude along global subduction faults based on geometrical features alone. Our results suggest that most subduction zones are capable of hosting great-sized earthquakes over much of their area. Many bathymetric features previously identified as barriers are indistinguishable from the surrounding megathrust from the perspective of slab curvature, meaning that they either do not play an important role in arresting earthquake rupture or that their influence on slab geometry at depth is not resolvable at the spatial scale of our subduction zone geometry models.

Geophysical Journal International

Lithology and evolution of the crust-mantle boundary region in the southwestern Basin and Range Province

Mantle and crustal xenoliths from volcanic rocks in the southwestern Basin and Range province and Colorado Plateau Transition Zone reveal histories of episodic magmatism and deformation that have profoundly influenced the crustal structure of this region. Seismic transects in this area show a strongly reflective Moho of generally low relief, which, in the area of modern transects, consists of a thin zone (<2 km thick) of short reflectors. The upper mantle is transparent and has a P n of 7.8–8.0 km/s similar to much of the western United States. A lower crustal zone, 2–13 km thick, has variable internal reflectivity and a relatively low velocity of 6.6–6.8 km/s. Upper mantle peridotite xenoliths show both ductile and brittle deformational features and have structures and compositions affected by magmatic intrusion; intrusions form complex dike systems and extensive zones of grain boundary infiltration in peridotite xenoliths. Whereas melt infiltration preceded and followed ductile deformation, brittle deformation, represented by closely spaced joint systems and faults, followed ductile deformation and is related to the youngest magmatic episodes. These structural characteristics and high uppermost mantle temperature (∼1000°C) may combine to explain the relatively low P n . Alternating layers of ductily deformed and undeformed peridotites, with or without igneous intrusions, may contribute to the reflectivity of the Moho. Lower crustal xenoliths are dominantly igneous-textured pyroxenites and mafic to intermediate gabbros identical to the dikes in peridotite xenoliths. The crustal xenoliths also commonly are jointed, and in addition many show partial melting and have abundant cavities that probably were filled with CO 2 -rich fluids. These rocks are interpreted as products of underplated magmas that were fed through the mantle dike systems and may represent the lowest crustal unit identified in the seismic records. The mafic compositions and high densities of the crustal xenoliths indicate that the low velocity of the lower crust may be caused in part by fracture systems, partial melts, and high temperatures. Garnet granulite xenoliths from a locality with no mantle peridotite xenoliths probably represent crust of the region before late Miocene extension. Felsic granulite xenoliths from two localities have velocities like those of the two lower crustal units identified seismically and could be present in the modern crust as unequilibrated remnants of old crust. The preferred model for the evolution of the lower lithosphere is one in which extension affects the upper mantle as well as the crust and is overlapped in time by multiple magmatic episodes. The earliest magmatic events preceded extension, and later events accompanied and followed extension.

Arizona, California, Nevada

Insights on the characteristics and sources of gas from an underground coal mine using compositional data analysis

Coal mine gas originates from the gas emission zone (GEZ) of the mine, as well as the longwall face and pillars. Gas emissions are controlled directly at the sources using horizontal or vertical boreholes drilled from surface or from the entries in advance of mining, or it is captured from the fractured and caved zones (gob) using ventholes during mining. The rest of the gas, especially that gas that originates from the longwall face and caved zone, mixes with the ventilation air and travels through bleeder and return entries before being exhausted to atmosphere from air shafts. Although the gas associated with mining mostly focuses on methane, the gas is not pure methane but is a mixture (where both the components and their quantities vary depending on the sampling location). Understanding the evolution of the composition of the gas from source to different sampling and evaluation points in mines using proper statistical analysis and interpretation methods can lead to better designed degasification and ventilation systems and the selection of the most adequate utilization method for generating energy from the gas. In this work, we present the results of compositional data analysis (CoDa) of gases sampled at different locations in a longwall mine operating in Pennsylvanian coal-bearing strata in the Northern Appalachian coal basin. Sampling locations were from accessible parts of the bleeder entries, returns, bleeder evaluation points and shafts within the mine, and also from gob gas ventholes (GGV) and coal degasification boreholes drilled in the panel areas. In addition, desorbed gas samples from seams that are important for the mine were included in the analyses for comparison. The compositional data analysis showed that the gas composition shifts based on the sampling location, putting in-mine and pure coalbed gases on the opposite ends. Removal of air from the samples did not change this observation suggesting that oxygen is depleted especially for in-mine samples, which is due to oxidation. Results also suggested that desorbed gas samples may not represent the composition of the coal seam gas.

International Journal of Coal Geology

Seismic tomography shows that upwelling beneath Iceland is confined to the upper mantle

We report the results of the highest-resolution teleseismic tomography study yet performed of the upper mantle beneath Iceland. The experiment used data gathered by the Iceland Hotspot Project, which operated a 35-station network of continuously recording, digital, broad-band seismometers over all of Iceland 1996-1998. The structure of the upper mantle was determined using the ACH damped least-squares method and involved 42 stations, 3159 P-wave, and 1338 S-wave arrival times, including the phases P, pP, sP, PP, SP, PcP, PKIKP, pPKIKP, S, sS, SS, SKS and Sdiff. Artefacts, both perceptual and parametric, were minimized by well-tested smoothing techniques involving layer thinning and offset-and-averaging. Resolution is good beneath most of Iceland from ??? 60 km depth to a maximum of ??? 450 km depth and beneath the Tjornes Fracture Zone and near-shore parts of the Reykjanes ridge. The results reveal a coherent, negative wave-speed anomaly with a diameter of 200-250 km and anomalies in P-wave speed, Vp, as strong as -2.7 per cent and in S-wave speed, Vs, as strong as -4.9 per cent. The anomaly extends from the surface to the limit of good resolution at ??? 450 km depth. In the upper ??? 250 km it is centred beneath the eastern part of the Middle Volcanic Zone, coincident with the centre of the ??? 100 mGal Bouguer gravity low over Iceland, and a lower crustal low-velocity zone identified by receiver functions. This is probably the true centre of the Iceland hotspot. In the upper ??? 200 km, the low-wave-speed body extends along the Reykjanes ridge but is sharply truncated beneath the Tjornes Fracture Zone. This suggests that material may flow unimpeded along the Reykjanes ridge from beneath Iceland but is blocked beneath the Tjornes Fracture Zone. The magnitudes of the Vp, Vs and Vp/Vs anomalies cannot be explained by elevated temperature alone, but favour a model of maximum temperature anomalies <200 K, along with up to ??? 2 per cent of partial melt in the depth range ??? 100-300 km beneath east-central Iceland. The anomalous body is approximately cylindrical in the top 250 km but tabular in shape at greater depth, elongated north-south and generally underlying the spreading plate boundary. Such a morphological change and its relationship to surface rift zones are predicted to occur in convective upwellings driven by basal heating, passive upwelling in response to plate separation and lateral temperature gradients. Although we cannot resolve structure deeper than ??? 450 km, and do not detect a bottom to the anomaly, these models suggest that it extends no deeper than the mantle transition zone. Such models thus suggest a shallow origin for the Iceland hotspot rather than a deep mantle plume, and imply that the hotspot has been located on the spreading ridge in the centre of the north Atlantic for its entire history, and is not fixed relative to other Atlantic hotspots. The results are consistent with recent, regional full-thickness mantle tomography and whole-mantle tomography images that show a strong, low-wave-speed anomaly beneath the Iceland region that is confined to the upper mantle and thus do not require a plume in the lower mantle. Seismic and geochemical observations that are interpreted as indicating a lower mantle, or core-mantle boundary origin for the North Atlantic Igneous Province and the Iceland hotspot should be re-examined to consider whether they are consistent with upper mantle processes.

Geophysical Journal International

Landsat time series analysis of fractional plant cover changes on abandoned energy development sites

Oil and natural gas development in the western United States has increased substantially in recent decades as technological advances like horizontal drilling and hydraulic fracturing have made extraction more commercially viable. Oil and gas pads are often developed for production, and then capped, reclaimed, and left to recover when no longer productive. Understanding the rates, controls, and degree of recovery of these reclaimed well sites to a state similar to pre-development conditions is critical for energy development and land management decision processes. Here we use a multi-decadal time series of satellite imagery (Landsat 5, 1984–2011) to assess vegetation regrowth on 365 abandoned well pads located across the Colorado Plateau in Utah, Colorado, and New Mexico. We developed high-frequency time series of the Soil-Adjusted Total Vegetation Index (SATVI) for each well pad using the Google Earth Engine cloud computing platform. BFAST time-series models were used to fit temporal trends, identifying when vegetation was cleared from the site and the magnitudes and rates of vegetation change after abandonment. The time series metrics are used to calculate the relative fractional vegetation cover (RFVC) of each pad, a measure of post-abandonment vegetation cover relative to pre-drilling condition. Mean and median RFVC were 36% (s.d. 33%) and 26%, respectively, five years after abandonment, with one third of well pads having RFVC greater than 50%. Statistical analyses suggest that much of the high vegetation cover is associated with weedy invasive annual species such as cheatgrass (Bromus tectorum) and Russian thistle (Salsola spp.). Climate conditions and the year of abandonment also play a role, with increased cover in later years associated with a wetter period. Non-linear change at many pads suggests longer recovery times than would be estimated by linear extrapolation. New techniques implemented here address a complex response of cover change to soils, management, and climate over time, and can be extended to the operational monitoring of energy development across large areas.

International Journal of Applied Earth Observation

Geologic mapping of Europa

Galileo data enable the major geological units, structures, and surface features to be identified on Europa. These include five primary units (plains, chaos, band, ridge, and crater materials) and their subunits, along with various tectonic structures such as faults. Plains units are the most widespread. Ridged plains material spans a wide range of geological ages, including the oldest recognizable features on Europa, and appears to represent a style of tectonic resurfacing, rather than cryovolcanism. Smooth plains material typically embays other terrains and units, possibly as a type of fluid emplacement, and is among the youngest material units observed. At global scales, plains are typically mapped as undifferentiated plains material, although in some areas differences can be discerned in the near infrared which might be related to differences in ice grain size. Chaos material is composed of plains and other preexisting materials that have been severely disrupted by inferred internal activity; chaos is characterized by blocks of icy material set in a hummocky matrix. Band material is arrayed in linear, curvilinear, wedge-shaped, or cuspate zones with contrasting albedo and surface textures with respect to the surrounding terrain. Bilateral symmetry observed in some bands and the relationships with the surrounding units suggest that band material forms by the lithosphere fracturing, spreading apart, and infilling with material derived from the subsurface. Ridge material is mapped as a unit on local and some regional maps but shown with symbols at global scales. Ridge material includes single ridges, doublet ridges, and ridge complexes. Ridge materials are considered to represent tectonic processes, possibly accompanied by the extrusion or intrusion of subsurface materials, such as diapirs. The tectonic processes might be related to tidal flexing of the icy lithosphere on diurnal or longer timescales. Crater materials include various interior (smooth central, rough inner, and annular massif) and exterior (continuous ejecta) subunits. Structural features and landforms are shown with conventional symbols. Type localities for the units are identified, along with suggestions for portraying the features on geological maps, including colors and letter abbreviations for material units. Implementing these suggestions by the planetary mapping community would facilitate comparisons of maps for different parts of Europa and contribute to an eventual global synthesis of its complex geology. On the basis of initial mapping results, a stratigraphic sequence is suggested in which ridged plains form the oldest unit on Europa, followed by development of band material and individual ridges. Band materials tend to be somewhat older than ridges, but in many areas the two units formed simultaneously. Similarly, the formation of most chaos follows the development of ridged plains; although chaos is among the youngest materials on Europa, some chaos units might have formed contemporaneously with ridged plains. Smooth plains generally embay all other units and are late-stage in the evolution of the surface. C1 craters are superposed on ridged plains but are crosscut by other materials, including bands and ridges. Most c2 craters postdate all other units, but a few c2 craters are cut by ridge material. C3 craters constitute the youngest recognizable material on Europa. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Computer simulations of large asteroid impacts into oceanic and continental sites--Preliminary results on atmospheric, cratering and ejecta dynamics

Computer simulations have been completed that describe passage of a 10-km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics caused by impact of the asteroid into both oceanic and continental sites. The asteroid was modeled as a spherical body moving vertically at 20 km/s with a kinetic energy of 2.6 × 10 30 ergs (6.2 × 10 7 Mt ). Detailed material modeling of the asteroid, ocean, crustal units, sedimentary unit, and mantle included effects of strength and fracturing, generic asteroid and rock properties, porosity, saturation, lithostatic stresses, and geothermal contributions, each selected to simulate impact and geologic conditions that were as realistic as possible. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock wave followed by a highly shock compressed and heated air mass. Rapid expansion of this shocked air created a large low-density region that also expanded away from the impact area. Shock temperatures in air reached ∼20, 000K near the surface of the uplifting crater rim and were as high as ∼2000K at more than 30 km range and 10 km altitude. Calculations to 30 s showed that the shock fronts in the air and in most of the expanding shocked air mass preceded the formation of the crater, ejecta, and rim uplift and did not interact with them. As cratering developed, uplifted rim and target material were ejected into the very low density, shock-heated air immediately above the forming crater, and complex interactions could be expected. Calculations of the impact events showed equally dramatic effects on the oceanic and continental targets through an interval of 120 s. Despite geologic differences in the targets, both cratering events developed comparable dynamic flow fields and by ∼29s had formed similar-sized transient craters ∼39km deep and ∼62km across. Transient-rim uplift of ocean and crust reached a maximum altitude of nearly 40 km at ∼30s and began to decay at velocities of 500 m/s to develop large-tsunami conditions. After ∼30s, strong gravitational rebound drove both craters toward broad flat-floored shapes. At 120 s, transient crater diameters were ∼80km (continental) and ∼105km (oceanic) and transient depths were ∼27km; crater floors consisting of melted and fragmented hot rock were rebounding rapidly upward. By 60 s, the continental crater had ejected ∼2 × 10 14 t, about twice the mass ejected from the oceanic crater. By 120 s, ∼70, 000km 3 (continental) and ∼90, 000km 3 (oceanic) target material were excavated (no mantle) and massive ejecta blankets were formed around the craters. We estimate that in excess of ∼70% of the ejecta would finally lie within ∼3 crater diameters of the impact, and the remaining ejecta (∼10 13 t), including the vaporized asteroid, would be ejected into the atmosphere to altitudes as high as the ionosphere. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to contribute substantial material to the atmosphere.

International Journal of Impact Engineering

India National Gas Hydrate Program Expedition-02: Operational and technical summary

The India National Gas Hydrate Program is being steered by the government of India's Ministry of Petroleum and Natural Gas (MoPNG) with participation of Directorate General of Hydrocarbons (DGH), Oil and Natural Gas Corporation Limited (ONGC), and the National Oil Companies and Research Institutes of India. The India National Gas Hydrate Program Expedition 01 (NGHP-01) established the presence of gas hydrate in the Krishna Godavari (KG) and Mahanadi Basins and in the offshore area of the Andaman Sea Basin. However, the gas hydrates discovered during NGHP-01 were mainly distributed as fracture-filling material in fine-grained clay-rich sediments. The India National Gas Hydrate Program Expedition 02 (NGHP-02) was carried out with an objective to discover gas hydrate in sand-rich sediment along the eastern offshore margin of India. ONGC planned and executed NGHP-02 on the behalf of the MoPNG. NGHP-02 started on March 3, 2015 and was completed on July 28, 2015 (total 147 days) using the Japanese scientific Drilling Vessel Chikyu (D/V Chikyu). During NGHP-02, 42 holes at 25 sites were drilled, cored, and/or surveyed with downhole logging tools. These sites were located in four areas along the eastern margin of India and formally named Area A (Mahanadi Basin, three sites), Area B (northern part of the KG-Basin, twelve sites), Area C (central part of the KG-Basin, six sites), and Area E (southern part to the KG-Basin, four sites). All 25 sites established during NGHP-02 were first drilled and logged with logging-while-drilling (LWD) tools and an additional 17 holes were then drilled and/or cored with conventional coring tools (HPCS/ESCS) or pressure coring tools (PCTB). Wireline logging was conducted in 10 holes and formation tests using a dual packer Modular Formation Dynamics Tester (MDT) tool were carried out in two holes. The onboard science team used the laboratory facilities on the D/V Chikyu to examine and analyse the physical properties, geochemistry, and sedimentology of all the cores collected during the expedition. Core samples were also analysed in additional post-expedition shore-based studies conducted in numerous domestic and international gas hydrate research laboratories. The NGHP-02 sediment cores were archived at the National Gas Hydrate Core Repository in Mumbai (India), which is associated with the ONGC Gas Hydrate Research and Technology Centre (GHRTC). The necessary data for characterizing the occurrence of gas hydrate, such as interstitial water chlorinities, core-derived gas chemistry, physical and sedimentological properties, thermal images of the recovered cores, pressure core and downhole measured logging data (LWD and/or conventional wireline log data), were obtained from most of the drill sites established during NGHP-02. Almost all the drill sites yielded evidence for the occurrence of gas hydrate; however, the inferred in situ concentration of gas hydrate varied substantially from site to site. For the most part, the interpretation of downhole logging data, core thermal images, interstitial water analyses, and pressure core images from the sites established during NGHP-02 indicate that the occurrence of concentrated gas hydrate is mostly associated with coarser grained (sand-rich) sediments. This paper presents the operational and technical summary of NGHP-02. NGHP-02 started on March 3, 2015 and was completed on July 28, 2015 (total 147 days) using the Japanese scientific Drilling Vessel Chikyu (D/V Chikyu). During NGHP-02, 42 holes at 25 sites were drilled, cored, and/or surveyed with downhole logging tools. These sites were located in four areas along the eastern margin of India and formally named Area A (Mahanadi Basin, three sites), Area B (northern part of the KG-Basin, twelve sites), Area C (central part of the KG-Basin, six sites), and Area E (southern part to the KG-Basin, four sites). All 25 sites established during NGHP-02 were first drilled and logged with logging-while-drilling (LWD) to

Journal of Marine and Petroleum Geology