Theoretical models of planetary system formation: mass vs. semi-major axis

Context. Planet formation models have been developed during the past years to try to reproduce what has been observed of both the solar system and the extrasolar planets. Some of these models have partially succeeded, but they focus on massive planets and, for the sake of simplicity, exclude planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of radial velocity, transit, and direct imaging surveys, seems to be even more pronounced for low-mass planets. These new observations require improving planet formation models, including new physics, and considering the formation of systems. Aims: In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc and show the effect, in terms of architecture and composition of this multiplicity. Methods: We used an N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. Results: We show that the masses and semi-major axes of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping. We find that the fraction of ejected planets increases from nearly 0 to 8% as we change the number of embryos we seed the system with from 2 to 20 planetary embryos. Moreover, our calculations show that, when considering planets more massive than ~5 M⊕, simulations with 10 or 20 planetary embryos statistically give the same results in terms of mass function and period distribution.

Core formation and mantle differentiation on Mars

Geochemical investigation of Martian meteorites (SNC meteorites) yields important constraints on the chemical and geodynamical evolution of Mars. These samples may not be representative of the whole of Mars; however, they provide constraints on the early differentiation processes on Mars. The bulk composition of Martian samples implies the presence of a metallic core that formed concurrently as the planet accreted. The strong depletion of highly siderophile elements in the Martian mantle is only possible if Mars had a large scale magma ocean early in its history allowing efficient separation of a metallic melt from molten silicate. The solidification of the magma ocean created chemical heterogeneities whose ancient origin is manifested in the heterogeneous 142Nd and 182W abundances observed in different meteorite groups derived from Mars. The isotope anomalies measured in SNC meteorites imply major chemical fractionation within the Martian mantle during the life time of the short-lived isotopes 146Sm and 182Hf. The Hf-W data are consistent with very rapid accretion of Mars within a few million years or, alternatively, a more protracted accretion history involving several large impacts and incomplete metal-silicate equilibration during core formation. In contrast to Earth early-formed chemical heterogeneities are still preserved on Mars, albeit slightly modified by mixing processes. The preservation of such ancient chemical differences is only possible if Mars did not undergo efficient whole mantle convection or vigorous plate tectonic style processes after the first few tens of millions of years of its history.

Petrogenesis of rift-related tephrites, phonolites and trachytes (Central European Volcanic Province, Rhon, FRG): Constraints from Sr, Nd, Pb and O isotopes

The volcanic rocks of the Rhön area (Central European Volcanic Province, Germany) belong to a moderately alkali basaltic suite that is associated with minor tephriphonolites, phonotephrites, tephrites, phonolites and trachytes. Based on isotope sytematics (87Sr/86Sr: 0.7033–0.7042; 143Nd/144Nd: 0.51279–0.51287; 206Pb/204Pb: 19.1–19.5), the inferred parental magmas formed by variable degrees of partial melting of a common asthenospheric mantle source (EAR: European Asthenospheric Reservoir of Cebriá and Wilson, 1995). Tephrites, tephriphonolites, phonotephrites, phonolites and trachytes show depletions and enrichments in some trace elements (Sr, Ba, Nb, Zr, Y) indicating that they were generated by broadly similar differentiation processes that were dominated by fractionation of olivine, clinopyroxene, amphibole, apatite and titaniferous magnetite ± plagioclase ± alkalifeldspar. The fractionated samples seem to have evolved by two distinct processes. One is characterized by pure fractional crystallization indicated by increasing Nb (and other incompatible trace element) concentrations at virtually constant 143Nd/144Nd ~ 0.51280 and 87Sr/86Sr ~ 0.7035. The other process involved an assimilation–fractional crystallization (AFC) process where moderate assimilation to crystallization rates produced evolved magmas characterized by higher Nb concentrations at slightly lower 143Nd/144Nd down to 0.51275. Literature data for some of the evolved rocks show more variable 87Sr/86Sr ranging from 0.7037 to 0.7089 at constant 143Nd/144Nd ~ 0.51280. These features may result from assimilation of upper crustal rocks by highly differentiated low-Sr (< 100 ppm Sr) lavas. However, based on the displacement of the differentiated rocks from this study towards lower 143Nd/144Nd ratios and modeled AFC processes in 143Nd/144Nd vs. 87Sr/86Sr and 207Pb/204Pb vs. 143Nd/144Nd space assimilation of lower crustal rocks seems more likely. The view that assimilation of lower crustal rocks played a role is confirmed by high-precision double-spike Pb isotope data that reveal higher 207Pb/204Pb ratios (15.62–15.63) in the differentiated rocks than in the primitive basanites (15.58–15.61). This is compatible with incorporation of radiogenic Pb from lower crustal xenoliths (207Pb/204Pb: 15.63–15.69) into the melt. However, 206Pb/204Pb ratios are similar for the differentiated rocks (19.13–19.35) and the primitive basanites (19.12–19.55) implying that assimilation involved an ancient crustal end member with a higher U/Pb ratio than the mantle source of the basanites. In addition, alteration-corrected δ18O values of the differentiated rocks range from c. 5 to 7‰ which is the same range as observed in the primitive alkaline rocks. This study confirms previous interpretations that highlighted the role of AFC processes in the evolution of alkaline volcanic rocks in the Rhön area of the Central European Volcanic Province.

