A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.     

Simulation of the influence of solar radiation variations on the global climate with an ocean-atmosphere general circulation model

 Two simulations with a global coupled ocean-atmosphere circulation model have been carried out to study the potential impact of solar variability on climate. The Hoyt and Schatten estimate of solar variability from 1700 to 1992 has been used to force the model. Results indicate that the near-surface temperature simulated by the model is dominated by the long periodic solar fluctuations (Gleissberg cycle), with global mean temperatures varying by about 0.5 K. Further results indicate that solar variability and an increase in greenhouse gases both induce to a first approximation a comparable pattern of surface temperature change, i.e., an increase of the land-sea contrast. However, the solar-induced warming pattern in annual means and summer is more centered over the subtropics, compared to a more uniform warming associated with the increase in greenhouse gases. The observed temperature rise over the most recent 30 and 100 years is larger than the trend in the solar forcing simulation during the same period, indicating a strong likelihood that, if the model forcing and response is realistic, other factors have contributed to the observed warming. Since the pattern of the recent observed warming agrees better with the greenhouse warming pattern than with the solar variability response, it is likely that one of these factors is the increase of the atmospheric greenhouse gas concentration. Received: 14 October 1996 / Accepted: 9 May 1997     

Dynamic interaction of bridge–train system under non‐uniform seismic ground motion

The probability that an earthquake occurs when a train is running over a bridge in earthquake‐prone regions is much higher than before, for high‐speed railway lines are rapidly developed to connect major cities worldwide. This paper presents a finite element method‐based framework for dynamic analysis of coupled bridge–train systems under non‐uniform seismic ground motion, in which rail–wheel interactions and possible separations between wheels and rails are taken into consideration. The governing equations of motion of the coupled bridge–train system are established in an absolute coordinate system. Without considering the decomposition of seismic responses into pseudo‐static and inertia‐dynamic components, the equations of motion of the coupled system are formed in terms of displacement seismic ground motions. The mode superposition method is applied to the bridge structure to make the problem manageable while the Newmark‐β method with an iterative computation scheme is used to find the best solution for the problem concerned. Eight high‐speed trains running over a multi‐span steel truss‐arch bridge subject to earthquakes are taken as a case study. The results from the case study demonstrate that the spatial variation of seismic ground motion affects dynamic responses of the bridge–train system. The ignorance of pseudo‐static component when using acceleration seismic ground motions as input may underestimate seismic responses of the bridge–train system. The probability of separation between wheels and rails becomes higher with increasing train speed. Copyright © 2011 John Wiley & Sons, Ltd.     

A cellular model of Holocene upland river basin and alluvial fan evolution

The CAESAR (Cellular Automaton Evolutionary Slope And River) model is used to simulate the Holocene development of a small upland catchment (4·2 km²) and the alluvial fan at its base. The model operates at a 3 m grid scale and simulates every flood over the last 9200 years, using a rainfall record reconstructed from peat bog wetness indices and land cover history derived from palynological sources. Model results show that the simulated catchment sediment discharge above the alluvial fan closely follows the climate signal, but with an increase in the amplitude of response after deforestation. The important effects of sediment storage and remobilization are shown, and findings suggest that soil creep rates may be an important control on long term (>1000 years) temperate catchment sediment yield. The simulated alluvial fan shows a complex and episodic behaviour, with frequent avulsions across the fan surface. However, there appears to be no clear link between fan response and climate or land use changes suggesting that Holocene alluvial fan dynamics may be the result of phases of sediment storage and remobilization, or instabilities and thresholds within the fan itself. Copyright © 2002 John Wiley & Sons, Ltd.     

Meso‐scale particle modeling of concrete deterioration caused by alkali‐aggregate reaction

A meso‐scale particle model is presented to simulate the expansion of concrete subjected to alkali‐aggregate reaction (AAR) and to analyze the AAR‐induced degradation of the mechanical properties. It is the first attempt to evaluate the deterioration mechanism due to AAR using the discrete‐element method. A three‐phase meso‐scale model for concrete composed of aggregates, mortar and the interface is established with the combination of a pre‐processing approach and the particle flow code, PFC2D. A homogeneous aggregate expansion approach is applied to model the AAR expansion. Uniaxial compression tests are conducted for the AAR‐affected concrete to examine the effects on the mechanical properties. Two specimens with different aggregate sizes are analyzed to consider the effects of aggregate size on AAR. The results show that the meso‐scale particle model is valid to predict the expansion and the internal micro‐cracking patterns caused by AAR. The two different specimens exhibit similar behavior. The Young's modulus and compressive strength are significantly reduced with the increase of AAR expansion. The shape of the stress–strain curves obtained from the compression tests clearly reflects the influence of internal micro‐cracks: an increased nonlinearity before the peak loading and a more gradual softening for more severely affected specimens. Similar macroscopic failure patterns of the specimens under compression are observed in terms of diagonal macroscopic cracks splitting the specimen into several triangular pieces, whereas localized micro‐cracks forming in slightly affected specimens are different from branching and diffusing cracks in severely affected ones, demonstrating different failure mechanisms. Copyright © 2012 John Wiley & Sons, Ltd.     

Recognition of geochemical footprints of mineral systems in the regolith at regional to continental scales

