Deep crustal structure of the Parana Basin from receiver functions and Rayleigh-wave dispersion: Evidence for a fragmented cratonic root

The deep crustal structure of the Parana Basin of southern Brazil is investigated by analyzing P- and PP-wave receiver functions at 17 Brazilian Lithosphere Seismic Project stations within the basin. The study area can be described as a typical Paleozoic intracratonic basin that hosts one of the largest Large Igneous Province of the world and makes a unique setting for investigating models of basin subsidence and their interaction with mantle plumes. Our study consists of (1) an analysis of the Moho interaction phases in the receiver functions to obtain the thickness and bulk Vp/Vs ratio of the basin`s underlying crust and (2) a joint inversion with Rayleigh-wave dispersion velocities from an independent tomographic study to delineate the detailed S-wave velocity variation with depth. The results of our analysis reveal that Moho depths and bulk Vp/Vs ratios (including sediments) vary between 41 and 48 km and between 1.70 and 1.76, respectively, with the largest values roughly coinciding with the basin`s axis, and that S-wave velocities in the lower crust are generally below 3.8 km/s. Select sites within the basin, however, show lower crustal S-wave velocities slightly above 3.9 km/s suggestive of underplated mafic material. We show that these observations are consistent with a fragmented cratonic root under the Parana basin that defined a zone of weakness for the initial Paleozoic subsidence of the basin and which allowed localized mafic underplating of the crust along the suture zones by Cenozoic magmatism.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Coda wave attenuation in the Parecis Basin, Amazon Craton, Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos GaA(0)chos, Amazon craton-Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q (c)) of coda waves modeled usingwhere Q (0) is the coda quality factor at frequency of 1 Hz and eta is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239-252, 2008), based in the temporal attenuation coefficient, chi(f) instead of Q(f), given by the equation for the calculation of the geometrical attenuation (gamma) and effective attenuation Q (c) values have been computed at central frequencies (and band) of 1.5 (1-2), 3.0 (2-4), 6.0 (4-8), 9.0 (6-12), 12 (8-16), and 18 (12-24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin for the surrounding shield and for the whole region of Porto dos GaA(0)chos Using the independent frequency model, we found: for the cratonic zone, gamma = 0.014 s (-aEuro parts per thousand 1), nu a parts per thousand 1.12; for the basin zone with sediments of similar to 500 m, gamma = 0.031 s (-aEuro parts per thousand 1), nu a parts per thousand 1.27; and for the Parecis basin with sediments of similar to 1,000 m, gamma = 0.047 s (-aEuro parts per thousand 1), nu a parts per thousand 1.42. Analysis of the attenuation factor (Q (c)) for different values of the geometrical spreading parameter (nu) indicated that an increase of nu generally causes an increase in Q (c), both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, (in the deepest part of the basin), than in the basement, (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.

Astrophysical S factors for fusion reactions involving C, O, Ne, and Mg isotopes

Using the Sao Paulo potential and the barrier penetration formalism we have calculated the astrophysical factor S(E) for 946 fusion reactions involving stable and neutron-rich isotopes of C, O, Ne, and Mg for center-of-mass energies E varying from 2 to approximate to 18-30 MeV (covering the range below and above the Coulomb barrier). We have parameterized the energy dependence, S(E), by an accurate universal 9-parameter analytic expression and present tables of fit parameters for all the reactions. We also discuss the reduced 3-parameter version of our fit which is highly accurate at energies below the Coulomb barrier, and outline the procedure for calculating the reaction rates. The results can be easily converted to thermonuclear or pycnonuclear reaction rates to simulate various nuclear burning phenomena, in particular, stellar burning at high temperatures and nucleosynthesis in high density environments. (C) 2010 Elsevier Inc. All rights reserved

The Itajai foreland basin: a tectono-sedimentary record of the Ediacaran period, Southern Brazil

The Itajai Basin located in the southern border of the Luis Alves Microplate is considered as a peripheral foreland basin related to the Dom Feliciano Belt. It presents an excellent record of the Ediacaran period, and its upper parts display the best Brazilian example of Precambrian turbiditic deposits. The basal succession of Itajai Group is represented by sandstones and conglomerates (BaA(0) Formation) deposited in alluvial and deltaic-fan systems. The marine upper sequences correspond to the Ribeiro Carvalho (channelized and non-channelized proximal silty-argillaceous rhythmic turbidites), Ribeiro Neisse (arkosic sandstones and siltites), and Ribeiro do Bode (distal silty turbidites) formations. The ApiA(0)na Formation felsic volcanic rocks crosscut the sedimentary succession. The Cambrian Subida leucosyenogranite represents the last felsic magmatic activity to affect the Itajai Basin. The Brusque Group and the Florianpolis Batholith are proposed as source areas for the sediments of the upper sequence. For the lower continental units the source areas are the Santa Catarina, So Miguel and CamboriA(0) complexes. The lack of any oceanic crust in the Itajai Basin suggests that the marine units were deposited in a restricted, internal sea. The sedimentation started around 600 Ma and ended before 560 Ma as indicated by the emplacement of rhyolitic domes. The Itajai Basin is temporally and tectonically correlated with the Camaqu Basin in Rio Grande do Sul and the Arroyo del Soldado/Piriapolis Basin in Uruguay. It also has several tectono-sedimentary characteristics in common with the African-equivalent Nama Basin.

High T/P evolution and metamorphic ages of the migmatitic basement of northern Sierras Pampeanas, Argentina: Characterization of a mid-crustal segment of the Famatinian belt

New petrologic, thermobarometric and U-Pb monazite geochronologic information allowed to resolve the metamorphic evolution of a high temperature mid-crustal segment of an ancient subduction-related orogen. The EI Portezuelo Metamorphic-Igneous Complex, in the northern Sierras Pampeanas, is mainly composed of migmatites that evolved from amphibolite to granulite metamorphic facies, reaching thermal peak conditions of 670-820 degrees C and 4.5-5.3 kbar. The petrographic study combined with conventional and pseudosection thermobarometry led to deducing a short prograde metamorphic evolution within migmatite blocks. The garnet-absent migmatites represent amphibolite-facies rocks, whereas the cordierite-garnet-K-feldspar-sillimanite migmatites represent higher metamorphic grade rocks. U-Pb geochronology on monazite grains within leucosome record the time of migmatization between approximate to 477 and 470 Ma. Thus, the El Portezuelo Metamorphic-Igneous Complex is an example of exhumed Early Ordovician anatectic middle crust of the Famatinian mobile belt. Homogeneous exposure of similar paleo-depths throughout the Famatinian back-arc and isobaric cooling paths suggest slow exhumation and consequent longstanding crustal residence at high temperatures. High thermal gradients uniformly distributed in the Famatinian back-arc can be explained by shallow convection of a low-viscosity asthenosphere promoted by subducting-slab dehydration. (C) 2011 Elsevier Ltd. All rights reserved.

Evaluation of petrogenetic models for intermediate and silicic plutonic rocks from the Sierra de Valle Fertil-La Huerta, Argentina: Petrologic constraints on the origin of igneous rocks in the Ordovician Famatinian-Puna paleoarc

The whole Valle Fertil-La Huerta section appears as a calc-alkaline plutonic suite typical of a destructive plate margin. New Sr and Nd isotopic whole-rock data and published whole-rock geochemistry suggest that the less-evolved intermediate (dioritic) rocks can be derived by magmatic differentiation, mainly by hornblende + plagioclase +/- Fe-Ti oxide fractional crystallization, from mafic (gabbroic) igneous precursors. Closed-system differentiation, however, cannot produce the typical intermediate (tonalitic) and silicic (granodioritic) plutonic rocks, which requires a preponderant contribution of crustal components. Intermediate and silicic plutonic rocks from Valle Fertil-La Huerta section have formed in a plate subduction setting where the thermal and material input of mantle-derived magmas promoted fusion of fertile metasedimentary rocks and favored mixing of gabbroic or dioritic magmas with crustal granitic melts. Magma mixing is observable in the field and evident in variations of chemical elemental parameters and isotopic ratios, revealing that hybridization coupled with fractionation of magmas took place in the crust. Consideration of the whole-rock geochemical and isotopic data in the context of the Famatinian-Puna magmatic belt as a whole demonstrates that the petrologic model postulated for the Sierra Valle Fertil-La Huerta section has the potential to explain the generation of plutonic and volcanic rocks across the Early Ordovician paleoarc from central and northwestern Argentina. As the petrologic model does not require the intervention of old Precambrian continental crust, the nature of the basement on which thick accretionary turbiditic sequences were deposited remains a puzzling aspect. Discussion in this paper provides insights into the nature of magmatic source rocks and mechanisms of magma generation in Cordilleran-type volcano-plutonic arcs of destructive plate margins. (C) 2010 Elsevier Ltd. All rights reserved.

Contrasting magmatic signatures in the Rairakhol and Koraput alkaline complexes, Eastern Ghats belt, India

The relation between alkaline magmatism and tectonism has been a contentious issue, particularly for the Precambrian continental regions. Alkaline complexes at the southwestern margin of Eastern Ghats belt, India, have been interpreted as rift-valley magmatism. However, those complexes occurring in granulite ensemble in the interior segments of the Eastern Ghats belt could not possibly be related to the rift-system, assumed for the western margin of the Eastern Ghats belt. Koraput complex was emplaced in a pull-apart structure, dominated by magmatic fabrics and geochemically similar to a fractionated alkaline complex, compatible with an alkalibasalt series. Rairakhol complex, on the other hand, shows dominantly solid-state deformation fabrics and geochemically similar to a fractionated calc-alkaline suite. Isotopic data for the Koraput complex indicate ca. 917 Ma alkaline magmatism from a depleted mantle source and postcrystalline thermal overprint at ca. 745 Ma, also recorded from sheared metapelitic country rocks. The calc-alkaline magmatism of the Rairakhol complex occurred around 938 Ma, from an enriched mantle source, closely following Grenvillian granulite facies imprint in the charnockitic country rocks.

Petrogenesis of Upper Paleozoic post-collisional peraluminous leucogranites, Sierra de Ancasti, northwest Argentina

The Santa Rosa and Sauce Guacho plutons are two post-collisional peraluminous Late Devonian to Early Carboniferous leucogranites that intruded the banded schists of the Ancasti Formation. The leucogranites are composed of microcline phenocrysts along with quartz, plagioclase, muscovite, biotite, ilmenite, tourmaline, apatite, monazite and zircon. Their geochemical composition is consistent with S-type granites and mineralogically they belong to MPG granites (muscovite-peraluminous granites). It is proposed that granite magma generation was related to shear zones that concentrated fluids in the metasedimentary crust during a collision or transcurrent tectonics. U-Pb analyses on monazite gave an age of 369.8 +/- 5.3 Ma, while Sm/Nd isotopic data yield epsilon(Nd(t)) values of -5.3 for Sauce Guacho and -5.7 for Santa Rosa indicating crustal provenance. Nd model ages between 1,544 and 1,571 Ma are within the range of magmatic rocks from the Lower Ordovician Famatinian Arc in the Central Sierras Pampeanas.

Thermochronology of central Ribeira Fold Belt, SE Brazil: Petrological and geochronological evidence for long-term high temperature maintenance during Western Gondwana amalgamation

The studied sector of the central Ribeira Fold Belt (SE Brazil) comprises metatexites, diatexites, charnockites and blastomylonites. This study integrates petrological and thermochronological data in order to constrain the thermotectonic and geodynamic evolution of this Neoproterozoic-Ordovician mobile belt during Western Gondwana amalgamation. New data indicate that after an earlier collision stage at similar to 610 Ma (zircon, U-Pb age), peak metamorphism and lower crust partial melting, coeval with the main regional high grade D(1) thrust deformation, occurred at 572-562 Ma (zircon, U-Pb ages). The overall average cooling rate was low (<5 degrees C/Ma) from 750 to 250 degrees C (at similar to 455 Ma; biotite-WR Rb-Sr age), but disparate cooling paths indicate differential uplift between distinct lithotypes: (a) metatexites and blastomylonites show a overall stable 3-5 degrees C/Ma cooling rate; (b) charnockites and associated rocks remained at T>650 degrees C during sub-horizontal D(2) shearing until similar to 510-470 Ma (garnet-WR Sm-Nd ages) (1-2 degrees C/Ma), being then rapidly exhumed/cooled (8-30 degrees C/Ma) during post-orogenic D(3) deformation with late granite emplacement at similar to 490 Ma (zircon, U-Pb age). Cooling rates based on garnet-biotite Fe-Mg diffusion are broadly consistent with the geochronological cooling rates: (a) metatexites were cooled faster at high temperatures (6 degrees C/Ma) and slowly at low temperatures (0.1 degrees C/Ma), decreasing cooling rates with time; (b) charnockites show low cooling rates (2 degrees C/Ma) near metamorphic peak conditions and high cooling rates (120 degrees C/Ma) at lower temperatures, increasing cooling rates during retrogression. The charnockite thermal evolution and the extensive production of granitoid melts in the area imply that high geothermal gradients were sustained fora long period of time (50-90 Ma). This thermal anomaly most likely reflects upwelling of asthenospheric mantle and magma underplating coupled with long-term generation of high HPE (heat producing elements) granitoids. These factors must have sustained elevated crustal geotherms for similar to 100 Ma, promoting widespread charnockite generation at middle to lower crustal levels. (C) 2010 Elsevier B.V. All rights reserved.

Sr and Nd isotope composition of the metamorphic, sedimentary and ultramafic xenoliths of Lanzarote (Canary Islands): Implications for magma sources

The lavas produced by the Timanfaya eruption of 1730-1736 (Lanzarote, Canary Islands) contain a great many sedimentary and metamorphic (metasedimentary), and mafic and ultramafic plutonic xenoliths. Among the metamorphosed carbonate rocks (calc-silicate rocks [CSRs]) are monomineral rocks with forsterite or wollastonite, as well as rocks containing olivine +/- orthopyroxene +/- clinopyroxene +/- plagioclase: their mineralogical compositions are identical to those of the mafic (gabbros) and ultramafic (dunite, wherlite and lherzolite) xenoliths. The (87)Sr/(16)Sr (around 0.703) and (143)Nd/(144)Nd (around 0.512) isotope ratios of the ultramafic and metasedimentary xenoliths are similar, while the (147)Sm/(144)Nd ratios show crustal values (0.13-0.16) in the ultramafic xenoliths and mantle values (0.18-0.25) in some CSRs. The apparent isotopic anomaly of the metamorphic xenoliths can be explained in terms of the heat source (basaltic intrusion) inducing strong isotopic exchange ((87)Sr/(86)Sr and (143)Nd/(144)Nd) between metasedimentary and basaltic rocks. Petrofabric analysis also showed a possible relationship between the ultramafic and metamorphic xenoliths. (C) 2009 Elsevier B.V. All rights reserved.

Relics of eclogite facies assemblages in the Ceara Central Domain, NW Borborema Province, NE Brazil: Implications for the assembly of West Gondwana

The Borborema Province, in the NE of Brazil, is a rather complex piece in the Brazil-Africa puzzle as it represents the junction of the Dahomeyide/Pharusian, Central African, Aracuai and Brasilia fold belts located between the West-African/Sao Luis, Congo/Sao Francisco and Amazonas craton. The correlation between the Dahomeyides from W-Africa (Ghana, Benin, Togo, and Mali) and the Borborema Province involves the Medio Coreau and Central Ceara domains. The inferred continuation of the main oceanic suture zone exposed in the Dahomeyides of W Africa is buried beneath the Phanerozoic Parnaiba Basin in Brazil (northwest of the Medio Coreau domain) where some high density gravity anomalies may represent hidden remnants of an oceanic suture. In addition to this major suture a narrow, nearly continuous strip composed of mainly mafic pods containing relics of eclogite-facies assemblages associated with partially migmatized granulite-facies metapelitic gneisses has been found further east in the NW Borborema Province. These high pressure mafic rocks, interpreted as retrograded eclogites, are located between the Transbrasiliano Lineament and the Santa Quiteria continental arc and comprise primitive to evolved arc-related rocks with either arc- or MORB-type imprints that can indicate either deep subduction of oceanic lithosphere or roots of continental and oceanic magmatic arcs. Average peak P-T conditions under eclogite-facies metamorphism (T=770 degrees C and P = 17.3 kbar) were estimated using garnet-clinopyroxene thermometry and Jd content in clinopyroxene. Transition to granulite-facies conditions, as well as later widespread re-equilibration under amphibolite facies, were registered both in the basic and the metapelitic rocks and suggest a clockwise P-T path characterized by an increase in temperature followed by strong decompression. A phenomenon possibly related to the exhumation of a highly thickened crust associated with the suturing of the Medio Coreau and Central Ceara domains, two distinct crustal blocks separated by the Transbrasiliano Lineament. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Metamorphic P-T evolution of granulites in the central Ribeira Fold Belt, SE Brazil

Pseudosections, geothermobarometric estimates and careful petrographic observations of gneissic migmatites and granulites from Neoproterozoic central Ribeira Fold Belt (SE Brazil) were performed in order to quantify the metamorphic P-T conditions during prograde and retrograde evolution of the Brasiliano Orogeny. Results establish a prograde metamorphic trajectory from amphibolite facies conditions to metamorphic peak (T = 850 +/- 50 A degrees C; P = 8 +/- 1 kbar) that promoted widespread dehydrationmelting of 30 to 40% of the gneisses and high-grade granitization. After the metamorphic peak, migmatites evolved with cooling and decompression to T a parts per thousand 500 A degrees C and P a parts per thousand 5 kbar coupled with aH2O increase, replacing the high-grade paragenesis plagioclase-quartz-K-feldspar-garnet by quartz-biotite-sillimanite-(muscovite). Cordierite absence, microtextural observations and P-T results constrain the migmatite metamorphic evolution in the pseudosections as a clockwise P-T path with retrograde cooling and decompression. High-temperature conditions further dehydrated the lower crust with biotite and amphibole-dehydration melting and granulite formation coupled with 10% melt generation. Granulites can thus be envisaged as middle to lower crust dehydrated restites. Granulites were slowly (nearly isobarically) cooled, followed by late exhumation/retrograde rapid decompression and cooling, reflecting a two step P-T path. This retrograde evolution, coupled with water influx, chemically reequilibrated the rocks from granulite to amphibolite/greenschist facies, promoting the replacement of the plagioclase-quartz-garnet-hypersthene peak assemblage by quartz-biotite- K-feldspar symplectites.

GEOCHEMISTRY AND U-PB ZIRCON AGE OF THE INTERNAL OSSA-MORENA ZONE OPHIOLITE SEQUENCES: A REMNANT OF RHEIC OCEAN IN SW IBERIA

The Early Paleozoic geodynamic evolution in SW Iberia is believed to have been dominated by the opening of the Rheic Ocean. The Rheic Ocean is generally accepted to have resulted from the drift of peri-Gondwanan terranes such as Avalonia from the northern margin of Gondwana during Late Cambrian-Early Ordovician times. The closure of the Rheic Ocean was the final result of a continent-continent collision between Gondwana and Laurussia that produced the Variscan orogen. The Ossa-Morena Zone is a peri-Gondwana terrane, which preserves spread fragments of ophiolites - the Internal Ossa-Morena Zones Ophiolite Sequences (IOMZOS). The final patchwork of the IOMZOS shows a complete oceanic lithospheric sequence with geochemical characteristics similar to the ocean-floor basalts, without any orogenic fingerprint and/or crustal contamination. The IOMZOS were obducted and imbricated with high pressure lithologies. Based on structural, petrological and whole-rock geochemical data, the authors argue that the IOMZOS represent fragments of the oceanic lithosphere from the Rheic Ocean. Zircon SHRIMP U-Pb geochronological data on metagabbros point to an age of ca. 480 Ma for IOMZOS, providing evidence of a well-developed ocean in SW Iberia during this period, reinforcing the interpretation of the Rheic Ocean as a wide ocean among the peri-Gondwanan terranes during Early Ordovician times.

Age and geochemistry of mantle peridotites and diorite dykes from the Baldissero body: Insights into the Paleozoic-Mesozoic evolution of the Southern Alps

Trace element and isotopic data obtained for mantle spinel Iherzolites and diorite dykes from the Baldissero massif (Ivrea-Verbano Zone, Western Italy) provide new, valuable constraints on the petrologic and geodynamic evolution of the Southern Alps in Paleozoic to Mesozoic times. Whole rock and mineral chemistry indicates that Baldissero Iherzolites can be regarded as refractory mantle residues following limited melt extraction. In particular, the Light Rare Earth Elements (LREE)-depleted and fractionated compositions of whole rock and clinopyroxene closely match modelling results for refractory residues after low degrees (similar to 4-5%) of near-fractional melting of depleted mantle, possibly under garnet-facies conditions. Following this, the peridotite sequence experienced subsolidus re-equilibration at lithospheric spinel-facies conditions and intrusion of several generations of dykes. However, Iherzolites far from dykes show very modest metasomatic changes, as evidenced by the crystallisation of accessory titanian pargasite and the occurrence of very slight enrichments in highly incompatible trace elements (e.g. Nb). The Re-Os data for Iherzolites far from the dykes yield a 376 Ma (Upper Devonian) model age that is considered to record a partial melting event related to the Variscan orogenic cycle s.l. Dioritic dykes cutting the mantle sequence have whole rock, clinopyroxene and plagioclase characterised by high radiogenic Nd and low radiogenic Sr, which point to a depleted to slightly enriched mantle source. Whole rock and mafic phases of diorites have high Mg# values that positively correlate with the incompatible trace element concentrations. The peridotite at the dyke contact is enriched in orthopyroxene, iron and incompatible trace elements with respect to the Iherzolites far from dykes. Numerical simulations indicate that the geochemical characteristics of the diorites can be explained by flow of a hydrous, silica-saturated melt accompanied by reaction with the ambient peridotite and fractional crystallisation. The composition of the more primitive melts calculated in equilibrium with the diorite minerals show tholeiitic to transitional affinity. Internal Sm-Nd, three-point isochrons obtained for two dykes suggest an Upper Triassic-Lower Jurassic emplacement age (from 204 31 to 198 29 Ma). Mesozoic igneous events are unknown in the southern Ivrea-Verbano Zone (IVZ), but the intrusion of hydrous melts, mostly silica-saturated, have been well documented in the Finero region, i.e. the northernmost part of IVZ and Triassic magmatism with calc-alkaline to shoshonitic affinity is abundant throughout the Central-Eastern Alps. The geochemical and chronological features of the Baldissero diorites shed new light on the geodynamic evolution of the Southern Alps before the opening of the Jurassic Tethys. (C) 2010 Elsevier B.V. All rights reserved.