Hydrodynamic response of subduction zones to seismic activity: A case study for the Costa Rica margin

Thermal anomalies in tectonically active areas are often attributed to sub-seafloor fluid circulation and faulting mechanisms, particularly in subduction zones where the largest thrust earthquakes occur. Postseismic fluid flow is enabled by the poroelastic response of the fault system to the earthquake's strain field, as well as by the rupturing of permeability barriers in the vicinity of the fault zone. We investigated the relative importance of these mechanisms on postseismic pore-pressure diffusion and advective heat transport in the subduction zone setting. A two-dimensional numerical fluid flow and heat transport model was developed for the Costa Rica subduction zone offshore of the Nicoya Peninsula. The flow and transport model was coupled with an earthquake strain model to quantify the effects of coseismic strain and permeability enhancement on fluid pressures and temperatures within the Costa Rica margin. Coseismic changes in pore pressure and postseismic pore-pressure diffusion were found to be sensitive to the compressibility of the porous medium, and patterns of pore-pressure recovery were more complex than that predicted by theoretical faulting models. Coseismic contraction and extension of the crust produced high fluid pressures close to the fault, while the inflow of fluid from depth increased fluid pressures several years following the simulated fault slip. Crustal deformation alone was not observed to perturb the temperature field. Laterally extensive permeability increases of two orders of magnitude along the décollement were required to produce small changes in heat flow. Local permeability changes in the upper slope region of least five orders of magnitude were necessary to noticeably affect heat flow. The results of the numerical simulations may help to refine conceptual faulting models and provide guidance for locating long-term hydrologic monitoring sites at Costa Rica and other subduction zones.     

Optimization strategies for sediment reduction practices on roads in steep,forested terrain

Many forested steeplands in the western United States display a legacy of disturbances due to timber harvest, mining or wildfires, for example. Such disturbances have caused accelerated hillslope erosion, leading to increased sedimentation in fish‐bearing streams. Several restoration techniques have been implemented to address these problems in mountain catchments, many of which involve the removal of abandoned roads and re‐establishing drainage networks across road prisms. With limited restoration funds to be applied across large catchments, land managers are faced with deciding which areas and problems should be treated first, and by which technique, in order to design the most effective and cost‐effective sediment reduction strategy. Currently most restoration is conducted on a site‐specific scale according to uniform treatment policies. To create catchment‐scale policies for restoration, we developed two optimization models – dynamic programming and genetic algorithms – to determine the most cost‐effective treatment level for roads and stream crossings in a pilot study basin with approximately 700 road segments and crossings. These models considered the trade‐offs between the cost and effectiveness of different restoration strategies to minimize the predicted erosion from all forest roads within a catchment, while meeting a specified budget constraint. The optimal sediment reduction strategies developed by these models performed much better than two strategies of uniform erosion control which are commonly applied to road erosion problems by land managers, with sediment savings increased by an additional 48 to 80 per cent. These optimization models can be used to formulate the most cost‐effective restoration policy for sediment reduction on a catchment scale. Thus, cost savings can be applied to further restoration work within the catchment. Nevertheless, the models are based on erosion rates measured on past restoration sites, and need to be updated as additional monitoring studies evaluate long‐term basin response to erosion control treatments. Copyright © 2006 John Wiley & Sons, Ltd.     

Repeated Aqueous Flooding from the Cerberus Fossae: Evidence for Very Recently Extant, Deep Groundwater on Mars

The geomorphology and topography of the Cerberus Plains region of Mars show three spatially and temporally distinct, young, aqueous flood channel systems. Flood geomorphology in each of these channels, as seen in Mars Orbiter Camera images, consists of streamlined forms, longitudinal lineations, and a single occurrence of transverse dunes, features similar to those in the flood-carved terrain of the Channeled Scabland in the northwestern United States. As additional geomorphic evidence of flooding, small cones (interpreted as phreatic) are found preferentially in the channels or at their distal ends. Glaciers, lava flows, and CO₂-charged density flows are each inconsistent with these geomorphic features. Mars Orbiter Laser Altimeter data show two of the three channel systems (Athabasca Valles and an unnamed northern channel system) emanating from the Cerberus Fossae; we suggest that the third channel system (Marte Vallis) also originated at the fissures. The discharges for two of the three systems (Athabasca Valles and Marte Vallis) have been estimated from surface topography to have been on the order of 10⁶ m³/s. Crater counts indicate that the channels are not only young (extreme Late Amazonian), but also were carved asynchronously. Geomorphic evidence suggests that two of the channels (Athabasca and Marte Valles) experienced more than one flood. Emanation from volcanotectonic fissures instead of chaotic terrain distinguishes these Cerberus Plains channels from the larger, older circum-Chryse channels. Groundwater must have collected in a liquid state prior to flood onset to flow at the estimated discharge rates. Lack of large-scale subsidence near the channels' origination points along the Cerberus Fossae indicates that this groundwater was at least several kilometers deep.     

