Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

Yucatán karst features and the size of Chicxulub crater

The buried Chicxulub impact structure is marked by a dramatic ring of sinkholes (called cenotes if containing water), and adjacent less prominent partial rings, which have been shown to coincide with maxima in horizontal gravity gradients and a topographic depression. These observations, along with the discreteness and spacing of the features, suggest a formation mechanism involving faulting in the outer slump zone of the crater, which would thus have a diameter of approximately 180 km.
An opposing view, based primarily on the interpretation of gravity data, is that (he crater is much larger than the cenote ring implies. Given the association of the known cenote ring with faults, we here examine northern Yucatan for similar rings in gravity, surface features and elevation, which we might expect to be associated with outer concentric faults in the case of a larger, possibly multiring, structure.
No such outer rings have been found, although definite patterns are seen in the distribution of karst features outside the crater rim. We explain these patterns as resulting mainly from deformation related to the block fault zone that parallels tbe shelf edge of eastern Yucatan.     

The use of GPR and VES in delineating a contamination plume in a landfill site: a case study in SE Brazil

This paper presents the results of the application of the Ground Penetrating Radar (GPR) method, or Georadar, in outlining a zone of contamination due to solid residues at the waste burial site of Rio Claro in the state of São Paulo, SE Brazil. A total of eight GPR profiles with 50- and 100-MHz antennae were surveyed. Six profiles were located within the landfill site and the remaining two were outside. The main objective of the GPR survey was to evaluate the side extension of contamination. A Vertical Electric Sounding (VES) survey was performed at four points within the site in order to investigate the vertical extent of the contamination plume and to define the bottom of the landfill. Two additional VESs were done outside the landfill with the purpose of determining the top of the ground water table and the geoelectric stratigraphy of the background. From the interpretation of the GPR profiles, it was possible to locate the top of the contamination plume and to infer that it was migrating laterally beyond the limits of the waste disposal site. This was observed along the profile situated close to the highway SP-127, which was about 20 m from the limit of the site. The signature of the contaminant appears as a discontinuous reflector that is believed to be a shallow ground water table. The discontinuity is marked by a shadow zone, which is characteristic of conductive contaminant residues. The contamination did not move far enough to reach a sugar cane plantation located at approximately 100 m from the border of the site. In the regions free from contamination, the ground water table was mapped at approximately 10 m of depth, and it was characterized by a strong and continuous reflector. The radar signal penetrated deep enough and enabled the identification of a second reflector at approximately 14 m deep, interpreted as the contact between the Rio Claro and the Corumbataí formations. The contact is marked by the presence of gravel characterized by ferruginous concretes, which cause the strong amplitude reflection in the GPR profile. Within the landfill site, the quantitative interpretation of the VES results showed the contamination zone. The base of the landfill varies between 11 and 15 m deep. Outside the landfill site, the VES results showed no indication of contamination and allowed the determination of the top of the ground water table and the contact between the Rio Claro and the Corumbataí formations. The results of GPR and VES showed a good agreement and the integrated interpretations were supported by local geology and information from several boreholes, about 17 m depth, on average. The bottom of the landfill reaches a maximum of 14.5 m depth.     

Sedimentology and palynology of the Calafate Formation (Maastrichtian), Austral Basin, Southern Patagonia, Argentina

The Calafate Formation crops out in south-western Santa Cruz Province, Argentina, and displays a stacking of asymmetrical coarsening–fining-upward cycles. These cycles are interpreted as the product of short-lived transgressive-regressive events in which the coarsening upward part represents sedimentary aggradation with a stable or decreasing sea level. Sedimentological and palynological analyses indicate nearshore marine conditions. Even though the existence of an estuary or incised valley cannot be determined, this is the most probable palaeogeographic model. Based on dinoflagellate cysts, the base of the section is considered to be not older than Maastrichtian. The presence of the oyster Ambigostrea clarae (Ihering) occurring together with the dinoflagellate cyst species Manumiella druggii (Stover) Bujak and Davies and Eisenackia circumtabulata Drugg in the middle part of the section indicates an age no older than late Maastrichtian. According to sedimentological data, deposits representing the Cretaceous–Palaeogene transition would have been eroded, which is confirmed by the presence of Grapnelispora loncochensis Papú. This megaspore is a consistent component of the Maastrichtian assemblages from Patagonia.     

Vegetation dynamics during the late Pleistocene in the Barreirinhas region, Maranhão State, northeastern Brazil, based on carbon isotopes in soil organic matter

The study place is in the Barreirinhas region, Maranhão State, northeastern Brazil. A vegetation transect of 78 km was studied among four vegetation types: Restinga (coastal vegetation), Cerrado (woody savanna), Cerradão (dense woody savanna), and Forest, as well as three forested sites around Lagoa do Caçó, located approximately 10 km of the transect. Soil profiles in this transect were sampled for δ¹³C analysis, as well as buried charcoal fragments were used for ¹⁴C dating. The data interpretation indicated that approximately between 15,000 and ∼9000 ¹⁴C yr B.P., arboreal vegetation prevailed in the whole transect, probably due to the presence of a humid climate. Approximately between ∼9000 and 4000-3000 ¹⁴C yr B.P., there was the expansion of the savanna, probably related to the presence of drier climate. From ∼4000-3000 ¹⁴C yr B.P. to the present, the results indicated an increase in the arboreal density in the area, due to the return to a more humid and probably similar climate to the present. The presence of buried charcoal fragments in several soil depths suggested the occurrence of palaeofires during the Holocene. The vegetation dynamic inferred in this study for northeastern Brazil is in agreement with the results obtained in areas of Amazon region, based on pollen analysis of lake sediments and carbon isotope analysis of soil organic matter (SOM), implying than similar climatic conditions have affected these areas during the late Pleistocene until the present.     

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.     

Crustal and mantle strengths in continental lithosphere: is the jelly sandwich model obsolete?

The relative importance of the contribution of the lower crust and of the lithospheric mantle to the total strength of the continental lithosphere is assessed systematically for realistic ranges of layer thickness, composition, and temperature. Results are presented as relative strength maps, giving the ratio of the lower crust to upper mantle contribution in terms of crustal thickness and surface heat flow. The lithosphere shows a “jelly sandwich” rheological layering for low surface heat flow, thin to average crustal thickness, and felsic or wet mafic lower crustal compositions. On the other hand, most of the total strength resides in the seismogenic crust in regions of high surface heat flow, crust of any thickness, and dry mafic lower crustal composition.     

Coastal neotectonics in Southern Central Andes: uplift and deformation of marine terraces in Northern Chile (27°S)

Neotectonic observations allow a new interpretation of the recent tectonic behaviour of the outer fore arc in the Caldera area, northern Chile (27°S). Two periods of deformation are distinguished, based on large-scale Neogene to Quaternary features of the westernmost part of the Coastal Cordillera: Late Miocene to Early Pliocene deformations, characterized by a weak NE–SW to E–W extension is followed by uppermost Pliocene NW–SE to E–W compression. The Middle Pleistocene to Recent time is characterized by vertical uplift and NW–SE extension. These deformations provide clear indications of the occurrence of moderate to large earthquakes. Microseismic observations, however, indicate a lack of shallow crustal seismicity in coastal zone. We propose that both long-term brittle deformation and uplift are linked to the subduction seismic cycle.