Multiple-scale hydrological and geophysical data integration through non-linear Bayesian sequential simulation

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the scale of a field site represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed downscaling procedure based on a non-linear Bayesian sequential simulation approach. The main objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity logged at collocated wells and surface resistivity measurements, which are available throughout the studied site. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariatekernel density function. Then a stochastic integration of low-resolution, large-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities is applied. The overall viability of this downscaling approach is tested and validated by comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure allows obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Structural characterization of Turtle Mountain anticline (Alberta, Canada) and impact on rock slope failure

his paper proposes a structural investigation of the Turtle Mountain anticline (Alberta, Canada) to better understand the role of the different tectonic features on the development of both local and large scale rock slope instabilities occurring in Turtle Mountain. The study area is investigated by combining remote methods with detailed field surveys. In particular, the benefit of Terrestrial Laser Scanning for ductile and brittle tectonic structure interpretations is illustrated. The proposed tectonic interpretation allows the characterization of the fracturing pattern, the fold geometry and the role of these tectonic features in rock slope instability development. Ten discontinuity sets are identified in the study area, their local variations permitting the differentiation of the study zone into 20 homogenous structural domains. The anticline is described as an eastern verging fold that displays considerable geometry differences along its axis and developed by both flexural slip and tangential longitudinal strain folding mechanisms. Moreover, the origins of the discontinuity sets are determined according to the tectonic phases affecting the region (pre-folding, folding, post-folding). The localization and interpretation of kinematics of the different instabilities revealed the importance of considering the discrete brittle planes of weakness, which largely control the kinematic release of the local instabilities, and also the rock mass damage induced by large tectonic structures (fold hinge, thrust).

CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes.

BACKGROUND: The Complete Arabidopsis Transcript MicroArray (CATMA) initiative combines the efforts of laboratories in eight European countries 1 to deliver gene-specific sequence tags (GSTs) for the Arabidopsis research community. The CATMA initiative offers the power and flexibility to regularly update the GST collection according to evolving knowledge about the gene repertoire. These GST amplicons can easily be reamplified and shared, subsets can be picked at will to print dedicated arrays, and the GSTs can be cloned and used for other functional studies. This ongoing initiative has already produced approximately 24,000 GSTs that have been made publicly available for spotted microarray printing and RNA interference. RESULTS: GSTs from the CATMA version 2 repertoire (CATMAv2, created in 2002) were mapped onto the gene models from two independent Arabidopsis nuclear genome annotation efforts, TIGR5 and PSB-EuGène, to consolidate a list of genes that were targeted by previously designed CATMA tags. A total of 9,027 gene models were not tagged by any amplified CATMAv2 GST, and 2,533 amplified GSTs were no longer predicted to tag an updated gene model. To validate the efficacy of GST mapping criteria and design rules, the predicted and experimentally observed hybridization characteristics associated to GST features were correlated in transcript profiling datasets obtained with the CATMAv2 microarray, confirming the reliability of this platform. To complete the CATMA repertoire, all 9,027 gene models for which no GST had yet been designed were processed with an adjusted version of the Specific Primer and Amplicon Design Software (SPADS). A total of 5,756 novel GSTs were designed and amplified by PCR from genomic DNA. Together with the pre-existing GST collection, this new addition constitutes the CATMAv3 repertoire. It comprises 30,343 unique amplified sequences that tag 24,202 and 23,009 protein-encoding nuclear gene models in the TAIR6 and EuGène genome annotations, respectively. To cover the remaining untagged genes, we identified 543 additional GSTs using less stringent design criteria and designed 990 sequence tags matching multiple members of gene families (Gene Family Tags or GFTs) to cover any remaining untagged genes. These latter 1,533 features constitute the CATMAv4 addition. CONCLUSION: To update the CATMA GST repertoire, we designed 7,289 additional sequence tags, bringing the total number of tagged TAIR6-annotated Arabidopsis nuclear protein-coding genes to 26,173. This resource is used both for the production of spotted microarrays and the large-scale cloning of hairpin RNA silencing vectors. All information about the resulting updated CATMA repertoire is available through the CATMA database http://www.catma.org.

Networks in Argentine agriculture: a multiple-case study approach

