ALBERT

Beschreibung: Directions recovered from palaeomagnetic samples are usually archived with some quantitative information about their precision, most often in the form of a so-called α 95 angle. Such angles are classically co-estimated with the recovered palaeomagnetic direction from a collection of samples providing individual estimates of this direction. In some instances, however, palaeomagnetic directions have to be inferred from a single sample in which case no α 95 angle can be recovered in this way. Fortunately, the progressive demagnetization techniques and principal component analysis universally used to recover directional information from single samples provide alternative measures of the error affecting the recovered direction, known as Maximum Angular Deviation (MAD) angles. These have so far only been considered as rough quality indicators. Here, however, we show that directions recovered in this way can be assumed to satisfy a Fisher distribution, and that the corresponding MAD angles can be rescaled into α 95 estimates by multiplying it by an appropriate factor, which only depends on the number of demagnetization steps used in the principal component analysis and on whether one relies on a standard or a so-called ‘anchored’ principal component analysis. These coefficients have been tabulated and practical recommendations for taking advantage of them outlined in the final section of the text. They provide simple means for users to produce much needed error bars on declination and inclination time series recovered from sedimentary long sequences.

Beschreibung: We investigate the effect of crop price and climate variables on rainfed corn and soybean yields and acreage in the United States using a large panel dataset for the 1977–2007 period. Instrumental variables are used to control for endogeneity of prices in yield and acreage regressions, while allowing for spatially auto-correlated errors. We find that an increase in corn price has a statistically significant positive impact on corn yield, but the effect of soybean price on soybean yields is not statistically significant. The estimated price elasticities of corn yield and acreage are 0.23 and 0.45, respectively. Of the increase in corn supply caused by an increase in corn price, we find that 33.8% is due to price-induced yield enhancement and 66.2% is due to price-induced acreage expansion. We also find that the impact of climate change on corn production ranges from $-$ 7% to $-$ 41% and on soybean ranges from $-$ 8% to $-$ 45%, depending on the climate change scenarios, time horizon, and global climate models used to predict climate change. We show that the aggregate net impact of omitting price variables is an overestimation of the effect of climate change on corn yield by up to 9% and on soybean yield by up to 15%.

Beschreibung: The downward continuation of the observed geomagnetic field from the Earth's surface to the core–mantle boundary (CMB) is complicated due to induction and diffusion processes in the electrically conducting Earth mantle, which modify the amplitudes and morphology of the geomagnetic field. Various methods have been developed to solve this problem, for example, the perturbation approach by Benton & Whaler, or the non-harmonic downward continuation by Ballani et al. In this paper, we present a new approach for determining the geomagnetic field at the CMB by reformulating the ill-posed, one-sided boundary-value problem with time-variable boundary-value function on the Earth's surface into an optimization problem for the boundary condition at the CMB. The reformulated well-posed problem is solved by a conjugate gradient technique using the adjoint gradient of a misfit. For this purpose, we formulate the geomagnetic adjoint-state equations for efficient computations of the misfit gradient. Beside the theoretical description of the new adjoint-state method (ASM), the first applications to a global geomagnetic field model are presented. The comparison with other methods demonstrates the capability of the new method to determine the geomagnetic field at the CMB and allows us to investigate the variability of the determined field with respect to the applied methods. This shows that it is necessary to apply the ASM when investigating the effect of the Earth's mantle conductivity because the difference between the results of approximate methods (harmonic downward continuation, perturbation approach) and the rigorous ASM are of the same order as the difference between the results of the ASM applied for different mantle conductivities.

Beschreibung: Recent first principles calculations of the Earth's outer core thermal and electrical conductivities have raised their values by a factor of three. This has significant implications for geodynamo operation, in particular, forcing the development of a stably stratified layer at the core–mantle boundary (CMB). This study seeks to test the hypothesis of a stably stratified layer in the uppermost core by analysing geomagnetic observations made by the CHAMP satellite. An inversion method is utilized that jointly solves for the time-dependent main field and the core surface flow, where we assume the temporal variability of the main field, its secular variation (SV), to be entirely due to advective motion within the liquid outer core. The results show that a large-scale pure toroidal flow, consistent with a stably stratified layer atop the outer core, is not compatible with the observed magnetic field during the CHAMP era. However, allowing just a small amount of poloidal flow leads to a model fitting the observations satisfactorily. As this poloidal flow component is large scale, within a predominantly toroidal, essentially tangentially geostrophic flow, it is compatible with a stably stratified upper outer core. Further, our assumption of little or no diffusive SV may not hold, and a small amount of SV generated locally by diffusion might lead to a large-scale pure toroidal flow providing an acceptable fit to the data.

Beschreibung: We present the first inversion of geomagnetic Sq data in a framework of 3-D conductivity models. This problem has been considered as immensely difficult due to the complex spatial structure of the Sq source which, in addition, varies with season and solar activity. Recently, we developed a 3-D electromagnetic (EM) inversion solution that allows one to work in a consistent manner with data that originates from sources, irrespective of their spatial complexity. In this paper, we apply our 3-D EM inversion scheme to Sq data collected during the Australian Wide Array of Geomagnetic Stations project. Within this project, three components of the geomagnetic field were recorded between 1989 November and 1990 December with the use of 53 portable vector magnetometers. The instruments were distributed over the Australian mainland with an average spacing of 275 km between sites. Inverting this unique—in a sense of its spatial regularity, density and long operational time—data set, we recovered the 3-D conductivity distribution beneath Australia at upper mantle depths (100–520 km). This depth range was justified in the paper from resolution studies using checkerboard tests. In addition, we performed extensive modelling to estimate quantitatively the influence of various factors on Sq signals, namely from hypothetical anomalies, inaccuracy in the source, ocean, and model discretization. As expected, the ocean (coastal) effect appeared to be the largest so that it has to be accounted for during 3-D inversion as accurately as possible. Our 3-D inversions—of data from either single or multiple days—revealed a strong offshore conductor near the south-east coast of Australia, which persists at all considered depths. Varying in details, this anomaly is remarkably robust irrespective of the considered day(s). We compared our results to those obtained from a different inversion scheme and an independent induction data set, and observed encouraging similarity. Combination of the two results suggests, that this conductor continues to the base of the mantle transition zone at 660 km. The nature of this anomaly is not fully understood but one possible explanation is that it is attributed to a reservoir responsible for three hotspots in the region.

Beschreibung: A detailed palaeomagnetic, rock-magnetic and palaeointensity study has been carried out on a Miocene volcanic sequence which consists of 39 consecutive lava flows recording a polarity transition in La Gomera (Canary Islands, Spain). In addition, new 40 Ar/ 39 Ar ages were obtained in two flows, yielding 9.63 ± 0.06 Ma in the lower and 9.72 ± 0.08 Ma in the upper part of the sequence. Palaeomagnetic results allowed determining a ChRM direction in all studied lavas: The 25 lowermost flows of the sequence display normal polarity directions and above, a sequence of 14 flows correspond to a transitional geomagnetic regime. If considered together with palaeomagnetic results from a previous study, which were obtained on the flows immediately overlying the upper part of the sequence analysed in this work, these results indicate that the reversal recorded in the Hermigua sequence corresponds to the normal to reverse C4Ar2n to C4Ar3r polarity transition. The lower-lying 25 normal-polarity flows yield a mean direction D  = 359.6°, I  = 42.4° ( α 95  = 5.1°; k  = 33) which agrees well with the expected values. Above, a sequence of 14 flows displays a more irregular directional behaviour, including several transitional directions, suggesting the occurrence of a precursor to the transition. The presence of a virtual geomagnetic pole (VGP) cluster in the western Atlantic Ocean observed in this study coincides with previous records of Miocene transitions. Interestingly, this region corresponds to a near-radial flux centre of the present-day non axial dipole field. Angular dispersion of VGPs calculated for the 25 lowermost normal polarity flows of the sequence shows a lower than expected result. Palaeointensity determinations were carried out using a Thellier type double heating method. 27 of the 48 analysed samples measured yielded successful results. Mean VDMs mean values range from 1.1 ± 0.5 to 8.8 ± 0.9 x 10 22 Am 2 . The intensity values decrease significantly on approaching the directional transitional zone, suggesting an earlier start of the polarity transition in the intensity record, typical of a decreasing axial dipole.

Beschreibung: We present a new grand mean palaeomagnetic pole (Plong: 222.1°, Plat: –64.0°, A 95 : 2.6°, N = 49) for the ca. 1110 Ma Umkondo large igneous province (LIP) of the Kalahari Craton. New palaeomagnetic data from 24 sills in Botswana and compiled reprocessed existing data are used to develop a palaeomagnetic pole as the Fisher mean of cooling unit virtual geomagnetic poles (VGPs). The mean and its associated uncertainty provide the best-constrained pole yet developed for the province. Comparing data from individual cooling units allows for evaluation of palaeosecular variation at this time in the Mesoproterozoic. The elongation of the population of VGPs is consistent with that predicted by the TK03.GAD model lending support to the dipolar nature of the field in the late Mesoproterozoic. In our new compilation, 4 of 59 (~7 per cent) of the igneous units have northerly declinations while the rest are south-directed indicating that a geomagnetic reversal occurred during magmatic activity. Interpreting which of these polarities corresponds with a normal or reversed geomagnetic field relative to other continents can constrain the relative orientations between cratons with time-equivalent data. This interpretation is particularly important in comparison to Laurentia as it bears on Kalahari's involvement and position in the supercontinent Rodinia. The dominance of south-directed declinations within the Umkondo Province was previously used to suggest that these directions are the same polarity as reversed directions from the early magmatic stage of the Keweenawan Midcontinent Rift of Laurentia. Two Umkondo sills with northerly declinations have U-Pb baddeleyite ages of ca. 1109 Ma that are temporally close to dated Midcontinent Rift units having reversed directions. Based on this comparison, and palaeomagnetic data from younger units in the Kalahari Craton, we favour the option in which the sites with northerly declinations from the Umkondo Province correspond to the reversed polarity directions from the early magmatic stage in the Midcontinent Rift. This interpretation allows for the Namaqua-Natal metamorphic belt of Kalahari to be a conjugate to the Grenville margin of North America and for Kalahari to have become conjoined with Laurentia within the supercontinent Rodinia subsequent to Umkondo LIP magmatic activity.

Beschreibung: The northeastern Caribbean island arc, which materializes the boundary between the North American and Caribbean plates, is particularly exposed to large earthquakes and tsunamis. The low level of preparedness of a large part of its population and the lack of risk reduction provisions in public policies in many countries of the region put their population and economy at high risk in case of large telluric events. Here, we investigate the impact of three possible earthquake scenarios, consistent with the regional seismotectonic setting, on northern Haiti through inundation by tsunami waves. These scenarios simulate the effect of a M 8.0 earthquake on the Septentrional strike-slip fault (possibly similar to the 1842 earthquake), a M 8.1 earthquake on the offshore thrust fault system north of Haiti, and an earthquake rupturing a large portion of the offshore thrust fault system north of Haiti and the Dominican Republic. We calculate run-up heights along the northern coast of Haiti, in particular in the densely populated Cap Haitien. We find that the rupture of the offshore North Hispaniola thrust fault could result in wave heights up to 10 m with inundation up to 4 km inland, with only 10–15 min between ground shaking and the first wave arrivals. The city of Cap Haitien is particularly exposed, with potential flooding of most of the city and its suburbs, including the international airport. We also find that the historical reports available for the 1842 earthquake, when compared to our simulations, favor a rupture of the North Hispaniola thrust fault, although much uncertainty remains. If the 1842 earthquake did not rupture the Septentional fault offshore Haiti, then it is currently capable of at least a M w 7.7 earthquake, significantly larger than previously thought. The simulations presented here provide a basis for developing conservative maps of run-up heights that can be transferred, with added factors of safety, into practical implementation for tsunami preparedness and protection.

Beschreibung: We present the results of palaeomagnetic analysis on Late Bronge Age pottery from Santorini carried out in order to estimate the thermal effect of the Minoan eruption on the pre-Minoan habitation level. A total of 170 specimens from 108 ceramic fragments have been studied. The ceramics were collected from the surface of the pre-Minoan palaeosol at six different sites, including also samples from the Akrotiri archaeological site. The deposition temperatures of the first pyroclastic products have been estimated by the maximum overlap of the re-heating temperature intervals given by the individual fragments at site level. A new statistical elaboration of the temperature data has also been proposed, calculating at 95 per cent of probability the re-heating temperatures at each site. The obtained results show that the precursor tephra layer and the first pumice fall of the eruption were hot enough to re-heat the underlying ceramics at temperatures 160–230 °C in the non-inhabited sites while the temperatures recorded inside the Akrotiri village are slightly lower, varying from 130 to 200 °C. The decrease of the temperatures registered in the human settlements suggests that there was some interaction between the buildings and the pumice fallout deposits while probably the buildings debris layer caused by the preceding and syn-eruption earthquakes has also contributed to the decrease of the recorded re-heating temperatures.

Beschreibung: We have developed an algorithm, which we call HexMT, for 3-D simulation and inversion of magnetotelluric (MT) responses using deformable hexahedral finite elements that permit incorporation of topography. Direct solvers parallelized on symmetric multiprocessor (SMP), single-chassis workstations with large RAM are used throughout, including the forward solution, parameter Jacobians and model parameter update. In Part I, the forward simulator and Jacobian calculations are presented. We use first-order edge elements to represent the secondary electric field ( E ), yielding accuracy O ( h ) for E and its curl (magnetic field). For very low frequencies or small material admittivities, the E -field requires divergence correction. With the help of Hodge decomposition, the correction may be applied in one step after the forward solution is calculated. This allows accurate E -field solutions in dielectric air. The system matrix factorization and source vector solutions are computed using the MKL PARDISO library, which shows good scalability through 24 processor cores. The factorized matrix is used to calculate the forward response as well as the Jacobians of electromagnetic (EM) field and MT responses using the reciprocity theorem. Comparison with other codes demonstrates accuracy of our forward calculations. We consider a popular conductive/resistive double brick structure, several synthetic topographic models and the natural topography of Mount Erebus in Antarctica. In particular, the ability of finite elements to represent smooth topographic slopes permits accurate simulation of refraction of EM waves normal to the slopes at high frequencies. Run-time tests of the parallelized algorithm indicate that for meshes as large as 176 x 176 x 70 elements, MT forward responses and Jacobians can be calculated in ~1.5 hr per frequency. Together with an efficient inversion parameter step described in Part II, MT inversion problems of 200–300 stations are computable with total run times of several days on such workstations.

