Terrestrial, planetary and interplanetary magnetic fields

A discussion meeting at Burlington House on 10 March 2000 featured eight speakers and an audience of about 50. Andy Smith reports.     

Catastrophe revisited – disastrous flow failures of mine and municipal solid waste

Catastrophic flow failures have occurred with alarming frequency in mine tailings dams and dumps of discards and other mine waste. In recent years, catastrophic flow failures have also occurred in dumps of municipal solid waste and even in what were intended to be carefully controlled and well-engineered landfills. Apart from the environmental devastation caused by these flows, they are also dangerous to human life and society. Examples include the Buffalo Creek disaster in the USA in 1972 (Vick, Planning, Design and Analysis of Tailings Dams, Wiley, 1983) that killed 118 people, made 4000 homeless and destroyed 50 million US dollars worth of property and facilities. The flow slide that occurred in the Umraniye– Hekimbasi refuse dump in Turkey in 1993 (Kocasoy and Curi, Waste Management and Research, 13, 1995, 305) killed 39 people, destroying their homes in the process. This paper will briefly review some of the more notable flow slides in waste materials, analyzing conditions necessary for a flow failure to occur and pointing to ways of preventing this type of failure by a combination of sound design and operating procedures.     

A two degree-of-freedom earthquake model with static/dynamic friction

Can a simple multi-block-spring model with total symmetry make interesting predictions for fault behaviour? Our model consists of a symmetric, slowly driven, two degree-of-freedom block-spring system with static/dynamic friction. The simple friction law and slow driving rate allow the state of this fourth order system to be described between slip events by a single variable, the difference in the stretch of the driving springs. This stretch difference measures the locked-in stress and is closely related to fault stress inhomogeneity. In general,smoothing is not observed. A spatially homogeneous stress state is found to almost always be unstable, in that the system tends toward an inhomogeneous state after many slip events. The system evolves either to a cycle that alternates between two types of earthquakes, or to a cycle with repeating but identical asymmetric earthquakes. One type of alternating earthquake solution is structurally unstable, which implies a great sensitivity to model perturbations. For this simple model, spatial asymmetry necessarily occurs, despite the symmetry in the model, thus suggesting that spatial structure in seismicity patterns may be a consequence of earthquake dynamics, not just fault heterogeneity.     

Modeling outflow from the Ernest Henry Fe oxide Cu–Au deposit: implications for ore genesis and exploration

The Ernest Henry Fe oxide Cu–Au (IOCG) deposit (>ca. 1.51 Ga) is hosted by breccia produced during the waning stages of an evolving hydrothermal system that formed a number of tens of metres to a kilometre scale, pre- and syn-ore alteration halos, although no demonstrable patterns have been attributed to fluids expelled through the outflow zones. However, the recognition of a population of hypersaline fluid inclusions representing the ‘spent’ fluids after Cu–Au deposition at Ernest Henry provides the basis to model the geochemical characteristics of the deposit's outflow zones. Geochemical modeling at 300 °C was undertaken at both high and low fluid/rock ratios via FLUSH models involving three host rock types: (1) granite, (2) calc–silicate rock, and (3) graphitic schist. In models run at high fluid/rock ratios, all rock types are essentially fluid-buffered, and produce an albite–quartz–hematite–barite-rich assemblage, although in low fluid–rock environments, the pH, redox, and geochemical character of the host rock exerts a greater influence on the mineralogy of the alteration assemblages (e.g., andradite, Fe–chlorite, and magnetite). Significant sulphide mineralization was predicted in graphitic schist where sphalerite occurred in both low- and high-porosity models, which indicates the possibility of an association between high-temperature IOCG mineralization and lower temperature base metal mineralization.Cooling experiments (from 300 to 100 °C) using the ‘spent fluids’ predict early high-T (300–200 °C) Na-, Ca-, Fe-, and Mn-rich, magnetite-bearing hydrothermal associations, whereas with cooling to below 200 °C, and with progressive fluid–rock interaction, the system produces rhodochrosite-bearing, hematite–quartz–muscovite–barite-rich assemblages. These results show that the radical geochemical and mineralogical changes associated with cooling and progressive fluid influx are likely to be accompanied by major transformations in the geophysical expression (e.g., spectral and magnetic character) of the alteration in the outflow zone, and highlight the potential link between magnetite- and hematite-bearing IOCG hydrothermal systems.     

Upscaling unsaturated hydraulic parameters for flow through heterogeneous anisotropic sediments

We compare two methods for determining the upscaled water characteristics and saturation-dependent anisotropy in unsaturated hydraulic conductivity from a field-scale injection test. In both approaches an effective medium approximation is used to reduce a porous medium of M textures to an equivalent homogenous medium. The first approach is a phenomenological approach based on homogenization and assumes that moisture-based Richards’ equation can be treated like the convective–dispersive equation (CDE). The gravity term, dK_z(θ)/d(θ), analogous to the vertical convective velocity in the CDE, is determined from the temporal evolution of the plume centroid along the vertical coordinate allowing calculation of an upscaled K_z(θ). As with the dispersion tensor in the CDE, the rate of change of the second spatial moment in 3D space is used to calculate the water diffusivity tensor, D(θ), from which an upscaled K(θ) is calculated. The second approach uses the combined parameter scale inverse technique (CPSIT). Parameter scaling is used first to reduce the number of parameters to be estimated by a factor M. Upscaled parameters are then optimized by inverse modeling to produce an upscaled K(θ) characterized by a pore tortuosity–connectivity tensor, L. Parameters for individual textures are finally determined from the optimized parameters by inverse scaling using scale factors determined a priori. Both methods produced upscaled K(θ) that showed evidence of saturation dependent anisotropy. Flow predictions with the STOMP simulator, parameterized with upscaled parameters, were compared with field observations. Predictions based on the homogenization method were able to capture the mean plume behavior but could not reproduce the asymmetry caused by heterogeneity and lateral spreading. The CPSIT method captured the effects of heterogeneity and anisotropy and reduced the mean squared residual by nearly 90% compared to local-scale and upscaled parameters from the homogenization method. The Pacific Northwest National Laboratory is operated for the US Department of Energy by Battelle under Contract DE-AC05-76RL01830.     

Mesozoic–Tertiary exhumation history of the Altai Mountains, northern Xinjiang, China: New constraints from apatite fission track data

This study uses apatite fission track (FT) analysis to constrain the exhumation history of bedrock samples collected from the Altai Mountains in northern Xinjiang, China. Samples were collected as transects across the main structures related to Palaeozoic crustal accretion events. FT results and modeling identify three stages in sample cooling history spanning the Mesozoic and Tertiary. Stage one records rapid cooling to the low temperature part of the fission track partial annealing zone circa 70 ± 10 °C. Stage two, records a period of relative stability with little if any cooling taking place between 75 and 25–20 Ma suggesting the Altai region had been reduced to an area of low relief. Support for this can be found in the adjacent Junngar Basin that received little if any sediment during this interval. Final stage cooling took place in the Miocene at an accelerated rate bringing the sampled rocks to the Earth's surface. This last stage, linked to the far field effects of the Himalayan collision, most likely generated the surface uplift and relief that define the present-day Altai Mountains.