Formation of high-grade ignimbrites Part II. A pyroclastic suspension current model with implications also for low-grade ignimbrites

 Analogue experiments in part I led to the conclusion that pyroclastic flows depositing very high-grade ignimbrite move as dilute suspension currents. In the thermo–fluid–dynamical model developed, the degree of cooling of expanded turbulent pyroclastic flows dynamically evolves in response to entrainment of air and mass loss to sedimentation. Initial conditions of the currents are derived from column-collapse modeling for magmas with an initial H₂O content of 1–3 wt.% erupting through circular vents and caldera ring-fissures. The flows spread either longitudinally or radially from source up to a runout distance that increases with higher mass flux but decreases with higher gas content, temperature, bottom slope and coarser initial grain size. Progressive dilution by entrainment and sedimentation causes pyroclastic currents to transform into buoyant ash plumes at the runout distance. The ash plumes reach stratospheric heights and distribute 30–80% of the erupted material as widespread co-ignimbrite ash. Pyroclastic suspension currents with initial mass fluxes of 10⁷-10¹² kg/s can spread for tens of kilometers with only limited cooling, although they move as supercritical, strongly entraining currents for the eruption conditions considered here. With increasing eruption mass flux, cooling during passage through the fountain diminishes while cooling during flow transport increases. The net effect is that eruption temperature exerts the prime control on emplacement temperature. Pyroclastic suspension currents can form welded ignimbrite across their entire extent if eruption temperature is T_o>1.3^.T_mw, the minimum welding temperature. High eruption rates, a large fraction of fine ash, and a ring-fissure vent favor the formation of extensive high-grade ignimbrite. For very hot eruptions producing sticky, partially molten pyroclasts, analysis of particle aggregation systematics shows that factors favoring longer runout also favor more efficient aggregation, which reduces runout. As a result, very high-grade ignimbrites cannot spread more than a few tens of kilometers from their source. In cooler pyroclastic currents, particles do not aggregate, and the sedimentation process may involve re-entrainment of particles, which potentially leads to more extensive cooling and longer runout; such effects, however, are only significant when net erosion of substrate occurs. Model results can be employed to estimate mass flux and duration of ignimbrite eruptions from measured ignimbrite masses and aspect ratios. The model also provides an alternative explanation of the observed decrease in H/Lratios with ignimbrite mass. Received: 10 May 1998 / Accepted: 21 October 1998     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

The scattering of Rossby waves from finite abrupt topography

The scattering of first mode linear baroclinic Rossby waves by a top-hat ridge in a continuously stratified ocean, with Brunt-Väisälä frequency that decays exponentially with depth below a surface mixed layer, is the subject of this study. A numerical mode matching technique is used to calculate the transmission coefficients for the propagating modes over the ridge. It is found that the scattered field depends crucially upon the stratification. For example, when the majority of the density variation is confined to a thin thermocline, corresponding to a small e-folding scale, gamma ^?1, for the Brunt-Väisälä frequency, a large amount of the incident wave energy is reflected by a small amplitude ridge. Appreciable energy conversion between the propagating barotropic and baroclinic modes takes place in this case. An asymptotic analysis for a small amplitude ridge is presented that confirms these numerical results. In the limit gamma ^?1→ 0, it is demonstrated that the scattered field in the continuously stratified ocean model differs markedly from the two-layer solution. The latter does not exhibit appreciable reflection of the incident wave energy for a small amplitude ridge. In conclusion, the application of a two-layer ocean model to describe Rossby wave scattering by ridges in place of a continuously stratified model cannot be recommended.     

An empirical method for estimating travel times for wet volcanic mass flows

 Travel times for wet volcanic mass flows (debris avalanches and lahars) can be forecast as a function of distance from source when the approximate flow rate (peak discharge near the source) can be estimated beforehand. The near-source flow rate is primarily a function of initial flow volume, which should be possible to estimate to an order of magnitude on the basis of geologic, geomorphic, and hydrologic factors at a particular volcano. Least-squares best fits to plots of flow-front travel time as a function of distance from source provide predictive second-degree polynomial equations with high coefficients of determination for four broad size classes of flow based on near-source flow rate: extremely large flows (>1 000 000 m³/s), very large flows (10 000–1 000 000 m³/s), large flows (1000–10 000 m³/s), and moderate flows (100–1000 m³/s). A strong nonlinear correlation that exists between initial total flow volume and flow rate for "instantaneously" generated debris flows can be used to estimate near-source flow rates in advance. Differences in geomorphic controlling factors among different flows in the data sets have relatively little effect on the strong nonlinear correlations between travel time and distance from source. Differences in flow type may be important, especially for extremely large flows, but this could not be evaluated here. At a given distance away from a volcano, travel times can vary by approximately an order of magnitude depending on flow rate. The method can provide emergency-management officials a means for estimating time windows for evacuation of communities located in hazard zones downstream from potentially hazardous volcanoes. Received: 5 June 1997 / Accepted: 2 February 1998     

