Lu–Hf geochronology and trace element distribution in garnet: Implications for uplift and exhumation of ultra-high pressure granulites in the Sudetes, SW Poland

Combining Lu–Hf garnet geochronology with in situ trace element analyses in garnet allowed us to gain new insight into the metamorphic evolution of UHP–UHT rocks in the Stary Gierałtów region, in the Polish Sudetes. Prograde garnet growth recorded by Rayleigh-type heavy REE (HREE) zoning in the felsic granulites indicates that the obtained 386.6 ± 4.9 Ma Lu–Hf age represents the time of garnet crystallization on a prograde UHP metamorphic path. The surrounding rocks were metamorphosed at the same time as indicated by 381.2 ± 6.7 Ma Sm–Nd garnet age obtained for the mid-crustal metapelites. The second metamorphic episode, which affected most of the lower crust in the Orlica–Śnieżnik Massif (OSM) occurred at ca. 340 Ma as determined by U–Pb zircon and Sm–Nd garnet dating of granulites in this and previous studies is interpreted as a high temperature event, which took place on a retrograde path.

Trace element distribution in garnets from the layered granulites showed significant differences in distribution of medium and HREE in garnets from mafic and felsic protoliths over the course of the metamorphic evolution. This had strong impact on the isotopic dating results and led to “decoupling” of the Sm–Nd and Lu–Hf clocks, which recorded timing of the two different metamorphic episodes separated by as much as 40 Ma. Moreover, the preservation of the HREE growth zonation profile in garnets from the felsic granulites whose minimum metamorphic temperature was established at 900 °C implies that the Lu–Hf system under relatively dry conditions does not undergo significant diffusional re-equilibration even at such extreme temperatures and therefore it sill provides the age of prograde garnet growth. Under hydrous conditions, at least some resetting will take place, as documented by the partially relaxed HREE zonation profile in the amphibolitised mafic granulite, which yielded a 10 Ma younger age. The HREE distribution study appeared to be a particularly valuable and essential tool, which allowed us to distinguish garnet growth from post-growth complexities and hence, provide improved age interpretation. Medium REE, on the other hand, did not show any obvious correlation with the isotopic signature of garnet.

Two distinct metamorphic episodes recorded in the Stary Gierałtów region show that buoyancy-driven uplift of UHP rocks can be arrested at the base of a continental crust if not supported by any additional force. In our case study, the UHP rocks would have never reached the surface if their uplift had not been resumed after a long pause under a different tectonic regime. The multistage, discontinuous uplift revealed by the UHP rocks of the OSM provides a new scenario for the exhumation of continental crust from mantle depths distinct from the fast-track exhumation histories recognized in UHP terranes elsewhere.     

Pits,rifts and slumps: the summit structure of Piton de la Fournaise

A clear model of structures and associated stress fields of a volcano can provide a framework in which to study and monitor activity. We propose a volcano-tectonic model for the dynamics of the summit of Piton de la Fournaise (La Reunion Island, Indian Ocean). The summit contains two main pit crater structures (Dolomieu and Bory), two active rift zones, and a slumping eastern sector, all of which contribute to the actual fracture system. Dolomieu has developed over 100 years by sudden large collapse events and subsequent smaller drops that include terrace formation. Small intra-pit collapse scars and eruptive fissures are located along the southern floor of Dolomieu. The western pit wall of Dolomieu has a superficial inward dipping normal fault boundary connected to a deeper ring fault system. Outside Dolomieu, an oval extension zone containing sub-parallel pit-related fractures extends to a maximum distance of 225 m from the pit. At the summit the main trend for eruptive fissures is N80°, normal to the north–south rift zone. The terraced structure of Dolomieu has been reproduced by analogue models with a roof to width ratio of approximately 1, suggesting an original magma chamber depth of about 1 km. Such a chamber may continue to act as a storage location today. The east flank has a convex–concave profile and is bounded by strike-slip fractures that define a gravity slump. This zone is bound to the north by strike-slip fractures that may delineate a shear zone. The southern reciprocal shear zone is probably marked by an alignment of large scoria cones and is hidden by recent aa lavas. The slump head intersects Dolomieu pit and may slide on a hydrothermally altered layer known to be located at a depth of around 300 m. Our model has the summit activity controlled by the pit crater collapse structure, not the rifts. The rifts become important on the mid-flanks of the cone, away from pit-related fractures. On the east flank the superficial structures are controlled by the slump. We suggest that during pit subsidence intra-pit eruptions may occur. During tumescence, however, the pit system may become blocked and a flank eruption is more likely. Intrusions along the rift may cause deformation that subsequently increases the slump’s potential to deform. Conversely, slumping may influence the east flank stress distribution and locally control intrusion direction. These predictions can be tested with monitoring data to validate the model and, eventually, improve monitoring.     