O/O mapping and hydrogeology of a short-lived (≈10 years) fumarolic (>500°C) meteoric–hydrothermal event in the upper part of the 0.76 Ma Bishop Tuff outflow sheet, California

¹⁸O/¹⁶O data from the 200-m-thick, 0.76 Ma Bishop Tuff outflow sheet provide evidence for a vigorous, short-lived (≈10 years), high-temperature, fumarolic meteoric–hydrothermal event. This is proved by: (1) the juxtaposition in the upper, partially welded Bishop Tuff of low-¹⁸O groundmass/glass (δ¹⁸O=−5 to +3) with coexisting quartz and feldspar phenocrysts having magmatic δ¹⁸O values (+8.7±0.3; +7.5±0.3); and (2) the fact that these kinds of ¹⁸O/¹⁶O signatures correlate very well with morphological features and mapped zones of fumarolic activity. Profiles of δ¹⁸O with depth in the Bishop Tuff within the fumarole area define a 40- to 50-m-thick, low-¹⁸O, stratigraphic zone that is sandwiched between the essentially unwelded near-surface portion of the tuff and an underlying, densely welded black tuff that displays magmatic ¹⁸O/¹⁶O values. Shallow-dipping columnar joints and other fumarolic features (i.e., subhorizontal tubular conduits and steep fissures) correlate very well with these pervasively devitrified, low-¹⁸O zones. The base of the low-¹⁸O zone is extremely sharp (3‰ per meter) and is located directly above the transition from partially welded tuff to densely welded black tuff. The observed average whole-rock ¹⁸O-depletions within this low-¹⁸O zone are about 6–7‰, requiring meteoric water/rock ratios in excess of 0.24 in mass units. Rainfall on the surface of the tuff would not have been high enough to supply this much H₂O in the short lifetime of fumarolic activity, suggesting that some recharge must have been from groundwater flow through the upper part of the tuff, above the sloping (1°–5°) top of the impermeable lower zone. This is compatible with the observation that the fumarolic areas roughly correlate with the preeruptive regional drainage pattern. Some of this recharge may in part have been from the lake that filled Long Valley caldera, which was dammed by the Bishop Tuff up to the level of this boundary between the partially and densely welded zones (≈7000 ft, the elevation of the highest Long Valley Lake shorelines). Gazis et al. had previously shown that the 2.8-Ma intracaldera Chegem Tuff from the Caucasus Mountains exhibits exactly the same kind of ¹⁸O-signature that we have correlated with fossil fumaroles in the Bishop Tuff outflow sheet. Although not recognized as such by McConnell et al.; ¹⁸O/¹⁶O data from drill-hole samples from the intracaldera Bishop Tuff in Long Valley also display this characteristic ¹⁸O signature (i.e., analogous δ¹⁸O-depth profiles, as well as low-¹⁸O groundmass coexisting with high-¹⁸O feldspar phenocrysts). This fumarolic ¹⁸O/¹⁶O signature is observed to much greater depths (≈650–750 m) in the intracaldera tuffs (≈1500 m thick) than it is in the ≈200-m-thick Bishop Tuff outflow sheet (≈80 m depth).     

An assessment of recent trophic changes in Windermere South Basin (England) based on diatom remains and fossil pigments

Cyanobacterial carotenoids and diatom remains have been analyzed in recent sediments from the Windermere South Basin (WSB) to study the trophic evolution experienced by the lake. Dates in the top 30 cm were specifically established through radionuclide (²¹⁰Pb and¹³⁷Cs) analyses. Diatom stratigraphy shows dominance of the centric diatomsCyclotella comensis andC. radiosa and several benthic taxa in the early postglacial. This indicates oligotrophy in the WSB during that period. This assemblage was replaced by another dominated by the diatomAsterionella formosa in the 1870's, as has been established from the²¹⁰Pb dating. From that date onwards, the lake underwent a progression towards eutrophy, indicated by the progressive increase inAulacoseira subarctica (c. 1930's),Fragilaria crotonensis (c. 1943), and more recently, of the centricsStephanodiscus parvus (c. 1971) andCyclotella meneghiniana (1988).Carotenoid stratigraphy reveals the differences between different sections of the core. Oscillaxanthin and myxoxanthophyll had very low records in the early and medium parts of the core, but increased fromc. 1950's, showing peaks atc. 1967, 1979 and 1987. Some of these peaks indicated a differential abundance ofOscillatoria, and are matched to those observed directly during the ongoing monitoring of the phytoplankton of the lake.The coincidence between the historic appearance of diatoms associated with nutrient-rich waters and the enhanced carotenoid occurrence suggest a common response to phosphorus enrichment, and that the progressive change towards eutrophy has been accentuated during the last twenty-five years.     

The influence of urbanization of sinkhole development in central Pennsylvania

