MODFLOW-Style parameters in underdetermined parameter estimation

In this article, we discuss the use of MODFLOW-Style parameters in the numerical codes MODFLOW_2005 and MODFLOW_2005-Adjoint for the definition of variables in the Layer Property Flow package. Parameters are a useful tool to represent aquifer properties in both codes and are the only option available in the adjoint version. Moreover, for overdetermined parameter estimation problems, the parameter approach for model input can make data input easier. We found that if each estimable parameter is defined by one parameter, the codes require a large computational effort and substantial gains in efficiency are achieved by removing logical comparison of character strings that represent the names and types of the parameters. An alternative formulation already available in the current implementation of the code can also alleviate the efficiency degradation due to character comparisons in the special case of distributed parameters defined through multiplication matrices. The authors also hope that lessons learned in analyzing the performance of the MODFLOW family codes will be enlightening to developers of other Fortran implementations of numerical codes.     

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.     

SVME: an ensemble of support vector machines for detecting oil spills from full resolution MODIS images

This paper addresses oil spill detection from remotely sensed optical images. In particular, it focuses on the automatic classification of regions of interest (ROIs) in two classes, namely oil spills or look-alikes. Candidate regions and the corresponding boundaries have been manually identified from full resolution Moderate Resolution Imaging Spectroradiometer images, related to the Mediterranean Sea over the years 2008 and 2009. Then, a set of features has been extracted from each ROI, allowing to formulate the oil spill detection problem as a two-class classification task on the provided regions (i.e. using a supervised learning strategy). Since ROI classification is challenging, some desired characteristics for the classification algorithm are first identified, such as accuracy, robustness, etc. Then, a solution (called SVME) is provided: it is based on an ensemble of incremental/decremental cost-oriented Support Vector Machines, aggregated with the Receiving Operating Characteristic (ROC) convex hull method in the ROC space. Such a solution addresses all the desired characteristics. Finally, the results obtained on the collected dataset are shown. The importance of this study is the devising of a powerful classification technique that may have an impact on optical oil spill detection from space, especially if fused with satellite synthetic aperture radar data. Moreover, it is shown how the proposed system can be used as a decision support tool, to help a junior operator in making more reliable detections.     

Genetic variability of the blue and red shrimp Aristeus antennatus in the Western Mediterranean Sea inferred by DNA microsatellite loci

Genetic variation at eight microsatellite loci was studied in nine populations of the blue and red shrimp Aristeus antennatus to investigate whether distinct stocks are present in the Western Mediterranean Sea. A high level of gene flow and no evidence of genetic partitioning were discovered. No significant variation was found (F_ST = 0.00673, P-value = 0.067) even when shrimps from exploited and those from deep-water unexploited grounds were compared. No evidence of reduction or expansion of population size in the recent past was found, as indicated by the bottleneck and interlocus g-tests. Our results are consistent with previous studies using mitochondrial gene methods and allozymes, indicating that, for this species, extensive pelagic larval dispersal and adult migration are probably responsible for the genetic homogeneity observed. In particular, due to a different bathymetric distribution of males and females, reported to be associated with different water masses and hence with possible differential dispersal capacity between sexes, the hypothesis of sex-biased dispersal was tested. Mean values of corrected assignment indices and mean relatedness values were higher for males, suggesting that females are the more widely dispersing sex. Molecular assessment of A. antennatus from the Western Mediterranean provides data of biological and evolutionary interest for the successful management of such a highly valuable fishery resource.     

Tectonic and climatic controls on historical landscape modifications: The avulsion of the lower Cecina River (Tuscany, central Italy)