Argentina is among the four largest producers of soybeans, sunflower, corn, and wheat, among other agricultural products. Institutional and policy changes during the 1990s fostered the development of Argentine agriculture and the introduction of innovative process and product technologies (no-till, agrochemicals, GMO, GPS) and new investments in modern, large-scale sunflower and soybean processing plants. In addition to technological changes, a "quiet revolution" occurred in the way agricultural production was carried out and organized: from self-production or ownership agriculture to a contract-based agriculture. The objective of this paper is to explore and describe the emergence of networks in the Argentine crop production sector. The paper presents and describes four cases that currently represent about 50% of total grain and oilseed production in Argentina: "informal hybrid form", "agricultural trust fund", "investor-oriented corporate structure", and "network of networks". In all cases, hybrid forms involve a group of actors linked by common objectives, mainly to gain scale, share resources, and improve the profitability of the business. Informal contracts seem to be the most common way of organizing the agriculture process, but using short-term contracts and sequential interfirm collaboration. Networks of networks involve long-term relationships and social development, and reciprocal interfirm collaboration. Agricultural trust fund and investor-oriented corporate structures have combined interfirm collaboration and medium-term relationships. These organizational forms are highly flexible and show a great capacity to adapt to challenges; they are competitive because they enjoy aligned incentives, flexibility, and adaptability.

Solving Large Location-Allocation problems by Clustering and Simulated Annealing

Globalization involves several facility location problems that need to be handled at large scale. Location Allocation (LA) is a combinatorial problem in which the distance among points in the data space matter. Precisely, taking advantage of the distance property of the domain we exploit the capability of clustering techniques to partition the data space in order to convert an initial large LA problem into several simpler LA problems. Particularly, our motivation problem involves a huge geographical area that can be partitioned under overall conditions. We present different types of clustering techniques and then we perform a cluster analysis over our dataset in order to partition it. After that, we solve the LA problem applying simulated annealing algorithm to the clustered and non-clustered data in order to work out how profitable is the clustering and which of the presented methods is the most suitable

Structural and metamorphic evolution of the northern Himachal Himalaya, NW India - (Spiti-eastern Lahul-Parvati valley traverse)

The Himalayan orogen is the result of the collision between the Indian and Asian continents that began 55-50 Ma ago, causing intracontinental thrusting and nappe formation. Detailed mapping as well as structural and microfabric analyses on a traverse from the Tethyan Himalaya southwestward through the High Himalayan Crystalline and the Main Central Thrust zone (MCT zone) to the Lesser Himalayan Sequence in the Spiti-eastern Lahul-Parvati valley area reveal eight main phases of deformation, a series of late stage phases and five stages of metamorphic crystallization. This sequence of events is integrated into a reconstruction of the tectonometamorphic evolution of the Himalayan orogen in northern Himachal Pradesh. The oldest phase D-1 is preserved as relies in the High Himalayan Crystalline. Its deformational conditions are poorly known, but the metamorphic evolution is well documented by a prograde metamorphism reaching peak conditions within the upper amphibolite facies. This indicates that D-1 was an important tectonometamorphic event including considerable crustal thickening. The structural, metamorphic and sedimentary record suggest that D-1 most probably represents an early stage of continental collision. The first event clearly attributed to the collision between India and Asia is documented by two converging nappe systems, the NE-verging Shikar Beh Nappe and the SW-verging north Himalayan nappes. The D-2 Shikar Beh Nappe is characterized by isoclinal folding and top-to-the NE shearing, representing the main deformation in the High Himalayan Crystalline. D-2 also caused the main metamorphism in the High Himalayan Crystalline that was of a Barrovian-type, reaching upper amphibolite facies peak conditions. The Shikar Beh Nappe is interpreted to have formed within the Indian crust SW of the subduction zone. Simultaneously with NE-directed nappe formation, incipient subduction of India below Asia caused stacking of the SW-verging north Himalayan Nappes, that were thrust from the northern edge of the subducted continent toward the front of the Shikar Beh Nappe. As a result, the SW-verging folds of the D-3 Main Fold Zone formed in the Tethyan Himalaya below the front of the north Himalayan nappes. D-3 represents the main deformation in the Tethyan Himalaya, associated with a greenschist facies metamorphism. Folding within the Main Fold Zone subsequently propagated toward SW into the High Himalayan Crystalline, where it overprinted the preexisting D-2 structures. After subduction at the base of the north Himalayan nappes, the subduction zone stepped to the base of the High Himalayan Crystalline, where D-3 folds were crosscut by SW-directed D-4 thrusting. During D-4, the Crystalline Nappe, comprising the Main Fold Zone and relies of the Shikar Beh Nappe was thrust toward SW over the Lesser Himalayan Sequence along the 4 to 5 kms thick Main Central Thrust zone. Thrusting was related to a retrograde greenschist facies overprint at the base of the Crystalline Nappe and to pro-grade greenschist facies conditions in the Lesser Himalayan Sequence. Simultaneously with thrusting at the base of the Crystalline Nappe, higher crustal levels were affected by NE-directed D-5 normal extensional shearing and by dextral strike-slip motion, indicating that the high-grade metamorphic Crystalline Nappe was extruded between the low-grade metamorphic Lesser Himalayan Sequence at the base and the north Himalayan nappes at the top. The upper boundary of the Crystalline Nappe is not clearly delimited and passes gradually into the low-grade rocks at the front of the north Himalayan nappes. Extrusion of the Crystalline Nappe was followed by the phase D-6, characterized by large-scale, upright to steeply inclined, NE-verging folds and by another series of normal and extensional structures D-7+D-8 that may be related to ongoing extrusion of the Crystalline Nappe. The late stage evolution is represented by the phases D-A and D-B that indicate shortening parallel to the axis of the mountain chain and by D-C that is interpreted to account for the formation of large-scale domes with NNW-SSE-trending axes, an example of which is exposed in the Larji-Kullu-Rampur tectonic window.

