Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Flexible social organization and high incidence of drifting in the sweat bee, Halictus scabiosae.

The very diverse social systems of sweat bees make them interesting models to study social evolution. Here we focus on the dispersal behaviour and social organization of Halictus scabiosae, a common yet poorly known species of Europe. By combining field observations and genetic data, we show that females have multiple reproductive strategies, which generates a large diversity in the social structure of nests. A detailed microsatellite analysis of 60 nests revealed that 55% of the nests contained the offspring of a single female, whereas the rest had more complex social structures, with three clear cases of multiple females reproducing in the same nest and frequent occurrence of unrelated individuals. Drifting among nests was surprisingly common, as 16% of the 122 nests in the overall sample and 44% of the nests with complex social structure contained females that had genotypes consistent with being full-sisters of females sampled in other nests of the population. Drifters originated from nests with an above-average productivity and were unrelated to their nestmates, suggesting that drifting might be a strategy to avoid competition among related females. The sex-specific comparison of genetic differentiation indicated that dispersal was male-biased, which would reinforce local resource competition among females. The pattern of genetic differentiation among populations was consistent with a dynamic process of patch colonization and extinction, as expected from the unstable, anthropogenic habitat of this species. Overall, our data show that H. scabiosae varies greatly in dispersal behaviour and social organization. The surprisingly high frequency of drifters echoes recent findings in wasps and bees, calling for further investigation of the adaptive basis of drifting in the social insects.

Dispersion of natural arsenic in the Malcantone watershed, Southern Switzerland: field evidence for repeated sorption-desorption and oxidation-reduction processes

In recent years, elevated arsenic concentrations have been found in waters and soils of many, countries, often resulting in a health threat for the local population. Switzerland is not an exception and this paper deals with the release and subsequent fate of arsenic in a 200-km(2) mountainous watershed, characterized by crystalline silicate rocks (gneisses, schists, amphibolites) that contain abundant As-bearing sulfide ore deposits, some of which have been mined for iron and gold in the past. Using analytical methods common for mineralogical, ground water and soil studies (XRD, XRF, XAS-XANES and -EXAFS, electron microprobe, extraction, ICP, AAS with hydride generator, ion chromatography), seven different field situations and related dispersion processes of natural arsenic have been studied: (1) release by rock weathering, (2) transport and deposition by water and ice; (3) release of As to the ground and surface water due to increasing pH; (4) accumulation in humic soil horizons; (5) remobilization by reduction in water-saturated soils and stagnant ground waters; (6) remobilization by using P-rich fertilizers or dung and (7) oxidation, precipitation and dilution in surface waters. Comparison of the results with experimental adsorption studies and speciation diagrams from the literature allows us to reconstruct and identify the typical behavior of arsenic in a natural environment under temperate climatic conditions. The main parameters identified are: (a) once liberated from the primary minerals, sorption processes on Fe-oxy-hydroxides dominate over Al-phases, such as Al-hydroxides or clay minerals and limit the As concentrations in the spring and well waters between 20 and 300 mug/l. (b) Precipitation as secondary minerals is limited to the weathering domain, where the As concentrations are still high and not yet too diluted by rain and soils waters. (c) Although neutral and alkaline pH conditions clearly increase the mobility of As, the main factor to mobilize As is a low redox potential (Eh close or below 0 mV), which favors the dissolution of the Fe-oxy-hydroxides on which the As is sorbed. (d) X-ray absorption spectroscopy (XAS) of As in water-logged humic forest soils indicates that the reduction to As III only occurs at the solid-water interface and that the solid contains As as As V (e) A and Bh horizons of humic cambisols can effectively capture As when As-rich waters flow through them. Complex spatial and temporal variation of the various parameters in a watershed results in repeated mobilization and immobilization of As, which continuously transports As from the upper to the lower part of a watershed and ultimately to the ocean. (C) 2004 Elsevier B.V. All rights reserved.

Comparative testing of a miniature diffusion size classifier to assess airborne ultrafine particles under field conditions

