Modelling the fine-structure of the geoid: Methods,data requirements and some results

The requirements for precise geoid models on local and regional scales have increased in recent years, primarily due to the ongoing developments in height determination by GPS on land, but also due to oceanographic requirements in using satellite altimetry for recovering dynamic sea-surface topography. Suitable methods for geoid computations from gravity data include Stokes integration, FFT methods, and least-squares collocation. Especially the FFT methods are efficient in handling large amounts of gravity data, and new variants of the methods taking earth curvature rigorously into account provide attractive methods for obtaining continental-scale, high-resolution geoid models. The accuracy of such models may be from 2–5 cm locally, to 50–100 cm on regional scales, depending on gravity data coverage, long wave-length gravity field errors, and datum problems. When approaching the cm-level geoid basic geoid definition questions (geoid or quasigeoid?) become very significant, especially in rugged areas. In the paper the geoid modelling methods and problems are reviewed, and some investigations on local data requirements for cm-level geoid prediction are presented. Some actual results are presented from Scandinavia, where a recent regional high-resolution geoid model yields apparent accuracies of 2–10 cm over GPS baselines of 50 to 2000 km.     

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.     

Sample return of interstellar matter (SARIM)

The scientific community has expressed strong interest to re-fly Stardust-like missions with improved instrumentation. We propose a new mission concept, SARIM, that collects interstellar and interplanetary dust particles and returns them to Earth. SARIM is optimised for the collection and discrimination of interstellar dust grains. Improved active dust collectors on-board allow us to perform in-situ determination of individual dust impacts and their impact location. This will provide important constraints for subsequent laboratory analysis. The SARIM spacecraft will be placed at the L2 libration point of the Sun–Earth system, outside the Earth’s debris belts and inside the solar-wind charging environment. SARIM is three-axes stabilised and collects interstellar grains between July and October when the relative encounter speeds with interstellar dust grains are lowest (4 to 20 km/s). During a 3-year dust collection period several hundred interstellar and several thousand interplanetary grains will be collected by a total sensitive area of 1 m². At the end of the collection phase seven collector modules are stored and sealed in a MIRKA-type sample return capsule. SARIM will return the capsule containing the stardust to Earth to allow for an extraction and investigation of interstellar samples by latest laboratory technologies.     

A statistical test on the local effects of spatially structured variance

Spatial variance is an important characteristic of spatial random variables. It describes local deviations from average global conditions and is thus a proxy for spatial heterogeneity. Investigating instability in spatial variance is a useful way of detecting spatial boundaries, analysing the internal structure of spatial clusters and revealing simultaneously acting geographic phenomena. Recently, a corresponding test statistic called ‘Local Spatial Heteroscedasticity’ (LOSH) has been proposed. This test allows locally heterogeneous regions to be mapped and investigated by comparing them with the global average mean deviation in a data set. While this test is useful in stationary conditions, its value is limited in a global heterogeneous state. There is a risk that local structures might be overlooked and wrong inferences drawn. In this paper, we introduce a test that takes account of global spatial heterogeneity in assessing local spatial effects. The proposed measure, which we call ‘Local Spatial Dispersion’ (LSD), adapts LOSH to local conditions by omitting global information beyond the range of the local neighbourhood and by keeping the related inferential procedure at a local level. Thereby, the local neighbourhoods might be small and cause small-sample issues. In the view of this, we recommend an empirical Bayesian technique to increase the data that is available for resampling by employing empirical prior knowledge. The usefulness of this approach is demonstrated by applying it to a Light Detection and Ranging-derived data set with height differences and by making a comparison with LOSH. Our results show that LSD is uncorrelated with non-spatial variance as well as local spatial autocorrelation. It thus discloses patterns that would be missed by LOSH or indicators of spatial autocorrelation. Furthermore, the empirical outcomes suggest that interpreting LOSH and LSD together is of greater value than interpreting each of the measures individually. In the given example, local interactions can be statistically detected between variance and spatial patterns in the presence of global structuring, and thus reveal details that might otherwise be overlooked.     

Assessment of long-range kinematic GPS positioning errors by comparison with airborne laser altimetry and satellite altimetry