Planet formation models: the interplay with the planetesimal disc

Context. According to the sequential accretion model (or core-nucleated accretion model), giant planet formation is based first on the formation of a solid core which, when massive enough, can gravitationally bind gas from the nebula to form the envelope. The most critical part of the model is the formation time of the core: to trigger the accretion of gas, the core has to grow up to several Earth masses before the gas component of the protoplanetary disc dissipates. Aims: We calculate planetary formation models including a detailed description of the dynamics of the planetesimal disc, taking into account both gas drag and excitation of forming planets. Methods: We computed the formation of planets, considering the oligarchic regime for the growth of the solid core. Embryos growing in the disc stir their neighbour planetesimals, exciting their relative velocities, which makes accretion more difficult. Here we introduce a more realistic treatment for the evolution of planetesimals' relative velocities, which directly impact on the formation timescale. For this, we computed the excitation state of planetesimals, as a result of stirring by forming planets, and gas-solid interactions. Results: We find that the formation of giant planets is favoured by the accretion of small planetesimals, as their random velocities are more easily damped by the gas drag of the nebula. Moreover, the capture radius of a protoplanet with a (tiny) envelope is also larger for small planetesimals. However, planets migrate as a result of disc-planet angular momentum exchange, with important consequences for their survival: due to the slow growth of a protoplanet in the oligarchic regime, rapid inward type I migration has important implications on intermediate-mass planets that have not yet started their runaway accretion phase of gas. Most of these planets are lost in the central star. Surviving planets have masses either below 10 M⊕ or above several Jupiter masses. Conclusions: To form giant planets before the dissipation of the disc, small planetesimals (~0.1 km) have to be the major contributors of the solid accretion process. However, the combination of oligarchic growth and fast inward migration leads to the absence of intermediate-mass planets. Other processes must therefore be at work to explain the population of extrasolar planets that are presently known.

Metamorphic and geochronogical study of the Triassic El Oro metamorphic complex, Ecuador: Implications for high-temperature metamorphism in a forearc zone