Beschreibung: Following the creation described in Part I of a deformable edge finite-element simulator for 3-D magnetotelluric (MT) responses using direct solvers, in Part II we develop an algorithm named HexMT for 3-D regularized inversion of MT data including topography. Direct solvers parallelized on large-RAM, symmetric multiprocessor (SMP) workstations are used also for the Gauss–Newton model update. By exploiting the data-space approach, the computational cost of the model update becomes much less in both time and computer memory than the cost of the forward simulation. In order to regularize using the second norm of the gradient, we factor the matrix related to the regularization term and apply its inverse to the Jacobian, which is done using the MKL PARDISO library. For dense matrix multiplication and factorization related to the model update, we use the PLASMA library which shows very good scalability across processor cores. A synthetic test inversion using a simple hill model shows that including topography can be important; in this case depression of the electric field by the hill can cause false conductors at depth or mask the presence of resistive structure. With a simple model of two buried bricks, a uniform spatial weighting for the norm of model smoothing recovered more accurate locations for the tomographic images compared to weightings which were a function of parameter Jacobians. We implement joint inversion for static distortion matrices tested using the Dublin secret model 2, for which we are able to reduce nRMS to ~1.1 while avoiding oscillatory convergence. Finally we test the code on field data by inverting full impedance and tipper MT responses collected around Mount St Helens in the Cascade volcanic chain. Among several prominent structures, the north–south trending, eruption-controlling shear zone is clearly imaged in the inversion.

Beschreibung: Palaeointensity variation is investigated for an inferred time period spanning from 2.34 to 1.96 Ma. Twenty-nine consecutive lava flows are sampled along cliffs 350 m high generated by normal faulting on the Dobi section of Afar depression, Ethiopia. Magnetostratigraphy and K–Ar measurements indicate a lava sequence of R–N–R–N geomagnetic field polarities in ascending order; the lower normal polarity is identified as the Réunion Subchron. Reliability of palaeomagnetic data is ascertained through careful thermal demagnetization and by the reversal test. The Tsunakawa–Shaw method yielded 70 successful palaeointensity results from 24 lava flows and gave 11 acceptable mean palaeointensities. Reliability in palaeointensity data is ascertained by the similar values obtained by the IZZI–Thellier method and thus 11 reliable mean values are obtained from our combined results. After the older reverse polarity with the field intensity of 19.6 ± 7.8 μT, an extremely low palaeointensity period with an average of 6.4 μT is shown to occur prior to the Réunion Subchron. During the Réunion Subchron, the dipole field strength is shown to have returned to an average of 19.5 μT, followed by second extreme low of 3.6 μT and rejuvenation with 17.1 ± 5.3 μT in the younger reverse polarity. This ‘W-shape’ palaeointensity variation is characterized by occurrences of two extremely weak fields lower than 8 μT prior to and during the Réunion Subchron and a relatively weak time-averaged field of approximately 15 μT. This feature is also found in sedimentary cores from the Ontong Java Plateau and the north Atlantic, indicative of a possibly global geomagnetic field phenomenon rather than a local effect on Ethiopia. Furthermore, we estimate a weak virtual axial dipole moment of 3.66 (±1.85) x 10 22 Am 2 during early stage of the Matuyama Chron (inferred time period of 2.34–1.96 Ma).

Beschreibung: Adopting a model with two half-spaces that consist of solid and porous materials, we numerically investigate the seismoelectric conversion at the solid–porous interface. First, the wave fields in a low-porosity two-layer model are compared with those in a homogeneous full-space model. The consistency of seismic waves is a validation of our program. We are interested in the quasi-coseismic electromagnetic (EM) signals recorded in the solid area near the interface because they seemingly accompany seismic waves. Then, further numerical simulations on an ordinary two-layer model are conducted. On the basis of time slice snapshots and theoretical analysis, we determine that quasi-coseismic EM signals are essentially non-coseismic EM fields, which include radiation and evanescent EM waves. Evanescent EM waves are induced by the seismic waves that arrive at the interface with the incident angle greater than the critical angle. These waves decay faster than radiation EM waves when moving away from the interface. In the porous layer, evanescent EM waves can hardly be recognized unless they are separated from coseismic EM signals. This finding can be the reason why evanescent EM waves have not been identified in previous seismoelectric studies. Awareness of the fact that seismoelectric conversion at an interface can generate evanescent and EM waves is likely to result in a comprehensive understanding and improved interpretation of the seismoelectric coupling phenomenon.

Beschreibung: We present a new technique for modelling the global lithospheric magnetic field at Earth's surface based on the estimation of equivalent potential field sources. As a demonstration we show an application to magnetic field measurements made by the CHAMP satellite during the period 2009–2010 when it was at its lowest altitude and solar activity was quiet. All three components of the vector field data are utilized at all available latitudes. Estimates of core and large-scale magnetospheric sources are removed from the measurements using the CHAOS-4 model. Quiet-time and night-side data selection criteria are also employed to minimize the influence of the ionospheric field. The model for the remaining lithospheric magnetic field consists of magnetic equivalent potential field sources (monopoles) arranged in an icosahedron grid at a depth of 100 km below the surface. The corresponding model parameters are estimated using an iteratively reweighted least-squares algorithm that includes model regularization (either quadratic or maximum entropy) and Huber weighting. Data error covariance matrices are implemented, accounting for the dependence of data variances on quasi-dipole latitude. The resulting equivalent source lithospheric field models show a degree correlation to MF7 greater than 0.7 out to spherical harmonic degree 100. Compared to the quadratic regularization approach, the entropy regularized model possesses notably lower power above degree 70 and a lower number of degrees of freedom despite fitting the observations to a very similar level. Advantages of our equivalent source method include its local nature, the possibility for regional grid refinement and the production of local power spectra, the ability to implement constraints and regularization depending on geographical position, and the ease of transforming the equivalent source values into spherical harmonics.

Beschreibung: Previous works revealed a close relationship between magnetic susceptibility (MS) and heavy metal (HM) contents originating from industrial sources. However, despite general statements on the usefulness of magnetic mapping, the benefit of this procedure for geochemistry was not quantified yet. We present a study on fly ash pollution in soil around a coal-burning power plant complex and simulate a stepwise approach of magnetic pre-screening and subsequent targeted sampling for chemical analysis. The aim of this study is not to discuss correlations between MS and HM, but to show that a combined stepwise magnetic-chemical approach is the most efficient way for outlining HM contamination. In order to provide quantitative evidence, we explored map similarities of spatial HM distributions based on magnetochemical data and chemical data only. We determined 3-D triangular planes defined by categorized HM values at the sampling coordinates and calculated the average dihedral angle of the normal vectors as a similarity result. The study shows that the ‘Targeted’ HM map (selection of 30 sites based magnetic pre-screening) has a higher similarity with the ‘True’ Pollution HM map (85 sites) than HM maps resulting from site selections (30 sites) without using magnetic pre-screening information.

Beschreibung: Archaeological baked clay remains provide valuable information about the geomagnetic field in historical past, but determination of the geomagnetic field characteristics, especially intensity, is often a difficult task. This study was undertaken to elucidate the reasons for unsuccessful intensity determination experiments obtained from two different Bulgarian archaeological sites (Nessebar – Early Byzantine period and Malenovo – Early Iron Age). With this aim, artificial clay samples were formed in the laboratory and investigated. The clay used for the artificial samples preparation differs according to its initial state. Nessebar clay was baked in the antiquity, but Malenovo clay was raw, taken from the clay deposit near the site. The obtained artificial samples were repeatedly heated eight times in known magnetic field to 700 °C. X-ray diffraction analyses and rock-magnetic experiments were performed to obtain information about the mineralogical content and magnetic properties of the initial and laboratory heated clays. Two different protocols were applied for the intensity determination—Coe version of Thellier and Thellier method and multispecimen parallel differential pTRM protocol. Various combinations of laboratory fields and mutual positions of the directions of laboratory field and carried thermoremanence were used in the performed Coe experiment. The obtained results indicate that the failure of this experiment is probably related to unfavourable grain sizes of the prevailing magnetic carriers combined with the chosen experimental conditions. The multispecimen parallel differential pTRM protocol in its original form gives excellent results for the artificial samples, but failed for the real samples (samples coming from previously studied kilns of Nessebar and Malenovo sites). Obviously the strong dependence of this method on the homogeneity of the used subsamples hinders its implementation in its original form for archaeomaterials. The latter are often heterogeneous due to variable heating conditions in the different parts of the archaeological structures. The study draws attention to the importance of multiple heating for the stabilization of grain size distribution in baked clay materials and the need of elucidation of this question.

Beschreibung: Magnetic minerals acquire a remanent magnetization when cooled or warmed in an applied field H through a phase transition. This paper reports the first observations of transition cooling and warming remanent magnetizations (TrCRM and TrWRM) in pyrrhotite (Fe 7 S 8 ) and hematite (αFe 2 O 3 ), with phase transitions at 32 K (Besnus transition) and 250 K (Morin transition). TrWRM was produced by warming in a 2 mT field from 10 to 300 K and TrCRM by cooling in 2 mT from 300 to 10 K. In both experiments, magnetization M was measured at 5 K intervals. The largest changes in M occurred over a 20 K interval flanking the Besnus or Morin transition but substantial changes also occurred away from any transition, presumably due to continuous changes in anisotropy. The triclinic low-temperature phase of pyrrhotite has high remanence efficiency compared to the room-temperature monoclinic phase. When H was zeroed at 300 K, only 35–40% of M survived as TrWRM, but in TrCRM experiments when H was zeroed at 10 K, practically 100% of M was preserved. After zero-field cooling to 10 K, 50–60% of TrWRM was retained but 〈15% of TrCRM survived zero-field warming to 300 K. In hematite, both TrWRM of the room-temperature spin-canted phase and TrCRM of the low-temperature defect phase exist, as shown by net remanence production following complete warming–cooling or cooling–warming cycles between 10 and 300 K. Details of remanence acquisition/demagnetization in crossing the Morin transition could not be resolved, however. Symmetries were observed for both minerals. Single-domain size pyrrhotite crystals had mirror-image field-off and field-on warming curves below 32 K, during TrWRM production/TrCRM loss at the 32 K transition, and up to 150 K. Submicron single-domain hematite had almost complete symmetry between field-off and field-on curves during cooling as well as warming, both below and above the Morin transition but not in the transition region itself. This is the first time that mirror-image symmetries, observed previously for pseudo-single-domain and multidomain magnetite, have been documented for mineral grains in a single-domain state.

Beschreibung: A computational method is given for solving the forward modelling problem for transient electromagnetic exploration. Its key features are the discretization of the quasi-static Maxwell's equations in space using the first-kind family of curl-conforming Nédélec elements combined with time integration using rational Krylov methods. We show how rational Krylov methods can also be used to solve the same problem in the frequency domain followed by a synthesis of the transient solution using the fast Hankel transform, and we argue that the pure time-domain solution is more efficient. We also propose a new surrogate optimization approach for selecting the pole parameters of the rational Krylov method which leads to convergence within an a priori determined number of iterations independent of mesh size and conductivity structure. These poles are repeated in a cyclic fashion, which, in combination with direct solvers for the discrete problem, results in significantly faster solution times than previously proposed schemes.

Beschreibung: This article analyses the sequence of changes in land used for milk production on dairy farms before, during and towards the abolition of milk quotas. Using a unique dataset comprising farm level data of the Netherlands between 1971 and 2011 we estimate two duration models, analysing the time period between increases and decreases in dairy land use. The impact of milk quota, socio-economic, farm income and economic-political variables on the likelihood of a farm changing its land use are assessed. Results show that changes are highly farm specific, but that quota abolition will lead to a more dynamic dairy sector.

Beschreibung: In 2010, a wide-angle seismic reflection/refraction profile was acquired along the Hecataeus Rise, an area of shallow seabed immediately south of Cyprus in the eastern Mediterranean. The profile crossed from the Hecataeus Rise, through the Cyprus Arc to the Levantine Basin beyond. Due to the short length of the profile and the corresponding lack of deep ray coverage, velocity modelling was complemented by gravity modelling to gain constraints on deep crustal structure. The resultant model reveals velocities for the Hecataeus Rise that show no evidence of shallow ophiolites like those seen on mainland Cyprus, and the velocities are not diagnostic of a unique crustal affinity. Low-velocity sediments make up at least 7 km of the upper structure of Hecataeus Rise and these sediments overlie a two-layered crust. From the gravity modelling, the combined sediments and crust of Hecataeus Rise appear to be thinner than the Eratosthenes Seamount block to the southwest. A high-velocity lower crustal block is modelled under the seaward edge of Hecataeus Rise and, based on the gravity modelling, is inferred to extend landwards beneath the Rise. Similar high-velocity blocks were identified on the southwestern edge of Hecataeus Rise along nearby refraction lines and were interpreted as remnant Tethyan oceanic crust, foundered in the Cyprus Arc, along which subduction has ceased in this area. Given the thin two-layered crust beneath a thick accumulation of sediments modelled for Hecataeus Rise, we interpret that Hecataeus Rise represents a collage of oceanic fragments, accreted together within the failed subduction zone. Outboard of the crust of Hecataeus Rise, a 5-km deep low-velocity basin, possibly an accretionary wedge, is imaged that appears to correspond with the Cyprus Arc deformation zone imaged on both coincident and along-strike seismic reflection lines. A similar and wider feature is observed on seismic refraction lines to the west and combined, these may be revealing an eastward tapering zone of crustal deformation. To the south of the profile, the Levantine Basin appears undisturbed by the collision to the north and exhibits a uniform and homogeneous velocity structure.

Beschreibung: The Rogaland Igneous Complex (RIC) in southern Norway intruded into post-Sveconorwegian granulite facies crust ~930 Ma. It includes three massif anorthosites, several small leuconorite bodies and the ~7 km thick norite-quartz mangerite layered Bjerkreim-Sokndal (BKS) intrusion. The intrusion consists of five rhythmic megaunits created by repeated magma influxes topped by a transition zone and more evolved mangerites and quartz mangerites. Over 70 palaeomagnetic sites have been collected in the BKS, sampling all the megacyclic subunits and overlying mangerites. Remanence within the BKS is held in hemo-ilmenite-only rocks (lower parts of the megacyclic units), mixed hemo-ilmenite and magnetite rocks (upper parts of the lower megacyclic units) and magnetite only rocks in the upper highest megacyclic unit and overlying mangerites. Due to the different oxides present magnetic susceptibility varies over four orders of magnitude with a bimodal distribution (mean susceptibility of 6.4 x 10 –3 SI for hemo-ilmenite rocks, and 8.7 x 10 –2 SI for magnetite rocks). NRM values do not show a strong bimodal distribution as many of the rocks lacking magnetite have hemo-ilmenite with strong lamellar magnetism; average NRM for the entire suite is 8.83 A m –1 . All sites within the cyclic part of the intrusion have stable remanence and produce well-clustered site means. Samples from the upper mangerite rocks, dominated by MD magnetite, are commonly unstable and not all sites provide acceptable data. Mean directions for 66 sites spanning the entire intrusion are I = –73.5°, D = 303.4°, with α 95 = 3.7° and k = 24. The resulting pole position is at 35.9°S and 217.9°E, with a palaeolatitude for this part of Baltica of –59.4°. Examination of the magnetic mineralogy combined with geochronology for RIC rocks and cooling rates for the region yields an age of magnetization of 916 Ma. Metamorphic country rocks yield similar directions at least 10 km from the contact, confirming the presence of a contact aureole around the intrusion. Comparison to other early Neoproterozoic palaeomagnetic poles from southern Scandinavia confirms the high southern latitude position of Baltica at this time, and combined with the few other ~900 Ma poles from Baltica affirm the presence of the Rodinia supercontinent at this time.

