Investigation of Small‐Scale Preferential Flow with a Forced‐Gradient Tracer Test

A new tracer experiment (referred to as MADE‐5) was conducted at the well‐known Macrodispersion Experiment (MADE) site to investigate the influence of small‐scale mass‐transfer and dispersion processes on well‐to‐well transport. The test was performed under dipole forced‐gradient flow conditions and concentrations were monitored in an extraction well and in two multilevel sampler (MLS) wells located at 6, 1.5, and 3.75 m from the source, respectively. The shape of the breakthrough curve (BTC) measured at the extraction well is strongly asymmetric showing a rapidly arriving peak and an extensive late‐time tail. The BTCs measured at seven different depths in the two MLSs are radically different from one another in terms of shape, arrival times, and magnitude of the concentration peaks. All of these characteristics indicate the presence of a complex network of preferential flow pathways controlling solute transport at the test site. Field‐experimental data were also used to evaluate two transport models: a stochastic advection‐dispersion model (ADM) based on conditional multivariate Gaussian realizations of the hydraulic conductivity field and a dual‐domain single‐rate (DDSR) mass‐transfer model based on a deterministic reconstruction of the aquifer heterogeneity. Unlike the stochastic ADM realizations, the DDSR accurately predicted the magnitude of the concentration peak and its arrival time (within a 1.5% error). For the multilevel BTCs between the injection and extraction wells, neither model reproduced the observed values, indicating that a high‐resolution characterization of the aquifer heterogeneity at the subdecimeter scale would be needed to fully capture 3D transport details.     

Lessons Learned from 25 Years of Research at the MADE Site

Field studies at well‐instrumented research sites have provided extensive data sets and important insights essential for development and testing of transport theories and mathematical models. This paper provides an overview of over 25 years of research and lessons learned at one of such field research sites on the Columbus Air Force Base in Mississippi, commonly known as the Macrodispersion Experiment (MADE) site. Since the mid‐1980s, field data from the MADE site have been used extensively by researchers around the world to explore complex contaminant transport phenomena in highly heterogeneous porous media. Results from field investigations and modeling analyses suggested that connected networks of small‐scale preferential flow paths and relative flow barriers exert dominant control on solute transport processes. The classical advection‐dispersion model was shown to inadequately represent plume‐scale transport, while the dual‐domain mass transfer model was found to reproduce the primary observed plume characteristics. The MADE site has served as a valuable natural observatory for contaminant transport studies where new observations have led to better understanding and improved models have sprung out analysis of new data.     

The response of braided planform configuration to flow variations,bed reworking and vegetation: the case of the Tagliamento River,Italy

Morphological features of braided rivers (bars, channels and pools) experience major changes in area, shape and spatial distribution as a response to (i) the pulsation of discharge during a flood and (ii) the bed evolution induced by floods. In this work, at‐a‐station relationships between water level and planform configuration were investigated on the Tagliamento River, a large gravel‐bed braided river in northeast Italy, over a 2‐year study period comprising three bankfull events and several small‐to‐medium floods. The analysis was performed on two 1‐km‐long reaches, characterized by different riparian vegetation cover. Ground‐based images with an hourly temporal resolution were acquired using software‐controlled, digital cameras. Bars, channels, pools and vegetated patches were manually digitized on more than 100 rectified images. Sequences of constant‐level images spanning the study period were used to quantify the impact of floods on the stability of at‐a‐station relationships and on the turnover rate of water bodies. The analysis shows that wetted area increased almost linearly with water level in both reaches. The average number of branches per cross‐section peaked at intermediate flow levels, increasing from 2 at low flow up to 6–7. The number of branches displayed the largest fluctuations over time, with significant changes produced also by moderate floods. Turnover rates were high in both reaches, with more than 30% of wetted areas at low flow converting into bare gravel in less than 2 months. Vegetation colonization was found to limit the mobility of the low flow channels over time by concentrating the flow in fewer, deeper anabranches. The number of channels per cross‐section was 30–40% less in the vegetated reach and the proportion of low flow water bodies in the same position after 12 months increased from 3% to 14%. Copyright © 2012 John Wiley & Sons, Ltd.     

The influence of local topography for wind direction on Mars: two examples of dune fields in crater basins

In this work, we perform an analysis of large dark dunes within Moreux Crater and Herschel Crater on Mars using High Resolution Imaging Science Experiment (HiRISE) and Context Camera (CTX) data sets. These data allow us to conduct a detailed analysis of dune morphology and slip faces, concluding that the studied dune fields are influenced by topographically‐controlled complex wind directions. Our morphological analysis reveals that inside Moreux Crater in particular, the topographic setting dominates the wind flow direction, leading to the development of a sand transport pathway encircling the central peak of the crater. The dune fields in Herschel Crater are also affected by winds controlled by variable topography as suggested by the presence of complex dunes and dune fields. Our analysis indicate that the studied dune systems is not the result of paleo‐wind regimes. Furthermore, we perform thermal inertia measurements using thermal emission spectrometer (TES) data, which indicate that the studied dune fields consist of medium sand 250–500 µm in diameter. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved 2-D attenuation analysis for Northern Italy using a merged dataset from selected regional seismic networks