In the forearc of the Andean active margin in southwest Ecuador, the El Oro metamorphic complex exhibits a well exposed tilted forearc section partially migmatized. We used Raman spectroscopy on carbonaceous matter (RSCM) thermometry and pseudosections coupled with mineralogical and textural studies to constrain the pressure–temperature (P–T) evolution of the El Oro metamorphic complex during Triassic times. Our results show that anatexis of the continental crust occurred by white-mica and biotite dehydration melting along a 10 km thick crustal domain (from 4.5 to 8 kbar) with increasing temperature from 650 to 700 °C. In the biotite dehydration melting zone, temperature was buffered at 750–820 °C in a 5 km thick layer. The estimated average thermal gradient during peak metamorphism is of 30 °C/km within the migmatitic domain can be partitioned into two apparent gradients parts. The upper part from surface to 7 km depth records a 40–45 °C/km gradient. The lower part records a quasi-adiabatic geotherm with a 10 °C/km gradient consistent with an isothermal melting zone. Migmatites U–Th–Pb geochronology yielded zircon and monazite ages of 229.3 ± 2.1 Ma and 224.5 ± 2.3 Ma, respectively. This thermal event generated S-type magmatism (the Marcabeli granitoid) and was immediately followed by underplating of the high-pressure low-temperature (HP-LT) Arenillas–Panupalí unit at 225.8 ± 1.8 Ma. The association of high-temperature low-pressure (HT-LP) migmatites with HP-LT unit constitutes a new example of a paired metamorphic belt along the South American margin. We propose that in addition to crustal thinning, underplating of the Piedras gabbroic unit before 230 Ma provided the heat source necessary to foster crustal anatexis. Furthermore, its MORB signature shows that the asthenosphere was involved as the source of the heat anomaly. S-type felsic magmatism is widespread during this time and suggests that a large-scale thermal anomaly affected a large part of the South American margin during the late Triassic. We propose that crustal anatexis is related to an anomaly that arose during subduction of the Panthalassa ocean under the South American margin. Slab verticalization or slab break-off can be invoked as the origin of the upwelling of the asthenosphere.

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

Paleozoic-aged brine and authigenic helium preserved in an Ordovician shale aquiclude

Consideration of the geosphere for isolation of nuclear waste has generated substantial interest in the origin, age, and movement of fl uids and gases in low-permeability rock formations. Here, we present profi les of isotopes, solutes, and helium in porewaters recovered from 860 m of Cambrian to Devonian strata on the eastern fl ank of the Michigan Basin. Of particular interest is a 240-m-thick, halite-mineralized, Ordovician shale and carbonate aquiclude, which hosts Br–-enriched, post-dolomitic brine (5.8 molal Cl) originating as evaporated Silurian seawater. Authigenic helium that has been accumulating in the aquiclude for more than 260 m.y. is found to be isolated from underlying allochthonous, 3He-enriched helium that originated from the rifted base of the Michigan Basin and the Canadian Shield. The Paleozoic age and immobility of the pore fl uids in this Ordovician aquiclude considerably strengthen the safety case for deep geological repositories, but also provide new insights into the origin of deep crustal brines and opportunities for research on other components of a preserved Paleozoic porewater system.

Neutron capture on Pt isotopes in iron meteorites and the Hf–W chronology of core formation in planetesimals

The short-lived 182Hf–182W isotope system can provide powerful constraints on the timescales of planetary core formation, but its application to iron meteorites is hampered by neutron capture reactions on W isotopes resulting from exposure to galactic cosmic rays. Here we show that Pt isotopes in magmatic iron meteorites are also affected by capture of (epi)thermal neutrons and that the Pt isotope variations are correlated with variations in 182W/184W. This makes Pt isotopes a sensitive neutron dosimeter for correcting cosmic ray-induced W isotope shifts. The pre-exposure 182W/184W derived from the Pt–W isotope correlations of the IID, IVA and IVB iron meteorites are higher than most previous estimates and are more radiogenic than the initial 182W/184W of Ca–Al-rich inclusions (CAI). The Hf–W model ages for core formation range from +1.6±1.0 million years (Ma; for the IVA irons) to +2.7±1.3 Ma after CAI formation (for the IID irons), indicating that there was a time gap of at least ∼1 Ma between CAI formation and metal segregation in the parent bodies of some iron meteorites. From the Hf–W ages a time limit of <1.5–2 Ma after CAI formation can be inferred for the accretion of the IID, IVA and IVB iron meteorite parent bodies, consistent with earlier conclusions that the accretion of differentiated planetesimals predated that of most chondrite parent bodies.

Ar-40/Ar-39 Evidence for Middle Proterozoic (1300-1500 Ma) Slow Cooling of the Southern Black Hills, South Dakota, Midcontinent, North America: Implications for Early Proterozoic P-T Evolution and Posttectonic Magmatism

