Approximation of interfacial properties in multiphase porous medium systems

We investigate the estimation of interfacial areas, curvatures, and common curve lengths in multiphase porous medium systems. Algorithms are developed to obtain estimates of these quantities based upon a variety of potential data sources and estimation approaches. The accuracy of the derived approximations are evaluated as a function of the data type and resolution of the data. The methods advanced improve upon standard approaches now in use and show excellent accuracy at resolutions on the order of five lattice points per minimum radius of curvature of the object being resolved. Finally, we suggest a promising class of extensions that could lead to further improvements in the accuracy of such methods.     

Investigating the ancient Martian magnetic field using microwaves

The new microwave palaeointensity technique has been used to investigate samples from the Martian meteorite Nakhla. This technique is a promising new way to obtain absolute palaeointensity information regarding the ancient Martian magnetic field as recorded by the Martian meteorites. Assuming that a part of the magnetic remanence is of thermal origin and originating on Mars the two samples studied yield estimates of 4 μT for the Martian magnetic field at 1.35 Ga.     

The physical volcanology of the 1600 eruption of Huaynaputina, southern Peru

Volcán Huaynaputina is a group of four vents located at 16°36'S, 70°51'W in southern Peru that produced one of the largest eruptions of historical times when ~11 km³ of magma was erupted during the period 19 February to 6 March 1600. The main eruptive vents are located at 4200 m within an erosion-modified amphitheater of a significantly older stratovolcano. The eruption proceeded in three stages. Stage I was an ~20-h sustained plinian eruption on 19-20 February that produced an extensive dacite pumice fall deposit (magma volume ~2.6 km³). Throughout medial-distal and distal parts of the dispersal area, a fine-grained plinian ashfall unit overlies the pumice fall deposit. This very widespread ash (magma volume ~6.2 km³) has been recognized in Antarctic ice cores. A short period of quiescence allowed local erosion of the uppermost stage-I deposits and was followed by renewed but intermittent explosive activity between 22 and 26 February (stage II). This activity resulted in intercalated pyroclastic flow and pumice fall deposits (~1 km³). The flow deposits are valley confined, whereas associated co-ignimbrite ash fall is found overlying the plinian ash deposit. Following another period of quiescence, vulcanian-type explosions of stage III commenced on 28 February and produced crudely bedded ash, lapilli, and bombs of dense dacite (~1 km³). Activity ceased on 6 March. Compositions erupted are predominantly high-K dacites with a phenocryst assemblage of plagioclase>hornblende>biotite>Fe-Ti oxides-apatite. Major elements are broadly similar in all three stages, but there are a few important differences. Stage-I pumice has less evolved glass compositions (~73% SiO₂), lower crystal contents (17-20%), lower density (1.0-1.3 g/cm³), and phase equilibria suggest higher temperature and volatile contents. Stage-II and stage-III juvenile clasts have more evolved glass (~76% SiO₂) compositions, higher crystal contents (25-35%), higher densities (up to 2.2 g/cm³), and lower temperature and volatile contents. All juvenile clasts show mineralogical evidence for thermal disequilibrium. Inflections on a plot of log thickness vs area^1/2 for the fall deposits suggest that the pumice fall and the plinian ash fall were dispersed under different conditions and may have been derived from different parts of the eruption column system. The ash appears to have been dispersed mainly from the uppermost parts of the umbrella cloud by upper-level winds, whereas the pumice fall may have been derived from the lower parts of the umbrella cloud and vertical part of the eruption column and transported by a lower-altitude wind field. Thickness half distances and clast half distances for the pumice fall deposit suggests a column neutral buoyancy height of 24-32 km and a total column height of 34-46 km. The estimated mass discharge rate for the ~20-h-long stage-I eruption is 2.4᎒⁸ kg/s and the volumetric discharge rate is ~3.6᎒⁵ m³/s. The pumice fall deposit has a dispersal index (Hildreth and Drake 1992) of 4.4, and its index of fragmentation is at least 89%, reflecting the dominant volume of fines produced. Of the 11 km³ total volume of dacite magma erupted in 1600, approximately 85% was evacuated during stage 1. The three main vents range in size from ~70 to ~400 m. Alignment of these vents and a late-stage dyke parallel to the NNW-SSE trend defined by older volcanics suggest that the eruption initiated along a fissure that developed along pre-existing weaknesses. During stage I this fissure evolved into a large flared vent, vent 2, with a diameter of approximately 400 m. This vent was active throughout stage II, at the end of which a dome was emplaced within it. During stage III this dome was eviscerated forming the youngest vent in the group, vent 3. A minor extra-amphitheater vent was produced during the final event of the eruptive sequence. Recharge may have induced magma to rise away from a deep zone of magma generation and storage. Subsequently, vesiculation in the rising magma batch, possibly enhanced by interaction with an ancient hydrothermal system, triggered and fueled the sustained Plinian eruption of stage I. A lower volatile content in the stage-II and stage-III magma led to transitional column behavior and pyroclastic flow generation in stage II. Continued magma uprise led to emplacement of a dome which was subsequently destroyed during stage III. No caldera collapse occurred because no shallow magma chamber developed beneath this volcano.     