Melt-filled hybrid fractures in the oceanic mantle: Melt enhanced deformation during along-axis flow beneath a propagating spreading ridge axis

Mid-ocean ridges represent important locations for understanding the interactions between deformation and melt production, transport, and emplacement. Melt transport through the mantle beneath mid-ocean ridges is closely associated with deformation. Currently recognized transport and emplacement processes at ridges include: 1) dikes and sills filling stress-controlled fractures, 2) porous flow in a divergent flow field, 3) self-organizing porous dunite channels, and 4) shear zones. Our recent observations from the sub-oceanic mantle beneath a propagating ridge axis in the Oman ophiolite show that gabbronorite and olivine gabbro dikes fill hybrid fractures that show both shear and extensional components of strain. The magnitudes of shear strain recorded by the dikes are significant and comparable to the longitudinal extensions across the dikes. We suggest that the hybrid dikes form from the interactions between shear deformation and pressurized melt in regions of along-axis flow at mid-ocean ridges. The displacement across the dikes is kinematically compatible with high temperature flow recorded by plastic fabrics in host peridotites. Field observations and mechanical considerations indicate that the dikes record conditions of higher stress and lower temperature than those recorded by the plastic flow fabrics. The features of hybrid dikes suggest formation during progressive deformation as conditions changed from penetrative plastic flow to strain localization along melt-filled fractures. The combined dataset indicates that the dikes are formed during along-axis flow away from regions of diapiric upwelling at propagating ridge segments. Hybrid dikes provide a potentially powerful kinematic indicator and strain recorder and define a previously unrecognized mechanism of melt migration. Our calculations show that hybrid dikes require less melt pressure to form than purely tensile dikes and thus may provide a mechanism to tap melt reservoirs that are under-pressurized with respect to lithostatic pressure.     

Tectonics of the Easter plate

We present a new model for the Easter plate in which rift propagation has resulted in the formation of a rigid plate between the propagating and dying ridges. We use the distribution of earthquakes, eleven new focal mechanisms, and existing bathymetric and magnetic data to describe the tectonics of this area. Both the Easter-Nazca and Easter-Pacific Euler poles are sufficiently close to the Easter plate to cause rapid changes in rates and directions of motion along the boundaries. The east and west boundaries are propagating and dying ridges; the southwest boundary is a slow-spreading ridge and the northern boundary is a complex zone of convergent and transform motion.The Easter plate may reflect the tectonics of rift propagation on a large scale, where rigid plate tectonics requires boundary reorientation. We use simple schematic models to illustrate the general features and processes which occur at plates resulting from large-scale rift propagation.     

Complex kame belt morphology,stratigraphy and architecture

The development of glacier karst at the margins of melting ice sheets produces complex glaciofluvial sediment-landform assemblages that provide information on ice sheet downwasting processes. We present the first combined geomorphological, sedimentological and geophysical investigation of the Brampton Kame Belt, an important glaciofluvial depositional zone at the centre of the last British-Irish Ice Sheet. Ground-penetrating radar (GPR) data allow the broad scale internal architecture of ridges (eskers) and flat-topped hills (ice-walled lake plains) to be determined at four sites. In combination with sediment exposures, these provide information on lateral and vertical variations in accretion styles, depositional boundaries, and grain size changes. Building on existing work on the subject, we propose a refined model for the formation of ice-walled lake plains resulting from the evolution and collapse of major drainage axes into lakes as stable glacier karst develops during deglaciation. The internal structure of esker ridges demonstrates variations in sedimentation that can be linked to differences in ridge morphologies across the kame belt. This includes low energy flow conditions and multiple accretion phases identified within large S-N oriented esker ridges; and fluctuating water pressures, hyperconcentrated flows, and significant deformation within a fragmented SW–NE oriented esker ridge. In combination with updated geomorphological mapping, this work allows us to identify two main styles of drainage within the kame belt: (1) major drainage axes aligned broadly S-N that extend through the entire kame belt and collapsed into a chain of ice-walled lakes; and (2) a series of smaller, fragmented SW–NE aligned esker ridges that represent ice-marginal drainage as the ice sheet receded south-eastwards up the Vale of Eden. Our study demonstrates the importance of integrated geomorphological, sedimentological and geophysical investigations in order to understand complex and polyphase glaciofluvial sediment-landform assemblages. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Rotating flow over long shallow ridges