Miniature diffusion size classifiers (miniDiSC) are novel handheld devices to measure ultrafine particles (UFP). UFP have been linked to the development of cardiovascular and pulmonary diseases; thus, detection and quantification of these particles are important for evaluating their potential health hazards. As part of the UFP exposure assessments of highwaymaintenance workers in western Switzerland, we compared a miniDiSC with a portable condensation particle counter (P-TRAK). In addition, we performed stationary measurements with a miniDiSC and a scanning mobility particle sizer (SMPS) at a site immediately adjacent to a highway. Measurements with miniDiSC and P-TRAK correlated well (correlation of r = 0.84) but average particle numbers of the miniDiSC were 30%âeuro"60% higher. This difference was significantly increased for mean particle diameters below 40 nm. The correlation between theminiDiSC and the SMPSduring stationary measurements was very high (r = 0.98) although particle numbers from the miniDiSC were 30% lower. Differences between the three devices were attributed to the different cutoff diameters for detection. Correction for this size dependent effect led to very similar results across all counters.We did not observe any significant influence of other particle characteristics. Our results suggest that the miniDiSC provides accurate particle number concentrations and geometric mean diameters at traffic-influenced sites, making it a useful tool for personal exposure assessment in such settings.

A High-Resolution 2D/3D Seismic Study of a Thrust Fault Zone in Lake Geneva, Switzerland

P130 A HIGH-RESOLUTION 2D/3D SEISMIC STUDY OF A THRUST FAULT ZONE IN LAKE GENEVA SWITZERLAND M. SCHEIDHAUER M. BERES D. DUPUY and F. MARILLIER Institute of Geophysics University of Lausanne 1015 Lausanne, Switzerland Summary A high-resolution three-dimensional (3D) seismic reflection survey has been conducted in Lake Geneva near the city of Lausanne Switzerland where the faulted molasse basement (Tertiary sandstones) is overlain by complex Quaternary sedimentary structures. Using a single 48-channel streamer an area of 1200 m x 600 m was surveyed in 10 days. With a 5-m shot spacing and a receiver spacing of 2.5 m in the inline direction and 7.5 m in the crossline direction, a 12-fold data coverage was achieved. A maximum penetration depth of ~150 m was achieved with a 15 cu. in. water gun operated at 140 bars. The multi-channel data allow the determination of an accurate velocity field for 3D processing, and they show particularly clean images of the fault zone and the overlying sediments in horizontal and vertical sections. In order to compare different sources, inline 55 was repeated with a 30/30 and a 15/15 cu. in. double-chamber air gun (Mini GI) operated at 100 and 80 bars, respectively. A maximum penetration depth of ~450 m was achieved with this source.

Thrusting, extension, and doming during the polyphase tectonometamorphic evolution of the High Himalayan Crystalline Zone in NW India

In the NW Himalaya of India, high-grade metamorphic rocks of the High Himalayan Crystalline Zone (HHCZ) are exposed as a 50 km large dome along the Miyar and Gianbul valleys. This Gianbul dome is cored by migmatitic paragneiss formed at peak conditions around 750 degreesC and 8 kbar, and symmetrically surrounded by sillimanite, kyanite +/- staurolite, garnet, biotite, and chlorite Barrovian mineral zones. Thermobarometric and structural investigations reveal that the Gianbul dome results from a polyphase tectono-metamorphic evolution. The first phase corresponds to the NE-directed thrusting of the Shikar Beh nappe, that is responsible for the Barrovian prograde metamorphic field gradient in the southern limb of the dome. In the northern limb of the dome, the Barrovian prograde metamorphism is the consequence of a second tectonic phase, associated with the SW-directed thrusting of the Nyimaling-Tsarap nappe. Following these crustal thickening events, exhumation and doming of the HHCZ high-grade rocks were controlled by extension along the north-dipping Zanskar Shear Zone, in the frontal part of the Nyimaling-Tsarap nappe, as well as by coeval to late extension along the south-dipping Khanjar Shear Zone, in the southern limb of the Gianbul dome. Rapid syn-convergence extension along both of these detachments induced a nearly isothermal decompression, resulting in a high-temperature/low-pressure metamorphic overprint, as well as enhanced partial melting. Such a rapid exhumation within a compressional orogenic context appears unlikely to be controlled solely by granitic diapirism. Alternatively, large-scale doming in the Himalaya could reflect a sub-vertical ductile extrusion of partially melted rocks.

Microquasars as high energy gamma-ray sources

Galactic microquasars are certainly one of the most recent additions to the field of high energy astrophysics and have attracted increasing interest over the last decade. However, the high energy part of the spectrum of microquasars is the most poorly known, mainly due the lack of sensitive instrumentation in the past. Microquasars are now primary targets for all of the observatories working in the X-ray and gamma-ray domains. They also appear as the possible counterparts for some of the unidentified sources of high-energy gamma-rays detected by the experiment EGRET on board the satellite COMPTON-GRO. This paper provides a general review of the main observational results obtained up to now as well as a summary of the scenarios for production of high-energy gamma-rays at the present moment.

