Woodleigh, Southern Carnarvon Basin, Western Australia: history of discovery, Late Devonian age, and geophysical and morphometric evidence for a 120 km-diameter impact structure

The discovery of the Woodleigh impact structure, first identified by R. P. lasky, bears a number of parallels with that of the Chlcxulub impact structure of K-T boundary age, underpinning complications inherent in the study of buried impact structures by geophysical techniques and drilling. Questions raised in connection with the diameter of the Woodleigh impact structure reflect uncertainties in criteria used to define original crater sizes in eroded and buried impact structures as well as limits on the geological controls at Woodleigh. The truncation of the regional Ajona - Wandagee gravity ridges by the outer aureole of the Woodleigh structure, a superposed arcuate magnetic anomaly along the eastern part of the structure, seismic-reflection data indicating a central > 37 km-diameter dome, correlation of fault patterns between Woodleigh and less-deeply eroded impact structures (Ries crater, Chesapeake Bay), and morphometric estimates all indicate a final diameter of 120 km. At Woodleigh, pre-hydrothermal shock-induced melting and diaplectic transformations are heavily masked by pervasive alteration of the shocked gneisses to montmorillonite-dominated clays, accounting for the high MgO and low K2O of cryptocrystalline components. The possible contamination of sub-crater levels of the Woodlelgh impact structure by meteoritic components, suggested by high Ni, Co, Cr, Ni/ Co and Ni/Cr ratios, requires further siderophile element analyses of vein materials. Although stratigraphic age constraints on the impact event are broad (post-Middle Devonian to pre-Early Jurassic) high-temperature (200-250 degrees C) pervasive hydrothermal activity dated by K-Ar isotopes of illite - smectite indicates an age of 359 +/- 4 Ma. To date neither Late Devonian crater fill, nor impact ejecta fallout units have been identified, although metallic meteoritic ablation spherules of a similar age have been found in the Conning Basin.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Submarine groundwater discharge and associated chemical input to a coastal sea

This paper presents a theoretical model of flow and chemical transport processes in subterranean estuaries (unconfined brackish groundwater aquifers at the ocean-land interface). The model shows that groundwater circulation and oscillating flow, caused by wave setup and tide, may constitute up to 96% of submarine groundwater discharge (SGWD) compared with 4% due to the net groundwater discharge. While these local flow processes do not change the total amount of land-derived chemical input to the ocean over a long period (e.g., yearly), they induce fluctuations of the chemical transfer rate as the aquifer undergoes saltwater intrusion. This may result in a substantial increase in chemical fluxes to the ocean over a short period (e.g., monthly and by a factor of 20 above the averaged level), imposing a possible threat to the marine environment. These results are essentially consistent with the experimental findings of Moore [1996] and have important implications for coastal resources management.

Groundwater waves in a coastal aquifer: A new governing equation including vertical effects and capillarity

Groundwater waves, that is, water table fluctuations, are a natural phenomenon in coastal aquifers. They represent an important part of the interaction between the ocean and aquifer and affect the mass exchange between them. This paper presents a new groundwater wave equation. Because it includes the effects of vertical flows and capillarity, the new equation is applicable to both intermediate-depth aquifers and high-frequency waves. Compared with the wave equation derived by Nielsen ed al. [1997], the present equation provides a closer representation of groundwater waves. In particular, it predicts high-frequency water table fluctuations as observed in the field. A validation of the new equation has been carried out by comparing the analytical solutions to it with predictions from direct simulations using the numerical model SUTRA. The effects of various physical parameters and their relative importance are also discussed.

Tidal fluctuations in a leaky confined aquifer: Dynamic effects of an overlying phreatic aquifer

Tidal fluctuations in a leaky confined coastal aquifer are damped significantly due to leakage into an overlying phreatic aquifer. Jiao and Tang [1999] presented an analytical solution to a simple model describing this phenomenon. Their solution assumes that the tidal fluctuations in the overlying phreatic aquifer are negligible (i.e,, a static phreatic aquifer), Here we examine dynamic effects of the overlying aquifer based on a new approximate analytical solution. The numerical results indicate that the dynamic effects can be significant for a relatively large leakage and a high transmissivity of the phreatic aquifer.