Implications of model uncertainty for the mapping of hillslope‐scale soil erosion predictions

An integrated modelling approach (MIRSED) which utilizes the process‐based soil erosion model WEPP (Water Erosion Prediction Project) is presented for the assessment of hillslope‐scale soil erosion at five sites throughout England and Wales. The methodology draws upon previous uncertainty analysis of the WEPP hillslope soil erosion model by the authors to qualify model results within an uncertainty framework. A method for incorporating model uncertainty from a range of sources is discussed as a first step towards using and learning from results produced through the GLUE (Generalized Likelihood Uncertainty Estimation) technique. Results are presented and compared to available observed data, which illustrate that levels of uncertainty are significant and must be taken into account if a meaningful understanding of output from models such as WEPP is to be achieved. Furthermore, the collection of quality, observed data is underlined for two reasons: as an essential tool in the development of soil erosion modelling and also to allow further constraint of model uncertainty. Copyright © 2001 John Wiley & Sons, Ltd.     

Yearly changes in the seasonal frequency and duration of short-term acid pulses in some Nova Scotia, Canada streams

 The seasonal frequency and duration of low pH events at three sites located in southwestern Nova Scotia were analyzed using measured and synthetic daily pH data. The basins varied in size from 0.3 to 300 km² and were subject to frequent snowmelt events in winter and spring, as well as occasional runoff events during summer and fall. Results showed that, in order to fill in missing data from periods where collection was interrupted, statistical approximations using discharge were not totally acceptable, as generated data consistently missed the extreme values measured. Despite a lack of totally accurate event pH estimates for periods where daily data were missing, analysis showed that in this region, low pH episodes can occur year-round including summer. The highest probabilities of low pH episodes nevertheless occurred in the winter and spring when snowmelts were frequent. Received: 12 November 1999 · Accepted: 18 July 2000     

Sources, distribution and storage of heavy metals in the Río Pilcomayo, Bolivia

The Río Pilcomayo rises on the Cerro Rico de Potosí precious metal-polymetallic tin deposits of Bolivia, and flows in a southeasterly direction for ca. 600 km to Bolivia's southern border with Argentina. Mining of the Potosí deposits has occurred continuously since 1545, generating large quantities of waste materials in the headwater of the basin. In addition, a tailings dam breach at the Porco mine in 1996 released an estimated 235 000 m³ of tailings and fluid into the upper reaches of the Río Pilaya, the largest tributary to the Pilcomayo.Concentrations of As, Sb, Cd, Cu, Pb, Hg, Ag, Tl and Zn in contemporary channel sediments upstream of the Pilaya confluence are significantly elevated above background values. Elevated levels appear to be associated with pyrite- and other sulphide mineral-bearing tailings materials transported more than 200 km downstream of the Potosí mines. Significant downstream declines in elemental concentrations occur within 15 km, and again between 150 and 200 km, from the mines. The initial decrease in concentrations is due to the rapid dilution of nearly pure tailings effluent released to the river from milling facilities near Potosí. The latter decrease results from a combination of geomorphic processes including the storage of sediment-borne metals within the channel bed and the influx of ‘clean’ sediment from several large tributaries. Downstream of the Pilaya confluence, concentrations of Cu, Pb, Hg and Zn are only slightly elevated above background values, and Ag, Cd, Sb and Tl cannot be distinguished from background levels. These data suggest that while the Porco tailings spill may have had a significant short-term impact on sediment and water quality along the lower reaches of the Río Pilcomayo, its longer-term impacts were limited. Metals stored and eroded from alluvial deposits of historical age in upstream reaches appear to be an important source of metals to the river today. An additional, and perhaps more significant source, is the release of tailings effluent to the river from modern milling operations. The transport of these contaminants downstream of Icla (203 km from Potosí) appears to be restricted by aggradational processes occurring in the vicinity of Puente Sucre. In addition, downstream of the confluence of the Río Pilaya, inputs of large amounts of ‘clean’ sediment have caused dilution of the metal contaminants. Data from other studies where similar geomorphic processes have occurred suggest that the metals in the upper Pilcomayo may eventually be moved downvalley as the aggradational processes are reversed and channel stabilisation occurs. Thus, the most significant impacts of metal contamination may not be realised in downstream areas for decades.