Mineral assemblages of the Francisco I. Madero Zn–Cu–Pb–(Ag) deposit, Zacatecas, Mexico: Implications for ore deposit genesis

The Francisco I. Madero deposit, central Mexico, occurs in the Mesozoic Guerrero Terrane, which hosts many ore deposits, both Cretaceous (volcanogenic massive sulfides) and Tertiary (epithermal and skarn deposits). It is hosted by a 600 m-thick calcareous-pelitic unit, of Lower Cretaceous age, crosscut by porphyritic dikes that strike NW–SE. A thick felsic volcanic Tertiary sequence, consisting of andesites and rhyolitic ignimbrites, unconformably overlies the Cretaceous series. At the base, the mineralization consists of several mantos developed within calcareous beds. They are dominantly composed of sphalerite, pyrrhotite and pyrite with minor chalcopyrite, arsenopyrite and galena. At the top of the orebody, there are calcic skarns formed through prograde and retrograde stages. The resulting mineral assemblages are rich in manganoan hedenbergite (Hd_75–28Di_40–4Jh_40–20), andraditic garnets (Adr_100–62Grs_38–0), epidote (Ep_95–36Czo_60–5Pie_8–0), chamosite, calcite and quartz. The temperature of ore deposition, estimated by chlorite and arsenopyrite geothermometry, ranges from 243° to 277 °C and from 300° to 340 °C, respectively. The pressure estimated from sphalerite geobarometry averages 2.1 kbar. This value corresponds to a moderately deep skarn and agrees with the high Cu content of the deposit. Paragenesis, P–T conditions and geological characteristics are compatible with a distal, dike-related, Zn skarn deposit. Its style of mineralization is similar to that of many high-temperature carbonate replacement skarn deposits in the Southern Cordillera.     

New thermochronologic constraints on the evolution of the Zaldívar porphyry copper deposit,Northern Chile

Life spans and thermal evolution of hydrothermal systems are of fundamental metallogenic importance. We were able to establish the chronology and cooling history of the Zaldívar porphyry copper deposit (Northern Chile) by applying a combination of different isotopic dating methods in minerals with different closure temperatures, including ⁴⁰Ar/³⁹Ar geochronology and zircon fission track thermochronology, together with fluid inclusion thermometry and previous published U–Pb zircon geochronology. The hydrothermal mineralization in the Zaldívar deposit is genetically related to the Llamo Porphyry unit. Samples of igneous biotites from this intrusion yielded ⁴⁰Ar/³⁹Ar plateau ages between 35.5 ± 0.7 and 37.7 ± 0.4 Ma defining a weighted average of 36.6 ± 0.5 Ma (2σ). In contrast, one sample from the Zaldívar porphyry, one from the andesites, and two from the Llamo porphyry yielded considerably younger fission track ages of approximately 29 Ma with a weighted mean for all ages of 29.1 ± 1.7 Ma (2σ). Thermal and compositional constraints for the hydrothermal system in the Zaldívar deposit from fluid inclusions thermometry show that at least three fluid types broadly characterize two main hydrothermal episodes during the evolution of the deposit. The main mineralization and alteration event is characterized by high temperature (above 320°C) hypersaline fluids (salinity between 30 and 56 wt.% NaCl equivalents) coexisting with low-density gas-rich inclusions (salinity less than 17 wt.% NaCl equivalents) that homogenizing into the gas phase at temperatures above 350°C. The second episode corresponds to a low-temperature event which is characterized by liquid-rich fluid inclusions that homogenize into the liquid phase at temperatures ranging from 200°C to 300°C with salinities lower than 10 wt.% NaCl equivalents. The ⁴⁰Ar/³⁹Ar data (36.6 ± 0.5 Ma, weighted average) obtained from igneous biotites represent the minimum age for the last high-temperature (above 300°C) hydrothermal pulse. When compared with previously published U–Pb ages (38.7 ± 1.3 Ma) in zircons from the Llamo porphyry, a close temporal relationship between crystallization of the parental intrusion and the thermal collapse of the last high-temperature hydrothermal event is evident. Cooling took place from approximately 800°C (crystallization of the intrusive complex defined by zircon U–Pb ages) to below 300 ± 50°C (biotite ⁴⁰Ar/³⁹Ar closure temperature) within approximately 1.5 m.y. Because the thermal annealing of fission tracks in zircons occurs at temperatures of 240 ± 30°, the zircon fission track (ZFT) ages of 29.1 ± 1.7 Ma (2σ) mark the end of the thermal activity in the Zaldívar area, specifically the time when the whole area cooled below this temperature, well after the collapse of the main hydrothermal event in the Zaldívar porphyry copper deposit. This cooling age roughly coincides with the age defined for the emplacement of dacitic dikes at 31 ± 2.8 Ma (2σ) (published K–Ar whole rock), 5 km south of the Zaldívar deposit, in the Escondida area. This late magmatic pulse probably is responsible for high heat flow in the Zaldívar deposit as late as 29 Ma. There is no evidence that the low temperature hydrothermal pulse recognized by fluid inclusion studies is related to this thermal event. The zircon fission track cooling ages are interpreted to be related to the time lag required for complete relaxation of the perturbation of the isotherms in the geothermal field imposed by the intrusion of magmatic bodies, with or without any association with low temperature hydrothermal activity.     

