Computer simulation studies of Crystallization under confinement conditions

Crystallization under confinement conditions is a very important process in geochemistry and geophysics. Computer simulations of fluids in nanometer scale pore spaces can provide a unique microscopic insight into the structure, dynamics and forces arising from the crystallization process. We discuss in this paper molecular dynamics computer simulations of crystallization in pores of nanometer dimensions. The crystallization pressure due to the freezing of a model of Argon in a nanopore is computed using molecular dynamics simulations. We also investigate the influence of pore geometry in determining the dynamics of confined fluids, as well as mass separation in binary mixtures. It turns out that the pore geometry reveals itself as an important variable, leading to 1) new mechanisms for fast diffusion in confined spaces, and 2) accumulation of solute in specific regions inside the pore.     

Genesis of monazite and Y zoning in garnet from the Black Hills, South Dakota

The paragenesis of monazite in metapelitic rocks from the contact aureole of the Harney Peak Granite, Black Hills, South Dakota, was investigated using zoning patterns of monazite and garnet, electron microprobe dating of monazite, bulk-rock compositions, and major phase mineral equilibria. The area is characterized by low-pressure and high-temperature metamorphism with metamorphic zones ranging from garnet to sillimanite zones. Garnet porphyroblasts containing euhedral Y annuli are observed from the garnet to sillimanite zones. Although major phase mineral equilibria predict resorption of garnet at the staurolite isograd and regrowth at the andalusite isograd, textural and mass balance analyses suggest that the formation of the Y annuli is not related to the resorption-and-regrowth of garnet having formed instead during garnet growth in the garnet zone. Monazite grains in Black Hills pelites were divided into two generations on the basis of zoning patterns of Y and U: monazite 1 with low-Y and -U and monazite 2 with high-Y and -U. Monazite 1 occurs in the garnet zone and persists into the sillimanite zone as cores shielded by monazite 2 which starts to form in the andalusite zone. Pelites containing garnet porphyroblasts with Y annuli and monazite 1 with patchy Th zoning are more calcic than those with garnet with no Y annuli and monazite with concentric Th zoning. Monazite 1 is attributed to breakdown of allanite in the garnet zone, additionally giving rise to the Y annuli observed in garnet. Monazite 2 grows in the andalusite zone, probably at the expense of garnet and monazite 1 in the andalusite and sillimanite zones. The ages of the two different generations of monazite are within the precision of chemical dating of electron microprobe. The electron microprobe ages of all monazites from the Black Hills show a single ca. 1713 Ma population, close to the intrusion age of the Harney Peak Granite (1715 Ma). This study demonstrates that Y zoning in garnet and monazite are critical to the interpretation of monazite petrogenesis and therefore monazite ages.     

Preservation/exhumation of ultrahigh-pressure subduction complexes

Ultrahigh-pressure (UHP) metamorphic terranes reflect subduction of continental crust to depths of 90–140 km in Phanerozoic contractional orogens. Rocks are intensely overprinted by lower pressure mineral assemblages; traces of relict UHP phases are preserved only under kinetically inhibiting circumstances. Most UHP complexes present in the upper crust are thin, imbricate sheets consisting chiefly of felsic units ± serpentinites; dense mafic and peridotitic rocks make up less than  10% of each exhumed subduction complex. Roundtrip prograde–retrograde P–T paths are completed in 10–20 Myr, and rates of ascent to mid-crustal levels approximate descent velocities. Late-stage domical uplifts typify many UHP complexes.

Sialic crust may be deeply subducted, reflecting profound underflow of an oceanic plate prior to collisional suturing. Exhumation involves decompression through the P–T stability fields of lower pressure metamorphic facies. Scattered UHP relics are retained in strong, refractory, watertight host minerals (e.g., zircon, pyroxene, garnet) typified by low rates of intracrystalline diffusion. Isolation of such inclusions from the recrystallizing rock matrix impedes back reaction. Thin-aspect ratio, ductile-deformed nappes are formed in the subduction zone; heat is conducted away from UHP complexes as they rise along the subduction channel. The low aggregate density of continental crust is much less than that of the mantle it displaces during underflow; its rapid ascent to mid-crustal levels is driven by buoyancy. Return to shallow levels does not require removal of the overlying mantle wedge. Late-stage underplating, structural contraction, tectonic aneurysms and/or plate shallowing convey mid-crustal UHP décollements surfaceward in domical uplifts where they are exposed by erosion. Unless these situations are mutually satisfied, UHP complexes are completely transformed to low-pressure assemblages, obliterating all evidence of profound subduction.     

Analysis of a rigid footing lying on three-layered soil using the finite difference method

A numerical procedure is described for the analysis of the vertical deformation and the stress distribution of the strip footings on layered soil media. Three layers of soil with different stiffness are considered with the middle soil layer the thinnest and most stiff layer. The soil media is discretized and using the theory of elasticity, the governing differential equations are obtained in terms of vertical and horizontal displacements. These equations along with appropriate boundary and continuity conditions are solved by using the finite difference method. The vertical and horizontal displacements, strains and stresses are found at various nodes in the soil media. Parametric studies are carried out to study the effect of the placement depth of the middle soil layer, the relative ratios of the moduli of deformation of the soil layers on the vertical displacement of the footing and the vertical stress distribution. These studies reveal that the middle thin but very stiff layer acts like a plate and redistributes the stresses on the lower soft soil layer uniformly. The displacement on the top and bottom of the middle soil layer is almost the same showing that the compression of the middle layer is negligible as it is very stiff.     

