A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.     

Remote measurement of river morphology via fusion of LiDAR topography and spectrally based bathymetry

This study developed and evaluated a hybrid approach to remote measurement of river morphology that combines LiDAR topography with spectrally based bathymetry. Comparison of filtered LiDAR point clouds with surveyed cross‐sections indicated that subtle features on low‐relief floodplains were accurately resolved by LiDAR but that submerged areas could not be detected due to strong absorption of near‐infrared laser pulses by water. The reduced number of returns made the active channel evident in a LiDAR point density map. A second dataset suggested that pulse intensity also could be used to discriminate land from water via a threshold‐based masking procedure. Fusion of LiDAR and optical data required accurate co‐registration of images to the LiDAR, and we developed an object‐oriented procedure for achieving this alignment. Information on flow depths was derived by correlating pixel values with field measurements of depth. Highly turbid conditions dictated a positive relation between green band radiance and flow depth and contributed to under‐prediction of pool depths. Water surface elevations extracted from the LiDAR along the water's edge were used to produce a continuous water surface that preserved along‐channel variations in slope. Subtracting local flow depths from this surface yielded estimates of the bed elevation that were then combined with LiDAR topography for exposed areas to create a composite representation of the riverine terrain. The accuracy of this terrain model was assessed via comparison with detailed field surveys. A map of elevation residuals showed that the greatest errors were associated with underestimation of pool depths and failure to capture cross‐stream differences in water surface elevation. Nevertheless, fusion of LiDAR and passive optical image data provided an efficient means of characterizing river morphology that would not have been possible if either dataset had been used in isolation. Copyright © 2011 John Wiley & Sons, Ltd.     

Ionospheric measurements of relative coronal brightness during the total solar eclipses of 11 August, 1999 and 9 July, 1945

Swept-frequency (1/10 MHz) ionosonde measurements were made at Helston, Cornwall (50°06N, 5°18W) during the total solar eclipse on August 11, 1999. Soundings were made every three minutes. We present a method for estimating the percentage of the ionising solar radiation which remains unobscured at any time during the eclipse by comparing the variation of the ionospheric E-layer with the behaviour of the layer during a control day. Application to the ionosonde date for 11 August, 1999, shows that the flux of solar ionising radiation fell to a minimum of 25±2% of the value before and after the eclipse. For comparison, the same technique was also applied to measurements made during the total solar eclipse of 9 July, 1945, at Sörmjöle (63°68N, 20°20E) and yielded a corresponding minimum of 16 ± 2%. Therefore the method can detect variations in the fraction of solar emissions that originate from the unobscured corona and chromosphere. We discuss the differences between these two eclipses in terms of the nature of the eclipse, short-term fluctuations, the sunspot cycle and the recently-discovered long-term change in the coronal magnetic field.     

An experimental study on non-linear progressive wave-induced dynamic stresses in seabed

In this study, the dynamic stresses within the seabed induced by non-linear progressive waves were explored through a series of hydraulic model tests on a movable bed within a wave flume. By comparing Stokes’ 2nd-order wave theory with the theory of wave-induced dynamic stresses within the seabed as proposed by Yamamoto et al. [1978. On the response of a poro-elastic bed to water waves. Journal of Fluid Mechanics 87 (1), 193–206.] and Hsu and Jeng [1994. Wane-induced soil response in an unsaturated anisotropic seabed of finite thickness. International Journal for Numerical and Analytical Methods in Geomechanics 18, 785–807], the experimental results show that the pressure on the seabed surface, the pore water pressure within the seabed as well as the vertical and the horizontal stresses are all smaller than their theoretical values. If we were to obtain the characteristics of seabed soil, the analytical solution of Hsu and Jeng [1994. Wane-induced soil response in an unsaturated anisotropic seabed of finite thickness. International Journal for Numerical and Analytical Methods in Geomechanics 18, 785–807] might agree to the simulation of the wave-induced effective stresses and shear stress in the sandy seabed. A different phase shift exists among all the three soil stresses. Their influences on the three dynamic stresses within seabed soil are important for seabed stability, and can be used in the verification of numerical models. In the whole, the non-linear progressive waves and the naturally deposited seabed are found to have a strong interaction, and the behavior of the induced dynamic stresses within the seabed is very complicated, and should be investigated integrally.     

Contribution to diffuse gamma-rays in the Galactic Centre region from unresolved millisecond pulsars

