New developments in eclipsing binary light curve modeling

The light curve modeling of binary stars has continued to evolve since its founding by Henry Norris Russell (see Russell and Merrill 1952 and citations therein) nearly a century ago, accelerated in the 1950s by Kopal's introduction of Roche geometry into models and by the development of synthetic light curve computer code in the 1970's. Improved physics and the use of more kinds of observational input are providing another round of important advances that promise to enlarge our knowledge of both binary stars and ensembles containing them. Here we discuss the newer horizons of light curve modeling and the steps being taken toward them.     

Vertical stress distribution due to an arbitrarily shaped foundation using triangulation and an analytical expression for a triangular loaded region

Using the basic Boussinesq's equation, the expression for the vertical stress distribution (σ_z) underneath any point on the ground surface due to a general triangular loaded region in a preferred orientation with a linearly varied loading has been successfully derived. When the triangle is not in a preferred orientation, a simple axis transformation is required and the expression will be equally applicable. Based on this expression, σ_z due to an arbitrarily shaped loaded foundation can simply be determined by first triangulating the loaded area and summing up the contributions from each generated triangular region. The procedures for triangulating and calculating the stress distribution can be simply automated through computer programs.     

Response of circular footings and anchor plates in non-homogeneous elastic soils

The axisymmetric elastic response of circular footings and anchor plates in a linearly non-homogeneous elastic soil is analysed. It is assumed that footings/anchors are flexible and subjected to axisymmetric vertical loads. The response of the footing/anchor is modelled by using the classical Poisson–Kirchhoff thin plate theory. A variational technique is used to analyse the interaction problem. A representation for the contact stress is established by using a fundamental solution corresponding to a unit vertical pressure acting over an annular region in the interior of the non-homogeneous soil. The fundamental solution can be derived by using rigorous analytical procedures. The influence of the footing flexibility and the degree of soil non-homogeneity on the displacements, bending moments and contact stresses of a surface footing is examined over a wide range of governing parameters. In the case of anchor plates the influence of depth of embedment, degree of soil non-homogeneity and anchor flexibility on the anchor displacement is investigated.     

Sources and sinks of sea salt in a newfoundland blanket bog

Over an oceanic peatland, the concentration of Na in fog averaged 38.1 mgl^?1 compared with 1.8 mgl^?1 in rain, resulting in a significant flux of mineral elements to the surface. Between 16 May and 20 June 1990 the average mass flux of Na to the bog surface by fog, rain, and dry deposition was 21.9, 10.4 and 7.0 mg m^?2 d^?1. There was little long-term storage of Na within the peatland system, where Na losses measured in stream runoff averaged 34.8 mg m² d^?1, and deep groundwater losses 4 mg m^?2 d^?1. Calcium and Mg were preferentially retained in the organic soil, whereas K was relatively mobile. Potassium tended to become concentrated in the unsaturated zone. Stream runoff had a consistently higher pH than groundwater, corresponding to higher Ca and Mg concentrations, which may have been from mineral sources in the headwater ponds. Otherwise, the stream water chemistry was closely related to groundwater in the upper layers of the peat deposit.     

Comparison of results from several AOGCMs for global and regional sea-level change 1900–2100

 Sea-level rise is an important aspect of climate change because of its impact on society and ecosystems. Here we present an intercomparison of results from ten coupled atmosphere-ocean general circulation models (AOGCMs) for sea-level changes simulated for the twentieth century and projected to occur during the twenty first century in experiments following scenario IS92a for greenhouse gases and sulphate aerosols. The model results suggest that the rate of sea-level rise due to thermal expansion of sea water has increased during the twentieth century, but the small set of tide gauges with long records might not be adequate to detect this acceleration. The rate of sea-level rise due to thermal expansion continues to increase throughout the twenty first century, and the projected total is consequently larger than in the twentieth century; for 1990–2090 it amounts to 0.20–0.37 m. This wide range results from systematic uncertainty in modelling of climate change and of heat uptake by the ocean. The AOGCMs agree that sea-level rise is expected to be geographically non-uniform, with some regions experiencing as much as twice the global average, and others practically zero, but they do not agree about the geographical pattern. The lack of agreement indicates that we cannot currently have confidence in projections of local sea-level changes, and reveals a need for detailed analysis and intercomparison in order to understand and reduce the disagreements. Received: 1 September 2000 / Accepted: 20 April 2001     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

Why Mt Etna?

The Etna volcano is located in an apparently anomalous position on the hinge zone of the Apennines subduction and its Na-alkaline geochemistry does not favour a magma source from the deep slab as indicated for the Aeolian K-alkaline magmatism. The steeper dip of the regional foreland monocline at the front of the Apennines in the Ionian Sea than in Sicily, implies a larger rollback of the subduction hinge in the Ionian Sea. Moreover, the lengthening of the Apennines arc needs extension parallel to the arc. Therefore, the larger southeastward subduction rollback of the Ionian lithosphere with respect to the Hyblean plateau in Sicily, should kinematically produce right-lateral transtension and a sort of vertical 'slab window' which might explain (i) the Plio-Pleistocene alkaline magmatism of eastern Sicily (e.g. the Etna volcano) and (ii) the late Pliocene to present right lateral transtensional tectonics and seismicity of eastern Sicily. The area of transfer of different dip and rollback occurs along the inherited Mesozoic passive continental margin between Sicily and the oceanic Ionian Sea, i.e. the Malta escarpment.     

Grain boundary diffusion in enstatite

We have investigated grain boundary diffusion rates in enstatite by heating single crystals of quartz packed in powdered San Carlos olivine (Mg_0.90Fe_0.10)₂SiO₄ at controlled oxygen fugacities in the range 10^?5.7 to 10^?8.7?atm and temperatures from 1350° to 1450?°C for times from 5 to 100?h at 1?atm total pressure. Following the experiments, the thickness of the coherent polycrystalline reaction rim of pyroxene that had formed between the quartz and olivine was measured using backscatter scanning imaging in the electron microprobe. Quantitative microprobe analysis indicated that the composition of this reaction phase is (Mg_0.92Fe_0.08)₂Si₂O₆. The rate of growth of the pyroxene increases with increasing temperature, is independent of the oxygen fugacity, and is consistent with a parabolic rate law, indicating that the growth rate is controlled by ionic diffusion through the pyroxene rim. Microstructural observations and platinum marker experiments suggest that the reaction phase is formed at the olivine-pyroxene interface, and is therefore controlled by the diffusion of silicon and oxygen. The parabolic rate constants determined from the experiments were analyzed in terms of the oxide activity gradient across the rim to yield mean effective diffusivities for the rate-limiting ionic species, assuming bulk transport through the pyroxene layer. These effective diffusivities are faster than the lattice diffusivities for the slowest species (silicon) calculated from creep experiments, but slower than measured lattice diffusivities for oxygen in enstatite. Thus, silicon grain boundary diffusion is most likely to be the rate-limiting process in the growth of the pyroxene rims. Also, as oxygen transport through the pyroxene rims must be faster than silicon transport, diffusion of oxygen along the grain boundaries must be faster than through the lattice. The grain boundary diffusivity for silicon in orthopyroxenite is then given by D¯gbSiδ=(3.3±3.0)×10^?9f0.0O₂e^?400±65/RT?m³s^?1, where the activation energy for diffusion is in kJ/mol, and δ is the grain boundary width in m. Calculated growth rates for enstatite under these conditions are significantly slower than predicted by an extrapolation from similar experiments performed at 1000?°C under high pressure (hydrous) conditions by Yund and Tullis (1992), perhaps due to water-enhancement of diffusion in their experiments.