3D structure of the Earth's crust beneath the northern part of the Bohemian Massif

The SUDETES 2003 wide-angle refraction/reflection experiment covered the area of the south-western Poland and the northern Bohemian Massif. The good quality data that were gathered combined with the data from previous experiments (POLONAISE'97, CELEBRATION 2000) allowed us to prepare a 3D seismic model of the crust and uppermost mantle for this area. We inverted travel times of both refracted and reflected P waves using the JIVE3D package. This allowed us to obtain a model of P-wave velocity distribution as well as the shape of major boundaries in the crust. We also present a detailed uncertainty analysis for both the boundary depths and the velocity field. In doing the uncertainty analysis we found an interesting, strong dependence between uncertainty and inversion scheme (order of used phases). We also compared the model with surface geology and found good correlation between velocity inhomogeneities in the uppermost crust (down to 2 km) and major geological units. The higher velocity lower crust (6.9–7.2 km/s) could result from remelting of the lower crust or magmatic underplating.     

Diurnal atmospheric forcing and temporal variations of the nutation amplitudes

A 29-year time-series of four-times-daily atmospheric effective angular momentum (EAM) estimates is used to study the atmospheric influence on nutation. The most important atmospheric contributions are found for the prograde annual (77 μas), retrograde annual (53 as), prograde semiannual (45 as), and for the constant offset of the pole (δψsinɛ₀=−86 as, δɛ=77 as). Among them only the prograde semiannual component is driven mostly by the wind term of the EAM function, while in all other cases the pressure term is dominant. These are nonnegligible quantities which should be taken into account in the new theory of nutation. Comparison with the VLBI corrections to the IAU 1980 nutation model taking into account the ocean tide contribution yields good agreement for the prograde annual and semiannual nutations. We also investigated time variability of the atmospheric contribution to the nutation amplitudes by performing the sliding-window least-squares analysis of both the atmospheric excitation and VLBI nutation data. Almost all detected variations of atmospheric origin can be attributed to the pressure term, the biggest being the in-phase annual prograde component (about 30 as) and the retrograde one (as much as 100200 as). These variations, if physical, limit the precision of classical modeling of nutation to the level of 0.1 mas. Comparison with the VLBI data shows significant correlation for the retrograde annual nutation after 1989, while for the prograde annual term there is a high correlation in shape but the size of the atmospherically driven variations is about three times less than deduced from the VLBI data. This discrepancy in size can be attributed either to inaccuracy of the theoretical transfer function or the frequency-dependent ocean response to the pressure variations. Our comparison also yields a considerably better agreement with the VLBI nutation data when using the EAM function without the IB correction for ocean response, which indicates that this correction is not adequate for nearly diurnal variations. Received: 10 September 1997 / Accepted: 5 March 1998     

Tourist farms in Lower Silesia, Poland

In view of the small size of most Polish farms, attempts have been made to encourage diversification into rural tourism. The paper examines the progress made in Lower Silesia which has good resources for domestic and international tourism. Thanks to the promotional work of the Agricultural Consultative Centres over a hundred farms are now providing services. They do not generate a large share of the farm income but business is expanding and some areas are becoming very popular. Thus a base has been established for future expansion. Further development depends on easier access to credit and government support for ecological agriculture. Better promotion is also needed to increase interest in rural tourism in Poland and attract more visitors from abroad. This revised version was published online in July 2006 with corrections to the Cover Date.     

Improvement of the geoid in local areas by satellite-to-satellite tracking

Summary The concept of satellite-to-satellite tracking measuring the relative velocity of two orbiting satellites spaced some hundreds kilometers on a close orbit, provides now possibilities for the investigation of the Earth’s gravity field. In the paper only medium and short wave length effects affecting the measured relative velocity have been considered. Collocation is used in such an analysis of local geoid improvement, because this method allows to combine heterogeneous data in a consistent way. Covariance functions relevant for the particular case of a circular equatorial orbit are given. Two kinds of observation equations have been formulated. The choice of observation equation with regard to satellites configuration is discussed. It is found that it is sufficient to have a limited number of satellite-to-satellite observations in a 7^o×7^o area around the estimation point with distances between profiles of about 1^o.5 and between the two satellites forming the pair of 200+350 km; the altitude of satellite-to-satellite observations should be as low as possible. The accuracy of the geoid determination strongly depends on the degree and order of the reference field used. An accuracy of about ±1 m can be achieved with an assumed reference field of (40,40). The influence of measuring errors is discussed and it is shown that only satellite-to-satellite observations with accuracy better then 0.1 mm/sec will give an improvement of the geoid. Finally, some results on the combination of low-low satellite-to-satellite tracking and terrestrial gravity data are given. The proposed method seems to be especially interesting for unsurveyed areas. Furthermore, it has the practical advantage that only a local coverage data is needed.     

Reconnaissance geochronology of basement rocks from the Amazonas and Maranhão Basins in Brazil

Rb/Sr whole rock as well as K/Ar whole rock and mineral age determinations have been carried out on drill core samples obtained from the crystalline basement underlying the Amazonas and Maranhão Basins in Brazil. The data confirm the existence of an ancient basement under the Upper and Middle Amazonas Basin. The primary age of the granitic rocks is estimated at about 1560 Ma which correlates with recent age data obtained on granitic rocks on both sides of the Amazonas syneclise. In some areas the K/Ar ages show a decrease probably due to the effect of the Nickerie thermal episode described previously in the Guayana Shield.In the basement of the interconnected Lower Amazonas and Maranhão Basins Late Precambrian (Caririan) metamorphism and igneous activity seems to have been pervasive. The data allow the continuation of the Brazilian age province known in northeastern Brazil under the Maranhão Basin, linking this province to the Brasilia-Paraguay orogenic belt of similar age.     

Improvement of the geoid in local areas by satellite gradiometry

Summary Satellite gradiometry is studied as a means to improve the geoid in local areas from a limited data coverage. Least-squares collocation is used for this purpose because it allows to combine heterogeneous data in a consistent way and to estimate the integrated effect of the attenuated spectrum. In this way accuracy studies can be performed in a general and reliable manner. It is shown that only three second-order gradients contribute significantly to the estimation of the geoidal undulations and that it is sufficient to have gradiometer data in a 5°×5° area around the estimation point. The accuracy of the geoid determination is strongly dependent on the degree and order of the reference field used. An accuracy of about ±1 m can be achieved with a reference field of (12, 12). There is an optimal satellite altitude for each reference field and this altitude may be higher than 300 km for a field of low degree and order. The influence of measuring errors is discussed and it is shown that only gradiometer data with accuracies better than ±0.05 E will give a significant improvement of the geoid. Finally, some results on the combination of satellite gradiometry and terrestrial gravity measurements are given. The proposed method seems to be well suited for local geoid determinations down to the meter range. It is especially interesting for unsurveyed and difficult areas because no terrestrial measurements are necessary. Furthermore, it has the practical advantage that only a local data coverage is needed.