Maize response to nitrogen fertilization timing in two tillage systems in a soil with high organic matter content

No-tillage systems, associated to black oat as preceding cover crop, have been increasingly adopted. This has motivated anticipated maize nitrogen fertilization, transferring it from the side-dress system at the stage when plants have five to six expanded leaves to when the preceding cover crop is eliminated or to maize sowing. This study was conducted to evaluate the effects of soil tillage system and timing of N fertilization on maize grain yield and agronomic efficiency of N applied to a soil with high organic matter content. A three-year field experiment was conducted in Lages, state of Santa Catarina, from 1999 onwards. Treatments were set up in a split plot arrangement. Two soil tillage systems were tested in the main plots: conventional tillage (CT) and no-tillage (NT). Six N management systems were assessed in the split-plots: S1 - control, without N application; S2 - all N (100 kg ha-1) applied at oat desiccation; S3 - all N applied at maize sowing; S4 - all N side-dressed when maize had five expanded leaves (V5 growth stage); S5 - 1/3 of N rate applied at maize sowing and 2/3 at V5; S6 - 2/3 of nitrogen rate applied at maize sowing and 1/3 at V5. Maize response to the time and form of splitting N was not affected by the soil tillage system. Grain yield ranged from 6.0 to 11.8 t ha-1. The anticipation of N application (S2 and S3) decreased grain yield in two of three years. In the rainiest early spring season (2000/2001) of the experiment, S4 promoted an yield advantage of 2.2 t ha-1 over S2 and S3. Application of total N rate before or at sowing decreased the number of kernels produced per ear in 2000/2001 and 2001/2002 and the number of ears produced per area in 2001/2002, resulting in reduced grain yield. The agronomic efficiency of applied N (kg grain increase/kg of N applied) ranged from 13.9 to 38.8 and was always higher in the S4 than in the S2 and S3 N systems. Short-term N immobilization did not reduce grain yield when no N was applied before or at maize sowing in a soil with high organic matter content, regardless of the soil tillage system.

The central Sesia-Lanzo zone (Western Italian Alps) - New field observations and lithostratigraphic subdivisions