Abstract

The flow of a rotating homogeneous, incompressible fluid past a long ridge is investigated. An analysis is presented for flows in which E ? 1, Ro ～ E^½, H/D ～ E⁰, h/D ～ E^½ and cosα ～ E⁰ where E is the Ekman number, Ro the Rossby number, H/D the fluid depth to ridge width ratio, h/D the ridge height to ridge width ratio and α the angle between the free stream flow and a line perpendicular to the ridge axis. The analysis includes effects of the nonlinear inertial terms. Particular examples of a ridge of triangular cross section and a sinusoidal topography are investigated in some detail. Experiments are presented for a triangular ridge which are in good agreement with the theory.     

Depositional processes and gas pore pressure in pyroclastic flows: an experimental perspective

The depositional processes and gas pore pressure in pyroclastic flows are investigated through scaled experiments on transient, initially fluidized granular flows. The flow structure consists of a sliding head whose basal velocity decreases backwards from the front velocity (U _f) until onset of deposition occurs, which marks transition to the flow body where the basal deposit grows continuously. The flows propagate in a fluid-inertial regime despite formation of the deposit. Their head generates underpressure proportional to U _f ² whereas their body generates overpressure whose values suggest that pore pressure diffuses during emplacement. Complementary experiments on defluidizing static columns prove that the concept of pore pressure diffusion is relevant for gas-particle mixtures and allow characterization of the diffusion timescale (t _d) as a function of the material properties. Initial material expansion increases the diffusion time compared with the nonexpanded state, suggesting that pore pressure is self-generated during compaction. Application to pyroclastic flows gives minimum diffusion timescales of seconds to tens of minutes, depending principally on the flow height and permeability. This study also helps to reconcile the concepts of en masse and progressive deposition of pyroclastic flow units or discrete pulses. Onset of deposition, whose causes deserve further investigation, is the most critical parameter for determining the structure of the deposits. Even if sedimentation is fundamentally continuous, it is proposed that late onset of deposition and rapid aggradation in relatively thin flows can generate deposits that are almost snapshots of the flow structure. In this context, deposition can be considered as occurring en masse, though not strictly instantaneously.     

Petrologic consequences of rift propagation on oceanic spreading ridges

The production of anomalously differentiated lava compositions at several mid-ocean spreading centers can be attributed to magmatic processes associated with propagating rifts. The degree of differentiation attained by magmas beneath oceanic spreading ridges depends mainly on the balance between cooling rate and the supply rate of new magma to shallow chambers. Low supply rates and moderate cooling rates allow advanced degrees of closed-system fractionation to occur. High supply rates result in open systems in which magma compositions are buffered by frequent replenishment with new hot magma. Propagating rift tips are a special class of ridge-transform intersection in which the balance between cooling and supply rates is conducive to the development of advanced degrees of differentiation over an expanded length of ridge. This balance is affected by the spreading rate, the propagation rate of the rift, the length of the bounding transform and proximity to hotspots. Maximum compositional variability and maximum degree of differentiation occur within 50 km of propagating rift tips and subsequently diminish with increasing distance. Rifts that propagate through plates in directions approximating their absolute motion relative to the lower mantle are characterized by the presence of anomalously differentiated lavas over longer ridge segments than are rifts that propagate against their absolute motion. Geochemical anomalies may persist, though changing in degree and extent, for several million years on ridge segments that stop propagating. The concept of “magnetic telechemistry” is generally supported by our study, but in the vicinity of hotspots, magnetic anomaly amplitude may be controlled more by bathymetric and/or thickened magnetic layer effects than by geochemistry.     

Pyroclastic flows and surges over water: an example from the 1883 Krakatau eruption

