Experimental partitioning of trace elements into schreibersite with applications to IIG iron meteorites

Some of the defining characteristics of the IIG iron meteorite group are their high bulk P contents and massive, coarse schreibersite, which have been calculated to make up roughly 11–14 wt% of each specimen. In this study, we produced two data sets to investigate the formation of schreibersites in IIG irons: measurements of trace elements in the IIG iron meteorite Twannberg and experimental determinations of trace element partitioning into schreibersite. The schreibersite‐bearing experiments were conducted with schreibersite in equilibrium with a P‐rich melt and with bulk Ni contents ranging from 0 to 40 wt%. The partitioning behavior for the 20 elements measured in this study did not vary with Ni content. Comparison of the Twannberg measurements with the experimental results required a correction factor to account for the fact that the experiments were conducted in a simplified system that did not contain a solid metal phase. Previously determined solid metal/P‐rich melt partition coefficients were applied to infer schreibersite/solid metal partitioning behavior from the experiments, and once this correction was applied, the two data sets showed broad similarities between the schreibersite/solid metal distribution of elements. However, there were also differences noted, in particular between the Ni and P contents of the solid metal relative to the schreibersite inferred from the experiments compared to that measured in the Twannberg sample. These differences support previous interpretations that subsolidus schreibersite evolution has strongly influenced the Ni and P content now present in the solid metal phase of IIG irons. Quantitative attempts to match the IIG solid metal composition to that of late‐stage IIAB irons through subsolidus schreibersite growth were not successful, but qualitatively, this study corroborates the striking similarities between the IIAB and IIG groups, which are highly suggestive of a possible genetic link between the groups as has been previously proposed.     

乌鲁木齐大气混合层厚度和稳定度与大气污染的关系

利用乌鲁木齐市4座10层100 m梯度气象塔2013年6月~2014年4月气象观测资料和7个环境监测站[WTBX]AQI[WTBZ]资料，计算并分析了大气混合层厚度和稳定度特征，探讨了大气混合层厚度和稳定度与污染的关系。结果表明：乌鲁木齐市混合层厚度夏季郊区高、城区低，冬季从南郊—城区—北郊随地势降低依次降低；夏季和冬季分别在1 559～1 772 m和526～1 156 m之间。地面至2 km以上每500 m高度间隔统计混合层厚度，500～1 000 m出现频率最多；月变化为6～9月基本在500 m以上，且每个高度区间其概率均超过10%，10月～次年2月1 500 m以上区间概率明显减小；日变化为中午13:00～16:00达到最高值，下午和傍晚迅速下降。白天较大的感热输送提供充足的热力条件，这也体现出白天以不稳定层结为主，夜间则以稳定层结为主。大气稳定度分类结果，夏季郊区和城区不稳定（A～C类）所占比例差不多，冬季北郊稳定（E、F类）所占比较最大、城区最弱。[WTBX]AQI指数冬季最大，从南郊—城区—北郊依次增大，这与采暖期污染物多、南郊比北郊地势高有利于扩散输送有关。总体来看，乌鲁木齐大气混合层厚度空间分布与气象要素、大气稳定度、地形等密切相关，对AQI[WTBZ]指数分布有重要影响，这对近地层大气污染状况预报有着重要的指导意义。     

Creating an isotopically similar Earth–Moon system with correct angular momentum from a giant impact

The giant impact hypothesis is the dominant theory explaining the formation of our Moon. However, the inability to produce an isotopically similar Earth–Moon system with correct angular momentum has cast a shadow on its validity. Computer-generated impacts have been successful in producing virtual systems that possess many of the observed physical properties. However, addressing the isotopic similarities between the Earth and Moon coupled with correct angular momentum has proven to be challenging. Equilibration and evection resonance have been proposed as means of reconciling the models. In the summer of 2013, the Royal Society called a meeting solely to discuss the formation of the Moon. In this meeting, evection resonance and equilibration were both questioned as viable means of removing the deficiencies from giant impact models. The main concerns were that models were multi-staged and too complex. We present here initial impact conditions that produce an isotopically similar Earth–Moon system with correct angular momentum. This is done in a single-staged simulation. The initial parameters are straightforward and the results evolve solely from the impact. This was accomplished by colliding two roughly half-Earth-sized impactors, rotating in approximately the same plane in a high-energy, off-centered impact, where both impactors spin into the collision.     

