Condition number estimates for combined potential boundary integral operators in acoustic scattering

We study the classical combined field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle, namely the indirect formulation due to Brakhage-Werner/Leis/Panic, and the direct formulation associated with the names of Burton and Miller. We obtain lower and upper bounds on the condition numbers for these formulations, emphasising dependence on the frequency, the geometry of the scatterer, and the coupling parameter. Of independent interest we also obtain upper and lower bounds on the norms of two oscillatory integral operators, namely the classical acoustic single- and double-layer potential operators.

Storm Tracking and eScience

A new 'storm-tracking approach' to analysing the prediction of storms by different forecast systems has recently been developed. This paper provides a brief illustration of the type of results/information that can be obtained using the approach. It also describes in detail how eScience methodologies have been used to help apply the storm-tracking approach to very large datasets

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Modelling habitat selection and distribution of the critically endangered Jerdon's courser Rhinoptilus bitorquatus in scrub jungle: an application of a new tracking method

1. Jerdon's courser Rhinoptilus bitorquatus is a nocturnally active cursorial bird that is only known to occur in a small area of scrub jungle in Andhra Pradesh, India, and is listed as critically endangered by the IUCN. Information on its habitat requirements is needed urgently to underpin conservation measures. We quantified the habitat features that correlated with the use of different areas of scrub jungle by Jerdon's coursers, and developed a model to map potentially suitable habitat over large areas from satellite imagery and facilitate the design of surveys of Jerdon's courser distribution. 2. We used 11 arrays of 5-m long tracking strips consisting of smoothed fine soil to detect the footprints of Jerdon's coursers, and measured tracking rates (tracking events per strip night). We counted the number of bushes and trees, and described other attributes of vegetation and substrate in a 10-m square plot centred on each strip. We obtained reflectance data from Landsat 7 satellite imagery for the pixel within which each strip lay. 3. We used logistic regression models to describe the relationship between tracking rate by Jerdon's coursers and characteristics of the habitat around the strips, using ground-based survey data and satellite imagery. 4. Jerdon's coursers were most likely to occur where the density of large (>2 m tall) bushes was in the range 300-700 ha(-1) and where the density of smaller bushes was less than 1000 ha(-1). This habitat was detectable using satellite imagery. 5. Synthesis and applications. The occurrence of Jerdon's courser is strongly correlated with the density of bushes and trees, and is in turn affected by grazing with domestic livestock, woodcutting and mechanical clearance of bushes to create pasture, orchards and farmland. It is likely that there is an optimal level of grazing and woodcutting that would maintain or create suitable conditions for the species. Knowledge of the species' distribution is incomplete and there is considerable pressure from human use of apparently suitable habitats. Hence, distribution mapping is a high conservation priority. A two-step procedure is proposed, involving the use of ground surveys of bush density to calibrate satellite image-based mapping of potential habitat. These maps could then be used to select priority areas for Jerdon's courser surveys. The use of tracking strips to study habitat selection and distribution has potential in studies of other scarce and secretive species.

Tracking larval insect movement within soil using high resolution X-ray microtomography

1. In contrast to above-ground insects, comparatively little is known about the behaviour of subterranean insects, due largely to the difficulty of studying them in situ. 2. The movement of newly hatched (neonate) clover root weevil (Sitona lepidus L. Coleoptera: Curculinidae) larvae was studied non-invasively using recently developed high resolution X-ray microtomography. 3. The movement and final position of S. lepidus larvae in the soil was reliably established using X-ray microtomography, when compared with larval positions that were determined by destructively sectioning the soil column. 4. Newly hatched S. lepidus larvae were seen to attack the root rhizobial nodules of their host plant, white clover (Trifolium repens L.). Sitona lepidus larvae travelled between 9 and 27 mm in 9 h at a mean speed of 1.8 mm h(-1). 5. Sitona lepidus larvae did not move through the soil in a linear manner, but changed trajectory in both the lateral and vertical planes.

Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks

This paper describes the measurements of the acoustic and petrophysical properties of two suites of low-shale sandstone samples from North Sea hydrocarbon reservoirs, under simulated reservoir conditions. The acoustic velocities and quality factors of the samples, saturated with different pore fluids (brine, dead oil and kerosene), were measured at a frequency of about 0.8 MHz and over a range of pressures from 5 MPa to 40 MPa. The compressional-wave velocity is strongly correlated with the shear-wave velocity in this suite of rocks. The ratio V-P/V-S varies significantly with change of both pore-fluid type and differential pressure, confirming the usefulness of this parameter for seismic monitoring of producing reservoirs. The results of quality factor measurements were compared with predictions from Biot-flow and squirt-flow loss mechanisms. The results suggested that the dominating loss in these samples is due to squirt-flow of fluid between the pores of various geometries. The contribution of the Biot-flow loss mechanism to the total loss is negligible. The compressional-wave quality factor was shown to be inversely correlated with rock permeability, suggesting the possibility of using attenuation as a permeability indicator tool in low-shale, high-porosity sandstone reservoirs.