Beschreibung: The Internal Sierras (IS) in the southern margin of the Western and Central Axial Zone (Southern Pyrenees) are affected by a syn-orogenic remagnetization that provides information to reconstruct deformation geometries at the time of acquisition of magnetization. Furthermore, the IS structure changes strike along its structural trend, from ~N120 to 130°E in the western and eastern margins to ~N070–090°E in the central part. Palaeomagnetic techniques have been used to (i) accurately define the timing of remagnetization with regard to deformation and (ii) determine if the along-strike trend variation in the IS was induced by deformation and thrust emplacement during the Pyrenean compression or, on the contrary, was the result of a primary orientation controlled by structures inherited from pre-orogenic times. From 23 new palaeomagnetic sites, collected in Upper Cretaceous marls and marly limestones, two meaningful and stable palaeomagnetic components were resolved, principally carried by magnetite: (1) a lower-temperature component (B) that unblocks between 200 °C and 325–400 °C and (2) a higher-temperature component (C) that has been successfully isolated by means of combined thermal (up to 400 °C) and AF demagnetization (generally up to 50–100 mT). The B component is a late remagnetization that post-dates folding and emplacement of basement thrust sheets in the IS (mainly the Gavarnie thrust). It supports small but statistically significant clockwise rotations in the western part of the IS (from +18 to +26°). These rotations can be attributed to the westwards shortening decrease in the thrust system below the Gavarnie unit that results from its along-strike structural change, with a higher number of basement thrusts to the east. The C component has been interpreted as an early remagnetization, based on the results of conglomerate and fold tests. This component predates basement thrusting and is diachronous across the study area: reverse and normal polarities dominate in the eastern and western margins of the IS, respectively. New and previous palaeomagnetic data point out that curvature in the IS is probably a primary feature and the along-strike change in their trend could be interpreted as the result of basement geometrical features inherited from Variscan, Late Variscan or Mesozoic times. A complex, multi-episodic remagnetization probably related to burial and deformation processes occurred during Eocene times.

Beschreibung: A comprehensive magnetic field model named CM5 has been derived from CHAMP, Ørsted and SAC-C satellite and observatory hourly-means data from 2000 August to 2013 January using the Swarm Level-2 Comprehensive Inversion (CI) algorithm. Swarm is a recently launched constellation of three satellites to map the Earth's magnetic field. The CI technique includes several interesting features such as the bias mitigation scheme known as Selective Infinite Variance Weighting (SIVW), a new treatment for attitude error in satellite vector measurements, and the inclusion of 3-D conductivity for ionospheric induction. SIVW has allowed for a much improved lithospheric field recovery over CM4 by exploiting CHAMP along-track difference data yielding resolution levels up to spherical harmonic degree 107, and has allowed for the successful extraction of the oceanic M 2 tidal magnetic field from quiet, nightside data. The 3-D induction now captures anomalous Solar-quiet features in coastal observatory daily records. CM5 provides a satisfactory, continuous description of the major magnetic fields in the near-Earth region over this time span, and its lithospheric, ionospheric and oceanic M 2 tidal constituents may be used as validation tools for future Swarm Level-2 products coming from the CI algorithm and other dedicated product algorithms.

Beschreibung: The observation of electromagnetic signals by ocean bottom electromagnetometers in association with the Tohoku tsunami of 2011 March 11 has raised the opportunity to re-examine the physics of motional induction due to oceanic long waves in the framework of 1-D magnetotellurics (MT). Although a propagating tsunami has a complex structure, the induced electromagnetic field can be simply approximated as a plane wave (though a simple thin-sheet approximation is not valid at higher frequencies). We found that the MT impedance due to a surface gravity wave (or the ‘motional impedance’) is influenced largely by the dispersion of the wave if the period is sufficiently short or the electrical conductivity of the seabed is low. The tipper due to the motional induction (or the ‘motional tipper’) and motional impedance are essentially identical if the underneath structure is 1-D. It would be possible to estimate the motional impedance and tipper from the observed ocean bottom electromagnetic field at the time of passing of a tsunami. The wave amplitude must be much greater than several tens of centimetres for the motional impedance and tipper estimation to be free from the effects of external sources. However, the obtained motional impedance and tipper will mostly represent the property of the wave and use of them may not be suitable to discuss the subseafloor conductivity structure.

Beschreibung: The lacustrine deposits in Heqing basin provide an excellent archive for long-term high-resolution palaeoclimate studies in the monsoon-dominated southeastern Tibetan Plateau region. In this study, we investigate the climatic significance of magnetic parameters for analysing the variability of the past climate. For this, we performed comprehensive time-series and statistical analysis on previously published proxy data from a 168-m-long drill core (Core-HQ) that spans ~900–30 ka, comprising mainly magnetic parameters and carbonate content (CC). Moreover, we investigated magnetic properties of modern soil in the catchment that predominantly formed on limestone bedrock. Key findings are: (1) modern soils and sediments of Core-HQ both contain a mixture of magnetite (Mt), maghemite (Mgh) and hematite (Ht), but magnetic concentration of the soils is one order higher; (2) a superparamagnetic (SP) fraction of Mt/Mgh dominates in the soils whereas in Core-HQ the SP contribution is generally very low; (3) a larger grain-size fraction of Mt/Mgh and Ht is also present in the soils. We explain variations of magnetic concentration and CC in Core-HQ by an increased wind transport of soil and a decreased surface water transport of carbonate and soil in less humid periods. Low-temperature oxidation of magnetite in the catchment is as another crucial process that reflects weathering conditions and is likely sensitive to humidity; the degree of LTO can be semi-quantified by the magnetic parameters ARM/SIRM and S -ratio. Combining CC, ARM/SIRM and S -ratio values, we derive a weathering intensity (WI). The WI index variation along Core-HQ shows strong fluctuations on a 100-kyr eccentricity scale in the lower part, especially during ~630–380 ka, followed by a long period (~320–80 ka) of persistently weaker weathering (drier?) conditions with low variability, and a rapid return of much stronger weathering (wetter?) conditions at ~80 ka. We suggest that a reduced influence of the Indian summer monsoon accounts for less moisture supply to the region and lower climatic variability.

Beschreibung: Magnetic mineralogy of crustal rocks has important implications for understanding continental crustal evolution and origin of regional magnetic anomalies. However, magnetic properties of the deep continental crust are still poorly understood. In this paper, measurements of density ( ), mass-specific magnetic susceptibility ( ), natural remanent magnetization (NRM) and magnetic hysteresis loops, temperature-dependent magnetic susceptibility ( – T ), chemical and mineral analyses were conducted on Archean gneiss samples from the Jixian petrophysical section in the Precambrian terrain, northeastern North China Craton, with the aim of refining understanding of magnetic phase transformations in the deep crustal rocks. Results show that density and rock magnetic properties change distinctly with metamorphic facies. The dominant magnetic mineral is magnetite, while a little hematite is present in a few samples. Together with geochemical and mineralogical compositions, it is inferred that progressive increase in metamorphic grade from east to west is the major cause for magnetic enhancement of the lower crust in the studied section. Therefore, we conclude that study of magnetic phases of deep crustal rocks can offer important insights into the history of high metamorphic grade terranes.

Beschreibung: The C -response is a conventional transfer function in global electromagnetic induction research and is traditionally determined from observations of magnetic variations in the vertical and horizontal components. Its interpretation relies on the assumption that the source of the variations is well approximated by a large-scale symmetric (magnetospheric) ring current, described by a single spherical harmonic. However, there is growing evidence for a more complex structure of this source. In this paper, we investigate the variability of C -responses due to sources different from the dominating large-scale symmetric ring current. We show that the effect is significant and persists at all periods. Describing the magnetospheric source by a single spherical harmonic coefficient thus injects substantial errors into the estimated responses. To overcome the problem, we introduce arrays of alternative transfer functions that relate the components of the magnetic variation to different spherical harmonic coefficients. These transfer functions can handle a complex spatial structure of the magnetospheric source. Compared to C -responses, we observe a significant increase in the coherencies relating input and output quantities of the new transfer functions, especially at high latitudes. This increases the usability of observatory magnetic data for the recovery of global 3-D mantle conductivity structure.

Beschreibung: Magnetic mineralogy of crustal rocks has important implications for understanding continental crustal evolution and origin of regional magnetic anomalies. However, magnetic properties of the deep continental crust are still poorly understood. In this paper, measurements of density ( ), mass-specific magnetic susceptibility ( ), natural remanent magnetization (NRM) and magnetic hysteresis loops, temperature-dependent magnetic susceptibility ( – T ), chemical and mineral analyses were conducted on Archean gneiss samples from the Jixian petrophysical section in the Precambrian terrain, northeastern North China Craton, with the aim of refining understanding of magnetic phase transformations in the deep crustal rocks. Results show that density and rock magnetic properties change distinctly with metamorphic facies. The dominant magnetic mineral is magnetite, while a little hematite is present in a few samples. Together with geochemical and mineralogical compositions, it is inferred that progressive increase in metamorphic grade from east to west is the major cause for magnetic enhancement of the lower crust in the studied section. Therefore, we conclude that study of magnetic phases of deep crustal rocks can offer important insights into the history of high metamorphic grade terranes.

Beschreibung: Surface nuclear magnetic resonance (surface-NMR) is a promising technique for exploring shallow subsurface aquifer structures. Surface-NMR can be applied in environments that are characterized as a 1-D layered Earth. The technique utilizes a single loop and is referred to as magnetic resonance sounding. The technique referred to as magnetic resonance tomography (MRT) allows complex 2-D aquifer structures to be explored. Currently, MRT requires multiple loops and a roll along measurement scheme, which causes long survey time. We propose a loop layout using an elongated transmitter and an in-loop receiver arrays (ETRA) to conduct a 2-D survey with just one measurement. We present a comprehensive comparison between the new layout and the common approaches based on sensitivity and resolution analyses and show synthetic and field data. The results show that ETRA generates subsurface images at sufficient resolution with significantly lower survey times than other loop layouts.

Beschreibung: We studied ionospheric responses to the 2012 April 11 M w 8.6 North Sumatra earthquake using total electron content (TEC) measurements with the regional Global Navigation Satellite System network. This earthquake ruptured the oceanic lithosphere off the Indian Ocean coast of North Sumatra, and is known as the largest strike-slip earthquake ever recorded. Coseismic ionospheric disturbances (CIDs) with rapid TEC enhancement of a few TEC units propagated northward with a speed of acoustic waves (~1 km s –1 ). Resonant atmospheric oscillation with a frequency ~4 mHz have been found as monochromatic oscillation of TEC lasting for an hour after the main shock and the largest aftershock. We compared CID amplitudes of 21 earthquakes world-wide with moment magnitudes ( M w ) 6.6–9.2. They roughly obeyed a law such that CID amplitude increases by two orders of magnitude for the M w increase of three. The 2012 North Sumatra earthquakes slightly deviated negatively from the trend possibly reflecting their strike-slip mechanisms, that is small vertical crustal movements for their magnitudes.

Beschreibung: Water theft carried out by manipulating water meters constrains volumetric pricing in semi-arid regions. Cooperative management can reduce theft and improve incentives for efficient water use by inducing peer monitoring. Using a theoretical model, we show that theft is more likely when prices are high, punishments are weak, and cooperatives are large. We also show how cooperative membership and punishment levels are determined endogenously by constraints on monitoring. We test the model on data from Tunisia for the years 2001–2003, relying on instruments that proxy for unobservable monitoring costs. The results confirm that well-designed incentives can reduce theft, and that constraints on monitoring costs affect institutional design.

Beschreibung: A large variety of subsidized crop insurance products are available to U.S. crop growers. Distinct and perhaps puzzling patterns in the choices of insurance products and coverage levels can be discerned. Where production conditions are better and yields are less risky then ( a ) higher insurance coverage levels are chosen; and ( b ) revenue insurance is preferred over yield insurance. Also, ( c ) the extent of preference for revenue insurance is stronger in more productive areas. Assuming, as many do, that growers seek to maximize subsidy transfers, point ( a ) can be explained by the interaction between yield technology and natural resource endowments. Points ( b ) and ( c ) can be explained by location in conjunction with the "natural hedge" and a contract design bias in how revenue insurance guarantees are computed. Empirical study of Risk Management Agency data on corn, soybean, and wheat yields, and insurance contract choices lend support to our model inferences.

Beschreibung: Perpetual conservation easements permanently remove the option to convert existing habitat to more intensive agricultural production. If existing habitat is at threat of conversion, removing the option to convert will reduce land values. In this article, we estimate the land value discount resulting from perpetual habitat conservation easements by using propensity score matching. We find that on the average eased parcel, land values fall by approximately $86 per acre for every acre of eased habitat. On average, our results suggest that landowners have been adequately compensated and conservation agencies have successfully secured habitat at risk of conversion.

Beschreibung: This paper investigates the spatial effects that the provision of environmental public goods have on residential location choices in a suburban context. Specifically, a spatial general equilibrium framework is developed to analyze the consequences of adopting an agri-environmental policy promoting the provision of positive farming externalities. We use a static monocentric model of an open city where agricultural bid-rents and agricultural amenities vary endogenously in space, and where the positive externalities associated with agricultural production are valued by households. Consistent with empirical evidence of the potential side effects that conservation policies may have in terms of urbanization patterns and land price changes, we show that under certain conditions implementing an agri-environmental policy may promote additional suburban development. Moreover, we demonstrate that the emergence of disconnected suburban areas may be significantly influenced by the location of land regulated by an agri-environmental policy. Finally, we discuss distributional aspects and show that while introducing an agri-environmental policy has a negative impact on most residential land value, it can have positive effects on farmland and residential land located within the regulated areas, suggesting the non-neutrality of such policies regarding the agents’ assets.

