Highly reducing conditions during core formation on Mercury: Implications for internal structure and the origin of a magnetic field

The high average density and low surface FeO content of the planet Mercury are shown to be consistent with very low oxygen fugacity during core segregation, in the range 3-6 log units below the iron-wüstite buffer. These low oxygen fugacities, and associated high metal content, are characteristic of high-iron enstatite (EH) and Bencubbinite (CB) chondrites, raising the possibility that such materials may have been important building blocks for this planet. With this idea in mind we have explored the internal structure of a Mercury sized planet of EH or CB bulk composition. Phase equilibria in the silicate mantle have been modeled using the thermodynamic calculator p-MELTS, and these simulations suggest that orthopyroxene will be the dominant mantle phase for both EH and CB compositions, with crystalline SiO₂ being an important minor phase at all pressures. Simulations for both compositions predict a plagioclase-bearing “crust” at low pressure, significant clinopyroxene also being calculated for the CB bulk composition. Concerning the core, comparison with recent high pressure and high temperature experiments relevant to the formation of enstatite meteorites, suggest that the core of Mercury may contain several wt.% silicon, in addition to sulfur. In light of the pressure of the core-mantle boundary on Mercury (∼7 GPa) and the pressure at which the immiscibility gap in the system Fe-S-Si closes (∼15 GPa) we suggest that Mercury’s core may have a complex shell structure comprising: (i) an outer layer of Fe-S liquid, poor in Si; (ii) a middle layer of Fe-Si liquid, poor in S; and (iii) an inner core of solid metal. The distribution of heat-producing elements between mantle and core, and within a layered core have been quantified. Available data for Th and K suggest that these elements will not enter the core in significant amounts. On the other hand, for the case of U both recently published metal/silicate partitioning data, as well as observations of U distribution in enstatite chondrites, suggest that this element behaves as a chalcophile element at low oxygen fugacity. Using these new data we predict that U will be concentrated in the outer layer of the mercurian core. Heat from the decay of U could thus act to maintain this part of Mercury’s core molten, potentially contributing to the origin of Mercury’s magnetic field. This result contrasts with the Earth where the radioactive decay of U represents a negligible contribution to core heating.     

Planetesimal formation by turbulent concentration

The formation of 1-1000 km diameter planetesimals from dust grains in a protoplanetary disk is a key step in planet formation. Conventional models for planetesimal formation involve pairwise sticking of dust grains, or the sedimentation of dust grains to a thin layer at the disk midplane followed by gravitational instability. Each of these mechanisms is likely to be frustrated if the disk is turbulent. Particles with stopping times comparable to the turnover time of the smallest eddies in a turbulent disk can become concentrated into dense clumps that may be the precursors of planetesimals. Such particles are roughly millimeter-sized for a typical protoplanetary disk. To survive to become planetesimals, clumps need to form in regions of low vorticity to avoid rotational breakup. In addition, clumps must have sufficient self gravity to avoid break up due to the ram pressure of the surrounding gas. Given these constraints, the rate of planetesimal formation can be estimated using a cascade model for the distribution of particle concentration and vorticity within eddies of various sizes in a turbulent disk. We estimate planetesimal formation rates and planetesimal diameters as a function of distance from a star for a range of protoplanetary disk parameters. For material with a solar composition, the dust-to-gas ratio is too low to allow efficient planetesimal formation, and most solid material will remain in small particles. Enhancement of the dust-to-gas ratio by 1-2 orders of magnitude, either vertically or radially, allows most solid material to be converted into planetesimals within the typical lifetime of a disk. Such dust-to-gas ratios may occur near the disk midplane as a result of vertical settling of short-lived clumps prior to clump breakup. Planetesimal formation rates are sensitive to the assumed size and rotational speed of the largest eddies in the disk, and formation rates increase substantially if the largest eddies rotate more slowly than the disk itself. Planetesimal formation becomes more efficient with increasing distance from the star unless the disk surface density profile has a slope of −1.5 or steeper as a function of distance. Planetesimal formation rates typically increase by an order-of-magnitude or more moving outward across the snow line for a solid surface density increase of a factor of 2. In all cases considered, the modal planetesimal size increases with roughly the square root of distance from the star. Typical modal diameters are 100 km and 400 km in the regions corresponding to the asteroid belt and Kuiper belt in the Solar System, respectively.     

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

Metallicity of the massive protoplanets around HR 8799 If formed by gravitational instability

The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their ‘pre-collapse’ stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be and 5.67×10³ years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions.     

A submillimetre survey of the kinematics of the Perseus molecular cloud – I. Data

