The NAIAD experiment for WIMP searches at Boulby mine and recent results

The NAIAD experiment (NaI Advanced Detector) for weakly interacting massive particle (WIMP) dark matter searches at Boulby mine (UK) is described. The detector consists of an array of encapsulated and unencapsulated NaI(Tl) crystals with high light yield. Six crystals are collecting data at present. Data accumulated by four of them (10.6 kg × year exposure) have been used to set upper limits on the WIMP–nucleon spin-independent and WIMP–proton spin-dependent cross-sections. Pulse shape analysis has been applied to discriminate between nuclear recoils, as may be caused by WIMP interactions, and electron recoils due to gamma background. Various calibrations of crystals are presented.     

Simulations of a giant hybrid air shower detector

We have used Monte Carlo simulations to investigate the capabilities of a giant air shower observatory designed to detect showers initiated by cosmic rays with energies exceeding 10¹⁹ eV. The observatory is to consist of an array of detectors that will characterise the air shower at ground level, and optical detectors to measure the fluorescence light emitted by the shower in the atmosphere. Using these detectors together in a ‘hybrid’ configuration, we find that precise geometrical reconstruction of the shower axis is possible, leading to excellent resolution in energy. The technique is also shown to provide very good reconstruction below 10¹⁹ eV, at energies where the ground array is not fully efficient.     

Late Holocene relative sea level rise and the Neoglacial history of the Greenland ice sheet

In West Greenland, early and mid Holocene relative sea level (RSL) fall was replaced by late Holocene RSL rise during the Neoglacial, after 4–3 cal. ka BP (thousand calibrated years before present). Here we present the results of an isolation basin RSL study completed near to the coastal town of Sisimiut, in central West Greenland. RSL fell from 14 m above sea level at 5.7 cal. ka BP to reach a lowstand of ?4.0 m at 2.3–1.2 cal. ka BP, before rising by an equivalent amount to present. Differences in the timing and magnitude of the RSL lowstand between this and other sites in West and South Greenland record the varied interplay of local and non‐Greenland RSL processes, notably the reloading of the Earth's crust caused by a Neoglacial expansion of the Greenland Ice Sheet (GIS) and the subsidence associated with the collapse of the Laurentide Ice Sheet forebulge. This means that the timing of the sea level lowstand cannot be used to infer directly when the GIS advanced during the Neoglacial. The rise in Late Holocene RSL is contrary to recently reported bedrock uplift in the Sisimiut area, based on repeat GPS surveys. This indicates that a belt of peripheral subsidence around the current ice sheet margin was more extensive in the late Holocene, and that there has been a switch from subsidence to uplift at some point in the last thousand years or so. Copyright © 2008 John Wiley & Sons, Ltd.     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.