Beschreibung: We develop a model to evaluate the profitability of controlling rodent damage by placing barn owl nesting boxes in agricultural areas. The model incorporates the spatial patterns of barn owl predation pressure on rodents, and the impact of this predation pressure on nesting choices and agricultural output. We apply the model to data collected in Israel and find the installation of nesting boxes profitable. While this finding indicates that economic policy instruments to enhance the adoption of this biological control method are redundant, it does support stricter regulations on rodent control using rodenticides.

Beschreibung: Broad-band magnetic susceptibility (MS) measurement, a novel magnetic method capable of quantifying a narrow grain size distribution (GSD) of superparamagnetic (SP) particles by measuring low-field MS at a number of frequency steps spanning four orders of magnitude, has been tested in a loess/palaeosol section at Luochuan in the Chinese Loess Plateau. The studied succession consists of sequences from the latest palaeosol unit (S0) to the upper part of the loess unit (L2), spanning the last glacial–interglacial cycle. Reconstructed GSDs consist of volume fractions on the order of 10 –24 m 3 , and the mean GSDs are modal but with distinctive skewness among the loess, the weakly developed palaeosols (weak palaeosols), and the mature palaeosols. This indicates that the mean volume of SP particles in this loess/palaeosol sequence tends to increase during the transition from loess -〉 weak palaeosol -〉 palaeosol, an indication of grain growth as pedogenesis progresses. Total frequency dependence, or TFD(per cent), the difference between 130 at the lowest (130 Hz) and 500k at the highest (500 kHz) frequencies normalized to 130 , is judged to be a more suitable index than previous frequency dependence parameters for the concentration of SP particles. TFD(per cent) has a strong correlation with 130 , showing a continuous ‘growth curve’ with the rate of increase being highest for the loess, moderate for the weak palaeosols, and saturated for the palaeosols. The characteristic curve suggests that smaller SP particles are preferentially formed in the earlier stage of pedogenesis rather than the later phase when even larger particles are formed in mature palaeosols. These results demonstrate that the broad-band MS measurement method will be useful for the quantitative assessment of magnetic nanoparticles in soils and sediments.

Beschreibung: Environmental magnetism can provide simple and fast methods to semi-quantify the contamination levels by heavy metals in urban areas from the relationships between the magnetic properties and the heavy metal concentrations. The aim of this study is to explore the interpretation of selected magnetic parameters as proxies for the evaluation of heavy metal contamination in urban dusts. Dust samples were collected from different districts of Nanjing, Southeast China. The magnetic properties of dusts were analysed and heavy metal (Fe, Al, As, Ba, Cd, Cr, Co, Cu, Hg, Mn, Mo, Ni, Pb, Sr, V and Zn) concentrations were measured. Magnetic minerals, mainly from anthropogenic activities, are dominated by coarse stable single domain-grained ferrimagnetic minerals. The – T curves indicated that the main magnetic minerals are magnetite and hematite. Magnetic properties exhibited similar spatial distribution and significant positive correlations with the concentrations of most metals and with the pollution load index. Magnetic properties had a strong link with the abundance of heavy metals derived from industrial emission and traffic activities, whereas there was a poor correlation with the concentrations of metals from multi-sources such as commercial/domestic sources and natural processes. Besides the low-frequency magnetic susceptibility ( LF ) and the saturation isothermal remanent magnetization, other magnetic parameters can be considered as efficient indicators in the assessment of heavy metal contamination. The linkage between the magnetic properties and heavy metal concentrations in street dusts depends on the sources and nature of their magnetic fraction. Environmental magnetism analyses may provide simple and rapid methods for the assessment of heavy metal contamination in urban street dusts.

Beschreibung: We present in this paper an alternative approach to modelling and inversion of surface nuclear magnetic resonance (sNMR) data that has numerous attractive characteristics. By considering the forward and inverse problem in the frequency-domain (FD) it is possible to heavily truncate the band-limited data set with no appreciable loss of information. Furthermore, it introduces a natural method of signal demodulating and sidesteps the rotating frame transformation, which can be difficult to apply to large volumes. By instead remaining in the laboratory frame, phase shifts due to off-resonance transmission and static variation in the Larmor frequencies become simpler to account for. For these reasons, we present the first practical scheme for comprehensively inverting the complex sNMR data set. The use of complex data allows for improved depth and decay constant resolution. However, it requires the knowledge of subsurface electrical conductivity, and by corollary the ability to accurately model those effects in the sNMR record. Finally, the complex FD signal is sensitive to dephasing of the sNMR record due to static magnetic field inhomogeneity. We illustrate the ability to detect dephasing in synthetic and field data and demonstrate an effective inversion which is able to solve for static dephasing effects. This results in better estimation of meaningful decay constants as well as reduced error in extrapolated initial signal amplitude. Both of these parameters are critical for using sNMR in applications requiring estimation of hydraulic parameters.

Beschreibung: We report a palaeomagnetic investigation of the last full geomagnetic field reversal, the Matuyama-Brunhes (M-B) transition, as preserved in a continuous sequence of exposed lacustrine sediments in the Apennines of Central Italy. The palaeomagnetic record provides the most direct evidence for the tempo of transitional field behaviour yet obtained for the M-B transition. 40 Ar/ 39 Ar dating of tephra layers bracketing the M-B transition provides high-accuracy age constraints and indicates a mean sediment accumulation rate of about 0.2 mm yr –1 during the transition. Two relative palaeointensity (RPI) minima are present in the M-B transition. During the terminus of the upper RPI minimum, a directional change of about 180 ° occurred at an extremely fast rate, estimated to be less than 2 ° per year, with no intermediate virtual geomagnetic poles (VGPs) documented during the transit from the southern to northern hemisphere. Thus, the entry into the Brunhes Normal Chron as represented by the palaeomagnetic directions and VGPs developed in a time interval comparable to the duration of an average human life, which is an order of magnitude more rapid than suggested by current models. The reported investigation therefore provides high-resolution integrated palaeomagnetic and radioisotopic data that document the fine details of the anatomy and tempo of the M-B transition in Central Italy that in turn are crucial for a better understanding of Earth's magnetic field, and for the development of more sophisticated models that are able to describe its global structure and behaviour.

Beschreibung: The anisotropy of magnetic susceptibility (AMS) of soft rocks was measured in order to distinguish between the effect of remote and local strain fields, determine the size of the related inelastic damage zone and resolve the fault-plane solutions of past earthquakes. The AMS fabrics were explored next to late Pleistocene syndepositional normal faults (total displacement up to ~3.5 m) that cross soft lacustrine rocks within the seismically active Dead Sea basin. ‘Deposition fabrics’ prevail meters away from the fault planes and are characterized by scattered maximum and intermediate principal AMS axes. ‘Deformation fabrics’ are detected up to tens of centimetres from the fault planes and are characterized by well-grouped AMS axes, in which one of the principal axes is parallel to the strike of the nearby fault. Variations in the AMS fabrics and magnetic lineations define the size of the inelastic damage zone around the faults. The results demonstrate that the deformation-driven magnetic fabrics and the associated inelastic damage zones are compatible with coseismic dynamic faulting and the effects of the local strain field during earthquakes. Most of the AMS fabrics show a conspicuous similarity to that of the fault-plane solutions, i.e., the principal AMS axes and instantaneous strain ellipsoids are coaxial. These results suggest a novel application of the AMS method for defining the shape and size of the damage zones surrounding dynamic faults and determining the full tensor of the local strain field.

Beschreibung: Very low frequency (VLF) electromagnetic waves that penetrate conductive magma-filled dykes generate secondary fields on the surface that can be used to invert for dyke properties. The model used for the interpretation calculates currents induced in a conductive strip by an inducing field that decays exponentially with depth due to the conductivity of the surrounding medium. The differential equations are integrated to give an inhomogeneous Fredholm equation of the second kind with a kernel consisting of a modified Bessel function of the second kind. Numerical methods are typically used to solve for the induced currents in the strip. In this paper, we apply a modified Galerkin–Chebyshev method, which involves separating the kernel into source and field spectra and integrating the source terms to obtain a matrix equation for the unknown coefficients. The incident wave is expressed as a Chebyshev series. The modified Bessel function is separated into a logarithmic singularity and a non-singular remainder, both of which are expanded in complex Chebyshev polynomials. The Chebyshev coefficients for the remainder are evaluated using a fast Fourier transform, while the logarithmic term and incident field have analytic series. The deconvolution then involves a matrix inversion. The results depend on the ratio of strip-size to skin-depth. For infinite skin-depth and a singular conductivity distribution given by $\tau _0 a/\sqrt{a^2 - z^2 }$ (where 0 is the conductance, a is the half-length and z the distance from the centre), Parker gives an analytic solution. We present a similar analytic series solution for the finite skin-depth case, where the size to skin depth ratio is small. Results are presented for different ratios of size to skin depth that can be compared with numerical solutions. We compare full-space and half-space solutions. A fit of the model to VLF data taken above a magma filled dykes in Hawaii and Mt Etna demonstrates that while properties such as depth to top, conductivity ratio, tilt and dip can be determined, the depth to bottom is indeterminate due to the exponential decay of the inducing field.

Beschreibung: Correlated data errors (noise) are common in magnetotelluric (MT) data, but MT inversions typically neglect error correlations without investigating the impact of this simplification on inversion results. This paper examines effects of neglecting frequency- and spatially correlated noise on MT inversion, based on a nonlinear Bayesian formulation which quantifies the uncertainties of inversion results in terms of marginal posterior probability densities and credibility intervals. To do so, data with frequency- and spatially correlated noise of differing degrees are generated for several layered (1-D) synthetic cases. Bayesian MT inversions are carried out for these data sets with and without accounting for error correlation (i.e. applying full and diagonal covariance matrices, respectively, in the inversion), and the results are compared. For cases with noise that is strongly correlated over frequency or space, parameter uncertainties estimated using the diagonal-covariance simplification (neglecting error correlations) are found to often be significantly underestimated compared with results computed using the full covariance matrix.

Beschreibung: Efficient water management in agriculture is becoming critical due to increasing environmental constraints and global food and bio-energy demands. Farmers may respond to increased water scarcity along three main adjustment margins: a move towards rain-fed agriculture (super-extensive margin) or towards less water-intensive crops (extensive margin), and a reduction in water intensity for irrigated crops (intensive margin). Using a positive mathematical programming model of regional supply calibrated to economic and agronomic information, we decompose the total effect of reduced water availability on these adjustment margins in Beauce, a productive cereal region that relies on a groundwater resource to meet its irrigation needs. For realistic water scarcity scenarios, 57 per cent of the total response is attributable to super-extensive margin adjustments. The extensive margin represents 28 per cent of the total response, while the intensive margin accounts for 15 per cent. Crop-level analysis reveals more subtle adaptation patterns.

Beschreibung: Water quality regulations in the United States apply almost exclusively to point sources. In impaired watersheds where both point and nonpoint sources contribute to pollution, the U.S. Environmental Protection Agency (EPA) is encouraging the use of point-nonpoint trading to reduce the cost of point sources to meet their permit requirement, and to encourage nonpoint sources to voluntarily contribute more towards meeting overall water quality goals. The EPA guidance encourages trading programs to set a nonpoint source eligibility baseline that extracts some "extra" abatement from nonpoint sources. Research has shown that setting an eligibility baseline that is substantially more stringent than current management could discourage nonpoint source participation and significantly hinder trading. In this paper we examine how choosing the eligibility baseline for agricultural sources affects the efficiency goal of trading (reducing costs to point sources), as well as how it affects the EPA goal of encouraging nonpoint abatement. Using data from the Chesapeake Bay Watershed we find that eligibility baselines set to encourage additional nonpoint source abatement reduce the supply of credits in a market; the more stringent the baseline, the fewer the trades and the smaller the overall abatement from nonpoint sources. A subsidy to farmers for reducing the cost of meeting a baseline encourages greater nonpoint source abatement, but may not benefit the trading market.

Beschreibung: We report a new method to infer continuous time-series of the declination, inclination and intensity of the magnetic field from archaeomagnetic data. Adopting a Bayesian perspective, we need to specify a priori knowledge about the time evolution of the magnetic field. It consists in a time correlation function that we choose to be compatible with present knowledge about the geomagnetic time spectra. The results are presented as distributions of possible values for the declination, inclination or intensity. We find that the methodology can be adapted to account for the age uncertainties of archaeological artefacts and we use Markov chain Monte Carlo to explore the possible dates of observations. We apply the method to intensity data sets from Mari, Syria and to intensity and directional data sets from Paris, France. Our reconstructions display more rapid variations than previous studies and we find that the possible values of geomagnetic field elements are not necessarily normally distributed. Another output of the model is better age estimates of archaeological artefacts.

Beschreibung: We present a rigorous method for interpolation of electric and magnetic fields close to an interface with a conductivity contrast. The method takes into account not only a well-known discontinuity in the normal electric field, but also discontinuity in all the normal derivatives of electric and magnetic tangential fields. The proposed method is applied to marine 3-D controlled-source electromagnetic modelling (CSEM) where sources and receivers are located close to the seafloor separating conductive seawater and resistive formation. For the finite-difference scheme based on the Yee grid, the new interpolation is demonstrated to be much more accurate than alternative methods (interpolation using nodes on one side of the interface or interpolation using nodes on both sides, but ignoring the derivative jumps). The rigorous interpolation can handle arbitrary orientation of interface with respect to the grid, which is demonstrated on a marine CSEM example with a dipping seafloor. The interpolation coefficients are computed by minimizing a misfit between values at the nearest nodes and linear expansions of the continuous field components in the coordinate system aligned with the interface. The proposed interpolation operators can handle either uniform or non-uniform grids and can be applied to interpolation for both sources and receivers.

Beschreibung: As a bioinvasion spreads across a landscape from its point of introduction, damages rise roughly with the square of the distance from the original invasion. It is thus generally beneficial, at the landscape scale, to apply eradication or containment controls early if not immediately upon discovery. However, an individual property owner only has incentives to consider the costs and benefits of control on his/her own property rather than potential landscape-scale damages. Bioinvasions will therefore generally be under-controlled in a landscape of independent owners operating under a laissez-faire system. A mechanism is thus needed to induce early cooperative contributions to control costs from beneficiaries who would, without them, be invaded later. We develop a spatially-explicit, integrated model of invasion spread and human behavior to examine how different degrees of spatial cooperation affect patterns of invasion spread and the total costs and damages imposed. We compare individual laissez-faire, cooperative control by adjacent neighbors, and cooperative control by groups including more distant but nearby neighbors. As expected, private laissez-faire control decisions tend to under-control the invasion relative to socially optimal control under most circumstances. But a reasonably high fraction of first best payoffs can be achieved with only a modest geographical reach of cooperation. We also find that less extensive cooperation is needed to control invasions whose costs and damages otherwise lead to the largest externalities (circumstances with costs that are relatively low compared with damages). This suggests that even small amounts of cooperation to control bioinvasions can provide large social benefits.

