Parallel distributed computing using Python     

Lisandro D. Dalcin  Rodrigo R. Paz Pablo A. Kler  Alejandro Cosimo 《Advances in water resources》2011,34(9):1124-1139

This work presents two software components aimed to relieve the costs of accessing high-performance parallel computing resources within a Python programming environment: MPI for Python and PETSc for Python.MPI for Python is a general-purpose Python package that provides bindings for the Message Passing Interface (MPI) standard using any back-end MPI implementation. Its facilities allow parallel Python programs to easily exploit multiple processors using the message passing paradigm. PETSc for Python provides access to the Portable, Extensible Toolkit for Scientific Computation (PETSc) libraries. Its facilities allow sequential and parallel Python applications to exploit state of the art algorithms and data structures readily available in PETSc for the solution of large-scale problems in science and engineering.MPI for Python and PETSc for Python are fully integrated to PETSc-FEM, an MPI and PETSc based parallel, multiphysics, finite elements code developed at CIMEC laboratory. This software infrastructure supports research activities related to simulation of fluid flows with applications ranging from the design of microfluidic devices for biochemical analysis to modeling of large-scale stream/aquifer interactions.  相似文献   

From back-analysis to run-out prediction: a case study in the Western Italian Alps   总被引：3，自引：1，他引：2  

Marina Pirulli  Alessio Colombo  Claudio Scavia 《Landslides》2011,8(2):159-170

Rock avalanches are complex phenomena that occur with a low frequency but which have a high destructive potential. As a consequence, the people who are responsible for the management of a territory are more and more interested in predicting the possible evolutions of well-known potential events. Tackling the above problems from a quantitative point of view, the RASH3D code, based on continuum mechanics concepts, has been here used to predict the evolution of a potential rock avalanche in the Western Italian Alps. A calibration-based approach, in which rheological parameters are constrained by systematic adjustment during trial-and-error back-analysis of past events similar to the landslide under investigation, is proposed to set rheological parameter values to be used for prediction purposes. The back-analysis of a 2?10⁶ m³ rock avalanche located in the Divedro Valley, close to the area of the potential event, has then been analysed using both a frictional and a Voellmy rheology. The characteristics of the slope and the dynamics of the event have made the frictional rheology more suitable to come to the correct simulation of the historical case. The back-analysis results have contributed not only in the selection of the rheological parameter values but also in the choice of the type of rheological law to use in the carried out forward-analyses.  相似文献   

Recent evolution of Lake Arreo,northern Spain: influences of land use change and climate     

Juan Pablo Corella  Adel El Amrani  Javier Sigró  Mario Morellón  Eugenio Rico  Blas Lorenzo Valero-Garcés 《Journal of Paleolimnology》2011,46(3):469-485

We present a high-resolution, multiproxy reconstruction of the depositional history of Lake Arreo, northern Spain, for the last 60 years. We conducted sedimentological, geochemical and diatom analyses in short cores and made a detailed comparison with regional instrumental climate data (1952–2007), limnological monitoring of the lake (1992–2008) and recent land use changes that affect the lake catchment. Chronology is based on “floating” discontinuous varve counts and ¹³⁷Cs and ¹⁴C dates. Four periods were identified in the Lake Arreo recent history: (1) prior to 1963, varved facies intercalated with fine turbidite deposits, and diatom assemblages dominated by Cyclotella taxa indicate predominantly meromictic conditions, (2) from 1964 to 1978, permanent anoxia persisted in bottom waters, as shown by similar facies and diatom assemblages as before, though detrital layers were coarser, (3) from 1979 to 1994, sediment delivery to the lake increased and laminated, clastic facies were deposited, and (4) from 1995 to 2008, dominance of massive facies and an increase in Fragilaria tenera and Achnanthes minutissima reflect relatively lower lake levels, less frequent bottom anoxia with more frequent water column mixing, similar to modern conditions. The period 1952–1979 was a time of meromixis and varved facies deposition, and was characterized by higher rainfall and less intense agricultural pressure in the watershed. There were two short humid periods (1992–1993 and 1996–1998) when monitoring data show more anoxic weeks per year and relatively higher lake levels. Increased cultivation of small landholdings in 1963, and particularly after 1979, caused a large increase in sediment delivery to the lake. The inferred lake evolution is in agreement with monitoring data that suggest a transition from dominantly meromictic conditions prior to 1993–1994 to a predominantly monomictic pattern of circulation since then, particularly after 2000. The synergistic effects of intensive water extraction for irrigation and lower rainfall since 1979, and particularly since 1994, brought the long period of meromictic conditions in Lake Arreo to an end. Water balance and sediment delivery to the lake are dominant factors that control the limnological and mixing conditions in Lake Arreo and they must be considered in management and restoration plans.  相似文献   

