A note on frame transformations with applications to geodetic datums

Rigorous equations in compact symbolic matrix notation are introduced to transform coordinates and velocities between ITRF frames and modern GPS-based geocentric geodetic datums. The theory is general but, after neglecting higher than second-order terms, it is shown that the equations revert to the formulation currently applied in most major continental datums. We discuss several examples: the North American Datum of 1983 (NAD83), the European Terrestrial Reference System of 1989 (ETRS89), the Geodetic Datum of Australia of 1994 (GDA94), and the South American Geocentric Reference System (SIRGAS). Electronic Publication     

Land use change under conditions of high population pressure: the case of Java

A long history of increases in population pressure in Java has caused agricultural land use to expand and intensify. More recent land use changes caused the conversion of prime agricultural land into residential and industrial area. Results of a dynamic, regional-scale, land use change model are presented, defining the spatial distribution of these land use changes. The model is based on multi-scale modelling of the relations between land use and socio-economic and biophysical determinants. Historical validation showed that the model can adequately simulate the pattern of land use change. Future patterns of land use change between 1994 and 2010 are simulated assuming further urbanization. The results suggest that most intensive land use changes will occur in Java's lowland areas.     

Tsunami Probability in the Caribbean Region

We calculated tsunami runup probability (in excess of 0.5 m) at coastal sites throughout the Caribbean region. We applied a Poissonian probability model because of the variety of uncorrelated tsunami sources in the region. Coastlines were discretized into 20 km by 20 km cells, and the mean tsunami runup rate was determined for each cell. The remarkable ~500-year empirical record compiled by O’Loughlin and Lander (2003) was used to calculate an empirical tsunami probability map, the first of three constructed for this study. However, it is unclear whether the 500-year record is complete, so we conducted a seismic moment-balance exercise using a finite-element model of the Caribbean-North American plate boundaries and the earthquake catalog, and found that moment could be balanced if the seismic coupling coefficient is c = 0.32. Modeled moment release was therefore used to generate synthetic earthquake sequences to calculate 50 tsunami runup scenarios for 500-year periods. We made a second probability map from numerically-calculated runup rates in each cell. Differences between the first two probability maps based on empirical and numerical-modeled rates suggest that each captured different aspects of tsunami generation; the empirical model may be deficient in primary plate-boundary events, whereas numerical model rates lack backarc fault and landslide sources. We thus prepared a third probability map using Bayesian likelihood functions derived from the empirical and numerical rate models and their attendant uncertainty to weight a range of rates at each 20 km by 20 km coastal cell. Our best-estimate map gives a range of 30-year runup probability from 0–30% regionally.     

The geology of kimberlite pipes of the Ekati property, Northwest Territories, Canada

This paper reviews key characteristics of kimberlites on the Ekati property, NWT, Canada. To date 150 kimberlites have been discovered on the property, five of which are mined for diamonds. The kimberlites intrude Archean basement of the central Slave craton. Numerous Proterozoic diabase dykes intrude the area. The Precambrian rocks are overlain by Quaternary glacial sediments. No Phanerozoic rocks are present. However, mudstone xenoliths and disaggregated sediment within the kimberlites indicate that late-Cretaceous and Tertiary cover (likely <200 m) was present at the time of emplacement. The Ekati kimberlites range in age from 45 to 75 Ma. They are mostly small pipe-like bodies (surface area mostly <3 ha but up to 20 ha) that typically extend to projected depths of 400–600 m below current surface. Pipe morphologies are strongly controlled by joints and faults. The kimberlites consist primarily of variably bedded volcaniclastic kimberlite (VK). This is dominated by juvenile constituents (olivine and lesser kimberlitic ash) and variable amounts of exotic sediment (primarily mud), with minor amounts of xenolithic wall-rock material (generally <5%). Kimberlite types include: mud-rich resedimented VK (mRVK); olivine-rich VK (oVK); sedimentary kimberlite; primary VK (PVK); tuffisitic kimberlite (TK) and magmatic kimberlite (MK). The presence and arrangement of these rock types varies widely. The majority of bodies are dominated by oVK and mRVK, but PVK is prominent in the lower portions of certain kimberlites. TK is rare. MK occurs primarily as precursor dykes but, in a few cases, forms pipe-filling intrusions. The internal geology of the kimberlites ranges from simple single-phase pipes (RVK or MK), to complex bodies with multiple, distinct units of VK. The latter include pipes infilled with steep, irregular VK blocks/wedges and at least one case in which the pipe is occupied by well-defined sub-horizontal VK phases, including a unique, 100-m-thick graded sequence. The whole-rock compositions of VK samples suggest significant loss of kimberlitic fines during eruption followed by variable dilution by surface sediment and concurrent incorporation of kimberlitic ash. Diamond distribution within the kimberlites reflects the amount and nature of mantle material sampled by individual kimberlite phases, but is modified considerably by eruption and depositional processes. The characteristics of the Ekati kimberlites are consistent with a two-stage emplacement process: (1) explosive eruption/s causing vent clearing followed by formation of a significant tephra rim/cone of highly fragmented, olivine-enriched juvenile material with varying amounts of kimberlitic ash and surface sediments (predominantly mud); and (2) infilling of the vent by direct deposition from the eruption column and/or resedimentation of crater rim materials. The presence of less fragmented, juvenile-rich PVK in the lower portions of certain pipes and the intrusion of large volumes of MK to shallow levels in some bodies suggest emplacement of relatively volatile-depleted, less explosive kimberlite in the later stages of pipe formation and/or filling. Explosive devolatilisation of CO₂-rich kimberlite magma is interpreted to have been the dominant eruption mechanism, but phreatomagmatism is thought to have played a role and, in certain cases, may have been dominant.     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