Beschreibung: Farmland conservation policies typically use zoning and differentiated taxes to prevent urban development of farmland, but little is known about the effectiveness of these policies. This study adds to current knowledge by examining the impact of British Columbia's Agricultural Land Reserve (ALR), established in 1973, which severely restricts subdivision and nonagricultural uses for more than 4.7 million hectares of farmland. To determine the extent to which the ALR preserves farmland by reducing or removing the development option, a multilevel hedonic pricing model is used to estimate the impact of land use, geographic, and zoning characteristics on farmland value near the capital city of Victoria on Vancouver Island. Using sales data from 1974 through 2008, the model demonstrates a changing ALR impact over time that varies considerably by improved and unimproved land types. In 2008, landowners paid 19% less for the typical improved farmland parcel within the ALR versus that outside it. This suggests that would-be developers expect permanency in the zoning law, and prefer non-ALR zoned land. However, ALR land that is unimproved has a premium of 55%, suggesting that this land is more valuable for agriculture than for development. Farmland located closer to the city or the commuting highway commands a premium if it has a residence on it, with a residence also explaining why smaller agricultural properties sell at higher prices. However, it appears that zoning by itself is insufficient to protect farmland; other policies likely need to be implemented in conjunction with zoning to protect agricultural land.

Beschreibung: In many parts of the world, natural vegetation has been cleared to allow agricultural production. To ensure a long-term flow of ecosystem services without compromising agricultural activities, restoring the environment requires a balance between public and private benefits and costs. Information about private benefits generated by environmental assets can be utilized to identify conservation opportunities on private lands, evaluate environmental projects, and design effective policy instruments. We use a spatio-temporal hedonic model to estimate the private benefits of native vegetation on rural properties in the state of Victoria, Australia. Specifically, we estimate the marginal value of native vegetation on private land and examine how it varies with the extent of vegetation on a property and across a range of property types and sizes. Private benefits of native vegetation are greater per unit area on small and medium-sized properties and smaller on large production-oriented farms. Native vegetation exhibits diminishing marginal benefits as its proportion of a property increases. The current extent of native vegetation cover is lower than the extent that would maximize the amenity value to many landowners. There is scope for improved targeting of investment in the study region by incorporating private benefits of environmental projects into environmental planning processes. Landowners with high marginal private benefits from revegetation would be more willing to participate in a revegetation program. Targeting these landowners would likely provide higher value for money because such projects could be implemented at lower public cost.

Beschreibung: According to field observations, electromagnetic (EM) signals accompanying seismic waves can be recorded. The orders of magnitude of observed coseismic electric and magnetic signals associated with earthquakes are usually around 1 to 10 1 mV km –1 and 10 –2 to 1 nT, respectively. In this paper, we carry out numerical simulation of coseismic EM signals associated with seismic waves due to electrokinetic effect and compare with field observations. The seismic source is represented by a finite fault measuring 15 x 9 km 2 with a max slip displacement 1.5 m, corresponding to a M w 5.9 earthquake. While using the EM surface boundary condition of continuous horizontal EM components, the magnetic signals only accompany the late-arriving S waves at receiver near the ground surface. This is obviously different from field observations. Thus, we adopt another EM surface boundary condition which assumes the ground surface carries surface charge. For the used half-space model, a surface-charge density magnitude | Q sc | in excess of 10 –4 C m –2 is sufficient to make horizontal magnetic components clearly show up at the whole time duration of seismic waves. When | Q sc | increases, the contribution of surface-charge density to coseismic EM signals becomes more and more dominant comparing with that of the seismically induced streaming-current. We estimate the Q sc expected at the Earth's surface might be a value between –5 x 10 –4 and –0.1 C m –2 by the comparison between numerical results and field observations. The vertical magnetic signals only accompany the late-arriving seismic waves, because they are theoretically only induced by SH wave. The field observation results of vertical magnetic signals may be resulted from the scattering effect or the seismic dynamo effect. We conclude that electrokinetic effect combined with surface-charge assumption is one possible generation mechanism of the observed coseismic EM signals.

Beschreibung: Fifteen sites of lava flows from Martinique Island (FWI) have been selected to document the geomagnetic field in the Caribbean area over the past 2.5 Ma and further constrain the time-averaged field during this period. Identical characteristic directions were isolated using both AF and thermal stepwise demagnetization techniques in all flows. Nine mean-site directions have a normal polarity, while three others are reversed. The mean geomagnetic pole position obtained after reducing all directions to the same polarity is indistinguishable from the present north geographic pole. The dispersion is at least 8° larger than the values derived from the time-averaged field models and remains unexplained otherwise than resulting from the relatively small number of directions. The other three flows are characterized by large deviations from the expected north–south direction. One lava flow dated at 1.69 Ma (±0.02 Ma) is likely associated with a transitional field during the Gilsà subchron. The lava flow dated at 770 ka (±11 ka) coincides with the age of the Brunhes-Matuyama geomagnetic reversal and is also coeval with another intermediate flow of the same age found at Guadeloupe Island. The 617 ka (±52 ka) old unit is characterized by reversed directions that are evidently not related to the last reversal, but with other reversed polarity and transitional lava flows of the same age recorded, respectively at Mexico and La Palma island. We infer that the presence of reversed directions with the same age at distinct localities confirms that a short episode of reversed polarity has occurred during this period.

Beschreibung: Presently, most loess/palaeosol magnetic susceptibility records are interpreted as following either the wind-vigour model or the pedogenic enhancement model. However redoxomorphic processes induced by waterlogging, often referred to gleying in the loess literature, are also known to alter loess deposits but their impact on loess/palaeosol magnetic susceptibility records has received little attention. The reported rock magnetic study aims to characterize the mineral magnetic response of loess to waterlogging-induced redoxomorphic processes, thus improving our understanding of mineral magnetic changes within loess deposits with respect to environmental and climate conditions. The Nussloch loess-palaeosol deposit (Rhine Valley, Germany) was targeted because it is one of the best-studied Pleniglacial deposits for Western Europe in which numerous tundra gley intervals have been identified. Moreover, a comprehensive high-resolution environmental magnetism study has never been undertaken for this site. Various rock magnetism experiments were conducted at both room and low temperatures to characterise the composition, concentration and relative magnetic grain size of the mineral magnetic assemblage. The relative changes in magnetic parameters within the investigated loess interval are primarily controlled by (1) varying concentrations of coarse-grained ferrimagnetic particles of detrital (aeolian) origin and (2) dissolution of fine-grained ferrimagnetic particles related to in situ post-depositional alteration promoted by waterlogging-induced redoxomorphic processes. Goethite is found to be ubiquitous throughout the studied interval and is argued to have both a primary (aeolian) and secondary ( in situ ) origin. We conclude, that redoxomorphic processes induced by waterlogging, if present, will hinder the interpretation of magnetic susceptibility variations within loess and palaeosol deposits following the expected relationships dictated by the wind-vigour and the pedogenic enhancement magnetism models.

Beschreibung: Palaeosecular variation (PSV), as estimated from the scatter of remanent magnetization directions or poles, can be used to shed light to processes in the geodynamo, and potentially, to model the growth of the solid inner core. To understand the temporal aspects of the geomagnetic field behavior in the far past, we have calculated the scatter of palaeomagnetic poles (S) from a set of 55 high-quality observations derived from a new Precambrian paleomagnetic data compilation. Our analysis supports former Phanerozoic and Precambrian analyses of PSV, which favour a lower reversal rate, a higher stability of the geodynamo and a substantially smaller size of the inner core prior to 1.5 billion years ago.

Beschreibung: A palaeomagnetic study has been carried out in 29 cores drilled at six different sites from the volcanic products of Lomo Negro eruption (El Hierro, Canary Islands, Spain). Systematic thermal and alternating field demagnetization of the samples’ natural remanent magnetization revealed a northward, stable palaeomagnetic direction similar in all the samples. Rock magnetic experiments indicate that this palaeomagnetic component is carried by a mixture of high-Ti and low-Ti titanomagnetite crystals typical of basaltic lithologies that have experienced a significant degree of oxyexsolution during subaerial cooling. The well constrained palaeomagnetic direction of Lomo Negro lavas was used to perform a palaeomagnetic dating of the volcanic event, using the SHA.DIF.14k global geomagnetic model restricted for the last 3000 yr. It can be unambiguously concluded that Lomo Negro eruption occurred well before the previously proposed date of 1793 AD, with three different age ranges being statistically possible during the last 3 ka: 115 BC-7 AD, 410–626 AD and 1499–1602 AD. The calibration of a previously published non-calibrated 14 C dating suggests a XVI c. date for Lomo Negro eruption. This conclusion leaves open the possibility that the seismic crisis occurred at El Hierro in 1793 AD was related to an intrusive magmatic event that either did not reach the surface or either culminated in an unregistered submarine eruption similar to the one occurred in 2011–2012 at the southern off-shore ridge of the island.

Beschreibung: We have performed an extended palaeomagnetic study of the Oslo Graben volcanics, compared to the study of half a century ago by van Everdingen, using modern techniques and a four times larger amount of sites, plus additional rock magnetic experiments. We conclude that the average direction ( D = 204.0, I = –37.9, k = 46.9, α 95 = 2.0) and associated palaeomagnetic pole ( = 48.3, = 155.5, K = 52.2, A 95 = 1.9) of the Krokskogen and Vestfold volcanics together are statistically identical to those of the earlier study. This gives confidence in the fact that older palaeomagnetic studies can be reliable and robust, even though methods have improved. Our larger number of samples, and better age constraints, enable us to separate the data into two major intervals: the younger, on average, Krokskogen area and the older Vestfold area. The results show firstly that palaeolatitudes are slightly higher than predicted by the latest apparent polar wander path (APWP) for Eurasia by Torsvik et al. These data support an early Permian Pangaea A configuration and do not necessitate a Pangaea B configuration. The larger data set also allows us to assess the distribution of the characteristic remanent magnetization directions of the Oslo Graben in terms of geomagnetic field behaviour, which were acquired during a long period of dominantly single polarity the Permo-Carboniferous Reversed Superchron (PCRS). The distributions show a significantly lower virtual geomagnetic pole (VGP) scatter at the observed (low) latitudes than expected from a compilation from lavas of the last 5 Myr. The data do however show excellent agreement with the scatter observed both during the Cretaceous Normal Superchron and the PCRS. A comparison of the directional distributions in terms of elongation is less discriminating, since the large errors in all cases allow a fit to the predicted elongation/inclination behaviour of the TK03.GAD model.

Beschreibung: We illustrate a new imaging method to estimate the depth to the sources of potential fields and the structural index. The method consists of applying the Depth from EXtreme Point (DEXP) transformation to the ratio (R) between two different-order partial derivatives of the field. While the scaling function of the potential field depends on the structural index, we show that the scaling function of R merely depends on the difference between the two used orders of differentiation. This allows three main features to be established for the DEXP transformation of R: (1) it is independent from the structural index; (2) the estimation of the source depths is fully automatic, simply consisting in the search of position of the extreme points of the DEXP image and (3) the structural index of each source is finally determined from the scaling function or the extreme points using the estimated depth. Besides the well-known characteristics of the DEXP transformation, such as high-resolution and stability, the DEXP transformation of R enjoys one more relevant feature: it can be applied to multisource cases, yielding simultaneously correct estimations of structural index and depth for each source in the same image. However, while the DEXP transformation is a linear transformation of the field, the DEXP transformation of R is non-linear, and a procedure is described to circumvent the non-linear effects. The method is tested with synthetic examples and the estimated source parameters show a good agreement with the true values. The method was applied also to real magnetic data from the Pima copper mine, Arizona, USA, Hamrawien area, Egypt and Cataldere, Bala district of Turkey. The results are consistent with the known information about the causative sources.

Beschreibung: Variational data assimilation (4DVar) is a powerful technique for tuning dynamic models to observations, in order not only to forecast future time evolution of the system, but to make inferences about quantities that are otherwise unconstrained by observation. We apply this technique, well-grounded in meteorology and oceanography, to the Earth's core where incompressible fluid motions in an electrically conducting medium are responsible for magnetic field generation. Our dynamic model's momentum equation neglects inertia such that the entire evolution depends only on the structure of the initial magnetic field; time evolution of the system is solely governed by the equation of magnetic induction. Nevertheless the dynamic system encompasses the effects of rotation, Lorentz forces and viscosity and aims to mimic a reasonable force-balance in the Earth's core. Building on the work of Li et al. , in order to optimize the data-fit subject to the dynamics, we further develop the mathematical structure of the adjoint equations of the system. We address the feasibility of recovering 3-D spatial properties of the system using only time-varying 2-D observations of different character. Using closed-loop testing, we demonstrate the retrievability of the initial state (and thus the entire trajectory) of the system over convective timescales, when sampling in regions in which magnetic induction dominates over diffusion. The results suggest the possibility of retrieving the entire trajectory of the dynamo system of the Earth using the 400-yr model of secular variation gufm1 .

Beschreibung: We investigate polarity reversals in the geodynamo using a rotating, convection-driven dynamo model. As the flow in rapidly rotating convection is dominated by columns aligned with the axis of rotation, the focus is on the dynamics of columnar vortices. By studying the growth of a seed magnetic field to a stable axial dipole field, we show that the magnetic field acts in ways that significantly enhance the relative helicity between cyclonic and anticyclonic vortices. This flow asymmetry is the hallmark of a dipolar dynamo. Strong buoyancy, on the other hand, offsets the effect of the magnetic field, establishing parity between positive and negative vortices. As the dipole field is deprived of the helicity required to support itself, the dynamo is pushed into a reversing state. This is a likely regime for polarity reversals in the Earth's core. The integral lengthscale at which buoyancy injects energy is not significantly different from the convective flow lengthscale, which implies that buoyancy does not feed vortices at the small scales where non-linear inertia is present. The lengthscale at which the Lorentz force acts in the reversing dynamo is small, which may allow the passive presence of non-linear inertia in the small scales.

