Long‐term final void salinity prediction for a post‐mining landscape in the Hunter Valley,New South Wales,Australia

Opencast mining alters surface and subsurface hydrology of a landscape both during and post‐mining. At mine closure, following opencast mining in mines with low overburden to coal ratios, a void is left in the final landform. This final void is the location of the active mine pit at closure. Voids are generally not infilled within the mines' lifetime, because of the prohibitive cost of earthwork operations, and they become post‐mining water bodies or pit lakes. Water quality is a significant issue for pit lakes. Groundwater within coal seams and associated rocks can be saline, depending on the nature of the strata and groundwater circulation patterns. This groundwater may be preferentially drawn to and collected in the final void. Surface runoff to the void will not only collect salts from rainfall and atmospheric fallout, but also from the ground surface and the weathering of fresh rock. As the void water level rises, its evaporative surface area increases, concentrating salts that are held in solution. This paper presents a study of the long term, water quality trends in a post‐mining final void in the Hunter Valley, New South Wales, Australia. This process is complex and occurs long term, and modelling offers the only method of evaluating water quality. Using available geochemical, climate and hydrogeological data as inputs into a mass‐balance model, water quality in the final void was found to increase rapidly in salinity through time (2452 to 8909 mg l⁻¹ over 500 years) as evaporation concentrates the salt in the void and regional groundwater containing high loads of salt continues to flow into the void. Copyright © 2004 John Wiley & Sons, Ltd.     

Observations of mine seismicity in the eastern Wasatch Plateau,Utah, U.S.A.: A possible case of implosional failure

In the summer of 1984, a three-dimensional, high-resolution microearthquake network was operated in the vicinity of two coal mines beneath Gentry Mountain in the eastern Wasatch Plateau, Utah. During a six-week period, approximately 3,000 seismic events were observed of which the majority were impulsive, higher frequency (>10 Hz), short duration (<2–3 sec) events probably associated with the caving of the roof from a longwall operation. In contrast, 234 of the largest located events appeared to occur predominantlybeneath the mines to a depth of 2 to 3 km consistent with previous studies. The magnitudes of these events ranged from less thanM _c 0 to 1.6. In addition to the unusual depths of these latter events, an anomalous aspect displayed by the events was an apparent dilatational focal mechanism suggesting a non-double-couple, possibly implosional source. Implosional events have been observed in other studies of mine seismicity; however, the generally inadequate instrumental coverage of the focal sphere has cast some doubt on the validity of such mechanisms. Previously suggested source mechanisms for such implosional events have included tensional failure through strata collapse, and a shear-implosional displacement mechanism. Shear failure must be involved in the failure process of the Gentry Mountain implosional events as evidenced by well-defined shear waves in the observed seismograms. Simultaneous monitoring in the East Mountain coal mining area to the south by the University of Utah revealed typical shear failure events mixed with implosional events. The observed double-couple, reverse focal mechanisms at East Mountain were similar to mechanisms determined in previous studies and a composite focal mechanism determined in this study for a sequence outside the mining areas. This suggested that the shear events within the mining areas are being influenced by the regional tectonic stress field. Thus in addition to the seismic events associated with caving of the roof from the longwall operation, there appear to be at least two other types of mining-induced seismic events occurring in the eastern Wasatch Plateau, both submine in origin: (1) events characterized by apparent non-double-couple possibly implosional focal mechanisms and well-defined shear waves; and (2) shear events, which are indistinguishable from tectonic earthquakes and may be considered mining triggered earthquakes. The small mining-induced stress changes that occur beyond a few hundred meters from the mine workings suggest both types of seismic events are occurring on critically stressed, pre-existing zones of weakness. Topography, overburden, method of mining, and mine configuration also appear to be significant factors influencing the occurrence of the implosional submine events.     

The rank evaluation of South Brazilian Gondwana coals on the basis of different chemical and physical parameters

Rank evaluation of South Brazilian Gondwana coals according to the American (American Society for Testing and Materials, ASTM) and German (Deutsche Industrie Normen, DIN) standard specifications or on the basis of Alpern's Universal Classification has too often led to conflicting results.When any attempt is made to determine the rank of Brazilian coals by using methods established for the Northern Hemisphere coals, two sorts of difficulties arise: namely, (1) the reflectence measurements do not indicate the real degree of coalification of the organic matter because of the liptinitic impregnations of the vitrinites, and (2) random reflectance values (R_ran, moisture, volatile-matter and carbon contents as well as the calorific value do not correlate when plotted on the ASTM and DIN classification schemes.In addition, when rank determinations are carried out by using fluorescence measurements and geochemical analyses the results are at variance in relation to those determined by the conventional chemical and physical parameters.It is suggested that the rank evaluation of South Brazilian coals is far more meaningful if carried out by reflectance measurements on vitrinites along with fluorescence measurements on sporinites, alginites and coal extracts.     

