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 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.     

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.     

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.     

On the use of the autocorrelation function: the constraint of using frequency band-limited signals for monitoring relative velocity changes

Correlations of seismic noise are commonly used to monitor temporal variations of relative seismic velocity in period ranges from 1 s up to 100 s. Of particular interest is the detection of small changes in the order of 0.01–0.1 % in propagation speeds. Measuring such small differences can, however, be significantly biased by temporal variations in the properties of the noise sources within the corresponding frequency band. Using synthetic data, we show that apparent relative velocity variations might appear only due to changes in the amplitude and frequency content caused by source variations. Removing such unwanted effects by applying narrow bandpass filters in the preprocessing restricts the high-resolution analysis of any signal due to Gabor’s uncertainty limit, i.e., the correlation function suffers a limited resolution to time delay estimates for small correlation times, low-frequency ranges, and in narrow frequency bands. Better understanding of spatiotemporal noise source properties and the theoretical limitations of time–frequency analysis is critical for accurate and reliable passive monitoring.     

Seismic fragility curves of as-built single-span masonry arch bridges

Masonry arch bridges are crucial elements in the railway transportation network throughout Europe. Although significant advances in seismic risk assessment of various bridge types have been made by developing fragility curves of generalized classes of structures, there are no comparable tools for masonry arch structures. In this context, this paper presents the construction of seismic fragility curves of single-span masonry bridges according to the limit analysis method. An iterative procedure is implemented to define the capacity curve of the equivalent single degree of freedom system through non-linear kinematic analysis. The process involves determination of the collapse mechanism, calculation of the limit load multiplier, and definition of the thrust line. The intrinsic variability of the seismic action is incorporated with the use of different sets of elastic spectra compatible with EC 8 Type-1 spectrum for various types of soil, with peak ground acceleration varying over the range 0.05–1.5 g. The fragility curves of the generalized classes of single-span masonry bridges are finally obtained from the effective ranges of the main geometric and material parameters affecting arch bridge capacity.     

Quartz as indicator mineral in the Central Swiss Alps: the quartz recrystallization isograd in the rock series of the northern Aar massif

Earlier observations on the deformation of rock-forming quartz during Alpine metamorphism that were made by Voll (1976) along a traverse following the Reuss valley (Aar massif, Central Swiss Alps) are still relevant but can be explained in more modern term. Voll found that a few kilometres south of the contact between the Aar massif and its sedimentary cover quartz has partly undergone a combination of dynamic recrystallization and nucleation along grain boundaries and fractures at the expense of the deformed parent grains. Further to the south, increasing grain growth of quartz under nearly static conditions is observed, with polygonal grain shapes typical of annealed microfabrics. Comparable stages of recrystallization could be found in several N–S traverses through the Aar massif so that “points of first occurrence” of newly formed quartz could be connected to define a quartz recrystallization isograd. This isograd extents over 90 km through the northern Aar massif, following a course parallel to the sanidine/microcline isograd but at a distance of 10–15 km further to the north. The evolution of the quartz microfabric is discussed with respect to recent estimates of the T-t path of the Aar massif.