Global photometric properties of Asteroid (4) Vesta observed with Dawn Framing Camera     

Jian-Yang Li  Lucille Le Corre  Stefan E. Schröder  Vishnu Reddy  Brett W. Denevi  Bonnie J. Buratti  Stefano Mottola  Martin Hoffmann  Pablo Gutierrez-Marques  Andreas Nathues  Christopher T. Russell  Carol A. Raymond 《Icarus》2013

Dawn spacecraft orbited Vesta for more than one year and collected a huge volume of multispectral, high-resolution data in the visible wavelengths with the Framing Camera. We present a detailed disk-integrated and disk-resolved photometric analysis using the Framing Camera images with the Minnaert model and the Hapke model, and report our results about the global photometric properties of Vesta. The photometric properties of Vesta show weak or no dependence on wavelengths, except for the albedo. At 554 nm, the global average geometric albedo of Vesta is 0.38 ± 0.04, and the Bond albedo range is 0.20 ± 0.02. The bolometric Bond albedo is 0.18 ± 0.01. The phase function of Vesta is similar to those of S-type asteroids. Vesta’s surface shows a single-peaked albedo distribution with a full-width-half-max ∼17% relative to the global average. This width is much smaller than the full range of albedos (from ∼0.55× to >2× global average) in localized bright and dark areas of a few tens of km in sizes, and is probably a consequence of significant regolith mixing on the global scale. Rheasilvia basin is ∼10% brighter than the global average. The phase reddening of Vesta measured from Dawn Framing Camera images is comparable or slightly stronger than that of Eros as measured by the Near Earth Asteroid Rendezvous mission, but weaker than previous measurements based on ground-based observations of Vesta and laboratory measurements of HED meteorites. The photometric behaviors of Vesta are best described by the Hapke model and the Akimov disk-function, when compared with the Minnaert model, Lommel–Seeliger model, and Lommel–Seeliger–Lambertian model. The traditional approach for photometric correction is validated for Vesta for >99% of its surface where reflectance is within ±30% of global average.  相似文献   

Weathering of ordinary chondrites from the Atacama Desert,Chile, by Mössbauer spectroscopy and synchrotron radiation X‐ray diffraction     

Pablo Munayco  Jimmy Munayco  Roberto R. de Avillez  Millarca Valenzuela  Pierre Rochette  Jérôme Gattacceca  Rosa B. Scorzelli 《Meteoritics & planetary science》2013,48(3):457-473

Some terrestrial areas have climatic and geomorphologic features that favor the preservation, and therefore, accumulation of meteorites. The Atacama Desert in Chile is among the most important of such areas, known as DCA. This desert is the driest on Earth, one of the most arid, uninhabitable localities with semiarid, arid, and hyper‐arid conditions. The meteorites studied here were collected from within the DCA of San Juan and Pampa de Mejillones, located, respectively, in the Central Depression and the Coastal Range of the Atacama Desert. ⁵⁷Fe Mössbauer spectroscopy was used for quantitative analysis of the degree of weathering of the meteorites, through the determination of the proportions of the various Fe‐bearing phases and in particular the amount of oxidized iron in terrestrial alteration products. The abundance of ferric ions in weathered chondrites can be related to specific precursor compositions and to the level of terrestrial weathering. The aim of the study was the identification, quantification, and differentiation of the weathering products in the ordinary chondrites found in the San Juan and the Pampa de Mejillones areas of the Atacama Desert. The ⁵⁷Fe Mössbauer spectroscopy study was complemented by synchrotron radiation X‐ray diffraction and magnetic susceptibility measurements. The results allow a clear differentiation of the rate of weathering in meteorite samples collected from the San Juan versus the Pampa de Mejillones areas of the Atacama Desert.  相似文献   