Pahoehoe and aa in Hawaii: volumetric flow rate controls the lava structure

The historical records of Kilauea and Mauna Loa volcanoes reveal that the rough-surfaced variety of basalt lava called aa forms when lava flows at a high volumetric rate (>5–10 m³/s), and the smooth-surfaced variety called pahoehoe forms at a low volumetric rate (<5–10 m³/s). This relationship is well illustrated by the 1983–1990 and 1969–1974 eruptions of Kilauea and the recent eruptions of Mauna Loa. It is also illustrated by the eruptions that produced the remarkable paired flows of Mauna Loa, in which aa formed during an initial short period of high discharge rate (associated with high fountaining) and was followed by the eruption of pahoehoe over a sustained period at a low discharge rate (with little or no fountaining). The finest examples of paired lava flows are those of 1859 and 1880–1881. We attribute aa formation to rapid and concentrated flow in open channels. There, rapid heat loss causes an increase in viscosity to a threshold value (that varies depending on the actual flow velocity) at which, when surface crust is torn by differential flow, the underlying lava is unable to move sufficiently fast to heal the tear. We attribute pahoehoe formation to the flowage of lava at a low volumetric rate, commonly in tubes that minimize heat loss. Flow units of pahoehoe are small (usually <1 m thick), move slowly, develop a chilled skin, and become virtually static before the viscosity has risen, to the threshold value. We infer that the high-discharge-rate eruptions that generate aa flows result from the rapid emptying of major or subsidiary magma chambers. Rapid near-surface vesiculation of gas-rich magma leads to eruptions with high discharge rates, high lava fountains, and fast-moving channelized flows. We also infer that long periods of sustained flow at a low discharge rate, which favor pahoehoe, result from the development of a free and unimpeded pathway from the deep plumbing system of the volcano and the separation of gases from the magma before eruption. Achievement of this condition requires one or more episodes of rapid magma excursion through the rift zone to establish a stable magma pathway.     

Enhancements of CO and O3 from burnings in sugar cane fields

The manual harvest of sugar cane requires the burning of its foliage. This burning has strongly increased in Brazil after the National Alcohol Program was started which substituted automobile gasoline engines for alcohol engines. Presently, the source strength per unit area of this rural pollution is comparable to the well-known biomass burning source in Amazonia. The observed concentrations of CO and O₃ in the rural area of the state of São Paulo during the 1988 burning season were twice as large as those reported from an aircraft experiment of 1985 for biomass burnings of the tropical rain forest. Results are reported from airplane measurements and from three fixed ground stations. Mixing ratios of ozone and carbon monoxide in the height range below 6 km are normally less than 40 and 100 ppbv, (parts per billion by volume), respectively, in the absence of burnings. A strong O₃ and CO layer was observed during the burning period with peak concentrations of 80 ppbv of ozone and 580 ppbv of CO at about 2 km. The concentrations of CH₄ and CO₂ were also large, 1756 ppbv and 409 ppmv, respectively, at 1500 m. During the dry season period of the experiment, the ground based O₃ average diurnal variations obtained at the rural sites were practically identical to the typical urban variation observed at São José dos Campos, with daytime ozone values between 45 and 60 ppbv. A second three-day airplane excursion to the surgar cane fields in the wet season of 1989 has produces results to be contrasted with the dry (burning) season of 1988 and 1989. Carbon monoxide concentrations were below 100 ppbv at all heights and ozone concentrations were around 30–40 ppbv. The maximum daytime concentrations at the ground station Bauru was 25 ppbv of O₃, and at Jaboticabal it was 35 ppbv of O₃, only one half of what was observed in the dry season.Universidade Estadual de São Paulo.     

Seasonal soil water storage changes beneath central Amazonian rainforest and pasture