Sudden drawdown and drainage of a horizontal aquifer

Drainage of a saturated horizontal aquifer following a sudden drawdown is reanalyzed using the Boussinesq equation. The effect of the finite length of the aquifer is considered in detail. An analytical approximation based on a superposition principle yields a very good estimate of the outflow when compared to accurate numerical solutions. An illustration of the new analytical approach to analyze basin-scale field data is used to demonstrate possible field applications of the new solution.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Infiltration from surface and buried point sources: The average wetting water content

The assumption in analytical solutions for flow from surface and buried point sources of an average water content, (θ) over bar, behind the wetting front is examined. Some recent work has shown that this assumption fitted some field data well. Here we calculated (θ) over bar using a steady state solution based on the work by Raats [1971] and an exponential dependence of the diffusivity upon the water content. This is compared with a constant value of (θ) over bar calculated from an assumption of a hydraulic conductivity at the wetting front of 1 mm day(-1) and the water content at saturation. This comparison was made for a wide range of soils. The constant (θ) over bar generally underestimated (θ) over bar at small wetted radii and overestimated (θ) over bar at large radii. The crossover point between under and overestimation changed with both soil properties and flow rate. The largest variance occurred for coarser texture soils at low-flow rates. At high-flow rates in finer-textured soils the use of a constant (θ) over bar results in underestimation of the time for the wetting front to reach a particular radius. The value of (θ) over bar is related to the time at which the wetting front reaches a given radius. In coarse-textured soils the use of a constant value of (θ) over bar can result in an error of the time when the wetting front reaches a particular radius, as large as 80% at low-flow rates and large radii.

Neogene and Quaternary rollback evolution of the Tyrrhenian Sea, the Apennines, and the Sicilian Maghrebides

Reconstruction of the evolution of the Tyrrhenian Sea shows that the major stage of rifting associated with the opening of this basin began at similar to10 Ma. It involved two episodes of back arc extension, which were induced by the rollback of a west dipping subducting slab. The first period of extension (10-6 Ma) was prominent in the northern Tyrrhenian Sea and in the western part of the southern Tyrrhenian Sea. The second period of extension, mainly affected the southern Tyrrhenian Sea, began in the latest Messinian (6-5 Ma) and has been accompanied by subduction rollback at rates of 60-100 km Myr(-1). Slab reconstruction, combined with paleomagnetic and paleogeographic constraints, indicates that in the central Apennines, the latest Messinian (6-5 Ma) arrival of a carbonate platform at the subduction zone impeded subduction and initiated a slab tear and major strike-slip faults. These processes resulted in the formation of a narrow subducting slab beneath the Ionian Sea that has undergone faster subduction rollback and induced extreme rates of back arc extension.

Investment prioritization based on broadscale spatial budgeting to meet downstream targets for suspended sediment loads

On the basis of a spatially distributed sediment budget across a large basin, costs of achieving certain sediment reduction targets in rivers were estimated. A range of investment prioritization scenarios were tested to identify the most cost-effective strategy to control suspended sediment loads. The scenarios were based on successively introducing more information from the sediment budget. The relationship between spatial heterogeneity of contributing sediment sources on cost effectiveness of prioritization was investigated. Cost effectiveness was shown to increase with sequential introduction of sediment budget terms. The solution which most decreased cost was achieved by including spatial information linking sediment sources to the downstream target location. This solution produced cost curves similar to those derived using a genetic algorithm formulation. Appropriate investment prioritization can offer large cost savings because the magnitude of the costs can vary by several times depending on what type of erosion source or sediment delivery mechanism is targeted. Target settings which only consider the erosion source rates can potentially result in spending more money than random management intervention for achieving downstream targets. Coherent spatial patterns of contributing sediment emerge from the budget model and its many inputs. The heterogeneity in these patterns can be summarized in a succinct form. This summary was shown to be consistent with the cost difference between local and regional prioritization for three of four test catchments. To explain the effect for the fourth catchment, the detail of the individual sediment sources needed to be taken into account.