A merged, high-quality waveform dataset from different seismic networks has been used to improve our understanding of lateral seismic attenuation for Northern Italy. In a previous study on the same region, Morasca et al. (Bull Seismol Soc Am 98:1936–1946, 2008) were able to resolve only a small area due to limited data coverage. For this reason, the interpretation of the attenuation anomalies was difficult given the complexity of the region and the poor resolution of the available data. In order to better understand the lateral changes in the crustal structure and thickness of this region, we selected 770 earthquakes recorded by 54 stations for a total of almost 16,000 waveforms derived from seismic networks operating totally or partially in Northern Italy. Direct S-wave and coda attenuation images were obtained using an amplitude ratio technique that eliminates source terms from the formulation. Both direct and early-coda amplitudes are used as input for the inversions, and the results are compared. Results were obtained for various frequency bands ranging between 0.3 and 25.0 Hz and in all cases show significant improvement with respect to the previous study since the resolved area has been extended and more crossing paths have been used to image smaller scale anomalies. Quality-factor estimates are consistent with the regional tectonic structure exhibiting a general trend of low attenuation under the Po Plain basin and higher values for the Western Alps and Northern Apennines. The interpretation of the results for the Eastern Alps is not simple, possibly because our resolution for this area is still not adequate to resolve small-scale structures.     

Compositionally zoned crystals and real-time degassing data reveal changes in magma transfer dynamics during the 2006 summit eruptive episodes of Mt. Etna

One of the major objectives of volcanology remains relating variations in surface monitoring signals to the magmatic processes at depth that cause these variations. We present a method that enables compositional and temporal information stored in zoning of minerals (olivine in this case) to be linked to observations of real-time degassing data. The integrated record may reveal details of the dynamics of gradual evolution of a plumbing system during eruption. We illustrate our approach using the 2006 summit eruptive episodes of Mt. Etna. We find that the history tracked by olivine crystals, and hence, most likely the magma pathways within the shallow plumbing system of Mt. Etna, differed considerably between the July and October eruptions. The compositional and temporal record preserved in the olivine zoning patterns reveal two mafic recharge events within months of each other (June and September 2006), and each of these magma supplies may have triggered the initiation of different eruptive cycles (July 14–24 and August 31–December 14). Correlation of these observations with gas monitoring data shows that the systematic rise of the CO₂/SO₂ gas values is associated with the gradual (pre- and syn-eruptive) supply of batches of gas-rich mafic magma into segments of Etna’s shallow plumbing system, where mixing with pre-existing and more evolved magma occurred.     

Motif Discovery on Seismic Amplitude Time Series: The Case Study of Mt Etna 2011 Eruptive Activity

Algorithms searching for similar patterns are widely used in seismology both when the waveforms of the events of interest are known and when there is no a priori-knowledge. Such methods usually make use of the cross-correlation coefficient as a measure of similarity; if there is no a-priori knowledge, they behave as brute-force searching algorithms. The disadvantage of these methods, preventing or limiting their application to very large datasets, is computational complexity. The Mueen–Keogh (MK) algorithm overcomes this limitation by means of two optimization techniques—the early abandoning concept and space indexing. Here, we apply the MK algorithm to amplitude time series retrieved from seismic signals recorded during episodic eruptive activity of Mt Etna in 2011. By adequately tuning the input to the MK algorithm we found eight motif groups characterized by distinct seismic amplitude trends, each related to a different phenomenon. In particular, we observed that earthquakes are accompanied by sharp increases and decreases in seismic amplitude whereas lava fountains are accompanied by slower changes. These results demonstrate that the MK algorithm, because of its particular features, may have wide applicability in seismology.     

Microlevelling procedures applied to regional aeromagnetic data: an example from the Transantarctic Mountains (Antarctica)