Simulating the Propagation of Infrasonic Waves and Estimating the Energy of the Chelyabinsk Meteoroid Explosion Observed on February 15, 2013

Results obtained from simulating the propagation of infrasonic waves from the Chelyabinsk meteoroid explosion observed on February 15, 2013, are given. The pseudodifferential parabolic equation (PDPE) method has been used for calculations. Data on infrasonic waves recorded at the IS31 station (Aktyubinsk, Kazakhstan), located 542.7 km from the likely location of the explosion, have been analyzed. Six infrasonic arrivals (isolated clearly defined pulse signals) were recorded. It is shown that the first “fast” arrival (F) corresponds to the propagation of infrasound in a surface acoustic waveguide. The rest of the arrivals (T1–T5) are thermospheric. The agreement between the results of calculations based on the PDPE method and experimental data is satisfactory. The energy E of the explosion has been estimated using two methods. One of these methods is based on the law of conservation of the acoustic pulse I, which is a product of the wave profile area S/2 of the signal under analysis and the distance to its source E _I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The other method is based on the relation between the energy of explosion and the dominant period T of recorded signal E _T [kt] = 1.02 × (T [s]²/σ)^3/2, where σ is the dimensionless distance determining the degree of nonlinear effects during the propagation of sound along ray trajectories. According to the data, the explosion energy E _I,T ranges from 1.87 to 32 kt TNT.     

Determination of Supermassive Black Hole Spins Based on the Standard Shakura–Sunyaev Accretion Disk Model and Polarimetric Observations

Based on spectropolarimetry for 47 type 1 active galactic nuclei observed with the 6-m BTA telescope, we have estimated the spins of the supermassive black holes at the centers of these galaxies. We have determined the spins based on the standard Shakura–Sunyaev accretion disk model. About 70% of the investigated active galactic nuclei are shown to have Kerr supermassive black holes with a dimensionless spin greater than 0.9.     

Critical Latitude in Tidal Dynamics Using the Kara Sea as an Example

It is well known that, within the linear nonviscous equations of tidal dynamics, the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components unlimitedly increase when approaching the critical latitude. It is also known that the linear equations of tidal dynamics with a constant and specified vertical eddy viscosity indicate the occurrence of significant tidal velocity shears in the near-bottom layer, which are responsible for increasing the baroclinic tidal energy dissipation, the turbulent kinetic energy, and the thickness of the bottom boundary layer. The first circumstance—the growth of the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components—is due to the elimination in the original equations of small terms, which are small everywhere except for the critical latitude zone. The second circumstance—the occurrence of significant tidal velocity shears—is due to the fact that internal tidal waves, which induce the dissipation of the baroclinic tidal energy and the diapycnal diffusion, are either not taken into account or described inadequately. It is suggested that diapycnal diffusion can lead to the degeneration (complete or partial) of tidal velocity shears, with all the ensuing consequences. The aforesaid is confirmed by simulation results obtained using the QUODDY-4 high-resolution three-dimensional finite-element hydrostatic model along the 66.25° E section, which passes in the Kara Sea across the critical latitude.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Factors influencing biodeposit production by the New Zealand freshwater mussel Echyridella menziesii

Laboratory experiments on the New Zealand freshwater mussel Echyridella menziesii were used to investigate the short-term effects (7–8 days) of food type on rates of biodeposition and benthic substrate respiration. Post-feeding biodeposition rates ranged from 0.34 to 1.52?mg?g^?1?h^?1 (mean?=?0.50?mg g^?1?h^?1) and were unaffected by the addition of toxin-producing Microcystis. Addition of suspended sediment (30?mg?L^?1) visibly altered substrate composition, and increased total and inorganic biodeposit production rates by 24–33% compared to mussels fed commercial phytoplankton stock. Biodeposition rates of mussels in lake bed substrates were 38% higher than those in silica sand for identical feeding regimes, suggesting that a significant proportion of material produced in this experiment could have been derived from feeding on organic matter in the lake bed sediments. Respiration rates were higher in treatments with Microcystis but were unaffected by the presence of mussels. This laboratory study suggests that biodeposition by E. menziesii is resilient to short-term exposure to Microcystis, and highlights the ability of mussels to alter benthic substrate composition by incorporating suspended sediment into substrates.     

