Radio Observations of Ions in Comet Hale-Bopp

We present radio observations of HCO⁺, H₃O⁺ and CO⁺ in Hale-Bopp. This revised version was published online in July 2006 with corrections to the Cover Date.     

The surface chemistry of multi-oxide silicates

The surface chemistry of natural wollastonite, diopside, enstatite, forsterite, and albite in aqueous solutions was characterized using both electrokinetic techniques and surface titrations performed for 20 min in batch reactors. Titrations performed in such reactors allow determination of both proton consumption and metal release from the mineral surface as a function of pH. The compositions, based on aqueous solution analysis, of all investigated surfaces vary dramatically with solution pH. Ca and Mg are preferentially released from the surfaces of all investigated divalent metal silicates at pH less than ∼8.5-10 but preferentially retained relative to silica at higher pH. As such, the surfaces of these minerals are Si-rich and divalent metal poor except in strongly alkaline solutions. The preferential removal of divalent cations from these surfaces is coupled to proton consumption. The number of protons consumed by the preferential removal of each divalent cation is pH independent but depends on the identity of the mineral; ∼1.5 protons are consumed by the preferential removal of each Ca atom from wollastonite, ∼3 protons are consumed by the preferential removal of each Mg or Ca atom from diopside or enstatite, and ∼4 protons are consumed by the preferential removal of each Mg from forsterite. These observations are interpreted to stem from the creation of additional ‘internal’ adsorption sites by the preferential removal of divalent metal cations which can be coupled to the condensation of partially detached Si. Similarly, Na and Al are preferentially removed from the albite surface at 2 > pH > 11; mass balance calculations suggest that three protons are consumed by the preferential removal of each Al atom from this surface over this entire pH range. Electrokinetic measurements on fresh mineral powders yield an isoelectric point (pH_IEP) 2.6, 4.4, 3.0, 4.5, and <1, for wollastonite, diopside, enstatite, forsterite, and albite, respectively, consistent with the predominance of SiO₂ in the surface layer of all of these multi-oxide silicates at acidic pH. Taken together, these observations suggest fundamental differences between the surface chemistry of simple versus multi-oxide minerals including (1) a dependency of the number and identity of multi-oxide silicate surface sites on the aqueous solution composition, and (2) the dominant role of metal-proton exchange reactions on the reactivity of multi-oxide mineral surfaces including their dissolution rate variation with aqueous solution composition.     

On the modelling of stationary planetary waves

The structure of the atmospheric stationary planetary waves is obtained by means of the quasi‐geostrophic linear model developed by Matsuno (1 970). To show the influence of the upper boundary condition on the structure of the waves, the latter is computed using a rigid top and a radiation condition. The modifications to the wave structure obtained when the upper boundary is lowered from a height of 65 km to 42.5 km, 32.5 km, and 22.5 km are examined. The effects of varying the vertical grid increment are studied through a comparison of the wave structures obtained with 6, 12, and 24 levels between 5 and 65 km.     

Mid-Holocene hemlock decline and diatom communities in van Nostrand Lake, Ontario, Canada

Terrestrial ecosystem disturbances inferred from the fossil hemlock pollen decline (ca. 4,800 BP) and recovery (ca. 3,800 BP) affected van Nostrand Lake, including the diatom communities. Ecological models suggest the lake responded by eutrophying, reflecting higher nutrient influx resulting from increased erosion. A decline in lake productivity followed as the forest vegetation recovered and erosion slowed. Lastly, as the forest switched from early and middle successional species to mature species, lake productivity increased as erosion increased nutrient inflow, especially phosphorus. However, this eutrophication response was delayed or buffered, perhaps due to wetland development surrounding the lake. The lake never fully returned to its initial, pre-hemlock decline state, but oscillated between more eutrophic and more mesotrophic or oligotrophic stages, possibly reflecting other disturbances in the catchment and climatic change.     

Spectral analysis of earthquake migration in South America

Earthquake migration along linear seismic belts is investigated by analyzing spacetime diagrams using spectral analysis. In order to sample the earthquakes in the space-time domain, they must first be convolved with a (sinx sint)/xt surface to obtain an unbiased and alias free twodimensional Fourier spectrum. Further enhancements are provided by selectively stacking patterns (a pattern is defined as the distribution of earthquakes in space and time before a particular earthquake), thereby reinforcing the similarities within the various patterns. With these techniques, it is possible to quantitatively estimate the migration rates (from theirspatial frequencies) and recurrence intervals (from theirtemporal frequencies) of large earthquakes in South America.Preliminary examination of the spectra for South America indicates that a low frequency peak occurs at approximately 2500 km and 27 years for earthquakes with magnitudes greater than 7.7. The results suggest a migration rate of approximately 95 km/yr from south to north and a recurrence interval of 27 years.     

Including prior information in the estimation of effective soil parameters in unsaturated zone modelling

In this paper we propose a methodology to include prior information in the estimation of effective soil parameters for modelling the soil moisture content in the unsaturated zone. Laboratory measurements on undisturbed soil cores were used to estimate the moisture retention curve and hydraulic conductivity curve parameters. The soil moisture content was measured at 25 locations along three transects and at three different depths (surface, 30 and 60 cm) on an 80×20 m hillslope for the year 2001. Soil cores were collected in 84 locations situated in three profile pits along the hillslope. For the estimation of the effective soil hydraulic parameters the joint probability distribution of measured parameter values was used as prior information. A two-horizon single column 1D MIKE SHE model based on Richards' equation was set-up for nine soil moisture measurement locations along the middle transect of the hillslope. The goal of the model is to simulate the soil moisture profile at each location. The shuffled complex evolution (SCE) algorithm has been applied to estimate effective model parameters using either wide parameter ranges, referred to as the ‘no-prior’ case, or the joint probability distribution of measured parameter values as prior information (‘prior’ case). When the prior information is incorporated in the SCE optimisation the goodness-of-fit of the model predictions is only slightly worse compared to when no-prior information is incorporated. However, the effective parameter estimates are more realistic when the prior information is incorporated. For both the no-prior and prior case the generalised likelihood uncertainty estimation procedure (GLUE) was subsequently used to estimate the uncertainty bounds (UB) on the model predictions. When incorporating the prior information more parameter sets were accepted for the estimation of the predictive uncertainty and the parameter values were more realistic. Moreover, UB better enclosed the observations. Thus, incorporating prior information in GLUE reduces the amount of model evaluations needed to obtain sufficient behavioural parameter sets. The results indicate the importance of prior information in the SCE and GLUE parameter estimation strategies.     

Seismology software: state of the practice

We analyzed the state of practice for software development in the seismology domain by comparing 30 software packages on four aspects: product, implementation, design, and process. We found room for improvement in most seismology software packages. The principal areas of concern include a lack of adequate requirements and design specification documents, a lack of test data to assess reliability, a lack of examples to get new users started, and a lack of technological tools to assist with managing the development process. To assist going forward, we provide recommendations for a document-driven development process that includes a problem statement, development plan, requirement specification, verification and validation (V&V) plan, design specification, code, V&V report, and a user manual. We also provide advice on tool use, including issue tracking, version control, code documentation, and testing tools.