On the Moment Method for the Modeling of Cloud Microphysics in Rarefied Turbulent Atmospheres: II. Stochastic Coagulation

The previously proposed (Rodin, 2002) method for calculating the microphysical properties of spatially inhomogeneous rarefied aerosol media with mixing using the lowest-order moments of the size distribution is generalized to particle coagulation. We show that when the problem is formulated in terms of moments, all of the solutions admitted by the stochastic coagulation equation lie within a narrow range whose boundaries can be determined by means of quadratic programming. We discuss the choice of an optimal solution within this range and compare the moment method with the results of our computations by the classical finite-difference method using a model of photochemical aerosols in Titan's atmosphere as an example. The moment method allows the efficiency of microphysical computations to be significantly increased by using precomputed low-dimension interpolation tables. It can be used to construct self-consistent models for the globular circulation of planetary atmospheres.     

Isobaric cooling and anti-clockwise P–T paths in the Variscan Odenwald Crystalline Complex, Germany

Mineral assemblages in Al₂O₃‐rich, SiO₂‐ and K₂O‐poor metapelitic rocks from the western Odenwald Crystalline Complex (Variscan Mid‐German Crystalline Rise, southern Germany) include corundum, spinel, cordierite, sillimanite, garnet and staurolite. Quartz is absent from almost all samples. Therefore, the applicability of conventional geothermobarometry is very limited or even impossible. Detailed petrographic investigation on selected samples permits inference of the sequence of appearance and disappearance of several mineral assemblages. The recognition of such partial re‐equilibration stages and their associated mineral assemblages, together with mineral stabilities predicted from KFMASH pseudosections, enables the determination of the pressure‐temperature (P–T) trajectories experienced by these rocks during the Variscan metamorphism. The rocks were metamorphosed under low‐P/high‐T conditions and underwent an anti‐clockwise P–T evolution. A pressure increase from about 2 kbar to 4 ± 0.5 kbar was accompanied by heating. Peak metamorphic conditions were reached at pressures of 4 ± 0.5 kbar and temperatures of at least 640 °C, probably even higher. The retrograde evolution is characterised by near‐isobaric cooling from ≥ 640 °C to approximately 550 °C. The rocks underwent the anti‐clockwise evolution in a subduction‐related magmatic arc setting. The close spatial association of the low‐P/high‐T rocks with recently discovered metabasic eclogites in the eastern part of the Odenwald Crystalline Complex may indicate a fossil paired metamorphic belt in the Central European Variscides.     

A new stereo‐analytical method for determination of removal blocks in discontinuous rock masses

The paper provides a new stereo‐analytical method, which is a combination of the stereographic method and analytical methods, to separate finite removable blocks from the infinite and tapered blocks in discontinuous rock masses. The methodology has applicability to both convex and concave blocks. Application of the methodology is illustrated through examples. Addition of this method to the existing block theory procedures available in the literature improves the capability of block theory in solving practical problems in rock engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Predicting the polarization of the microwave background from the WMAP temperature maps

In this paper, we study how to predict the polarization of the Cosmic Microwave Background (CMB) using knowledge of only the temperature (intensity) and the cross-correlation between temperature and polarization. We derive a “Wiener prediction” method and apply it to the Wilkinson Microwave Anisotropy Probe (WMAP) all-sky CMB temperature maps and to the MAXIMA field.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.