Pyroclastic deposits from the 1883 eruption of Krakatau are described from areas northeast of the volcano on the islands of Sebesi, Sebuku, and Lagoendi, and the southeast coast of Sumatra. Massive and poorly stratified units formed predominantly from pyroclastic flows and surges that traveled over the sea for distances up to 80 km. Granulometric and lithologic characteristics of the deposits indicate that they represent the complement of proximal subaerial and submarine pyroclastic flow deposits laid down on and close to the Krakatau islands. The distal deposits exhibit a decrease in sorting coefficient, median grain size, and thickness with increasing distance from Krakatau. Crystal fractionation is consistent with the distal facies being derived from the upper part of gravitationally segregated pyroclastic flows in which the relative amount of crystal enrichment and abundance of dense lithic clasts diminished upwards. The deposits are correlated to a major pyroclastic flow phase that occurred on the morning of 27 August at approximately 10 a.m. Energetic flows spread out away from the volcano at speeds in excess of 100 km/h and traveled up to 80 km from source. The flows retained temperatures high enough to burn victims on the SW coast of Sumatra. Historical accounts from ships in the Sunda Straits constrain the area affected by the flows to a minimum of 4x10³ km². At the distal edge of this area the flows were relatively dilute and turbulent, yet carried enough material to deposit several tens of centimeters of tephra. The great mobility of the Krakatau flows from the 10 a.m. activity may be the result of enhanced runout over the sea. It is proposed that the generation of steam at the flow/seawater interface may have led to a reduction in the sedimentation of particles and consequently a delay in the time before the flows ceased lateral motion and became buoyantly convective. The buoyant distal edge of these ash-and steam-laden clouds lifted off into the atmosphere, leading to cooling, condensation, and mud rain.     

On the inadequacies of Long's model for steady two-dimensional stratified flows

Nonlinear analysis of two-dimensional steady flows with density stratification in the presence of gravity is considered. Inadequacies of Long's model for steady stratified flow over topography are explored. These include occurrence of closed streamline regions and waves propagating upstream. The usual requirements in Long's model of constant dynamic pressure and constant vertical density gradient in the upstream condition are believed to be the cause of these inadequacies. In this article, we consider a relaxation of these requirements, and also provide a systematic framework to accomplish this. As illustrations of this generalized formulation, exact solutions are given for the following two special flow configurations: the stratified flow over a barrier in an infinite channel; the stratified flow due to a line sink in an infinite channel. These solutions exhibit again closed-streamline regions as well as waves propagating upstream. The persistence of these inadequacies in the generalized Long's model appears to indicate that they are not quite consequences of the assumptions of constant dynamic pressure and constant vertical density gradient in Long's model, contrary to previous belief.

On the other hand, solutions admitted by the generalized Long's model show that departures from Long's model become small as the flow becomes more and more supercritical. They provide a nonlinear mechanism for the generation of columnar disturbances upstream of the obstacle and lead in subcritical flows to qualitatively different streamline topological patterns involving saddle points, which may describe the lee-wave-breaking process in subcritical flows and could serve as seats of turbulence in real flows. The occurrences of upstream disturbances in the presence of lee-wave-breaking activity described by the present solution are in accord with the experiments of Long (Long, R.R., “Some aspects of the flow of stratified fluids, Part 3. Continuous density gradients”, Tellus 7, 341--357 (1955)) and Davis (Davis, R.E., “The two-dimensional flow of a stratified fluid over an obstacle”, J. Fluid Mech. 36, 127–143 ()).     

Late-Pleistocene flank collapse triggered by dome growth at Tacaná volcano, México-Guatemala, and its relationship to the regional stress regime

During late Pleistocene time, the extrusion of an andesitic dome at the summit of Tacaná volcano caused the collapse of its northwestern flank. The stratocone collapse was nearly parallel to the σ _min stress direction suggesting that failure was controlled by the regional stress field. The event produced a debris avalanche that was channelized in the San Rafael River and moved 8 km downstream. The deposit covered a minimum area of 4 km², had a volume of 0.8 ± 0.5 km³, with an H/L (vertical drop to horizontal transport distance ratio) of ~0.35, defining a degree of mobility that is atypical for volcanic debris avalanches. The flank failure undermined the summit dome leading to its collapse and the generation of a series of block-and-ash flows that were emplaced in quick succession and covered the avalanche surface. The collapse event left a 600-m-wide summit amphitheatre with a 30-degree opening to the northwest, and >200 m thick debris that blocked the San Rafael River. Remobilization of this material produced debris flows that eroded the primary deposits and cascaded into the Coatán River. After the collapse, the activity of Tacaná continued with the emission of the Agua Zarca lava flow dated at 10 ± 6 ka (⁴⁰Ar/³⁹Ar), and pyroclastic surges dated at 10,610 + 330/−315 yr BP (¹⁴C), which provide a minimum age for the collapse event. During the Holocene, Tacaná has been very active producing explosive and effusive eruptions that ended with the extrusion of two summit domes that today occupy the amphitheatre. The 1950 and 1986 phreatic outbursts occurred along the Pleistocene collapse scar. Currently ~300,000 inhabitants live within a 35 km radius of Tacaná, and could conceivably be impacted by future events of similar magnitude.     