A wavenumber independent boundary element method for an acoustic scattering problem

In this paper we consider the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data, a problem which models, for example, outdoor sound propagation over inhomogeneous. at terrain. To achieve good approximation at high frequencies with a relatively low number of degrees of freedom, we propose a novel Galerkin boundary element method, using a graded mesh with smaller elements adjacent to discontinuities in impedance and a special set of basis functions so that, on each element, the approximation space contains polynomials ( of degree.) multiplied by traces of plane waves on the boundary. We prove stability and convergence and show that the error in computing the total acoustic field is O( N-(v+1) log(1/2) N), where the number of degrees of freedom is proportional to N logN. This error estimate is independent of the wavenumber, and thus the number of degrees of freedom required to achieve a prescribed level of accuracy does not increase as the wavenumber tends to infinity.

Acoustic scattering by mildly rough unbounded surfaces in three dimensions

For a nonlocally perturbed half- space we consider the scattering of time-harmonic acoustic waves. A second kind boundary integral equation formulation is proposed for the sound-soft case, based on a standard ansatz as a combined single-and double-layer potential but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half- space Green's function. Due to the unboundedness of the surface, the integral operators are noncompact. In contrast to the two-dimensional case, the integral operators are also strongly singular, due to the slow decay at infinity of the fundamental solution of the three-dimensional Helmholtz equation. In the case when the surface is sufficiently smooth ( Lyapunov) we show that the integral operators are nevertheless bounded as operators on L-2(Gamma) and on L-2(Gamma G) boolean AND BC(Gamma) and that the operators depend continuously in norm on the wave number and on G. We further show that for mild roughness, i.e., a surface G which does not differ too much from a plane, the boundary integral equation is uniquely solvable in the space L-2(Gamma) boolean AND BC(Gamma) and the scattering problem has a unique solution which satisfies a limiting absorption principle in the case of real wave number.

A high-wavenumber boundary-element method for an acoustic scattering problem

In this paper we show stability and convergence for a novel Galerkin boundary element method approach to the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data. This problem models, for example, outdoor sound propagation over inhomogeneous flat terrain. To achieve a good approximation with a relatively low number of degrees of freedom we employ a graded mesh with smaller elements adjacent to discontinuities in impedance, and a special set of basis functions for the Galerkin method so that, on each element, the approximation space consists of polynomials (of degree $\nu$) multiplied by traces of plane waves on the boundary. In the case where the impedance is constant outside an interval $[a,b]$, which only requires the discretization of $[a,b]$, we show theoretically and experimentally that the $L_2$ error in computing the acoustic field on $[a,b]$ is ${\cal O}(\log^{\nu+3/2}|k(b-a)| M^{-(\nu+1)})$, where $M$ is the number of degrees of freedom and $k$ is the wavenumber. This indicates that the proposed method is especially commendable for large intervals or a high wavenumber. In a final section we sketch how the same methodology extends to more general scattering problems.

Investigating the practicality of using radio tracking to determine home range and movements of Picathartidae

Radiotelemetry is an important tool used to aid the understanding and conservation of cryptic and rare birds. The two bird species of the family Picathartidae are little-known, secretive, forest-dwelling birds endemic to western and central Africa. In 2005, we conducted a radio-tracking trial of Grey-necked Picathartes Picathartes oreas in the Mbam Minkom Mountain Forest, southern Cameroon, using neck collar (two birds) and tail-mounted (four birds) transmitters to investigate the practicality of radio-tracking Picathartidae. Three birds with tail-mounted transmitters were successfully tracked with the fourth, though not relocated for radio tracking, resighted the following breeding season. Two of these were breeding birds that continued to provision young during radio tracking. One neck-collared bird was found dead three days after transmitter attachment and the other neither relocated nor resighted. As mortality in one bird was potentially caused by the neck collar transmitter we recommend tail-mounted transmitters in future radio-tracking studies of Picathartidae. Home ranges, shown using minimum convex polygon and kernel estimation methods, were generally small (<0.5 km(2)) and centred around breeding sites. A minimum of 60 fixes were found to be sufficient for home range estimation.