Calibration of cyclic constitutive models by oscillating functions

Most structures are subjected to more cyclic loads during their life time than static loads. These cyclic action could be a result of either natural or man-made activities and may lead to soil failure. In order to understand the response of the foundation and its interaction with these complex cyclic loadings, various researchers have over the years developed different constitutive models. Although a lot of research is being carried out on these relatively new models, little or no details exist in literature about the model-based identification of the cyclic constitutive parameters which to a large extent govern the quality of the model output. This could be attributed to the difficulties and complexities of the inverse modeling of such complex phenomena. A variety of optimisation strategies are available for the solution of the sum of least-squares problems as usually done in the field of model calibration. However, for the back analysis (calibration) of the soil response to oscillatory load functions, this article gives insight into the model calibration challenges and also puts forward a method for the inverse modeling of cyclic loaded foundation response such that high-quality solutions are obtained with minimum computational effort.     

Detection of iodine monoxide radicals in the marine boundary layer using laser induced fluorescence spectroscopy

A Laser Induced Fluorescence (LIF) instrument has been developed to detect iodine monoxide (IO) radicals in the atmosphere. An all solid-state Nd:YAG pumped Ti:Sapphire laser operating at approximately 445 nm was used to excite the (2,0) band of the IO A²Π_3/2 ← X²Π_3/2 electronic transition, with off-resonance fluorescence in the (2,5) band detected at 521 nm. The sensitivity of the instrument was determined by calibration. IO (between 10 and 150 pptV) was generated following the 184.9 nm photolysis of N₂O/CF₃I/N₂ mixtures with O₃ actinometry used to determine the photolysis flux. The detection limit was determined to be 0.3 pptV for a 300 s integration period, with an uncertainty of 23% (1σ). The instrument was deployed in August/September 2006 during the RHaMBLe (Reactive Halogens in the Marine Boundary Layer) campaign in Roscoff, France. Located on a small jetty, a few metres from the water’s edge at high tide, the instrument measured significant levels of IO on 11 days, with a maximum of 27.6 ± 3.2 pptV observed on one day (averaged over 10 s) representing the highest IO mixing ratio recorded in the marine boundary layer to date. IO displayed a clear diurnal profile with a maximum at low tide during the daytime. These results represent the first point measurements of IO in the atmosphere by LIF.     

Crustal contaminants in the Permian Oslo Rift, South Norway: constraints from Precambrian geochemistry

During the late Paleozoic Oslo rifting event, the SW part of the Baltic Shield was penetrated by mantle-derived magmas from a depleted lithospheric or sublithospheric source. Along the way to their final emplacement, these magmas may have interacted with a heterogeneous continental crust, consisting of a mosaic of continental terranes, each with its unique composition and internal crustal history. Information on radiogenic isotope ratios and trace element distributions in the Precambrian terranes surrounding the rift can be used to define characteristic crustal components. These components may be used as endmembers in petrogenetic modelling of the Oslo Rift magmatic system. Based on available data, six endmember components can be identified, and (semi) quantitatively characterized in terms of Sr, Nd and Pb isotopes and selected trace elements. Data on the distribution of rock-types along the rift flanks allow estimates to be made of the relative importance of the components in different parts of the rift. Combining these data with petrological information may allow a realistic understanding of crust–magma interaction in the Oslo Rift magmatic system.