Ar-40/Ar-39 total gas and plateau dates from muscovite and biotite in the southern Black Hills, South Dakota, provide evidence for a period of Middle Proterozoic slow cooling. Early Proterozoic (1600-1650 Ma) mica dates were obtained from metasedimentary rocks located in a synformal structure between the Harney Peak and Bear Mountain domes and also south of Bear Mountain. Metamorphic rocks from the dome areas and undeformed samples of the similar to 1710 Ma Harney Peak Granite (HPG) yield Middle Proterozoic mica dates (similar to 1270-1500 Ma). Two samples collected between the synform and Bear Mountain dome yield intermediate total gas mica dates of similar to 1550 Ma. We suggest two end-member interpretations to explain the map pattern of cooling ages: (1) subhorizontal slow cooling of an area which exhibits variation in mica Ar retention intervals or (2) mild folding of a Middle Proterozoic (similar to 1500 Ma) similar to 300 degrees C isotherm. According to the second interpretation, the preservation of older dates between the domes may reflect reactivation of a preexisting synformal structure (and downwarping of relatively cold rocks) during a period of approximately east-west contraction and slow uplift during the Middle Proterozoic. The mica data, together with hornblende data from the Black Hills published elsewhere, indicate that the ambient country-rock temperature at the 3-4 kbar depth of emplacement of the HPG was between 350 degrees C and 500 degrees C, suggesting that the average upper crustal geothermal gradient was 25 degrees-40 degrees C/km prior to intrusion. The thermochronologic data suggest HPG emplacement was followed by a similar to 200 m.y. period of stability and tectonic quiescence with little uplift. We propose that crust thickened during the Early Proterozoic was uplifted and erosionally(?) thinned prior to similar to 1710 Ma and that the HPG magma was emplaced into isostatically stable crust of relatively normal thickness. We speculate that uplift and crustal thinning prior to HPG intrusion was the result of differential thinning of the subcrustal lithosphere beneath the Black Hills. If so, this process would have also caused an increase in mantle heat flux across the Moho and triggered vapor-absent melting of biotite to produce the HPG magma. This scenario for posttectonic granite generation is supported, in part, by the fact that in the whole of the Black Hills, the HPG is spatially associated with the deepest exposed Early Proterozoic country rock.

An Ice-Core Proxy for Antarctic Circumpolar Zonal Wind Intensity

Using US National Centers for Environmental Prediction/US National Center for Atmospheric Research re-analysis data, we investigate the relationships between crustal ion (nssCa(2+)) concentrations from three West Antarctic ice cores, namely, Siple Dome (SD), ITASE00-1 (IT001) and ITASE01-5 (IT015), and primary components of the climate system, namely, air pressure/geopotential height, zonal (u) and meridional (v) wind strength. Linear correlation analyses between nssCa(2+) concentrations and both air-pressure and wind fields for the period of overlap between records indicate that the SD nssCa(2+) variation is positively correlated with spring circumpolar zonal wind, while IT001 nssCa(2+) has a positive correlation with circumpolar zonal wind throughout the year (r > 0.3, p < 0.01). Intensified Southern Westerlies circulation is conducive to transport of more crustal aerosols to both sites. Further correlation analyses between nssCa(2+) concentrations from SD and IT001 and atmospheric circulation suggest that the high inland plateau (represented by core IT001) is largely influenced by transport from the upper troposphere. IT015 nssCa(2+) is negatively correlated with westerly wind in October and November, suggesting that stronger westerly circulation may weaken the transport of crustal species to IT015. Correlations of nssCa(2+) from the three ice cores with the Antarctic Oscillation index are consistent with results developed from the wind-field investigation. In addition, calibration between nssCa(2+) concentration and the multivariate El Nino-Southern Oscillation (ENSO) index shows that crustal species transport to IT001 is enhanced during strong ENSO events.

Thermal and Barometric Constraints on the Intrusive and Unroofing History of the Black Mountains: Implications for Timing, Initial Dip, and Kinematics of Detachment Faulting in the Death-Valley Region, California