We present submillimetre observations of the   J = 3 → 2  rotational transition of ¹²CO, ¹³CO and C¹⁸O across over 600 arcmin² of the Perseus molecular cloud, undertaken with the Heterodyne Array Receiver Programme (HARP), a new array spectrograph on the James Clerk Maxwell Telescope. The data encompass four regions of the cloud, containing the largest clusters of dust continuum condensations: NGC 1333, IC348, L1448 and L1455. A new procedure to remove striping artefacts from the raw HARP data is introduced. We compare the maps to those of the dust continuum emission mapped with the Submillimetre Common-User Bolometer Array (SCUBA; Hatchell et al.) and the positions of starless and protostellar cores (Hatchell et al.). No straightforward correlation is found between the masses of each region derived from the HARP CO and SCUBA data, underlining the care that must be exercised when comparing masses of the same object derived from different tracers. From the ¹³CO/C¹⁸O line ratio the relative abundance of the two species  ([¹³CO]/[C¹⁸O]∼ 7)  and their opacities (typically τ is 0.02–0.22 and 0.15–1.52 for the C¹⁸O and ¹³CO gas, respectively) are calculated. C¹⁸O is optically thin nearly everywhere, increasing in opacity towards star-forming cores but not beyond  τ₁₈∼ 0.9  . Assuming the ¹²CO gas is optically thick, we compute its excitation temperature, T _ex (around 8–30 K), which has little correlation with estimates of the dust temperature.     

Risks of shifting baselines highlighted by anecdotal accounts of New Zealand's snapper (Pagrus auratus) fishery

Anecdotal data sources may constitute an important component of the information available about an exploited species, as record keeping may not have occurred until after exploitation began. Here, we aimed to fill any gaps in the exploitative history of the sparid snapper (Pagrus auratus), using social and historical research methods. Social research consisted of interviews with recreational fishers, focusing on the most and largest snapper they had caught. In addition, the diary‐logs of two recreational fishers were analysed. Historical research consisted of investigation of old books, photos, archives and unpublished sources unconventional to fishery science. Interviews with fishers demonstrated no or weak trends in snapper abundance or size, and were likely impeded by a lack of ability to detect change in a fish stock that may still be considered abundant. The fishers’ perception of change, however, largely reflected recent experiences (last c. 10 years), when biomass is understood to have increased, and mostly did not consider experiences before the 1980s. Alternatively, diary‐logs of fisher catch rates produced a pattern that matched formal stock assessments of snapper biomass, suggesting declines in abundance up until the 1990s and an increase in biomass after that time. Historical research, although more qualitative, had the ability to investigate periods where formal records were not kept and described a fishery vastly different from the current one. Snapper were easily caught, in great abundance and in unusual locations. Localised depletion of snapper was first noticed in the early 20th century, despite spectacular catches of snapper occurring after that time. Snapper behaviour was also likely different, with visual sightings of snapper by onlookers a common occurrence. Although predictions from stock assessment models are consistent with that of the anecdotes listed here (i.e., high biomass in the past), these anecdotes are valuable as they explain lost biomass in a perspective meaningful to all. This perspective may be valuable for managers trying to consider the non‐financial value of a shared fishery but, if unrecognised, represents a shifting baseline.     

鲁西新泰孟家屯2.7Ga变质沉积岩与黑云斜长片麻岩锆石Hf同位素特征

对出露于山东新泰孟家屯2.7Ga的孟家屯岩组中石榴石石英岩(M08)、含十字石石榴石黑云母石英片岩(D242-Y2)和黑云斜长片麻岩(M03)(TTG质片麻岩)进行锆石Lu-Hf同位素分析。石榴石石英岩锆石核部176Lu/177Hf、176Hf/177Hf变化范围为0.001730～0.002546、0.281249～0.281360,锆石变质边部176Lu/177Hf、176Hf/177Hf变化范围为0.000123～0.002070、0.281241～0.281318;含十字石石榴石黑云母石英片岩锆石核部176Lu/177Hf、176Hf/177Hf变化范围为0.001334～0.002169、0.281226～0.281324,锆石变质边部176Lu/177Hf、176Hf/177Hf变化范围为0.000445～0.001375、0.281227～0.281309;黑云斜长片麻岩锆石176Lu/177Hf、176Hf/177Hf变化范围为0.000315～0.000847、0.281186～0.281265。孟家屯岩组石榴石英岩、含十字石石榴石黑去母石英片岩中碎屑(岩浆)锆石和黑云斜长片麻岩中岩浆锆石的εHf(t)分别为3.75～7.26、2.31～7.26和3.21～6.27,同时TDM1与其U-Pb年龄非常接近,表明它们起源于新生地壳的部分熔融。结合区域研究资料认为,2.7Ga是鲁西重要的一期构造岩浆热事件,也可能是华北克拉通重要的构造岩浆活动时期。变质沉积岩中许多变质增生锆石相对于核部岩浆锆石低176Lu/177Hf、高176Hf/177Hf,显示变质作用过程中石榴石的存在对锆石的Lu-Hf同位素体系有很大影响。     

闽西北地区混合花岗岩类型划分及其与金成矿关系

根据闽西北混合花岗岩产出的大地构造位置、成岩物质来源及其与周围地质体的关系等研究,认为闽西北地区混合花岗岩总体上可分为区域型混合花岗岩和边缘型混合花岗岩2种类型。通过2类混合花岗岩的空间展布和特征对比,并对其形成机理及其与金成矿关系进行了初步探讨。