Debris flow susceptibility mapping at a regional scale along the National Road N7, Argentina

A first assessment of debris flow susceptibility at a large scale was performed along the National Road N7, Argentina. Numerous catchments are prone to debris flows and likely to endanger the road-users. A 1:50,000 susceptibility map was created. The use of a DEM (grid 30 m) associated to three complementary criteria (slope, contributing area, curvature) allowed the identification of potential source areas. The debris flow spreading was estimated using a process- and GISbased model (Flow-R) based on basic probabilistic and energy calculations. The best-fit values for the coefficient of friction and the mass-to-drag ratio of the PCM model were found to be ? = 0.02 and M/D = 180 and the resulting propagation on one of the calibration site was validated using the Coulomb friction model. The results are realistic and will be useful to determine which areas need to be prioritized for detailed studies.

Biodiversity and global change in Mediterranean benthic communities: perspectives over large spatial and long temporal scales

Many terrestrial and marine systems are experiencing accelerating decline due to the effects of global change. This situation has raised concern about the consequences of biodiversity losses for ecosystem function, ecosystem service provision, and human well-being. Coastal marine habitats are a main focus of attention because they harbour a high biological diversity, are among the most productive systems of the world and present high anthropogenic interaction levels. The accelerating degradation of many terrestrial and marine systems highlights the urgent need to evaluate the consequence of biodiversity loss. Because marine biodiversity is a dynamic entity and this study was interested global change impacts, this study focused on benthic biodiversity trends over large spatial and long temporal scales. The main aim of this project was to investigate the current extent of biodiversity of the high diverse benthic coralligenous community in the Mediterranean Sea, detect its changes, and predict its future changes over broad spatial and long temporal scales. These marine communities are characterized by structural species with low growth rates and long life spans; therefore they are considered particularly sensitive to disturbances. For this purpose, this project analyzed permanent photographic plots over time at four locations in the NW Mediterranean Sea. The spatial scale of this study provided information on the level of species similarity between these locations, thus offering a solid background on the amount of large scale variability in coralligenous communities; whereas the temporal scale was fundamental to determine the natural variability in order to discriminate between changes observed due to natural factors and those related to the impact of disturbances (e.g. mass mortality events related to positive thermal temperatures, extreme catastrophic events). This study directly addressed the challenging task of analyzing quantitative biodiversity data of these high diverse marine benthic communities. Overall, the scientific knowledge gained with this research project will improve our understanding in the function of marine ecosystems and their trajectories related to global change.

SGP-1: prediction and validation of homologous genes based on sequence alignments

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of SGP-1 depends little on species-specific properties such as codon usage or the nucleotide distribution. SGP-1 may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

Comparative phylogeography of mutualists and the effect of the host on the genetic structure of its partners.

Whether or not species participating in specialized and obligate interactions display similar and simultaneous demographic variations at the intraspecific level remains an open question in phylogeography. In the present study, we used the mutualistic nursery pollination occurring between the European globeflower Trollius europaeus and its specialized pollinators in the genus Chiastocheta as a case study. Explicitly, we investigated if the phylogeographies of the pollinating flies are significantly different from the expectation under a scenario of plant-insect congruence. Based on a large-scale sampling, we first used mitochondrial data to infer the phylogeographical histories of each fly species. Then, we defined phylogeographical scenarios of congruence with the plant history, and used maximum likelihood and Bayesian approaches to test for plant-insect phylogeographical congruence for the three Chiastocheta species. We show that the phylogeographical histories of the three fly species differ. Only Chiastocheta lophota and Chiastocheta dentifera display strong spatial genetic structures, which do not appear to be statistically different from those expected under scenarios of phylogeographical congruence with the plant. The results of the present study indicate that the fly species responded in independent and different ways to shared evolutionary forces, displaying varying levels of congruence with the plant genetic structure

Local-scale spatial variation in diversity of social wasps in an Amazonian rain forest in Caxiuanã, Pará, Brazil (Hymenoptera, Vespidae, Polistinae)