Trace-element study and age dating of zircon from chromitites of the Bushveld Complex (South Africa)     

Marina Yudovskaya  Judith Kinnaird  Anthony J. Naldrett  Nickolay Rodionov  Anton Antonov  Sergey Simakin  Dmitry Kuzmin 《Mineralogy and Petrology》2013,107(6):915-942

The layered Bushveld Complex hosts a number of chromitite layers, which were found to contain significant amounts of zircon grains compared with adjacent silicate rocks. Cathodoluminescent-dark, partially metamict cores and transparent rims of composite zircon grains were analyzed for trace elements with SIMS and LA-ICPMS techniques. The cores are enriched in REE, Y, Th and U and are characterized by distinctly flatter REE patterns in contrast to those of the rims and transparent homogenous crystals. Zircon from the different stratigraphic units has specific Th/U ratios, the highest of which (1.5–4) occurs in a Merensky Reef zircon core. The Ti content of Bushveld zircon ranges from 12 to 52 ppm correlating to a crystallization temperature range of 760–930 °C. The geochemical characteristics of the first zircon generation are consistent with its high-temperature crystallization as the first major U, Th and REE acceptor from a highly-evolved residue of the high-Mg basalt magma, whereas the rims and coreless crystals have crystallized from percolating intercumulus liquid of new influx of the same magma. U-Pb SHRIMP dating of zircon cores and rims does not reveal a distinguishable difference between their ages indicating the absence of inherited zircon. Concordia ages of 2,051?±?9 Ma (2σ, MSWD?=?0.1) and 2,056?±?5 Ma (2σ, MSWD?=?0.05) for zircons from the Merensky Reef and the Upper Platreef located equally near the top of the Critical Zone are in agreement with published ages for the Merensky Reef. Zircon from the deeper-seated Lower Group, Middle Group and Lower Platreef chromitites yields younger concordia ages that may reflect prolonged late-stage volatile activity.  相似文献   

Web-GIS tool for the management of rural land markets     

Juan Porta  Jorge Parapar  Paula García  Gracia Fernández  Juan Touriño  Ramón Doallo  Francisco Ónega  Inés Santé  Pablo Díaz  David Miranda  Rafael Crecente 《Earth Science Informatics》2013,6(4):209-226

Land abandonment and stagnation of rural markets in the last few years have become one of the main concerns of rural administrations. The use of Web and GIS (Geographic Information System) technologies can help to mitigate the effects of these problems. This paper pro-poses a novel Web-GIS tool with spatial capabilities for the dynamization of rural land markets by encouraging the transfer of land from owners to farmers through the leasing of plots. The system, based on open source software, offers information about the properties, their environment and their owners. It uses standards for handling the geographic information and for communicating with external data sources. This system was used as the basis for the development of SITEGAL, the tool for the management of the Land Bank of Galicia (www.bantegal.com/sitegal). SITEGAL has been operational since 2007 obtaining benefits for both administration and users (farmers and land owners), and promoting the e-Government.  相似文献   

Coseismic Three-Dimensional Displacements Determined Using SAR Data: Theory and an Application Test   总被引：1，自引：0，他引：1  

Pablo J. González  José Fernández  Antonio G. Camacho 《Pure and Applied Geophysics》2009,166(8-9):1403-1424

We present a methodology to construct three-dimensional deformation maps using different parameters that can be extracted using SAR data. We apply the methodology using ENVISAT SAR data before and after the December 26^th, 2003 M_w 6.6 Bam earthquake in Iran to determine spatial quasi-continuous three-dimensional coseismic deformation maps. Two near vertical deformation offset components are computed using Envisat ASAR differential interferometry (DInSAR), while two horizontal deformation offset components are measured using sub-pixel correlation techniques applied to ASAR amplitude images. Using the presented methodology, we combine four unevenly precise independent projections of surface coseismic deformation to obtain the full three-dimensional coseismic deformation field caused by an earthquake. The full 3-D coseismic displacement vector is modeled using elastic deformation models constraining details of the fault geometry and slip distribution at depth. Results from the inversion are interpreted in the framework of the Iranian present-day tectonism. Full exploitation of dense 3-D coseismic surface deformation using SAR data, even for moderate earthquakes, should facilitate distinguishing between different interpretations of the mechanical properties of seismically active areas and within the inherent ambiguity of the geophysical inversion solutions.  相似文献   