The karsted limestone valleys of central Pennsylvania contain two populations of sinkholes. Solution sinkholes occur in the Champlainian limestone units along the margins of the valleys. Solution sinkholes are permanent parts of the landscape and, although a nuisance to construction, do not present other problems. The second population is the suffosional or soil-piping sinkholes These occur on all carbonate rock units including the Beekmantown and Gatesburg dolomites that comprise the two principal carbonate aquifers in the valley. Suffosional sinkholes are the principal land-use hazard. Suffosional sinkholes are transient phonomena. They occur naturally but are exacerbated by runoff modifications that accompany urbanization Suffosional sinkholes are typically 1.5–2.5 m in diameter depending on soil thickness and soil type. The vertical transport of soil to form the void space and soil arch that are the precursors to sinkhole collapse is through solutionally widened fractures and cross-joints and less often through large vertical openings in the bedrock. The limited solution development on the dolomite bedrock combined with soil thickness, seldom greater than 2 m, limits the size of the sinkholes. All aspects of suffosional sinkhole development are shallow processes: transport, piping, void and arch formation, and subsequent collapse take place usually less than 10 m below the land surface Factors exacerbating sinkhole development include pavement, street, and roof runoff which accelerates soil transport Such seemingly minor activities as replacing high grass and brush with mowed grass is observed to accelerate sinkhole development. Dewatering of the aquifer is not a major factor in this region     

The effect of a shallow low viscosity zone on the apparent compensation of mid-plate swells

A simple model for mid-plate swells is that of convection in a fluid which has a low viscosity layer lying between a rigid bed and a constant viscosity region. Finite element calculations have been used to determine the effects of the viscosity contrast, the layer thickness and the Rayleigh number on the flow and on the perceived compensation mechanism for the resulting topographic swell. As the viscosity decreases in the low viscosity zone, the effective local Rayleigh number for the top boundary layer of the convecting cell increases. Also, because the lower viscosity facilitates greater velocities in the low viscosity zone, the low viscosity layer produces proportionally greater horizontal flow near the conducting lid, causing the base of the conducting lid to appear like a free boundary. The change in the local Rayleigh number and in the effective boundary condition both cause the top boundary layer to thin. Through a Green's function analysis, we have found that the low viscosity zone damps the response of the surface topography to the temperature anomalies at depth, whereas it causes the gravity and geoid response functions to change sign at depth counteracting the positive contributions from the shallower temperature variations. By increasing the viscosity contrast, the conbined effects of the thinning of the boundary layer and the behaviour of the response functions allow the apparent depth of compensation to become arbitrarily small. Therefore, shallow depths of compensation cannot be used to argue against dynamic support of mid-plate swells. Furthermore, we compared the distribution of the effective compensating densities, which is used to obtain the geoid, to that of Pratt compensation, which is often used to calculate the depth of compensation from geoid and topography data for mid-plate swells. For all of our calculations including those with no low viscosity layer, the effective gravitational mass distribution is more complex than assumed in simple Pratt models, so that the Pratt models are not an appropriate gauge of the compensation mechanism.     

The origin of Samoa: new evidence from Sr, Nd, and Pb isotopes

We report Sr, Nd and Pb isotope ratios and parent and daughter element concentrations in 34 volcanic rocks from Samoa. The highly undersaturated post-erosional volcanics, which have erupted in Recent to Historic time along a 250-km-long fissure, have isotopic compositions that define fields distinct from those of the tholeiitic to alkalic lavas of the older Samoan shield volcanoes. Most shield lavas have²⁰⁶Pb/²⁰⁴Pb of 18.9–19.4,⁸⁷Sr/⁸⁶Sr of 0.7045–0.7055 and⁸⁷Sr/⁸⁶Sr (to 0.7075). In general, isotopic compositions of the shield lavas are similar to those of the Marquesas and Society Islands. Post-erosional samples have lower²⁰⁶Pb/²⁰⁴Pb and¹⁴³Nd/¹⁴⁴Nd and higher⁸⁷Sr/⁸⁶Sr than most shield lavas. The most striking feature of the post-erosional data is a negative correlation between²⁰⁷Pb/²⁰⁴Pb and²⁰⁶Pb/²⁰⁴Pb. This suggests that post-erosional lavas are derived from mixtures of the shield source and a high-²⁰⁷Pb/²⁰⁴Pb,⁸⁷Sr/⁸⁶Sr, low-²⁰⁶Pb/²⁰⁴Pb and¹⁴³Nd/¹⁴⁴Nd post-erosional source which may contain recycled ancient sediment. This enriched mantle domain may also underlie the Ontong-Java and Manihiki Plateaus to the north and west. Although both the Samoan shield and post-erosional lavas show chemical characteristics often associated with mantle plumes, only the shield volcanism can plausibly be related to a plume. The post-erosional eruptions appear to be the result of flexure and rifting due to plate bending at the northern termination of the Tonga Trench.