Beschreibung: A new methodology to estimate magnetotelluric (MT) tensor relationships, called Estimation of Local transfer-functIons by Combining Interstation Transfer-functions (ELICIT), is proposed whereby the MT tensor relationships of the local site are derived using only interstation transfer functions. The MT impedance tensor and the geomagnetic transfer function at the local site are characterised by combining interstation tensor relationships between electric and magnetic fields at the local site with the horizontal magnetic fields acquired at a neighbouring site. The main property of the proposed method is that the employed interstation transfer functions are independently constrained, without the need to acquire the electric and the magnetic fields at the local site simultaneously to recover the local MT tensor relationships. Due to this property, the ELICIT method offers new possibilities for MT data acquisition and processing, providing significant improvements when the magnetic time-series at the local site are affected by local noise or are truncated. Error analysis shows that, even when magnetic fields are truncated, the quality of the results obtained following the ELICIT method are similar to those we would obtain if the magnetic fields had not been truncated. Another important property is that different neighbouring sites can be used to recover the tensor relationships at the local site. Averaging of results obtained using different neighbouring sites can be performed to improve the statistics. For our example data, when the ELICIT method is used to improve the statistics, errors of the estimates for periods between 1000 and 20 000 s periodicities are clearly reduced. All interstation transfer functions are calculated doing remote reference and the bootstrap method is used to compute the errors, when necessary. Long period magnetotelluric data acquired in the Pyrenees and in the Atlas Mountains in Morocco, and magnetic data provided by Fürstenfeldbruck magnetic observatory have been used to test the proposed ELICIT method, with positive results. Due to the lack of requirement that the electric and the magnetic fields of the local site be acquired simultaneously, the proposed method also offers new possibilities for MT data acquisition, optimizing the available instrumentation.

Beschreibung: The robust statistical model of a Gaussian core contaminated by outlying data that underlies robust estimation of the magnetotelluric (MT) response function has been re-examined. The residuals from robust estimators are systematically long tailed compared to a distribution based on the Gaussian, and hence are inconsistent with the robust model. Instead, MT data are pervasively described by the alpha stable distribution family whose variance and sometimes mean are undefined. A maximum likelihood estimator (MLE) that exploits the stable nature of MT data is formulated, and its two-stage implementation in which stable parameters are first fit to the data and then the MT responses are solved for is described. The MLE is shown to be inherently robust, but differs from the conventional robust estimator because it is based on a model derived from the data, while robust estimators are ad hoc , being based on the robust model that is inconsistent with actual data. Propriety versus impropriety of the complex MT response was investigated, and a likelihood ratio test for propriety and its null distribution was established. The Cramér-Rao lower bounds for the covariance matrix of proper and improper MT responses were specified. The MLE was applied to exemplar long period and broad-band data sets from South Africa. Both are shown to be significantly stably distributed using the Kolmogorov–Smirnov goodness of fit and Ansari-Bradley non-parametric dispersion tests. Impropriety of the MT responses at both sites is pervasive, hence the improper Cramér-Rao bound was used to estimate the MLE covariance. The MLE is shown to be nearly unbiased and well described by a Gaussian distribution based on bootstrap simulation. The MLE was compared to a conventional robust estimator, establishing that the standard errors of the former are systematically smaller than for the latter and that the standardized differences between them exhibit excursions that are both too frequent and too large to be described by a Gaussian model. This is ascribed to pervasive bias of the robust estimator that is to some degree obscured by their systematically large confidence bounds. Finally, a series of topics for further investigation is proposed.

Beschreibung: Broad-band and long period magnetotelluric (MT) data were acquired at 39 stations along five NNW-SSE profiles crossing the Iapetus Suture Zone (ISZ) in Ireland. Regional strike analyses indicate that the vast majority of the MT data is consistent with an assumption of a 2-D geo-electric strike direction. Strike is N52°E for the three easternmost profiles and N75°E for the two westernmost profiles; these directions correlate well with the observed predominant geological strike of the study region. 2-D inversions of the galvanic distortion-corrected TE and TM mode data from each profile are shown and discussed. As mapped geological variations between the neighbouring profiles suggest a heterogeneous subsurface, it is important to verify the robustness of the presence and geometries of prominent conductivity anomalies by employing 3-D forward and inverse modelling. A high conductivity layer (resistivity of 1–10 m), found at middle to lower crustal depths and presumed to be indicative of metamorphosed graphitic sediments rich in sulphides deposited during the convergence of the Laurentian and Avalonian continents, essentially constitutes the electrical signature of the ISZ. Shallow conductors observed are probably due to black shales that were widely deposited within the sedimentary accretionary wedge during Ordovician time. We interpret the moderately low resistivity at shallow depths from west to east across Ireland as indicative of an increase in maturity of the black shales in the easterly direction. From our conductivity models the southern extent of the ISZ is inferred to lie between the Navan Silvermines Fault and the Navan Tipperary Line, and shows clear resistivity contrast along all the profiles at the southern MT stations. The change in resistivity deduced from the 2-D models is spatially related to the composition of Lower Palaeozoic Ordovician, Silurian, Devonian and Carboniferous rocks. At upper mantle depths of about 60 km, a high conductivity block below the central MT stations is found to lie within the accretionary wedge of the Iapetus suture, and the location of the conductive anomaly corroborates reasonably well with the inferred spreading head of the putative Iceland plume-related magmatic intrusion. The low resistivity upper crust beneath the ISZ is indeed rich in Ordovician rocks with black shale content in the eastern as well as the central part; the western part is largely underlain by a highly resistive block of volcanic and metamorphosed rocks forming crystalline basement.

Beschreibung: This paper presents new directional archeomagnetic data from nine Meso-/Neolithic fireplaces, sampled in a cave shelter, at Arconciel, in western Switzerland. Rock magnetic measurements indicate a homogenous magnetic mineralogy in all fireplaces, with magnetite as the main magnetic carrier. The remanent magnetization is stable and generally shows one characteristic directional component. Nine new directions, which were obtained from Arconciel, are combined with 356 other archeomagnetic data from a circular area with a radius of 700 km around this site, to obtain a penalized least square spline fit for the past 9000 yr. We found in general good agreement with other local compilations, such as the Balkan curve, the regional SCHA.DIF.8k model and with lake sediments from UK, Fennoscandia and Switzerland. Nevertheless, a time lag of several centuries is observed for a declination maximum between the archeomagnetic spline fit and the other European data records around 5900 BC. This time lag is also observed in the Swiss lake sediment record; therefore we interpret this shift as a local feature of the Earth's magnetic field.

Beschreibung: In Christian religion the sunrise is of great symbolic importance. Therefore, many churches constructed in the Middle Ages point towards geographic East. Although ‘easting’ of churches actually refers to the alignment towards the azimuth of sunrise on the individual churches’ patron's day, deviation of nave alignment from the geographic East direction is often assumed to be caused by the use of magnetic compasses. Therefore, the church alignment could provide information about historical magnetic declination. We investigate 124 churches in Lower Austria and 68 in northern Germany to clarify this question as well as the ‘easting’ hypothesis. Church orientations are determined from georeferenced satellite images. Metadata such as the construction year, possible reconstructions and the church patron are gathered to determine the date when current church direction was appointed, and to perform sunrise calculations. However, due to uncertainties of construction years and the declining importance of orientation tradition after the 15th century several churches are excluded from the study. Thus, 32 churches with reliable metadata remain for evaluation in each region. The analysis reveals a preferred alignment of naves towards geographic East in Lower Austria and northern Germany. The construction and alignment of churches was often affected by the pre-existing buildings and streets or topography and natural surroundings. Therefore, deviations from geographic East are more likely caused by town or landscape. The mean deviations from magnetic East in both regions are large compared to the mean deviations from geographic East and sunrise azimuths. Hence, the use of compasses cannot be confirmed. Despite a few churches indicating orientation according to their patron's day, a general trend cannot be observed in the data reviewed.

Beschreibung: A new finite difference solution to the global geomagnetic induction problem is developed and tested, based on a modified Lorenz gauge of the magnetic vector and electric scalar potentials and implementing a novel, overset ‘Yin-Yang’ grid that avoids unnecessary mesh refinement at the geographic poles. Previously used in whole-earth mantle convection models, the overset grid is built from a pair of partially overlapping mid-latitude latitude–longitude (lat/lon) grids, one of which is rotated with respect to the other for complete coverage of the sphere. Because of this symmetry, only one set of finite difference templates is required for global discretization of the governing Maxwell equations, a redundancy that is exploited for computational efficiency and multithreaded parallelization. Comparisons between solutions obtained by the proposed method show excellent agreement with those obtained by independent integral equation methods for 1-D, 2-D and 3-D problem geometries. The computational footprint of the method is minimized through a (non-symmetric) matrix-free BiCG-STAB iterative solver which computes finite difference matrix coefficients ‘on the fly’ as needed, rather than pulling stored values from memory. Scaling of the matrix-free BiCG-STAB algorithm with problem size shows behaviour similar to that seen with the (symmetric) QMR algorithm used in the Cartesian case from which the present algorithm is based. The proposed method may therefore provide a competitive addition to the existing body of global-scale geomagnetic induction modelling algorithms, allowing for resource-efficient forward modelling as the kernel for large-scale computing such as inversion of geomagnetic response functions, computational hypothesis testing and parametric studies of mantle geodynamics and physiochemical state.

Beschreibung: We present a robust and scalable solver for time-harmonic Maxwell's equations for problems with large conductivity contrasts, wide range of frequencies, stretched grids and locally refined meshes. The solver is part of the fully distributed adaptive 3-D electromagnetic modelling scheme which employs the finite element method and unstructured non-conforming hexahedral meshes for spatial discretization using the open-source software deal.II. We use the complex-valued electric field formulation and split it into two real-valued equations for which we utilize an optimal block-diagonal pre-conditioner. Application of this pre-conditioner requires the solution of two smaller real-valued symmetric problems. We solve them by using either a direct solver or the conjugate gradient method pre-conditioned with the recently introduced auxiliary space technique. The auxiliary space pre-conditioner reformulates the original problem in form of several simpler ones, which are then solved using highly efficient algebraic multigrid methods. In this paper, we consider the magnetotelluric case and verify our numerical scheme by using COMMEMI 3-D models. Afterwards, we run a series of numerical experiments and demonstrate that the solver converges in a small number of iterations for a wide frequency range and variable problem sizes. The number of iterations is independent of the problem size, but exhibits a mild dependency on frequency. To test the stability of the method on locally refined meshes, we have implemented a residual-based a posteriori error estimator and compared it with uniform mesh refinement for problems up to 200 million unknowns. We test the scalability of the most time consuming parts of our code and show that they fulfill the strong scaling assumption as long as each MPI process possesses enough degrees of freedom to alleviate communication overburden. Finally, we refer back to a direct solver-based pre-conditioner and analyse its complexity in time. The results show that for multiple right-hand sides the direct solver-based pre-conditioner can still be faster for problems of medium size. On the other hand, it also shows non-linear growth in memory, whereas the auxiliary space method increases only linearly.

Beschreibung: The configuration of the Pangea supercontinent has been a topic of intense debate for almost half a century, a controversy that stems from discrepancies between the geology-based Pangea-A and the palaeomagnetically based Pangea-B. Recent palaeomagnetic compilations aimed at resolving this controversy have identified the poor quality of palaeomagnetic data from Gondwana for Permian times as a major obstacle. Specifically, the vast majority of Gondwanan poles come from sedimentary rocks that are prone to biases from compaction or are poorly dated. Here, we present a new palaeomagnetic pole for cratonic South America based on impact melts from the 254.7 ± 2.5 Ma Araguainha impact structure (AIS). The granite basement, the impact-generated melt sheet and veins were sampled at 28 sites (169 specimens) and provided a reliable palaeomagnetic record similar to that of volcanic rocks. Alternating field and thermal demagnetization indicate a stable characteristic remanent magnetization carried by both magnetite and haematite. All sites but one show a single palaeomagnetic direction of normal polarity with a mean direction of Dec = 357.4°; Inc = –38.9°; N  = 28; k  = 62.35; α 95  = 3.5°, yielding a palaeomagnetic pole (AIS) at Lat = –84.2; Lon = 326.6; K  = 83.5; A 95  = 3.6°; S B  = 9.6°. The new pole provides a firm constraint on the position of Gondwana which is consistent with the Pangea A configuration.

Beschreibung: Reconstructions of the Holocene geomagnetic field and how it varies on millennial timescales are important for understanding processes in the core but may also be used to study long-term solar-terrestrial relationships and as relative dating tools for geological and archaeological archives. Here, we present a new family of spherical harmonic geomagnetic field models spanning the past 9000 yr based on magnetic field directions and intensity stored in archaeological artefacts, igneous rocks and sediment records. A new modelling strategy introduces alternative data treatments with a focus on extracting more information from sedimentary data. To reduce the influence of a few individual records all sedimentary data are resampled in 50-yr bins, which also means that more weight is given to archaeomagnetic data during the inversion. The sedimentary declination data are treated as relative values and adjusted iteratively based on prior information. Finally, an alternative way of treating the sediment data chronologies has enabled us to both assess the likely range of age uncertainties, often up to and possibly exceeding 500 yr and adjust the timescale of each record based on comparisons with predictions from a preliminary model. As a result of the data adjustments, power has been shifted from quadrupole and octupole to higher degrees compared with previous Holocene geomagnetic field models. We find evidence for dominantly westward drift of northern high latitude high intensity flux patches at the core mantle boundary for the last 4000 yr. The new models also show intermittent occurrence of reversed flux at the edge of or inside the inner core tangent cylinder, possibly originating from the equator.

Beschreibung: Several recent studies have used palaeomagnetic estimates of the virtual axial dipole moment to construct a quantitative stochastic model for fluctuations and reversals in the Earth's dipole field. We investigate the physical significance of the terms in a standard stochastic (Langevin) model using output from a numerical geodynamo model. The first term, known as the drift term, characterizes the slow adjustment of the dipole field toward a time-averaged state. We find that the timescale for this slow adjustment is set by the magnetic decay time of dipole fluctuations. These fluctuations are typically be represented by the first few decay modes. The second term is often called the noise term because it characterizes the influence of short-period convective fluctuations in the core. We establish a connection between the noise term and the rms variation in magnetic induction. Applying these results to the palaeomagnetic field suggests that the rms variation in dipole generation exceeds the mean rate of generation. Such large fluctuations may be necessary to permit magnetic reversals. Palaeomagnetic estimates of the drift term favour a high electrical conductivity in the core. A lower bound on electrical conductivity is 0.6 x 10 6 S m –1 . Similarly, we establish an upper bound on turbulent magnetic diffusivity (0.8 m 2 s –1 ), although realistic estimates may be much less.

Beschreibung: We demonstrate that first differences of polar orbiting satellite magnetic data in the along-track direction can be used to obtain high resolution models of the lithospheric field. Along-track differences approximate the north–south magnetic field gradients for non-polar latitudes. In a test case, using 2 yr of low altitude data from the CHAMP satellite, we show that use of along-track differences of vector field data results in an enhanced recovery of the small scale lithospheric field, compared to the use of the vector field data themselves. We show that the along-track technique performs especially well in the estimation of near zonal spherical harmonic coefficients. Moreover, lithospheric field models determined using along-track differences are found to be less sensitive to the presence of unmodelled external field contributions and problems associated with the polar gap are ameliorated. Experiments in modelling the Earth's lithospheric magnetic field with along-track differences are presented here as a proof of concept. We anticipate that use of such along-track differences in combination with east–west field differences, as are now provided by the Swarm satellite constellation, will be important in building the next generation of lithospheric field models.