Long-range airborne laser altimetry and laser scanning (LIDAR) or airborne gravity surveys in, for example, polar or oceanic areas require airborne kinematic GPS baselines of many hundreds of kilometers in length. In such instances, with the complications of ionospheric biases, it can be a real challenge for traditional differential kinematic GPS software to obtain reasonable solutions. In this paper, we will describe attempts to validate an implementation of the precise point positioning (PPP) technique on an aircraft without the use of a local GPS reference station. We will compare PPP solutions with other conventional GPS solutions, as well as with independent data by comparison of airborne laser data with “ground truth” heights. The comparisons involve two flights: A July 5, 2003, airborne laser flight line across the North Atlantic from Iceland to Scotland, and a May 24, 2004, flight in an area of the Arctic Ocean north of Greenland, near-coincident in time and space with the ICESat satellite laser altimeter. Both of these flights were more than 800 km long. Comparisons between different GPS methods and four different software packages do not suggest a clear preference for any one, with the heights generally showing decimeter-level agreement. For the comparison with the independent ICESat- and LIDAR-derived “ground truth” of ocean or sea-ice heights, the statistics of comparison show a typical fit of around 10 cm RMS in the North Atlantic, and 30 cm in the sea-ice region north of Greenland. Part of the latter 30 cm error is likely due to errors in the airborne LIDAR measurement and calibration, as well as errors in the “ground truth” ocean surfaces due to drifting sea-ice. Nevertheless, the potential of the PPP method for generating 10 cm level kinematic height positioning over long baselines is illustrated.     

Hydrologic Dynamics of a Subtropical Estuary Using Geochemical Tracers,Celestún,Yucatan, Mexico

Oxygen isotopes and strontium concentrations were used as geochemical tracers to discern the sources of water to Celestún Lagoon, a small subtropical estuary on the western side of the Yucatán Peninsula of Mexico. Celestún Lagoon is underlain by karstified limestone with numerous locations where groundwater is observed discharging directly to the lagoon. In this study, samples of groundwater, lagoon surface water, and seawater (SW) were collected in April 2008 and June 2009 and analyzed for salinity, stable isotopes of oxygen, and strontium (Sr²⁺) concentrations. These geochemical tracers were used in two tertiary mixing models to calculate the relative ratio inputs of fresh groundwater, brackish groundwater, and SW to the lagoon. Two sources of groundwater were found to contribute to the surface water in the lagoon; one fresh and the other brackish with an average salinity of 19 psu. The fresh groundwater had an oxygen isotopic signature (δ¹⁸O) and strontium concentration (Sr²⁺) of δ¹⁸O?=?-3.30‰ and Sr²⁺?=?0.03 mmol/l, respectively. The brackish groundwater observed in the northern end of the lagoon add a dissimilar oxygen isotopic signature and Sr²⁺ concentration of δ¹⁸O?=?3.01‰ and Sr²⁺?=?0.12 mmol/l, respectively. Local SW had an isotopic oxygen signature and Sr²⁺ concentration between the two fresher sources (δ¹⁸O?=?1.40‰, Sr²⁺?=?0.09 mmol/l). The lagoonwide results of the two tracer mixing models (δ¹⁸O and Sr²⁺) agreed well (within 5 %) and indicated a ratio of brackish groundwater–fresh groundwater– SW of 31 %–26 %–43 % (±5 %) for the Sr²⁺ model and 35 %–25 %–40 % (± 5 %) for the δ¹⁸O model. Brackish groundwater is dominant in the northern portion of the lagoon, while SW dominates the southern portion. Fresh groundwater discharge is a significant contributor of water along the entire eastern boundary of the lagoon where mangrove forests are the dominant vegetation.     

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.     

Normal moveout in focus1

Improving the accuracy of NMO corrections and of the corresponding interval velocities entails implementing a better approximation than the formula used since the beginning of seismic processing. The exact equations are not practical as they include many unknowns. The approximate expression has only two unknowns, the reflection time and the rms velocity, but becomes inaccurate for large apertures of the recording system and heterogeneous vertical velocities. Several methods of improving the accuracy have been considered, but the gains do not compensate for the dramatic increase in computing time. Two alternative equations are proposed: the first containing two parameters, the reflection time and the focusing time, is not valid for apertures much greater than is the standard formula, but has a much faster computing time and does not stretch the far traces; the other, containing three parameters, the reflection time, like focusing time and the tuning velocity, retains high frequencies for apertures about twice those allowed by the standard equation. Its computing time can be kept within the same limits. NMO equations, old and new, are designed strictly for horizontal layering, but remain reliable as long as the rays travel through the same layers in both the down and up directions. An equation, similar to Dix's formula, is given to compute the interval velocities. The entire scheme can be automated to produce interval-velocity sections without manual picking.     

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.