Geochemical studies of planetary accretion and evolution have invoked various degrees of collisional erosion to explain differences in bulk composition between planets and chondrites. Here we undertake a full, dynamical evaluation of 'crustal stripping' during accretion and its key geochemical consequences. We present smoothed particle hydrodynamics simulations of collisions between differentiated rocky planetesimals and planetary embryos. We find that the crust is preferentially lost relative to the mantle during impacts, and we have developed a scaling law that approximates the mass of crust that remains in the largest remnant. Using this scaling law and a recent set of N-body simulations, we have estimated the maximum effect of crustal stripping on incompatible element abundances during the accretion of planetary embryos. We find that on average one third of the initial crust is stripped from embryos as they accrete, which leads to a reduction of ~20% in the budgets of the heat producing elements if the stripped crust does not reaccrete. Erosion of crusts can lead to non-chondritic ratios of incompatible elements, but the magnitude of this effect depends sensitively on the details of the crust-forming melting process. The Lu/Hf system is fractionated for a wide range of crustal formation scenarios. Using eucrites (the products of planetesimal silicate melting, thought to represent the crust of Vesta) as a guide to the Lu/Hf of planetesimal crust partially lost during accretion, we predict the Earth could evolve to a superchondritic 176-Hf/177-Hf (3-5 parts per ten thousand) at present day. Such values are in keeping with compositional estimates of the bulk Earth. Stripping of planetary crusts during accretion can lead to detectable changes in bulk composition of lithophile elements, but the fractionation is relatively subtle, and sensitive to the efficiency of reaccretion.

The Earth appears non-chondritic in its abundances of refractory lithophile elements, posing a significant problem for our understanding of its formation and evolution. It has been suggested that this non-chondritic composition may be explained by collisional erosion of differentiated planetesimals of originally chondritic composition. In this work, we present N-body simulations of terrestrial planet formation that track the growth of planetary embryos from planetesimals. We simulate evolution through the runaway and oligarchic growth phases under the Grand Tack model and in the absence of giant planets. These simulations include a state-of-the-art collision model which allows multiple collision outcomes, such as accretion, erosion, and bouncing events, that enables tracking of the evolving core mass fraction of accreting planetesimals. We show that the embryos grown during this intermediate stage of planet formation exhibit a range of core mass fractions, and that with significant dynamical excitation, enough mantle can be stripped from growing embryos to account for the Earth's non-chondritic Fe/Mg ratio. We also find that there is a large diversity in the composition of remnant planetesimals, with both iron-rich and silicate-rich fragments produced via collisions.

Several lines of evidence indicate a non-chondritic composition for Bulk Earth. If Earth formed from the accretion of chondritic material, its non-chondritic composition, in particular the super-chondritic 142Nd/144Nd and low Mg/Fe ratios, might be explained by the collisional erosion of differentiated planetesimals during its formation. In this work we use an N-body code, that includes a state-of-the-art collision model, to follow the formation of protoplanets, similar to proto-Earth, from differentiated planetesimals (> 100 km) up to isolation mass (> 0.16 M_Earth). Collisions between differentiated bodies have the potential to change the core-mantle ratio of the accreted protoplanets. We show that sufficient mantle material can be stripped from the colliding bodies during runaway and oligarchic growth, such that the final protoplanets could have Mg/Fe and Si/Fe ratios similar to that of bulk Earth, but only if Earth is an extreme case and the core is assumed to contain 10% silicon by mass. This may indicate an important role for collisional differentiation during the giant impact phase if Earth formed from chondritic material.

Reverberation mapping (RM) is an important technique in studies of active galactic nuclei (AGN). The key idea of RM is to measure the time lag $\tau$ between variations in the continuum emission from the accretion disc and subsequent response of the broad line region (BLR). The measurement of $\tau$ is typically used to estimate the physical size of the BLR and is combined with other measurements to estimate the black hole mass $M_{\rm BH}$. A major difficulty with RM campaigns is the large amount of data needed to measure $\tau$. Recently, Fine et al (2012) introduced a new approach to RM where the BLR light curve is sparsely sampled, but this is counteracted by observing a large sample of AGN, rather than a single system. The results are combined to infer properties of the sample of AGN. In this letter we implement this method using a hierarchical Bayesian model and contrast this with the results from the previous stacked cross-correlation technique. We find that our inferences are more precise and allow for more straightforward interpretation than the stacked cross-correlation results.