Beschreibung: We develop an singular value decomposition-based compression of the Green's function matrix of an electromagnetic integral equation forward solver for global geomagnetic induction, on top of an fast Fourier transform reduction of the system to a block-diagonal form. With this approach, the memory usage and CPU time of Krylov subspace iterative solutions are significantly reduced at a very small cost of accuracy, making the accelerated forward solver well suited for 3-D inversions as well as forward simulations with multiple sources.

Beschreibung: The Fisher distribution is central to palaeomagnetism but presents several problems when used to characterize geomagnetic field directions as observed in sequences of volcanic rocks. First, it introduces a shallowing effect when used to define the mean of any group of directional unit vectors. This is problematic because it can suggest the presence of persistent non-axial dipole components when none are present. More importantly, it fails to capture the observed ‘long tail’ in distributions of both directions and associated virtual geomagnetic poles in terms of angular distance from a central direction. To achieve a good fit to data, it therefore requires the introduction of a second distribution (and therefore the estimation of additional parameters) or the arbitrary removal of data. Here we present a new distribution to describe palaeomagnetic directions and demonstrate that it overcomes both of these problems, generating robust indicators of both the central direction (or pole position) and the spread of palaeomagnetic data as defined by unit vectors. Starting from the assumption that poles (or directions) have an expected colatitude, rather than a mean location, we derive the spherical exponential distribution. We demonstrate that this new distribution provides a good fit to palaeomagnetic data sets from seven large igneous provinces between 15 and 65 Ma and also those produced by numerical dynamo models. We also use it to derive a new shape parameter which may be used as a diagnostic tool for testing goodness of fit of models to data and use this to argue for a shift in geomagnetic behaviour between 5 and 15 Ma. Furthermore, we point out that this new statistic can be used to determine the most appropriate distribution to be used when constructing confidence limits for poles.

Beschreibung: Haematite pigment is a common constituent of sedimentary rocks, but its contribution to the natural remanent magnetization of rocks is poorly understood. Here, we describe magnetic properties of two distinct pigment types that produce a characteristic decorative ‘print stone’ found in the ~2.5 Ga Mount McRae Shale Formation, Hamersley Province, Western Australia. Distinct magnetic remanence directions observed in the Print Stone can be correlated to each pigment type. By comparison with the Australian apparent polar wander path, the remanence carried by uniformly distributed pigment can be dated to ~15–25 Ma, while two age options, the Mesoproterozoic (~1.5 Ga) or the middle Carboniferous (~320–310 Ma), are permissible for the remanence carried by the pigment responsible for the distinctive ‘newsprint’ pattern. Magnetic properties and demagnetization characteristics of the different pigment types overlap significantly, and thus are not predictive of the dominant remanence carrier. Magnetic characteristics of the uniformly distributed pigment vary significantly on short spatial scales. Strong local control on pigment formation raises the possibility that a primary remanent magnetization may survive locally in pockets within sedimentary red bed formations.

Beschreibung: Remanent magnetizations of magnetites between single-domain (SD) threshold size (0.1 μm) and 20 μm have SD-like intensities and coercivities. This paper shows for the first time that magnetite's induced magnetization also has pseudo-single-domain behaviour. The first part of the paper reports temperature-dependent initial susceptibility data, k 0 ( T ), of sized magnetites and assesses their granulometric potential. The second part transforms coercive force data, H c ( T ), for the same magnetites into simulated k 0 ( T ) curves. The third part considers k 0 ( T ) results of coarse-grained mafic rocks as candidate sources of deep-seated magnetic anomalies. High-temperature susceptibility k 0 measured with a Kappabridge for eight fractions of crushed natural magnetites (median sizes of 0.6, 1, 3, 6, 9, 14, 110 and 135 μm) shows a progressive increase in the height of the Hopkinson peak below the Curie point as grain size decreases. The trend is systematic and has granulometric potential in the 1–14 μm range. Self-demagnetization should produce almost flat k 0 ( T ) in grains larger than SD size but experimentally, well-defined Hopkinson peaks are not limited to the finest grains. 1-μm magnetites have a peak 1.5 times k 0 at 20 °C and 14-μm grains have a peak of 1.25. Only 110 and 135 μm grains have T -independent k 0 . Using an empirical relationship between coercive force H c and k 0 , H c ( T ) data for the sized magnetites were used to simulate k 0 ( T ) results. A hump in the k 0 heating curve around 250 °C was traced to annealing out internal strains, evident in H c data measured in first heatings. For sizes ≤6 μm, observed Hopkinson peaks were smaller than predicted, possibly because of a previously unrecognized grain-size dependence of the empirical constant relating H c and k 0 . Two crystalline rocks, a gabbro and a diabase, combine SD-like Hopkinson peaks and multidomain (MD) flat ramps in their k 0 ( T ) data. In the diabase, a Hopkinson peak is prominent in separated plagioclase grains containing submicron magnetite, but is masked in whole-rock data. The gabbro has a clear superposition of SD and MD k 0 ( T ) functions in its whole-rock data, with a Hopkinson peak of 1.35. If oceanic layer-3 gabbros have similar susceptibility enhancement above 500 °C, they could be more important magnetic anomaly sources than room-temperature k 0 measurements on dredged or fault-uplifted samples would suggest.

Beschreibung: We examine the effects of energy prices on groundwater extraction using an econometric model of a farmer's irrigation water pumping decision that accounts for both the intensive and extensive margins. Our results show that energy prices have an effect on both types of margins. Increasing energy prices would affect crop selection decisions, crop acreage allocation decisions, and farmers’ demand for water. Our estimated total marginal effect, which sums the effects on the intensive and extensive margins, suggests that a $1 per million btu increase in the energy price would decrease water extraction by an individual farmer by 5.89 acre-feet per year, a decrease of 3.6 percent of the average annual extraction rate. Our estimated elasticity of water extraction with respect to energy price is –0.26.

Beschreibung: Recent increases in farm real estate values in the United States have increased farm equity. By exploiting periods of high and low appreciation that caused various increases in wealth for farmers owning various shares of their farmland, we examine whether U.S. grain farmers expanded their acres harvested or acres owned in response to an increase in their land wealth. We find that land wealth had little effect on farm size. However, for similarly-sized farms, a larger ownership share (10 percentage points) led to an increase in the growth of land owned (2 percentage points). Because older farmers own more of the land that they farm, greater land appreciation slows the rate at which younger farmers acquire land relative to older farmers.

Beschreibung: We investigate the influence of different thermal and velocity boundary conditions on numerical geodynamo models. We concentrate on the implications for magnetic field morphology, heat transport scaling laws, force balances and generation mechanisms. The field morphology most strongly depends on the local Rossby number, but there is some variation in the dipolarity of the field with boundary condition. Scaling laws also depend on the boundary conditions, but a diffusivity-free scaling is a good first order approximation for all our dipolar models. Our multipolar models, however, obey different scaling laws from dipolar models implying a different force balance in these models. We find that our dipolar models have a stronger degree of Lorentz–Coriolis balance compared to our multipolar models which have a stronger degree of Lorentz-inertial balance.The models with a stronger Lorentz–Coriolis dominance can be generated by either α, α 2 or α 2 mechanisms whereas the models with a stronger Lorentz-inertial balance are all α 2 dynamos. These results imply that some caution is necessary when extrapolating results from dynamo models to Earth-like parameters since the choice of boundary conditions can have important effects.

Beschreibung: Under economic structural changes, some households intensify farming while others reduce it, resulting in significant changes in landholding. This paper studies the link between such changes in landholding and household well-being during a period of rapid economic transformation in Vietnam. Using a rural household panel data set, we find a U-shaped relationship between landholding and well-being: both accumulating crop land and moving out of farming are associated with higher household income and expenditure. Notably, these relationships are greater in communes that are less developed, suggesting that the benefits of structural transformation may diminish at higher levels of development.

Beschreibung: Lake sediments are excellent sources of palaeoenvironmental and palaeoclimatic information because they provide continuous and high-resolution records. South America is of particular interest because it is the only landmass that stretches southward into the Pacific and Atlantic Oceans towards Antarctica. The aim of this study is to explore the relationship of magnetic parameters with elemental and palaeobiological data of Laguna Potrok Aike to develop a model of lake-level changes and related hydrological and climatic fluctuations. Magnetic measurements were performed on subsamples from 15 500 cal. BP to the present, and associated rock magnetic parameters were calculated to infer magnetic mineralogy, concentration and grain size. According to the model, parameters dependent on magnetic concentration and grain size are directly related to lake-level changes. During dry periods, the remanent coercivity displays high values, whereas the proportion of magnetite is relatively low. Low percentages of greigite are observed, indicating that the water of the lake was stratified at least four times during the studied period, at approximately 10 300, 8900, 8500 and 8300 cal. BP. The preservation of greigite by inhibiting its complete transformation into pyrite is associated with a rapid burial that occurs with high sedimentation rates. Thermal stratification could be caused by a slight cooling in the area triggered by a weakening of the Westerlies and/or low activity of the sun, sum to the effect of a flood of melt water in the North Atlantic.

Beschreibung: Measurement of magnetic vector or tensor quantities, namely of field or field gradient, delivers more details of the underlying geological setting in geomagnetic prospection than a scalar measurement of a single component or of the scalar total magnetic intensity. Currently, highest measurement resolutions are achievable with superconducting quantum interference device (SQUID)-based systems. Due to technological limitations, it is necessary to suppress the parasitic magnetic field response from the SQUID gradiometer signals, which are a superposition of one tensor component and all three orthogonal magnetic field components. This in turn requires an accurate estimation of the local magnetic field. Such a measurement can itself be achieved via three additional orthogonal SQUID reference magnetometers. It is the calibration of such a SQUID reference vector magnetometer system that is the subject of this paper. A number of vector magnetometer calibration methods are described in the literature. We present two methods that we have implemented and compared, for their suitability of rapid data processing and integration into a full tensor magnetic gradiometry, SQUID-based, system. We conclude that the calibration routines must necessarily model fabrication misalignments, field offset and scale factors, and include comparison with a reference magnetic field. In order to enable fast processing on site, the software must be able to function as a stand-alone toolbox.

Beschreibung: Motional induction in the ocean by tides has long been observed by both land and satellite measurements of magnetic fields. While these signals are weak (~10 nT) when compared to the main magnetic field, their persistent nature makes them important for consideration during geomagnetic field modelling. Previous studies have reported several discrepancies between observations and numerical predictions of the tidal magnetic signals and those studies were inconclusive of the source of the error. We address this issue by (1) analysing magnetometer data from ocean-bottom stations, where the low-noise and high-signal environment is most suitable for detecting the weak tidal magnetic signals, (2) by numerically predicting the magnetic field with a spatial resolution that is 16 times higher than the previous studies and (3) by using four different models of upper-mantle conductivity. We use vector magnetic data from six ocean-bottom electromagnetic (OBEM) stations located in the Northwestern Pacific Ocean. The OBEM tidal amplitudes were derived using an iteratively re-weighted least-squares (IRLS) method and by limiting the analysis of lunar semidiurnal (M2), lunar elliptic semidinurnal (N2) and diurnal (O1) tidal modes to the night-time. Using a 3-D electromagnetic induction solver and the TPX07.2 tidal model, we predict the tidal magnetic signal. We use earth models with non-uniform oceans and four 1-D mantle sections underneath taken from Kuvshinov and Olsen, Shimizu et al. and Baba et al. to compare the effect of upper-mantle conductivity. We find that in general, the predictions and observations match within 10–70 per cent across all the stations for each of the tidal modes. The median normalized percent difference (NPD) between observed and predicted amplitudes for the tidal modes M2, N2 and O1 were 15 per cent, 47 per cent and 98 per cent, respectively, for all the stations and models. At the majority of stations, and for each of the tidal modes, the higher resolution (0.25° 0.25°) modelling gave amplitudes consistently closer to the observations than the lower resolution (1° 1°) modelling. The difference in lithospheric resistance east and west of the Izu–Bonin trench system seems to be affecting the model response and observations in the O1 tidal mode. This response is not seen in the M2 and N2 modes, thereby indicating that the O1 mode is more sensitive to lithospheric resistance.

Beschreibung: Magnetotelluric observations over one or 2-D structures that are distorted by 3-D electro-galvanic effects, need to be corrected for their proper conversion into reliable images of the subsurface electrical resistivity distribution. One of the most widely used approaches for correcting the data is the Groom–Bailey decomposition of the impedance tensor in terms of the unknown parameters of strike, twist and shear, along with the also unknown 2-D impedances. The standard approach for recovering the 2-D impedances is to solve numerically for all the unknowns as a non-linear inverse problem. In this work, we pose the recovery of the undistorted impedances in terms of a quadratic equation whose solutions filter out the distortions after a final tune for the appropriate shear parameter. The formula relies on two known invariants of the impedance tensor, the series and parallel invariants that have special immunity to some of the distorting parameters. Compared with the standard numerical recovery, the analytical solution provides more accurate results; they are clear-cut for both amplitudes and phases. The recovered amplitudes of the impedance are independent of strike, twist and shear, and the phases, in addition, of static effects. The recovery formula is transparent to random noise; hence the data preserve the original uncertainties. The applications to individual soundings of the COPROD2S1 data set and to a profile over a seismogenic region in Baja California, illustrate the effectiveness of the approach.

Beschreibung: Great progress has been made in the numerical simulation of planetary dynamos, though these numerical experiments still operate in a regime very far from the planets. For example, it seems unlikely that viscous forces are at all significant in planetary interiors, yet some of the simulations display a significant dependence on viscosity, and indeed in some of the simulations the dynamo mechanism is itself viscously driven, taking the form of helical Ekman pumping within columnar convection rolls. Given the similarity of the external magnetic fields observed in the terrestrial planets and gas giants, and the extremely small value of the Ekman number in all such cases, it seems natural to suppose that the underlying dynamo mechanism in these planets is simple, robust, independent of viscosity and insensitive to mechanical boundary conditions. A key step to identifying this mechanism is to determine the source of helicity in planetary cores, which itself should be robust, independent of viscosity and insensitive to boundary conditions. In this paper, we explore the possibility that the helicity in the core of the Earth arises from the spontaneous emission of inertial waves, driven by the equatorial heat flux in the outer core. We also ask if a similar mechanism might operate in other planets, and perhaps act to supplement the helicity driven by Ekman pumping in the (viscous) numerical simulations. We demonstrate that such waves do indeed produce the required helicity distribution outside the tangent cylinder. Moreover, we show that these waves inevitably propagate along the axis of the columnar vortices, and indeed they are the very mechanism by which the columnar vortices form in the first place and the means by which the columns subsequently evolve. We also calculate the emf induced by such axially propagating inertial waves and show that, in principle, this emf is sufficient to support a self-sustaining dynamo of the α 2 type. Finally, we derive the scaling laws for this kind of inertial-wave dynamo. We compare these predictions with the (imperfect) simulations, and also with what little we know about the Earth's core. The numerical experiments fall into two categories; the slowly rotating simulations which cannot sustain inertial waves at the small scales and the rapidly rotating (planet-like) ones which can. Our scaling laws are consistent with the latter class of simulations, and also with what we know about the Earth.