The formation of high-grade ignimbrites, I: Experiments on high- and low-concentration transport systems containing sticky particles

 High-grade ignimbrites are thought to be deposited by pyroclastic flows at temperatures exceeding minimum welding temperature or even solidus temperature. Corresponding pyroclastic-flow particles range from plastic to partially liquid and are able to aggregate or coalesce. This contrasts with particles in pyroclastic flows producing unwelded ignimbrite, which are capable of elastic grain interactions. The low aspect ratio and great areal extent of high-grade ignimbrites requires transport in a particulate state either by (a) high-concentration mass flow facilitated by fluidizing gas reducing internal friction, or by (b) expanded turbulent flow of low but downward increasing concentration. This paper presents experiments designed to investigate the effects of plastic to liquid particles on these two contrasting transport mechanisms. Gas fluidization experiments using polyethyleneglycole (PEG) powders heated above minimum sintering (T_ms) and melting (T_m) temperatures cover a wide range of fluidization velocities (U_mf>U_a>0.6·U_t) but are always in the bubbly fluidization regime similar to fluidized ignimbrite ash, where particle volume concentration outside the bubbles is high (≈10^–1). When the powders reach a critical temperature T_m≥T≥T_ms, defluidization by catastrophic particle aggregation immediately commences in both stationary and laterally moving fluidized beds as well as in experiments using mixtures of high- and low-T_m (≥30 wt.%) PEG powders, when T≥T_ms of the lower-T_m powder. This indicates that extended particulate transport at T≥T_ms is not possible at such high particle concentrations. In the turbulent flow experiments, liquid sprays of molten PEG or water, vertically injected into a high-Re (>10⁴) horizontal air flow, form a low-concentration (10^–5 to 10^–4) turbulent suspension current. Proximal formation of partially coalesced aggregates, which settle faster than individual particles, causes the measured downstream decay of sedimentation rate to be steeper than predicted by theory of single solid-particle sedimentation from turbulent suspensions. As particles become finer downstream and coalescence efficiency decreases in response to cooling, more distally formed aggregates become too small and rare to modify sedimentation-rate decay from that of suspension flows containing solid particles. The key difference between the two transport systems is particle concentration, C. Since particle collision rate R_coll∝C², collision rates in fluidized beds are so high that all particles immediately aggregate when coalescence efficiency (1≥E_coal≥0) is larger than 10^-3. Low-concentration suspensions, on the other hand, require much higher values of E_coal for significant aggregation to occur. Dilute pyroclastic flows will have higher particle volume fractions (≈10^–3) than the experimental currents, but then viscous pyroclasts should have lower coalescence efficiencies than PEG droplets. Experimental results thus support an expanded turbulent transport mechanism of pyroclastic flows generating extensive high-grade ignimbrite sheets. Received: 28 August 1996 / Accepted: 3 December 1997     

Plume capture by divergent plate motions: implications for the distribution of hotspots,geochemistry of mid-ocean ridge basalts,and estimates of the heat flux at the core–mantle boundary

The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (10⁶≤Ra≤10⁹) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe^1/2 and when λ→1, q∝(Peλ)^1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.     

Numerical modeling of the effects of water flow,sediment transport and vegetation growth on the spatiotemporal patterning of the ridge and slough landscape of the Everglades wetland