Understanding the character of Australia's extensive regolith cover is crucial to the continuing success of mineral exploration. We hypothesise that the regolith contains geochemical fingerprints of processes related to the development and preservation of mineral systems at a range of scales. We test this hypothesis by analysing the composition of surface sediments within greenfield regional-scale (southern Thomson Orogen) and continental-scale (Australia) study areas. In the southern Thomson Orogen area, the first principal component (PC1) derived in our study [Ca, Sr, Cu, Mg, Au and Mo at one end; rare earth elements (REEs) and Th at the other] is very similar to the empirical vector used by a local company (enrichment in Sr, Ca and Au concomitant with depletion in REEs) to successfully site exploration drill holes for Cu–Au mineralisation. Mapping of the spatial distribution of PC1 in the region reveals several areas of elevated values and possible mineralisation potential. One of the strongest targets in the PC1 map is located between Brewarrina and Bourke in northern New South Wales. Here, exploration drilling has intersected porphyry Cu–Au mineralisation with up to 1 wt% Cu, 0.1 g/t Au, and 717 ppm Zn. The analysis of a comparable geochemical dataset at the continental scale yields a compositionally similar PC1 (Ca, Sr, Mg, Cu, Au and Mo at one end; REEs and Th at the other) to that of the regional study. Mapping PC1 at the continental scale shows patterns that (1) are spatially compatible with the regional study and (2) reveal several geological regions of elevated values, possibly suggesting an enhanced potential for porphyry Cu–Au mineralisation. These include well-endowed mineral provinces such as the Curnamona and Capricorn regions, but also some greenfield regions such as the Albany-Fraser/western Eucla, western Murray and Eromanga geological regions. We conclude that the geochemical composition of Australia's regolith may hold critical information pertaining to mineralisation within/beneath it.     

The molecular composition of ambers

Bulk (elemental composition, IR, CP/MAS ¹³C NMR) and molecular (GC-MS) analyses have been performed on a series of ambers and resins derived from different locations (Dominican Republic, Philippines, Canada, Israel, New Zealand, Chile) having diverse botanical affinities (Araucariaceae, Hymenaea) and variable age (from Holocene to Early Cretaceous). No major differences have been observed from the elemental composition and the spectroscopic data; however, the molecular analyses of the solvent extractable fraction show that a specific mixture of components is present in each sample. These are mainly diterpenoid products that in general are also found abundantly in the higher plants from which the ambers and resins originate. Nevertheless, a direct relationship between major terpenoid constituents in fossil resins and precursor plant materials can only be established for the younger samples.Irrespective of the geographical or botanical origin of the ambers and resins, several common age-dependent molecular transformation trends can be recognized: (1) progressive loss of olefinic bonds (especially those located in exocyclic positions), (2) decrease of functionalized products, and (3) increasing proportion of aromatized components. However, even in the samples of older age (Cretaceous) the degree of aromatization is very low when compared with that of other higher-plant related materials such as fossilized woods or low rank coals. This indicates that maturation must involve essentially olefin polymerization processes instead of extensive aromatization.     

Effects of sediment pulses on bed relief in bar‐pool channels

To further develop prediction of the range of morphological adjustments associated with sediment pulses in bar‐pool channels, we analyze channel bed topographic data collected prior to and following the removal of two dams in Oregon: Marmot Dam on the Sandy River and Brownsville Dam on the Calapooia River. We hypothesize that, in gravel‐bed, bar‐pool channels, the response of bed relief to sand and gravel sediment pulses is a function of initial relief and pulse magnitude. Modest increases in sediment supply to initially low‐relief, sediment‐poor cross‐sections will increase bed relief and variance of bed relief via bar deposition. Modest increases in sediment supply to initially high‐relief cross‐sections, characteristic of alternate bar morphology, will result in decreased bed relief and variance of relief via deposition in bar‐adjacent pools. These hypothesized adjustments are measured in terms of bed relief, which we define as the difference in elevation between the pool‐bottom and bar‐top. We evaluate how relief varies with sediment thickness, where both relief and mean sediment thickness at a cross‐section are normalized by the 90th percentile of observed relief values within a reach prior to a sediment pulse. Field measurements generally supported the stated hypotheses, demonstrating how introduction of a sediment pulse to low‐relief reaches can increase mean and variance of relief, while introduction to high‐relief reaches can decrease the mean and variance of bed relief, at least temporarily. In general, at both sites, the degree of impact increased with the thickness of sediment delivered to the cross‐section. Results thus suggest that the analysis is a useful step for understanding the morphological effects of sediment pulses introduced to gravel‐bed, bar‐pool channels. Copyright © 2015 John Wiley & Sons, Ltd.     

A new, instantaneous aeolian sand trap design for field use

A new aeolian sediment trap is described which can give up to 1 Hz measurement frequency in field conditions. The trap adopts a circular, horizontal trap design with a load cell connection to give continuous, unobtrusive trap measurement of sediment flux. Simultaneous velocity recording is carried out using an anemometer. Trap construction costs are approximately £200. Initial results in field conditions using a direct comparison of wind velocity data, sampled at an equivalent frequency, have given a first-order relationship between sediment flux and velocity. The trap enables simultaneous sampling of wind velocity and sediment flux at a sufficiently short interval to enable investigation of sediment transport dynamics under a variety of field conditions.