Beschreibung: The lack of robust water markets makes it difficult to value irrigation water. Because water rights are appurtenant to land, it is possible to infer the value of water from observed differences in the market price of land. We use panel data on repeat farmland sales in California's San Joaquin Valley to estimate a hedonic regression equation with parcel fixed effects. This controls for sources of omitted variables bias and allows us to recover the value of irrigation water to landowners in our sample. We show that a more traditional cross-sectional regression results in an artificially low value of irrigation water.

Beschreibung: The main sources of magnetic minerals in soils unaffected by anthropogenic pollution are iron oxides and hydroxides derived from parent materials through soil formation processes. Soil magnetic minerals can be used as indicators of environmental factors including soil forming processes, degree of pedogenesis, weathering processes and biological activities. In this study measurements of magnetic susceptibility are used to detect the presence and the concentration of soil magnetic minerals in topsoil and bulk samples in a small cultivated field, which forms a hydrological unit that can be considered to be representative of the rainfed agroecosystems of Mediterranean mountain environments. Additional magnetic studies such as isothermal remanent magnetization (IRM), anhysteretic remanent magnetization (ARM) and thermomagnetic measurements are used to identify and characterize the magnetic mineralogy of soil minerals. The objectives were to analyse the spatial variability of the magnetic parameters to assess whether topographic factors, soil redistribution processes, and soil properties such as soil texture, organic matter and carbonate contents analysed in this study, are related to the spatial distribution pattern of magnetic properties. The medians of mass specific magnetic susceptibility at low frequency ( lf ) were 36.0 and 31.1 10 –8 m 3 kg –1 in bulk and topsoil samples respectively. High correlation coefficients were found between the lf in topsoil and bulk core samples ( r  = 0.951, p  〈 0.01). In addition, volumetric magnetic susceptibility was measured in situ in the field ( is ) and values varied from 13.3 to 64.0 10 –5 SI. High correlation coefficients were found between lf in topsoil measured in the laboratory and volumetric magnetic susceptibility field measurements ( r  = 0.894, p  〈 0.01). The results obtained from magnetic studies such as IRM, ARM and thermomagnetic measurements show the presence of magnetite, which is the predominant magnetic carrier, and hematite. The predominance of superparamagnetic minerals in upper soil layers suggests enrichment in pedogenic minerals. The finer soil particles, the organic matter content and the magnetic susceptibility values are statistically correlated and their spatial variability is related to similar physical processes. Runoff redistributes soil components including magnetic minerals and exports fine particles out the field. This research contributed to further knowledge on the application of soil magnetic properties to derive useful information on soil processes in Mediterranean cultivated soils.

Beschreibung: State of the art numerical models of the Geodynamo are still performed in a parameter regime extremely remote from the values relevant to the physics of the Earth's core. In order to establish a connection between dynamo modelling and the geophysical motivation, it is necessary to use scaling laws. Such scaling laws establish the dependence of essential quantities (such as the magnetic field strength) on measured or controlled quantities. They allow for a direct confrontation of advanced models with geophysical constraints. We combine a numerical approach, based on a multiple linear regression method in the form of power laws, applied to a database of 102 direct numerical simulations (courtesy of U. Christensen), and a physical approach, based on energetics and forces balances. We show that previous empirical scaling laws for the magnetic field strength essentially reflect the statistical balance between energy production and dissipation for saturated dynamos. Such power based scaling laws are thus necessarily valid for any dynamo in statistical equilibrium and applicable to any numerical model, irrespectively of the dynamo mechanism. We show that direct numerical fits can provide contradictory results owing to biases in the parameters space covered in the numerics and to the role of a priori hypothesis on the fraction of ohmic dissipation. We introduce predictive scaling laws, that is relations involving input parameters of the governing equations only. We guide our reasoning on physical considerations. We show that our predictive scaling laws can properly describe the numerical database and reflect the dominant forces balance at work in these numerical simulations. We highlight the dependence of the magnetic field strength on the rotation rate. Finally, our results stress that available numerical models operate in a viscous dynamical regime, which is not relevant to the Earth's core.

Beschreibung: High-quality palaeomagnetic data for the early Carboniferous of Central Asia are scarce and the palaeogeographic evolution of this area prior to final amalgamation of the region east of the Ural mountains is still rather obscure. Here, we present palaeomagnetic data for early Carboniferous deposits from two areas in the Kyrgyz North Tianshan (NTS). Detailed rock-magnetic analysis indicates the presence of magnetite and haematite as magnetic carriers in these red sediments. In the Kazakh basin section (KEL), we identify a high-temperature component (HTC) of magnetization during stepwise thermal demagnetization at temperatures of up to ~680 °C yielding a site mean direction of D  = 176.2°, I  = –36.4°, k  = 57.4 and α 95  = 8.9° after tilt correction. Two HTCs of magnetization were identified in samples from the Sonkul Basin (DUN) with maximum blocking temperatures of ~600 °C (magnetite) and ~680 °C (haematite). The magnetite component was also identified with alternating field demagnetization. The resulting site mean directions for these two components identified in 16 and 14 sites, respectively, are D  = 149.3°, I  = –50.3°, k  = 73.6 and α 95  = 4.3° for the magnetite and D  = 139.6°, I  = –35.1°, k  = 71.6 and α 95  = 4.7° for the haematite component. All three mean directions show a significant increase of the precision parameter k after tilt correction indicating acquisition of the high-temperature magnetization prior to the main folding event in the Jurassic. We explain the difference of the two components of DUN by a process of inclination bias due to compaction to which the platy haematite particles are more susceptible. Applying the elongation-inclination (E/I) method to directional data from over 100 individual samples from location DUN results in a negligible correction for the magnetite component (〈5°), whereas the inclination of the haematite component corrects from –35.0° to –50.3° ( f  = 0.6, error interval –41.4° to –57.9°), which is then equal to the uncorrected magnetite inclination. The small number of samples from section KEL does not allow application of the E/I technique and inclination correction based on high field anisotropy of isothermal remanent magnetization was applied, yielding a corrected inclination of –75.2° ± 4°. Assuming comparable degrees of compaction for both study areas and applying the flattening factor obtained in DUN on samples from KEL, however, would result in comparable inclinations. The identification of inclination shallowing at both sections indicates that the age of magnetization is close to the deposition age. Assuming a reversed polarity of the directions from both areas results in palaeolatitudes of ~30°N for section DUN and ~60°N for the anisotropy-based correction of section KEL. The large difference, however, is geologically very unlikely. The inclination of the magnetite component of DUN (unaffected by inclination shallowing) favours a palaeoposition of ~30°N. This is supported by the inclination shallowing corrected haematite component of DUN yielding a comparable inclination. Therefore, our results indicate that the NTS domain was situated at ~30°N in the early Carboniferous. Furthermore, the NTS zone was probably not connected to Baltica or Siberia prior to the late Palaeozoic.

Beschreibung: An exploratory 3-D model of the electrical conductivity structure of the Australian continent is presented. The model is derived from the inversion of vertical magnetic-field transfer functions from the Australia-wide Array of Geomagnetic Stations. Crustal conductivity anomalies evident in the model are consistent with those previously mapped by independent magnetometer array studies and new electrical structures are suggested in the upper mantle. The model represents a seamless continent-scale basis for further models likely to be derived from subsequent studies. The model reveals three upper-mantle enhanced-conductivity anomalies beneath Archaean cratonic regions and two upper-mantle anomalies beneath Phanerozoic terranes in eastern Australia. Two of these anomalies have been investigated by recent magnetotelluric (MT) surveys, one in the Yilgarn Craton–Officer Basin–Musgrave Block the other in the Gawler Craton region, and are consistent with the MT results. Across much of central Australia enhanced conductivity at depths of 50–100 km is observed in the model. This region corresponds well with a recognized seismic velocity gradient at 75–100 km. Conductivity differences are also observed beneath Archaean cratons in Western Australia. The Pilbara Craton is represented as an enhanced conductivity anomaly at about 100 km, corresponding well with the lower-velocity anomaly evident in surface wave tomography models. The Yilgarn Craton is imaged as a low-conductivity body, with conductivity two orders of magnitude lower than the Pilbara Craton, continuing to greater depths.

Beschreibung: This paper presents new directional archeomagnetic data from nine Meso-/Neolithic fireplaces, sampled in a cave shelter, at Arconciel, in western Switzerland. Rock magnetic measurements indicate a homogenous magnetic mineralogy in all fireplaces, with magnetite as the main magnetic carrier. The remanent magnetization is stable and generally shows one characteristic directional component. Nine new directions, which were obtained from Arconciel, are combined with 356 other archeomagnetic data from a circular area with a radius of 700 km around this site, to obtain a penalized least square spline fit for the past 9000 yr. We found in general good agreement with other local compilations, such as the Balkan curve, the regional SCHA.DIF.8k model and with lake sediments from UK, Fennoscandia and Switzerland. Nevertheless, a time lag of several centuries is observed for a declination maximum between the archeomagnetic spline fit and the other European data records around 5900 BC. This time lag is also observed in the Swiss lake sediment record; therefore we interpret this shift as a local feature of the Earth's magnetic field.

Beschreibung: We present a new approach for 2-D magnetotelluric forward numerical modelling in contrast to traditional numerical methods like finite elements or finite differences. The method used for solving the partial differential equations is based on a mesh-free technique which does not need an underlaying mesh or grid. We use the Meshless Local Petrov–Galerkin (MLPG) method in combination with radial basis functions to simulate the response of a given conductivity model to a plane-wave source. We compare the mesh-free solution with known simulation programs and simple analytical solutions. Furthermore, we discuss the new magnetotelluric modelling method in terms of implementation and stability. First, we study the convergence and discretization errors of the new method with a simple half-space conductivity model. Then we compare our mesh-free simulation results with simple 2-D conductivity models with the results of a well-known finite element program. In the end, we provide a smooth conductivity model calculated with the mesh-free approach. The modelling results, even with randomly distributed nodes, are in a good agreement with those obtained by the finite element method.

Beschreibung: This article derives new implications for the land allocation and production decisions of profit-maximizing farm-firms where production of different crops is non-joint but subject to a constraint on the total land area. These implications are observable and thus subject to empirical scrutiny. An estimable model of crop production, land allocation, and input-use decisions is derived that permits joint production and enables the implications of non-joint but land-constrained production to be tested. This may improve econometric estimates of cross-price elasticities of supply by linking models of land use and production decisions, and allowing non-jointness to be imposed where appropriate.

Beschreibung: There is broad debate about including agriculture in greenhouse gas (GHG) reduction efforts such as the European Emissions Trading Scheme. Since most agricultural GHG emissions originate from non-point sources, they cannot be directly measured, and therefore have to be derived by calculation schemes (indicators). We designed five such GHG indicators for dairy farms and analyzed the trade-offs between their feasibility, measurement accuracy, and level of induced abatement costs. Analyses of induced abatement costs and calculation accuracy are based on emission reduction simulations with a highly-detailed single-farm optimization model. Feasibility is discussed in a qualitative manner. Our results indicate that the trade-offs depend on both farm characteristics and the targeted reduction level. In particular, the advantages of detailed indicators decrease for higher abatement levels. Only the least feasible indicator led to abatement costs that would result in emission efforts at given prices in the European Emissions Trading Scheme, although with a rather small potential. Our results thus suggest little potential for including dairy production into market-based reduction policies.

Beschreibung: Many agricultural support payments are based on past production with restrictions on how land may currently be used. When support payments to field crops are analyzed in a static framework, they do not directly impact current production decisions. However, over time, as relative profits change, these payments affect current output. The payments may keep land in less profitable production of program crops through restrictions prohibiting potentially more profitable endeavors such as cultivating fruits and vegetables. These payments have the potential to lead to production and trade distortions similar in magnitude to the distortions associated with direct production subsidies.

Beschreibung: Voluntary approaches have been used in a variety of contexts and for a variety of purposes in agriculture, including voluntary conservation programs and product labeling. This paper provides an overview of some of the general principles that emerge from the literature on voluntary approaches and their application in agriculture. The literature suggests that, to be effective, voluntary approaches must provide sufficiently strong participation incentives to a targeted population, clearly identify standards for behavior or performance that ensure additionality and avoid slippage, and monitor outcomes. Thus, reliance on voluntary approaches in agriculture is likely to be effective only if there is sufficient market demand for certain product characteristics, significant public funds are committed to pay for voluntary actions, or the political will exists to impose regulations if voluntary approaches fail.

Beschreibung: This paper analyses the impact of the recent decision by the European Union to ‘decouple’ agricultural support payments from agricultural production on Irish farmers' land market decisions. The land market participation decisions of Irish farmers are modelled using a dynamic probit model, while the extent of participation decisions is modelled using a dynamic tobit model. Decoupling does not appear to have significantly altered farmers' land market decisions. One likely explanation for this is the cross-compliance obligation for farmers to maintain land in a state fit for agricultural production in order to receive their full payments.

Beschreibung: We investigate farm size inequality in France using agricultural censuses and farm structure surveys at the NUTS3 level (‘départements’) during the period 1970–2007. Using calculated Gini coefficients, we show that farm size inequality has not systematically increased in France. An econometric analysis of the determinants of farm size inequality reveals that policy measures significantly affected farm size inequality, with most of the measures considered decreasing it. Empirical results suggest that the main contributor was the activity of the SAFER (Société d'Aménagement Foncier et d'Etablissement Rural), a specific feature of the French farm structural policy aimed at regulating rural land management. Besides, this research highlights the great complexity of the dynamics underlying the evolution of farm size distribution.

Beschreibung: This paper analyses implementation policies of environmental quota trade, with the Flemish nutrient production rights as an example. Implementation policies concern the transaction quantity, quota reduction and prevention of speculation. They are analysed with a static and a dynamic multi-agent quota trade model. The static model with discrete non-auctioned quota trade shows that the obligation for quota sellers to entirely stop their production stimulates structural change. The dynamic model version indicates that a flat rate reduction on traded quota and measures taken to prevent speculation combined with too low penalties for overuse stimulate the total production.