Integration of geomorphology, stratigraphy, sedimentology and morphotectonics in the analysis of the lower Cecina River reach, coastal Tuscany, reveals an undocumented historical channel avulsion. Geomorphological evidence and radiocarbon dating support that, from the Last Glacial Maximum until the end of the 16th century, the Cecina River flowed north of the present course and formed a well-developed cuspate delta. Two concurrent factors, active tectonics as a preparing factor and discharge regime as an activation factor, are thus inferred to have favored the avulsion of Cecina River. Fragmentary archaeological and historical records indicate that the late Holocene Cecina River plain was virtually unpopulated until the latest 16th century. This seems the main reason why high-magnitude hydrological events and prominent river channel avulsions were not reported in historical chronicles. From this perspective, geomorphological data may provide important knowledge and understanding of recent dynamics of environmental change when historical record is lacking or missing.     

Effects of CO<Subscript>2</Subscript> flushing on crystal textures and compositions: experimental evidence from recent K-trachybasalts erupted at Mt. Etna

Changes in magmatic assemblages and crystal stability as a response of CO₂-flushing in basaltic systems have rarely been directly addressed experimentally, making the role of CO₂ in magma dynamics still controversial and object of scientific debate. We conducted a series of experiments to understand the response of magmas from Etna volcano to CO₂ flushing. We performed a first experiment at 300 MPa to synthesize a starting material composed of crystals of some hundreds of µm and melt pools. This material is representative of an initial magmatic assemblage composed of plagioclase, clinopyroxene and a water-undersaturated melt with 1.6 wt% H₂O. In a second step, the initial assemblage was equilibrated at 300 and 100 MPa with fluids having different XCO ₂ ^fl (CO₂/(H₂O + CO₂)). At low XCO ₂ ^fl (< 0.2 to 0.4), plagioclase is completely dissolved and clinopyroxene show dissolution textures. For relatively high XCO ₂ ^fl (0.9 at 300 MPa), the flushing of a CO₂-rich fluid phase leads to an increase of the amount of clinopyroxene and a decrease of the abundance of plagioclase at 300 MPa. This decrease of plagioclase proportion is associated with a change in An content. Our experiments demonstrate that flushing basaltic systems with fluids may drastically affect crystal textures and phase equilibria depending on proportions of H₂O and CO₂ in the fluid phase. Since texture and crystal proportions are among the most important parameters governing the rheology of magmas, fluid flushing will also influence magma ascent to the Earth’s surface. The experimental results open new perspectives to decipher the textural and compositional record of minerals observed in volcanic rocks from Mt. Etna, and at the same time offer the basis for interpreting the information preserved in minerals from other basaltic volcanoes erupting magmas enriched in CO₂.     

Modeling ground deformations of Panarea volcano hydrothermal/geothermal system (Aeolian Islands, Italy) from GPS data

Panarea volcano (Aeolian Islands, Italy) was considered extinct until November 3, 2002, when a submarine gas eruption began in the area of the islets of Lisca Bianca, Bottaro, Lisca Nera, Dattilo, and Panarelli, about 2.5 km east of Panarea Island. The gas eruption decreased to a state of low degassing by July 2003. Before 2002, the activity of Panarea volcano was characterized by mild degassing of hydrothermal fluid. The compositions of the 2002 gases and their isotopic signatures suggested that the emissions originated from a hydrothermal/geothermal reservoir fed by magmatic fluids. We investigate crustal deformation of Panarea volcano using the global positioning system (GPS) velocity field obtained by the combination of continuous and episodic site observations of the Panarea GPS network in the time span 1995–2007. We present a combined model of Okada sources, which explains the GPS results acquired in the area from December 2002. The kinematics of Panarea volcano show two distinct active crustal domains characterized by different styles of horizontal deformation, supported also by volcanological and structural evidence. Subsidence on order of several millimeters/year is affecting the entire Panarea volcano, and a shortening of 10⁻⁶ year⁻¹ has been estimated in the Islets area. Our model reveals that the degassing intensity and distribution are strongly influenced by geophysical-geochemical changes within the hydrothermal/geothermal system. These variations may be triggered by changes in the regional stress field as suggested by the geophysical and volcanological events which occurred in 2002 in the Southern Tyrrhenian area.