Unroofing of the Black Mountains, Death Valley, California, has resulted in the exposure of 1.7 Ga crystalline basement, late Precambrian amphibolite facies metasedimentary rocks, and a Tertiary magmatic complex. The Ar-40/Ar-39 cooling ages, obtained from samples collected across the entire length of the range (>55 km), combined with geobarometric results from synextensional intrusions, provide time-depth constraints on the Miocene intrusive history and extensional unroofing of the Black Mountains. Data from the southeastern Black Mountains and adjacent Greenwater Range suggest unroofing from shallow depths between 9 and 10 Ma. To the northwest in the crystalline core of the range, biotite plateau ages from approximately 13 to 6.8 Ma from rocks making up the Death Valley turtlebacks indicate a midcrustal residence (with temperatures >300-degrees-C) prior to extensional unroofing. Biotite Ar-40/Ar-39 ages from both Precambrian basement and Tertiary plutons reveal a diachronous cooling pattern of decreasing ages toward the northwest, subparallel to the regional extension direction. Diachronous cooling was accompanied by dike intrusion which also decreases in age toward the northwest. The cooling age pattern and geobarometric constraints in crystalline rocks of the Black Mountains suggest denudation of 10-15 km along a northwest directed detachment system, consistent with regional reconstructions of Tertiary extension and with unroofing of a northwest deepening crustal section. Mica cooling ages that deviate from the northwest younging trend are consistent with northwestward transport of rocks initially at shallower crustal levels onto deeper levels along splays of the detachment. The well-known Amargosa chaos and perhaps the Badwater turtleback are examples of this "splaying" process. Considering the current distance of the structurally deepest samples away from moderately to steeply east tilted Tertiary strata in the southeastern Black Mountains, these data indicate an average initial dip of the detachment system of the order of 20-degrees, similar to that determined for detachment faults in west central Arizona and southeastern California. Beginning with an initially listric geometry, a pattern of footwall unroofing accompanied by dike intrusion progress northwestward. This pattern may be explained by a model where migration of footwall flexures occur below a scoop-shaped banging wall block. One consequence of this model is that gently dipping ductile fabrics developed in the middle crust steepen in the upper crust during unloading. This process resolves the low initial dips obtained here with mapping which suggests transport of the upper plate on moderately to steeply dipping surfaces in the middle and upper crust.

Landscape Evolution of the Dry Valleys, Transantarctic Mountains: Tectonic Implications

There are different views about the amount and timing of surface uplift in the Transantarctic Mountains and the geophysical mechanisms involved. Our new interpretation of the landscape evolution and tectonic history of the Dry Valleys area of the Transantarctic Mountains is based on geomorphic mapping of an area of 10,000 km(2). The landforms are dated mainly by their association with volcanic ashes and glaciomarine deposits and this permits a reconstruction of the stages and timing of landscape evolution. Following a lowering of base level about 55 m.y. ago, there was a phase of rapid denudation associated with planation and escarpment retreat, probably under semiarid conditions. Eventually, downcutting by rivers, aided in places by glaciers, graded valleys to near present sea level. The main valleys were flooded by the sea in the Miocene during a phase of subsidence before experiencing a final stage of modest upwarping near the coast. There has been remarkably little landform change under the stable, cold, polar conditions of the last 15 m.y. It is difficult to explain the Sirius Group deposits, which occur at high elevations in the area, if they are Pliocene in age. Overall, denudation may have removed a wedge of rock with a thickness of over 4 km at the coast declining to 1 km at a point 75 km inland, which is in good agreement with the results of existing apatite fission track analyses. It is suggested that denudation reflects the differences in base level caused by high elevation at the time of extension due to underplating and the subsequent role of thermal uplift and flexural isostasy. Most crustal uplift (2-4 km) is inferred to have occurred in the early Cenozoic with 400 m of subsidence in the Miocene followed by 300 m of uplift in the Pliocene.

Chemical Composition of Fresh Snow from Glaciar Marinelli, Tierra Del Fuego, Chile