Diffuse gamma-rays in the Galactic Centre region have been studied. We propose that there exists a population of millisecond pulsars in the Galactic Centre, which emit GeV gamma-rays through synchrotron-curvature radiation as predicted by outer gap models. These GeV gamma-rays from unresolved millisecond pulsars probably contribute to the diffuse gamma-ray spectrum detected by EGRET which displays a break at a few GeV. We have used a Monte Carlo method to obtain simulated samples of millisecond pulsars in the Galactic Centre region covered by EGRET  (∼ 15)  according to the different period and magnetic field distributions from observed millisecond pulsars in the Galactic field and globular clusters, and superposed their synchrotron-curvature spectra to derive the total GeV flux. Our simulated results suggest that there probably exist about 6000 unresolved millisecond pulsars in the region of angular resolution of EGRET, the emissions of which could contribute significantly to the observed diffuse gamma-rays in the Galactic Centre.     

GIS-based hazard mapping and zonation of debris flows in Xiaojiang Basin,southwestern China

Debris flow sites were identified at 140 locations in the Xiaojiang Basin in Yunnan province, southwestern China. Their spatial distribution and catchment characteristics are described in detail on the basis of previous research, air photo interpretation, field investigation and mapping using Geographical Information Systems (GIS). Using a statistical approach, a quantitative model of hazards assessment and zonation was developed through synthesis analysis of basin areas, gradients, and the relative reliefs of these debris flow sites. In terms of debris flow hazard assessment, areas within the Xiaojiang Basin can be classified as severe, heavy and light hazard regions.     

Fe and O isotope composition of meteorite fusion crusts: Possible natural analogues to chondrule formation?

Meteorite fusion crust formation is a brief event in a high‐temperature (2000–12,000 K) and high‐pressure (2–5 MPa) regime. We studied fusion crusts and bulk samples of 10 ordinary chondrite falls and 10 ordinary chondrite finds. The fusion crusts show a typical layering and most contain vesicles. All fusion crusts are enriched in heavy Fe isotopes, with δ⁵⁶Fe values up to +0.35‰ relative to the solar system mean. On average, the δ⁵⁶Fe of fusion crusts from finds is +0.23‰, which is 0.08‰ higher than the average from falls (+0.15‰). Higher δ⁵⁶Fe in fusion crusts of finds correlate with bulk chondrite enrichments in mobile elements such as Ba and Sr. The δ⁵⁶Fe signature of meteorite fusion crusts was produced by two processes (1) evaporation during atmospheric entry and (2) terrestrial weathering. Fusion crusts have either the same or higher δ¹⁸O (0.9–1.5‰) than their host chondrites, and the same is true for Δ¹⁷O. The differences in bulk chondrite and fusion crust oxygen isotope composition are explained by exchange of oxygen between the molten surface of the meteorites with the atmosphere and weathering. Meteorite fusion crust formation is qualitatively similar to conditions of chondrule formation. Therefore, fusion crusts may, at least to some extent, serve as a natural analogue to chondrule formation processes. Meteorite fusion crust and chondrules exhibit a similar extent of Fe isotope fractionation, supporting the idea that the Fe isotope signature of chondrules was established in a high‐pressure environment that prevented large isotope fractionations. The exchange of O between a chondrule melt and an ¹⁶O‐poor nebula as the cause for the observed nonmass dependent O isotope compositions in chondrules is supported by the same process, although to a much lower extent, in meteorite fusion crusts.     

Estimating contamination potential at waste-disposal sites using a natural tracer

Chloride is a conservative, natural tracer found in precipitation, soil water, and groundwater. The chloride mass-balance approach, long used to estimate groundwater recharge, also provides a downward flux of moisture and solute at sites where there is a potential for groundwater contamination. The flux is obtained by dividing the product of the mean annual precipitation and total annual chloride input (via precipitation and dust) by the mean soil-water chloride content. Chlorideversusdepth profiles can also be used to determine optimum depth of waste burial to minimize deterioration of waste containers. The method has been applied to three sites in arid alluvial-basin settings in New Mexico, U.S.A.: a proposed landfill, a battery recycling plant, and a hazardous-waste disposal facility. It is concluded that the method is reliable, economical, and practical. Furthermore, it can be applied at any stage in the development of a site. The chloride method should apply in any recharge area where the base of the root zone is separated from the water table by at least 3 m or so and chloride in soil water comes only from precipitation and dust.