A numerical model has been developed to simulate the spatiotemporal patterning of the ridge and slough landscape in wetlands, characterized by crests (ridges) and valleys (sloughs) that are elongated parallel to the direction of water flow. The model formulation consists of governing equations for integrated surface water and groundwater flow, sediment transport, and soil accretion, as well as litter production by vegetation growth. The model simulations show how the spatial pattern self-organizes over time with the generation of ridges and sloughs through sediment deposition and erosion driven by the water flow field. The spatial and temporal distributions of the water depth, flow rates and sediment transport processes are caused by differential flow due to vegetation and topography heterogeneities. The model was parameterized with values that are representative of the Everglades wetland in the southern portion of the Florida peninsula in the USA. Model simulation sensitivity was tested with respect to numerical grid size, lateral vegetation growth and the rate of litter production. The characteristic wavelengths of the pattern in the directions along and perpendicular to flow that are simulated with this model develop over time into ridge and slough shapes that resemble field observations. Also, the simulated elevation differences between the ridges and sloughs are of the same order of those typically found in the field. The width of ridges and sloughs was found to be controlled by a lateral vegetation growth distance parameter in a simplified formulation of vegetation growth, which complements earlier modeling results in which a differential peat accretion mechanism alone did not reproduce observations of ridge and slough lateral wavelengths. The results of this work suggest that ridge and slough patterning occurs as a result of vegetation's ability to grow laterally, enhancing sediment deposition in ridge areas, balanced by increased sediment erosion in slough areas to satisfy flow continuity. The interplay between sediment transport, water flow and vegetation and soil dynamic processes needs to be explored further through detailed field experiments, using a model formulation such as the one developed in this work to guide data collection and interpretation. This should be one of the focus areas of future investigations of pattern formation and stability in ridge and slough areas.     

Interactions between boulders and aeolian ridges on Mars

Here we present new observations of two different interactions between aeolian ridges and boulder fields on Mars that provide insight into past wind conditions. First, an analysis of ridge and boulder interactions at two test sites in Proctor Crater and an additional site ~430 km to the northeast shows that ridge geometry can be affected by changes in surface roughness elements generated by boulder fields. Second, a detailed examination of some of the boulder fields found that individual boulders can generate multi-armed ‘wakes’ that have no clear proxy on Earth. The ridge/boulder dynamics suggest that transverse aeolian ridges acted as roughness during their development, indicating that they formed at a length scale greater than wind ripples. The boulder wakes seem to represent an unusual interaction between flow separation and pre-exiting ridges; why this pattern is not observed on Earth remains uncertain.     

Fragmentation, transport and deposition of a low-grade ignimbrite: The Citlaltépetl Ignimbrite, Eastern México

 The Citlaltépetl Ignimbrite records one of the largest explosive events during the Holocene activity of Citlaltépetl Volcano (Pico de Orizaba). Multiple pyroclastic flow units, a fall deposit, and some lahar units were emplaced between 8500–9000 y B.P. as a result of repetitive but discrete explosive events. The whole ignimbrite resulted from discrete fluctuations in eruptive intensity that decreased with time. The initial pyroclastic flow pulse was by far the most violent and widespread event, and its deposits show conspicuous variations in structure and texture that could be associated with different mechanisms of transport and emplacement. Subpopulation Sequential Fragmentation Transport (SFT) analyses were carried out in order to determine the physical mechanisms that selectively concentrate or remove particles in the moving flows. We suggest that lateral and temporal changes in the flow rheology, in which fluidization, yield strength, entrainment of atmospheric air, and sedimentation played a dominant role in flow propagation and emplacement, may imprint a unique signature in the grain-size spectra. The lowermost unit of the Citlaltépetl Ignimbrite can be envisaged by a model in which progressive aggradation near the vent became replaced by en masse emplacement farther outward. Received: 24 March 1998 / Accepted: 9 October 1998     

The impact of ditch blocking on fluvial carbon export from a UK blanket bog

We investigated the effects of ditch blocking on fluvial carbon concentrations and fluxes at a 5‐year, replicated, control‐intervention field experiment on a blanket peatland in North Wales, UK. The site was hydrologically instrumented, and run‐off via open and blocked ditches was analysed for dissolved organic carbon (DOC), particulate organic carbon, dissolved carbon dioxide, and dissolved methane. DOC was also analysed in peat porewater and overland flow. The hillslope experiment was embedded within a paired control‐intervention catchment study, with 3 years of preblocking and 6 years of postblocking data. Results from the hillslope showed large reductions in discharge via blocked ditches, with water partly redirected into hillslope surface and subsurface flows, and partly into remaining open ditches. We observed no impacts of ditch blocking on DOC, particulate organic carbon, dissolved carbon dioxide or methane in ditch waters, DOC in porewaters or overland flow, or stream water DOC at the paired catchment scale. Similar DOC concentrations in ditch water, overland flow, and porewater suggest that diverting flow from the ditch network to surface or subsurface flow had a limited impact on concentrations or fluxes of DOC entering the stream network. The subdued response of fluvial carbon to ditch blocking in our study may be attributable to the relatively low susceptibility of blanket peatlands to drainage, or to physical alterations of the peat since drainage. We conclude that ditch blocking cannot be always be expected to deliver reductions in fluvial carbon loss, or improvements in the quality of drinking water supplies.