Estimating average shock pressures recorded by impactite samples based on universal stage investigations of planar deformation features in quartz—Sources of error and recommendations

Planar deformation features (PDFs) in quartz are the most widely used indicator of shock metamorphism in terrestrial rocks. They can also be used for estimating average shock pressures that quartz-bearing rocks have been subjected to. Here we report on a number of observations and problems that we have encountered when performing universal stage measurements and crystallographically indexing of PDF orientations in quartz. These include a comparison between manual and automated methods of indexing PDFs, an evaluation of the new stereographic projection template, and observations regarding the PDF statistics related to the c-axis position and rhombohedral plane symmetry. We further discuss the implications that our findings have for shock barometry studies. Our study shows that the currently used stereographic projection template for indexing PDFs in quartz might induce an overestimation of rhombohedral planes with low Miller–Bravais indices. We suggest, based on a comparison of different shock barometry methods, that a unified method of assigning shock pressures to samples based on PDFs in quartz is necessary to allow comparison of data sets. This method needs to take into account not only the average number of PDF sets/grain but also the number of high Miller–Bravais index planes, both of which are important factors according to our study. Finally, we present a suggestion for such a method (which is valid for nonporous quartz-bearing rock types), which consists of assigning quartz grains into types (A–E) based on the PDF orientation pattern, and then calculation of a mean shock pressure for each sample.     

Analysis of the exoplanet containing system Kepler-13

We have applied the close binary system analysis program WinFitter, with its physically detailed fitting function, to an intensive study of the complex multiple system Kepler-13 using photometry data from all 13 short cadence quarters downloaded from the NASA Exoplanet Archive (NEA) (http://exoplanetarchive.ipac.caltech.edu). The data-point error of our normalized, phase-sequenced and binned (380 points per bin: 0.00025 phase interval) flux values, at 14 ppm, allows the model’s specification for the mean reference flux level of the system to a precision better than 1 ppm. Our photometrically derived values for the mass and radius of KOI13.01 are \(6.8\pm0.6~\mbox{M}_{\mathrm{J}}\) and \(1.44\pm0.04~\mbox{R}_{\mathrm{J}}\). The star has a radius of \(1.67\pm0.05~\mbox{R}_{\odot}\). Our modelling sets the mean of the orbital inclination \(i\) at \(94.35\pm0.14^{\circ}\), with the star’s mean precession angle \(\phi_{p}\)—\(49.1\pm5.0^{\circ}\) and obliquity \(\theta_{o}\)\(67.9 \pm 3.0^{\circ}\), though there are known ambiguities about the sense in which such angles are measured.Our findings did not confirm secular variation in the transit modelling parameters greater than their full correlated errors, as argued by previous authors, when each quarter’s data was best-fitted with a determinable parameter set without prejudice. However, if we accept that most of the parameters remain the same for each transit, then we could confirm a small but steady diminution in the cosine of the orbital inclination over the 17 quarter timespan. This is accompanied by a slight increase of the star’s precession angle (less negative), but with no significant change in the obliquity of its spin axis. There are suggestions of a history of strong dynamical interaction with a highly distorted planet rotating in a 3:2 resonance with its revolution, together with a tidal lag of \(\sim30~\mbox{deg}\). The mean precessional period is derived to be about 1000 y, but at the present time the motion of the star’s rotation axis appears to be supporting the gravitational torque, rather than providing the balance against it that would be expected over long periods of time.The planet has a small but detectable backwarming effect on the star, which helps to explain the difference in brightness just after transit and just before occultation eclipses. In assessing these findings it is recognized that sources of uncertainty remain, notably with possible inherent micropulsational effects, variations from other components of the multiple star, stellar activity, differential rotation and the neglect of higher order terms (than \(r_{1}^{5}\)) in the fitting function, where \(r_{1}\) is the ratio of the radius of the star to the mean orbital separation of planet and host star.