Evaporation and infiltration were compared for tropical rainforest and pasture, near to Manaus, Brazil from October 1990 to February 1992 using measurements of soil water storage over a depth of 2 m. The soil is a clayey oxisol of low water available capacity. In both of the dry seasons studied, the maximum change in soil water storage in the forest was 154 mm and in the pasture it was 131 and 112 mm. Similar behaviour of the soil water reservoir below forest and pasture in the wet season implied that differences in evaporation and drainage were small. In the dry season, soil water storage behaviour in the upper metre of the soil was similar but there were marked differences in the second metre. The pasture took up little water from below 1.5 m but the forest appeared to utilise all of the available water in the 2 m profile in both seasons.

The water balance of the 2 m profile showed that the pasture evaporation rate was equal to that of the forest until storage had decreased 80 mm from the maximum. There was then a decline in pasture evaporation rate to 1.2 mm day⁻¹ as the storage decreased by a further 50 mm. In contrast, the forest uptake rate remained above 3.5 mm day⁻¹ until storage had decreased 140 mm from the maximum (within 15 mm of the extraction limit), before declining abruptly to less than 1.5 mm day⁻¹. There was strong evidence that the forest was able to abstract water from depths greater than 3.6 m.

Spatial variability of soil water storage was significantly greater beneath the pasture than beneath the forest, particularly following rainfall events in the dry season. This was largely the result of redistribution of rainfall as local surface runoff. There was no evidence of redistribution or runoff in the forest.     

A model for adjustment of differential gravity measurements with simultaneous gravimeter calibration

 A mathematical model is proposed for adjustment of differential or relative gravity measurements, involving simultaneously instrumental readings, coefficients of the calibration function, and gravity values of selected base stations. Tests were performed with LaCoste and Romberg model G gravimeter measurements for a set of base stations located along a north–south line with 1750 mGal gravity range. This line was linked to nine control stations, where absolute gravity values had been determined by the free-fall method, with an accuracy better than 10 μGal. The model shows good consistence and stability. Results show the possibility of improving the calibration functions of gravimeters, as well as a better estimation of the gravity values, due to the flexibility admitted to the values of the calibration coefficients. Received: 15 November 1999 / Accepted: 31 October 2000     

Pitfalls in Nonlinear Inversion

— We discuss and illustrate graphically with simple 2-D problems, four common pitfalls in geophysical nonlinear inversion. The first one establishes that the Lagrange multiplier, used to incorporate a priori information in the geophysical inverse problem, should be the largest value still compatible with a reasonable data fitting. This procedure should be used only when the interpreter is sure about the importance of the a priori information used to stabilize the inverse problem relative to the geophysical observations. Because this is rarely the case, the user should use the smallest Lagrange multiplier still producing stable solutions. The second pitfall is an attempt to automatically estimate the Lagrange multiplier by decreasing it along the iterative process used to solve the nonlinear optimization problem. Consequently, at the last iteration, the Lagrange multiplier may be so small that the problem may become ill-posed and any computed solution in this case is meaningless. The third pitfall is related to the incorporation of a priori information by a technique known as “Jumping.” This formulation, from the viewpoint of the class of Acceptable Gradient Methods, is incomplete and may lead to a premature halt in the iteration, and, consequently, to solutions far from the true one. Finally, the fourth pitfall is an inadequate convergence criterion which stops the iteration when the data misfit drops just below the noise level, irrespective of the fact that the functional to be minimized may not have attained its minimum. This means that the a priori information has not been completely incorporated, so that this stopping criterion partially neutralizes the effect of the stabilizing functional, and opens the possibility of obtaining unstable, meaningless estimates.     

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.     

Two Brazilian archaeological sites investigated by GPR: Serrano and Morro Grande

The application of the ground penetrating radar (GPR) at two archaeological sites, Serrano and Morro Grande, situated in Araruama County, Rio de Janeiro, Brazil, aids the study of a prehistoric indigenous culture, associated with the “Tupinambá” that inhabited the region during prehistoric times.The archaeological remains of the study area are mainly characterized by pottery artifacts for several uses, including funerary urns, which were buried within layers of sand and clay. Several profiles were acquired using a RAMAC system, with a 200 and 400 MHz frequency antennae. At the Serrano site, the profiles were acquired around some partially exposed pottery shards, due to sand exploitation. The resultant profiles provided a response model to guide the interpretation of new profiles acquired at other sites in the area, which present similar characteristics.The results showed the great importance of the dielectric permittivity contrast which exists between the targets and the host media, in order for possibly significant features to be identified in radar data.