Influence of Outcrop Scale Fractures on the Effective Stiffness of Fault Damage Zone Rocks     

W. Ashley Griffith  Pablo F. Sanz  David D. Pollard 《Pure and Applied Geophysics》2009,166(10-11):1595-1627

We combine detailed mapping and microstructural analyses of small fault zones in granodiorite with numerical mechanical models to estimate the effect of mesoscopic (outcrop-scale) damage zone fractures on the effective stiffness of the fault zone rocks. The Bear Creek fault zones were active at depths between 4 and 15 km and localize mesoscopic off-fault damage into tabular zones between two subparallel boundary faults, producing a fracture-induced material contrast across the boundary faults with softer rocks between the boundary faults and intact granodiorite outside the boundary faults. Using digitized fault zone fracture maps as the modeled fault geometries, we conduct nonlinear uniaxial compression tests using a novel finite-element method code as the experimental “laboratory” apparatus. Map measurements show that the fault zones have high nondimensional facture densities (>1), and damage zone fractures anastamose and intersect, making existing analytical effective medium models inadequate for estimation of the effective elastic properties. Numerical experiments show that the damage zone is strongly anisotropic and the bulk response of the fault zone is strain-weakening. Normal strains as small as 2% can induce a reduction of the overall stiffness of up to 75%. Fracture-induced effective stiffness changes are large enough to locally be greater than intact modulus changes across the fault due to juxtaposition of rocks of different lithologies; therefore mesoscopic fracturing is as important as rock type when considering material or bimaterial effects on earthquake mechanics. These results have important implications for earthquake rupture mechanics models, because mesoscopic damage zone fractures can cause a material contrast across the faults as large as any lithology-based material contrast at seismogenic depths, and the effective moduli can be highly variable during a single rupture event.  相似文献   

Evolution of the late Pleistocene Mojanda–Fuya Fuya volcanic complex (Ecuador), by progressive adakitic involvement in mantle magma sources     

Claude Robin  Jean-Philippe Eissen  Pablo Samaniego  Hervé Martin  Minard Hall  Joseph Cotten 《Bulletin of Volcanology》2009,71(3):233-258

The Mojanda–Fuya Fuya Volcanic Complex consists of two nearby volcanoes, Mojanda and Fuya Fuya. The older one, Mojanda volcano (0.6 to 0.2 Ma), was first constructed by andesites and high-silica andesites forming a large stratovolcano (Lower Mojanda). This edifice was capped by a basaltic andesite and andesitic cone (Upper Mojanda), which collapsed later to form a 3-km-wide summit caldera, after large phreatomagmatic eruptions. The Lower Fuya Fuya edifice was constructed by the extrusion of viscous Si-rich andesitic lavas and dacitic domes, and the emission of a thick sequence of pyroclastic-flow and fallout deposits which include two voluminous rhyolitic layers. An intermediate construction phase at Fuya Fuya is represented by a mainly effusive cone, andesitic in composition (San Bartolo edifice), the construction of which was interrupted by a major sector collapse in the Late Pleistocene. Finally, a complex of thick siliceous lavas and domes was emplaced within the avalanche amphitheatre, forming the Upper Fuya Fuya volcanic centre. This paper shows that the general evolution from an effusive to an explosive eruptive style is related to a progressive adakitic contribution to the magma source. Although all the rocks of the complex are included in the medium-K field of continental arcs, the Fuya Fuya suite (61–75 wt.% SiO₂) shows depletion in Y and HREE and high Sr/Y and La/Yb values, compared to the less silicic Mojanda suite (55–66.5 wt.% SiO₂). The Mojanda calc-alkaline suite was generated by partial melting of an adakite-metasomatised mantle source that left a residue with 2% garnet, followed by fractional crystallization of dominant plagioclase + pyroxene + olivine at shallow, intra-crustal depths. For Fuya Fuya, geochemical and mineralogical data suggest either (1) partial melting of a similar metasomatised mantle with more garnet in the residue (4%), followed by fractional crystallization involving plagioclase, amphibole and pyroxene, or (2) mixing of mafic mantle-derived magma from the Mojanda suite and slab melts, followed by the same fractional crystallization process.  相似文献