An extensive study of the central part of the Sesia Lanzo Zone has been undertaken to identify pre-Alpine protoliths and to reconstruct the lithologic and tectonic setting of this part of the Western Alps. Three main complexes have been defined: 1) the Polymetamorphic Basement Complex, corresponding to the lower unit of the Sesia Lanzo Zone after COMPAGNONI et al. (1977), is further subdivided into the three following units: a) an Internal Unit characterized by eo-Alpine high pressure (HP) assemblages (DAL PIAZ et al., 1972) (Eclogitic Micaschists); b) an Intermediate Unit where HP parageneses are partially re-equilibrated under greenschist conditions and c) an External Unit where the main foliation is defined by a greenschist paragenesis (Gneiss Minuti auct.). 2) the Monometamorphic Cover Complex, subdivided into the followings: a) the Bonze Unit, composed of sheared metagabbros, eclogitized metabasalts with MORB geochemical affinity and related metasediments (micaschists, quartzites and Mn-cherts) and b) the Scalaro Unit, containing predominantly metasediments of supposed Permo-Triassic age (yellow dolomitic marbles, calcschists and conglomeratic limestones, micaschists and quartzites with thin levels of basic rocks with within plate basalts [WPB] geochimical affinity). Multiple lithostratigraphic sequences for the Monometamorphic Cover Complex are proposed. The contact between the Bonze and Scalaro Units is defined by repetitions of dolomitic marbles and metabasalts; the ages of the metasediments have been assigned solely by analogy with other sediments of the Western Alps, due to the absence of fossils. The Monometamorphic Cover Complex can be considered as the autochthonous cover of the Sesia Lanzo Zone because of the primary contacts with the basement and because of the presence of preAlpine HT basement blocks in the cover sequences. 3) The pre-Alpine high temperature (HT) Basement Complex (or `'Seconda Zona Diorito-Kinzigitica''), comprises HT Hercynian rocks like kinzigites, amphibolites, granulites and calcite marbles; this Complex is always located between the Internal and the External Units and can be followed continuously for several kilometers south of the Gressoney Valley to the Orco Valley. A schematic evolution for the Sesia Lanzo Zone is proposed; based on available data together with new geochronological data, this study shows that the internal and external parts of the polymetamorphic basement of the Sesia Zone experienced different cooling histories .

Oscillatory magnetic tweezers based on ferromagnetic beads and simple coaxial coils

We report the design and validation of simple magnetic tweezers for oscillating ferromagnetic beads in the piconewton and nanometer scales. The system is based on a single pair of coaxial coils operating in two sequential modes: permanent magnetization of the beads through a large and brief pulse of magnetic field and generation of magnetic gradients to produce uniaxial oscillatory forces. By using this two step method, the magnetic moment of the beads remains constant during measurements. Therefore, the applied force can be computed and varies linearly with the driving signal. No feedback control is required to produce well defined force oscillations over a wide bandwidth. The design of the coils was optimized to obtain high magnetic fields (280 mT) and gradients (2 T/m) with high homogeneity (5% variation) within the sample. The magnetic tweezers were implemented in an inverted optical microscope with a videomicroscopy-based multiparticle tracking system. The apparatus was validated with 4.5 ¿m magnetite beads obtaining forces up to ~2 pN and subnanometer resolution. The applicability of the device includes microrheology of biopolymer and cell cytoplasm, molecular mechanics, and mechanotransduction in living cells.

Residual Stress Measurement on a MEMS Structure With High-Spatial Resolution

A new approach to the local measurement of residual stress in microstructures is described in this paper. The presented technique takes advantage of the combined milling-imaging features of a focused ion beam (FIB) equipment to scale down the widely known hole drilling method. This method consists of drilling a small hole in a solid with inherent residual stresses and measuring the strains/displacements caused by the local stress release, that takes place around the hole. In the presented case, the displacements caused by the milling are determined by applying digital image correlation (DIC) techniques to high resolution micrographs taken before and after the milling process. The residual stress value is then obtained by fitting the measured displacements to the analytical solution of the displacement fields. The feasibility of this approach has been demonstrated on a micromachined silicon nitride membrane showing that this method has high potential for applications in the field of mechanical characterization of micro/nanoelectromechanical systems.

Persistence of a biocontrol Pseudomonas inoculant as high populations of culturable and non-culturable cells in 200-cm-deep soil profiles

Little is known about the ecology of soil inoculants used for pathogen biocontrol, biofertilization and bioremediation under field conditions. We investigated the persistence and the physiological states of soil-inoculated Pseudomonas protegens (previously Pseudomonas fluorescens) CHA0 (108 CFU g−1 surface soil) in different soil microbial habitats in a planted ley (Medicago sativa L.) and an uncovered field plot. At 72 days, colony counts of the inoculant were low in surface soil (uncovered plot) and earthworm guts (ley plot), whereas soil above the plow pan (uncovered plot), and the rhizosphere and worm burrows present until 1.2 m depth (ley plot) were survival hot spots (105-106 CFU g−1 soil). Interestingly, strain CHA0 was also detected in the subsoil of both plots, at 102-105 CFU g−1 soil between 1.8 and 2 m depth. However, non-cultured CHA0 cells were also evidenced based on immunofluorescence microscopy. Kogure's direct viable counts of nutrient-responsive cells showed that many more CHA0 cells were in a viable but non-culturable (VBNC) or a non-responsive (dormant) state than in a culturable state, and the proportion of cells in those non-cultured states depended on soil microbial habitat. At the most, cells in a VBNC state amounted to 34% (above the plow pan) and those in a dormant state to 89% (in bulk soil between 0.6 and 2 m) of all CHA0 cells. The results indicate that field-released Pseudomonas inoculants may persist at high cell numbers, even in deeper soil layers, and display a combination of different physiological states whose prevalence fluctuates according to soil microbial habitats.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

Entropy change and magnetocaloric effect in Gd5(SixGe1-x)4

Isothermal magnetization curves up to 23 T have been measured in Gd5Si1.8Ge2.2. We show that the values of the entropy change at the first-order magnetostructural transition, obtained from the Clausius-Clapeyron equation and the Maxwell relation, are coincident, provided the Maxwell relation is evaluated only within the transition region and the maximum applied field is high enough to complete the transition. These values are also in agreement with the entropy change obtained from differential scanning calorimetry. We also show that a simple phenomenological model based on the temperature and field dependence of the magnetization accounts for these results.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.