The extensive application of digital enhancement and filtering as a powerful tool for aeromagnetic interpretation, not only of high resolution but also of regional data, requires an improved levelling. Two microlevelling techniques were thus compared in order to find an effective but relatively simple procedure to remove, or at least to reduce, residual magnetic errors remaining after standard levelling processes. This study was carried out on regional aeromagnetic data recently acquired at high magnetic latitudes along the Transantarctic Mountains in Antarctica, where it is particularly critical to remove time-dependent magnetic variations. Two-dimensional FFT filters applied to the gridded data, namely the Butterworth and a directional cosine filter, proved to be more effective than previously proposed one-dimensional space-domain filters in the reduction of the 'residual corrugation' not removed by statistical levelling. Tectonic interpretation of trends detected in the total field magnetic anomaly map and in the 3D analytic signal improved after application of frequency-domain microlevelling. However, we also show that when interpreting microlevelled data, two factors must be considered: (i) the possible presence of real geological trends aligned along the flight lines; (ii) modifications in the results yielded by depth estimates of magnetic sources due to the FFT filters applied during the microlevelling procedure. Such changes were seen both in the well-established 2D FFT method, based on the slope of the energy spectrum, and in the more recent 3D Euler deconvolution technique. Overall our results indicate that microlevelling could profitably be applied to older gridded aeromagnetic data sets in Antarctica, thus improving the accuracy and geological significance of future regional magnetic compilations, as already seen in other continents.     

Volcanological evolution of the Rivi–Capo Volcanic Complex at Salina,Aeolian Islands: magma storage processes and ascent dynamics

Critical to understanding explosive eruptions is establishing how accurately representative pyroclasts are of processes during magma vesiculation and fragmentation. Here, we present data on densities, and vesicle size and number characteristics, for representative pyroclasts from six silicic eruptions of contrasting size and style from Raoul volcano (Kermadec arc). We use these data to evaluate histories of bubble nucleation, coalescence, and growth in explosive eruptions and to provide comparisons with pumiceous dome carapace material. Density/vesicularity distributions show a scarcity of pyroclasts with ～65–75 % vesicularity; however, pyroclasts closest to this vesicularity range have the highest bubble number density (BND) values regardless of eruptive intensity or style. Clasts with vesicularities greater than this 65–75 % “pivotal” vesicularity range have decreasing BNDs with increasing vesicularities, interpreted to reflect continuing bubble growth and coalescence. Clasts with vesicularities less than the pivotal range have BNDs that decrease with decreasing vesicularity and preserve textures indicative of processes such as stalling and open system degassing prior to vesiculation in a microlite-rich magma, or vesiculation during slow ascent of degassing magma. Bubble size distributions (BSDs) and BNDs show variations consistent with 65–75 % representing the vesicularity at which vesiculating magma is most likely to undergo fragmentation, consistent with the closest packing of spheres. We consider that the observed vesicularity range may reflect the development of permeability in the magma through shearing as it flows through the conduit. These processes can act in concert with multiple nucleation events, generating a situation of heterogeneous bubble populations that permit some regions of the magma to expand and bubbles to coalesce with other regions in which permeable networks are formed. Fragmentation preserves the range in vesicularity seen as well as any post-fragmentation/pre-quenching expansion which may have occurred. We demonstrate that differing density pyroclasts from a single eruption interval can have widely varying BND values corresponding to the degree of bubble maturation that has occurred. The modal density clasts (the usual targets for vesicularity studies) have likely undergone some degree of bubble maturation and are therefore may not be representative of the magma at the onset of fragmentation.     

Benthic respiration and stoichiometry of regenerated nutrients in lake sediments with Dreissena polymorpha

The effect of patchy colonies of the invasive zebra mussel (Dreissena polymorpha) on sedimentary processes was investigated in a mesotrophic lake (Plateliai Lake, Lithuania). Benthic fluxes of O₂, TCO₂, CH₄, Mn²⁺, Fe²⁺, N₂, the inorganic forms of N, Si and P and dissolved organic C and N were quantified by dark incubations of sediments cores, with and without D. polymorpha. Individual mussels also were incubated for metabolic measurements. Sediments with D. polymorpha had significantly higher O₂ and TCO₂ fluxes and displayed higher rates of denitrification. The presence of mussels also resulted in higher regeneration of P and N (mostly as ammonium) while the effect on Si was not significant. However, likely due to the low zebra mussel biomass (57.2 ± 25.3 g_SFDW m^?2), biodeposition has not changed the ratio between anaerobic and total respiration. Methane and reduced metals fluxes were in fact similar in the presence and absence of D. polymorpha. Incubations of mussels without sediments confirmed that bivalve metabolism was the main driver of benthic respiration and nutrient recycling. Nitrate production suggested the presence of nitrifiers associated with the molluscs. The main outcome of this study was that zebra mussels alter the quantity and the stoichiometry of nutrients regenerated by the benthic compartment. The enhancement of nitrogen loss via denitrification, by a factor of 1.5, was much less pronounced than the increase in ammonium recycling rate, stimulated by a factor of 33. Negligible PO₄ ^3? fluxes in bare sediments (?3.4 ± 6.8 μmol m^?2 h^?1) increased in the presence of mussels and considerable amounts of this nutrient (69.6 ± 29.4 μmol m^?2 h^?1) were mobilized to the water column. Further research should address other nutrient sources to the lake to verify whether altered rates and stoichiometry of benthic regeneration can affect primary producer community composition and activity.