Computing Fresnel integrals via modified trapezium rules

In this paper we propose methods for computing Fresnel integrals based on truncated trapezium rule approximations to integrals on the real line, these trapezium rules modified to take into account poles of the integrand near the real axis. Our starting point is a method for computation of the error function of complex argument due to Matta and Reichel (J Math Phys 34:298–307, 1956) and Hunter and Regan (Math Comp 26:539–541, 1972). We construct approximations which we prove are exponentially convergent as a function of N , the number of quadrature points, obtaining explicit error bounds which show that accuracies of 10−15 uniformly on the real line are achieved with N=12 , this confirmed by computations. The approximations we obtain are attractive, additionally, in that they maintain small relative errors for small and large argument, are analytic on the real axis (echoing the analyticity of the Fresnel integrals), and are straightforward to implement.

Efficient evaluation of highly oscillatory acoustic scattering surface integrals

We consider the approximation of some highly oscillatory weakly singular surface integrals, arising from boundary integral methods with smooth global basis functions for solving problems of high frequency acoustic scattering by three-dimensional convex obstacles, described globally in spherical coordinates. As the frequency of the incident wave increases, the performance of standard quadrature schemes deteriorates. Naive application of asymptotic schemes also fails due to the weak singularity. We propose here a new scheme based on a combination of an asymptotic approach and exact treatment of singularities in an appropriate coordinate system. For the case of a spherical scatterer we demonstrate via error analysis and numerical results that, provided the observation point is sufficiently far from the shadow boundary, a high level of accuracy can be achieved with a minimal computational cost.

Mineral and geochemical characterization of a leptic aluandic soil and a thapto aluandic-ferralsol developed on trachytes in Mount Bambouto (Cameroon volcanic line)

Mineral and geochemical investigations were carried out on soil samples and fresh rock (trachytes) from two selected soil profiles (TM profile on leptic aluandic soils and TL profile on thapto aluandic-ferralsols) from Mount Bambouto to better understand geochemical processes and mineral paragenesis involved in the development of soils in this environment. In TM profile, the hydrated halloysites and goethite occur in the weathered saprolite boulders of BC horizon while dehydrated halloysite, gibbsite and goethite dominate the soils matrices of BC and A horizons. In TL profile, the dehydrated halloysites and goethite are the most abundant secondary minerals in the weathered saprolites of C and BC horizons while gibbsite, hematite and kaolinite occur in the soil matrices of BC, B and A horizons. The highest gibbsite content is in the platy nodules of B horizon. In both soil profiles, organo-metal complexes (most likely of AI and Fe) are present in the surface A horizon. Geochemically, between the fresh rock and the weathered saprolites in both soils, SiO2, K2O, CaO, Na2O and MgO contents decrease strongly while Fe2O3 and Al2O3 tend to accumulate. The molar ratio of SiO2/Al2O3 (Ki) and the sum of Ca, Mg, K and Na ions (TRB) also decreases abruptly between fresh rocks and the weathered saprolites, but increases significantly at the soil surface. The TM profile shows intense Al enrichment whereas the TL profile highlights enrichment in both AI and Fe as the weathering progresses upwards. Both soil profiles are enriched in Ni, Cu, Ba and Co and depleted in U, Th, Ta, Hf, Y, Sr, Pb, Zr and Zn relative to fresh rock. They also show a relatively low fractionation of the rare earth elements (REE: La, Nd, Sm, Eu, Tb, Yb and Lu), except for Ce which tends to be enriched in soils compared to CI chondrite. All these results give evidence of intense hydrolysis at soil deep in Mount Bambouto resulting in the formation of halloysite which progressively transforms into gibbsite and/or dehydrated halloysite. At the soil surface, the prominent pedogenetic process refers to andosolization with formation of organo-metal complexes. In TL profile, the presence of kaolinite in soil matrices BC and B horizons is consistent with ferralitization at soil deep. In conclusion, soil forming processes in Mount Bambouto are strongly influenced by local climate: (i) in the upper mountain (>2000 m), the fresh, misty and humid climate favors andosolization; whereas (ii) in the middle lands (1700-2000 m) with a relatively dry climate, both andosolization at the soil surface and ferralitization at soil deep act together. (C) 2009 Elsevier B.V. All rights reserved.

Assessment of the consistency of near-infrared water vapor line intensities using high-spectral-resolution ground-based Fourier transform measurements of solar radiation

Calculations of the absorption of solar radiation by atmospheric gases, and water vapor in particular, are dependent on the quality of databases of spectral line parameters. There has been increasing scrutiny of databases such as HITRAN in recent years, but this has mostly been performed on a band-by-band basis. We report nine high-spectral-resolution (0.03 cm(-1)) measurements of the solar radiation reaching the surface in southern England over the wave number range 2000 to 12,500 cm(-1) (0.8 to 5 mm) that allow a unique assessment of the consistency of the spectral line databases over this entire spectral region. The data are assessed in terms of the modeled water vapor column that is required to bring calculations and observations into agreement; for an entirely consistent database, this water vapor column should be constant with frequency. For the HITRAN01 database, the spread in water vapor column is about 11%, with distinct shifts between different spectral regions. The HITRAN04 database is in significantly better agreement (about 5% spread) in the completely updated 3000 to 8000 cm(-1) spectral region, but inconsistencies between individual spectral regions remain: for example, in the 8000 to 9500 cm(-1) spectral region, the results indicate an 18% (+/- 1%) underestimate in line intensities with respect to the 3000 to 8000 cm(-1) region. These measurements also indicate the impact of isotopic fractionation of water vapor in the 2500 to 2900 cm(-1) range, where HDO lines dominate over the lines of the most abundant isotope of H2O.

The Klein-Gordon equation on the half line: a Riemann-Hilbert approach

We solve an initial-boundary problem for the Klein-Gordon equation on the half line using the Riemann-Hilbert approach to solving linear boundary value problems advocated by Fokas. The approach we present can be also used to solve more complicated boundary value problems for this equation, such as problems posed on time-dependent domains. Furthermore, it can be extended to treat integrable nonlinearisations of the Klein-Gordon equation. In this respect, we briefly discuss how our results could motivate a novel treatment of the sine-Gordon equation.