Polistine wasps are important in Neotropical ecosystems due to their ubiquity and diversity. Inventories have not adequately considered spatial attributes of collected specimens. Spatial data on biodiversity are important for study and mitigation of anthropogenic impacts over natural ecosystems and for protecting species. We described and analyzed local-scale spatial patterns of collecting records of wasp species, as well as spatial variation of diversity descriptors in a 2500-hectare area of an Amazon forest in Brazil. Rare species comprised the largest fraction of the fauna. Close range spatial effects were detected for most of the more common species, with clustering of presence-data at short distances. Larger spatial lag effects could also be identified in some species, constituting probably cases of exogenous autocorrelation and candidates for explanations based on environmental factors. In a few cases, significant or near significant correlations were found between five species (of Agelaia, Angiopolybia, and Mischocyttarus) and three studied environmental variables: distance to nearest stream, terrain altitude, and the type of forest canopy. However, association between these factors and biodiversity variables were generally low. When used as predictors of polistine richness in a linear multiple regression, only the coefficient for the forest canopy variable resulted significant. Some level of prediction of wasp diversity variables can be attained based on environmental variables, especially vegetation structure. Large-scale landscape and regional studies should be scheduled to address this issue.

Pre-Mesozoic Alpine basements - Their place in the European Paleozoic framework

Prior to their Alpine overprinting, most of the pre-Mesozoic basement areas in Alpine orogenic structures shared a complex evolution, starting with Neoproterozoic sediments that are thought to have received detrital input from both West and East Gondwanan cratonic sources. A subsequent Neoproterozoic-Cambrian active margin setting at the Gondwana margin was followed by a Cambrian-Ordovician rifting period, including an Ordovician cordillera-like active margin setting. During the Late Ordovician and Silurian periods, the future Alpine domains recorded crustal extension along the Gondwana margin, announcing the future opening of the Paleotethys oceanic domain. Most areas then underwent Variscan orogenic events, including continental subduction and collisions with Avalonian-type basement areas along Laurussia and the juxtaposition and the duplication of terrane assemblages during strike slip, accompanied by contemporaneous crustal shortening and the subduction of Paleotethys under Laurussia. Thereafter, the final Pangea assemblage underwent Triassic and Jurassic extension, followed by Tertiary shortening, and leading to the buildup of the Alpine mountain chain. Recent plate-tectonic reconstructions place the Alpine domains in their supposed initial Cambrian-Ordovician positions in the eastern part of the Gondwana margin, where a stronger interference with the Chinese blocks is proposed, at least from the Ordovician onward. For the Visean time of the Variscan continental collision, the distinction of the former tectonic lower-plate situation is traceable but becomes blurred through the subsequent oblique subduction of Paleotethys under Laurussia accompanied by large-scale strike slip. Since the Pennsylvanian, this global collisional scenario has been replaced by subsequent and ongoing shortening and strike slip under rising geothermal conditions, and all of this occurred before all these puzzle elements underwent the complex Alpine reorganization.

A Small Baseline DIFSAR Approach for Investigating Deformations on Full Resolution SAR Interferograms

This paper presents a differential synthetic apertureradar (SAR) interferometry (DIFSAR) approach for investigatingdeformation phenomena on full-resolution DIFSAR interferograms.In particular, our algorithm extends the capabilityof the small-baseline subset (SBAS) technique that relies onsmall-baseline DIFSAR interferograms only and is mainly focusedon investigating large-scale deformations with spatial resolutionsof about 100 100 m. The proposed technique is implemented byusing two different sets of data generated at low (multilook data)and full (single-look data) spatial resolution, respectively. Theformer is used to identify and estimate, via the conventional SBAStechnique, large spatial scale deformation patterns, topographicerrors in the available digital elevation model, and possibleatmospheric phase artifacts; the latter allows us to detect, onthe full-resolution residual phase components, structures highlycoherent over time (buildings, rocks, lava, structures, etc.), as wellas their height and displacements. In particular, the estimation ofthe temporal evolution of these local deformations is easily implementedby applying the singular value decomposition technique.The proposed algorithm has been tested with data acquired by theEuropean Remote Sensing satellites relative to the Campania area(Italy) and validated by using geodetic measurements.

Development of Stage-Discharge Relations for Ungaged Bridge Waterways in Western Iowa, TR-567, 2010

Stage-discharge relations constitute a viable, alternative technique for estimating accurately flow for ungaged sites. In this research, we have utilized pressure transducers and Large Scale Particle Image Velocimetry techniques to develop stage-discharge relations at eleven sites in the Hungry Canyon Area (HCA) of southwestern Iowa under different hydrologic conditions. We have employed these data to calibrate and verify an established hydrologic model and then we have used this model to provide a stage discharge relation for different hydrologic conditions (i.e. rating curves). The benefits of the project are numerous including that the discharge data will be used for a number of purposes, including operational decision making in the HCA about the design of water-control and conveyance structures, input for hydraulic and hydrologic models, and calculation of sediment and other water quality constituents transport and “loads”, and for decision making. This project has also pointed out the difficulties in measuring flows in ungaged streams with ice jams, steep banks, erodible beds, and floating debris.

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.