An approximate projection of availability of the fresh groundwater resources in the South West district of NCT Delhi, India: a case study

The future availability and sustainability of fresh groundwater resources in the South West district of the national capital territory (NCT) Delhi, India, have been projected. Due to a rapid decline in groundwater level and quality, the district has been required by the Government of India to regulate development of groundwater resources. Shallow groundwater is mostly saline and water resources in the area are limited. The methodology applied here involves microzonation of the district in terms of thickness of fresh groundwater and then quantification of present and future availability of freshwater in different freshwater zones, including tentative timescales. The calculation method has been aided by data on historic trends in water level at representative groundwater monitoring stations, located either in fresh groundwater zones or near to them. It is estimated that the presently available 481 million m³ of resources will be reduced to 374 million m³ by year 2007 and to 303 million m³ by the year 2012, and by the year 2022 the district will have only 176 million m³ of available fresh groundwater resources.     

On the Importance of Geolithological Features for the Estimate of the Site Response: The Case of Catania Metropolitan Area (Italy)

The evaluation of seismic site response in the urban area of Catania was tackled by selecting test areas having peculiar lithological and structural features, potentially favourable to large local amplifications of ground motion. The two selected areas are located in the historical downtown and in the northern part of Catania where the presence of a fault is evident. Site response was evaluated using spectral ratio technique taking the horizontal- to-vertical component ratio of ambient noise. Inferences from microtremor measurements are compared with results from synthetic accelerograms and response spectra computed at all drillings available for this area. Such method is particularly suitable in urban areas where the nature of the outcropping geological units is masked by city growth and anthropic intervention on the surface geology. The microtremor H/V spectral ratios evaluated at soft sites located within the downtown profile tend to be smaller than that usually reported in the literature for such soils. A tendency for amplifications to peaks near 2 Hz is observed only in some sites located on recent alluvial deposits. Evidences for amplifications of site effects (frequency range 4–8 Hz) were observed in the sampling sites located on the fault, with a rapid decrease of spectral amplitude just a few tenths of metres away from the discontinuity. Numerical simulations evidenced the importance of geolithological features at depth levels even greater than 20–30 m. Besides this, the results strongly confirm the importance of the subsurface geological conditions, in the estimate of seismic hazard at urban scale.     

Comparative Seismic Risk Studies for German Earthquake Regions on the Basis of the European Macroseismic Scale EMS-98

Within the framework of recent research projects, basic tools for GIS-based seismic risk assessment technologies were developed and applied to the building stock and regional particularities of German earthquake regions. Two study areas are investigated, being comparable by the level of seismic hazard and the hazard-consistent scenario events (related to mean return periods of 475, 2475 and 10000 years). Significant differences exist with respect to the number of inhabitants, the grade and extent of urbanisation, the quality and quantity of building inventory: the case study of Schmölln in Eastern Thuringia seems to be representative for the majority of smaller towns in Germany, the case study of Cologne (Köln) stands for larger cities. Due to the similarities of hazard and scenario intensities, the considerable differences do not only require proper decisions concerning the appropriate methods and acceptable efforts, they enable conclusions about future research strategies and needs for disaster reduction management. Not least important, results can sharpen the focus of public interest. Seismic risk maps are prepared for different scenario intensities recognising the scatter and uncertainties of site-dependent ground motion and also of the applied vulnerability functions. The paper illustrates the impact of model assumptions and the step-wise refinements of input variables like site conditions, building stock or vulnerability functions on the distribution of expected building damage within the study areas. Furthermore, and in contrast to common research strategies, results support the conclusion that in the case of stronger earthquakes the damage will be of higher concentration within smaller cities like Schmölln due to the site-amplification potential and/or the increased vulnerability of the building stock. The extent of damage will be pronounced by the large number of masonry buildings for which lower vulnerability classes have to be assigned. Due to the effect of deep sedimentary layers and the composition of building types, the urban centre of Cologne will be less affected by an earthquake of comparable intensity.     

Modeling of Beams on Reinforced Granular Beds

The paper describes a mechanical model for estimating the flexural response of a strip footing, supporting a column (imposing a concentrated load), resting on a compacted granular bed overlying a reinforcement layer for example, geogrids, geomats etc. below which lies a loose soil deposit. The footing is idealized as a beam and the reinforcing element is assumed to have finite bending stiffness and negligible frictional resistance. The upper and lower soil layers are idealized by a series of linear and discrete springs (Winkler springs) of different stiffness values. To find the response of such a model the governing differential equations have been derived and expressed in a nondimensional form. A closed form analytical solution of the same has been obtained subjected to appropriate boundary conditions. Using the present approach the resulting solution for a degenerated case of a long beam is found to be identical to the same of Hetenyi (1946, Beams on elastic foundations, University of Michigan press, Ann Arbor, MI). Parametric studies reveal that the ratio of flexural rigidity of upper and lower beam and the ratio of stiffness of the upper and lower soil layers affect significantly the response of the foundation.     

River bank erosion control by soil nailing

A study has been done for analysing soil nailed cuts with circular type wedge failure by friction circle method. Various parameters such as nail length, nail diameter, nail inclination, wall inclination and angle of internal friction of soil have been considered to determine the factor of safety of nailed open cuts. The study shows that for cohesionless soil nailed cut, factor of safety increase with increase of parameters like angle of internal friction of soil, length of nail (L) versus height of cut (H) ratio, cohesion of soil and nail inclination (upto 15°) with horizontal. The study revealed that nails grouted with cement perform better than driven nails. A case study further confirms the analytical findings. Received 7 October