Large geographic and temporal extensions of the Río de la Plata Craton,South America,and its metacratonic eastern margin

Integration of existing isotopic and geological data allows a reconsideration of the distribution and age of the Río de la Plata Craton within South America. The reinterpretation increases the area of the craton to about 2,400,000 km² with implications for the tectonic map of South America and for global reconstruction of palaeocontinents. Four areas previously considered as separate cratons (Luís Alves, Curitiba, Tebicuary, and Paranapanema) are interpreted as part of the same Río de la Plata Craton. The craton is organized into six provinces and domains: Buenos Aires–Piedra Alta, Taquarembó, Tebicuary, Luís Alves, Encantadas, and Nico Pérez. The term ‘Transplatense’ is proposed to replace ‘Trans-Amazonian’ for Rhyacian events that occurred within the Río de la Plata Craton. The craton is formed not only by dominant Rhyacian rocks and local Archaean rocks, but also by Statherian and Mesoproterozoic rocks. The domains are all partially to totally covered by Phanerozoic basins (Paraná, Chacoparanense, Claromecó, Salado, Balcarce, and Colorado) which makes their investigation difficult. The Ediacaran–Cambrian collisions of the Brasilian orogen generated tectonic mixtures of orogenic and cratonic zones. This is more evident in the eastern margin of the craton, which behaved as a metacraton.     

Relationship between global mean sea-level and global mean temperature in a climate simulation of the past millennium

The possibility of using global mean near-surface temperature, its rate of change or the global mean ocean heat-flux as predictors to statistically estimate the change of global mean sea-level is explored in the context of a long climate simulation of the past millennium with the climate model ECHO-G. Such relationships have recently been proposed to by-pass the difficulty of estimating future sea-level changes based on simulations with coarse-resolution climate models. It is found that, in this simulation, a simple linear relationship between mean temperature and the rate of change of sea level does not exist. A regression parameter linking both variables, and estimated in sliding 120-year windows, varies widely along the simulation and, in some periods, even attains negative values. The ocean heat-flux and the rate-of-change of mean temperature seem to better capture the rate-of-change of sea level due to thermal expansion.     

Modelling and migration with orthogonal isochron rays

For increasing time values, isochrons can be regarded as expanding wavefronts and their perpendicular lines as the associated orthogonal isochron rays. The speed of the isochron movement depends on the medium velocity and the source-receiver position. We introduce the term equivalent-velocity to refer to the speed of isochron movement. In the particular case of zero-offset data, the equivalent velocity is half of the medium velocity. We use the concepts of orthogonal isochron-rays and equivalent velocity to extend the application of the exploding reflector model to non-zero offset imaging problems. In particular, we employ these concepts to extend the use of zero-offset wave-equation algorithms for modelling and imaging common-offset sections. In our imaging approach, the common-offset migration is implemented as a trace-by-trace algorithm in three steps: equivalent velocity computation, data conditioning for zero-offset migration and zero-offset wave-equation migration. We apply this methodology for modelling and imaging synthetic common-offset sections using two kinds of algorithms: finite-difference and split-step wavefield extrapolation. We also illustrate the isochron-ray imaging methodology with a field-data example and compare the results with conventional common-offset Kirchhoff migration. This methodology is attractive because it permits depth migration of common-offset sections or just pieces of that by using wave-equation algorithms, it extends the use of robust zero-offset algorithms, it presents favourable features for parallel processing, it permits the creation of hybrid migration algorithms and it is appropriate for migration velocity analysis.     

Daily hydrological modeling in the Amazon basin using TRMM rainfall estimates

Rainfall measurements by conventional raingauges provide relatively accurate estimates at a few points of a region. The actual rainfield can be approximated by interpolating the available raingauge data to the remaining of the area of interest. In places with relatively low gauge density such interpolated rainfields will be very rough estimates of the actual events. This is especially true for tropical regions where most rainfall has a convective origin with high spatial variability at the daily level. Estimates of rainfall by remote sensing can be very useful in regions such as the Amazon basin, where raingauge density is very low and rainfall highly variable. This paper evaluates the rainfall estimates of the Tropical Rainfall Measuring Mission (TRMM) satellite over the Tapajós river basin, a major tributary of the Amazon. Three-hour TRMM rainfall estimates were aggregated to daily values and were compared with catch of ground-level precipitation gauges on a daily basis after interpolating both data to a regular grid. Both daily TRMM and raingauge-interpolated rainfields were then used as input to a large-scale hydrological model for the whole basin; the calculated hydrographs were then compared to observations at several streamgauges along the river Tapajos and its main tributaries. Results of the rainfield comparisons showed that satellite estimates can be a practical tool for identifying damaged or aberrant raingauges at a basin-wide scale. Results of the hydrological modeling showed that TRMM-based calculated hydrographs are comparable with those obtained using raingauge data.