A fresh-snow sampling campaign was conducted during the late austral summer of 2006 in the accumulation zone of Glaciar Marinelli, located in the Cordillera Darwin, Tierra del Fuego, Chile. Snow samples were analyzed for stable isotopes (delta(18)O, major soluble ions (Na', K', Ca, Mg, a NO(3)(-), SO(4)(2-), MS(-)) and major and trace elements (Na, Mg, Al, Fe, Ca, Sr, Cd, Cs, Ba, La, Ce, Pr, Dy, Ho, Er, Bi, U, As, Ti, V, Cr, Mn). The dominance of marine chemistry resembles that in studies from Patagonian glaciers. Snow chemistry was dominantly loaded by marine species (Cl(-), Na(+) and ssSO(4)(2-)), while contributions of crustal species (e.g. Al and Fe) were very low. Empirical orthogonal function analysis suggests two possible dust sources, one represented by Al and Fe and the other by La, Ce and Pr. Enrichment-factor calculations suggest the majority of elements are within average upper-crustal ratios, but major enrichments of Bi and Cd (hundreds of times) suggest possible anthropogenic sources. Linear correlation of delta(18)O and barometric pressure (r = 0.60, p < 0.007) suggests a potential 'amount effect' relationship between depleted delta(18)O ratios and stronger storm conditions (e.g. greater precipitation). The snow-chemistry records from Glaciar Marinelli are the first measured in Tierra del Fuego, the southernmost glaciated region outside Antarctica.

Tertiary Minette and Melanephelinite Dikes, Wasatch Plateau, Utah - Records of Mantle Heterogeneities and Changing Tectonics

A swarm of minette and melanephelinite dikes is exposed over 2500 km2 in and near the Wasatch Plateau, central Utah, along the western margin of the Colorado Plateaus in the transition zone with the Basin and Range province. To date, 110 vertical dikes in 25 dike sets have been recognized. Strikes shift from about N80-degrees-W for 24 Ma dikes, to about N60-degrees-W for 18 Ma, to due north for 8-7 m.y. These orientations are consistent with a shift from east-west Oligocene compression associated with subduction to east-west late Miocene crustal extension. Minettes are the most common rock type; mica-rich minette and mica-bearing melanephelinite occurs in 24 Ma dikes, whereas more ordinary minette is found in 8-7 Ma dikes. One melanephelinite dike is 18 Ma. These mafic alkaline rocks are transitional to one another in modal and major element composition but have distinctive trace element patterns and isotopic compositions; they appear to have crystallized from primitive magmas. Major, trace element, and Nd-Sr isotopic data indicate that melanephelinite, which has similarities to ocean island basalt, was derived from small degree melts of mantle with a chondritic Sm/Nd ratio probably located in the asthenosphere, but it is difficult to rule out a lithospheric source. In contrast, mica-bearing rocks (mica melanephelinite and both types of minette) are more potassic and have trace element patterns with strong Nb-Ta depletions and Sr-Nd isotopic compositions caused by involvement with a component from heterogeneously enriched lithospheric mantle with long-term enrichment of Rb or light rare earth elements (REE) (epsilon Nd as low as - 15 in minette). Light REE enrichment must have occurred anciently in the mid-Proterozoic when the lithosphere was formed and is not a result of Cenozoic subduction processes. After about 25 Ma, foundering of the subducting Farallon plate may have triggered upwelling of warm asthenospheric mantle to the base of the lithosphere. Melanephelinite magma may have separated from the asthenosphere and, while rising through the lithosphere, provided heat for lithospheric magma generation. Varying degrees of interaction between melanephelinite and small potassic melt fractions derived from the lithospheric mantle can explain the gradational character of the melanephelinite to minette suite.

Continental degassing of 4He by surficial discharge of deep groundwater

Radiogenic He is produced by the decay of uranium and thorium in the Earth’s mantle and crust. From here, it is degassed to the atmosphere and eventually escapes to space. Assuming that all of the 4He produced is degassed, about 70% of the total He degassed from Earth comes from the continental crust. However, the outgoing flux of crustal He has not been directly measured at the Earth’s surface and the migration pathways are poorly understood. Here we present measurements of helium isotopes and the long-lived cosmogenic radio-isotope Kr in the deep, continental-scale Guarani aquifer in Brazil and show that crustal He reaches the atmosphere primarily by the surficial discharge of deep groundwater. We estimate that He in Guarani groundwater discharge accounts for about 20% of the assumed global flux from continental crust, and that other large aquifers may account for about 33%. Old groundwater ages suggest that He in the Guarani aquifer accumulates over half- to one-millionyear timescales. We conclude that He degassing from the continents is regulated by groundwater discharge, rather than episodic tectonic events, and suggest that the assumed steady state between crustal production and degassing of He, and its resulting atmospheric residence time, should be re-examined.