Effects of scattering and dust grain size on the temperature structure of protoplanetary discs: a three-layer approach

The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimetre observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, and for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations, including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size ≳10 μm) even without scattering. A three-layer interpretation is required for grain size ≲10 μm. When the grain size is 0.1–10 μm, the numerical simulations show that the isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size ≳10 μm (i.e. grey scattering), the numerical simulations show that the isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the 'green-house effect' due to the scattering of the cold diffuse radiation in the interior.     

Transformation textures,mechanisms of formation of high‐pressure minerals in shock melt veins of L6 chondrites,and pressure‐temperature conditions of the shock events

Abstract— The high‐pressure polymorphs of olivine, pyroxene, and plagioclase in or adjacent to shock melt veins (SMVs) in two L6 chondrites (Sahara 98222 and Yamato 74445) were investigated to clarify the related transformation mechanisms and to estimate the pressure‐temperature conditions of the shock events. Wadsleyite and jadeite were identified in Sahara 98222. Wadsleyite, ringwoodite, majorite, akimotoite, jadeite, and lingunite (NaAlSi₃O₈‐hollandite) were identified in Yamato 74445. Wadsleyite nucleated along the grain boundaries and fractures of original olivine. The nucleation and growth of ringwoodite occurred along the grain boundaries of original olivine, and as intracrystalline ringwoodite lamellae within original olivine. The nucleation and growth of majorite took place along the grain boundaries or fractures in original enstatite. Jadeite‐containing assemblages have complicated textures containing “particle‐like,” “stringer‐like,” and “polycrystalline‐like” phases. Coexistence of lingunite and jadeite‐containing assemblages shows a vein‐like texture. We discuss these transformation mechanisms based on our textural observations and chemical composition analyses. The shock pressure and temperature conditions in the SMVs of these meteorites were also estimated based on the mineral assemblages in the SMVs and in comparison with static high‐pressure experimental results as follows: 13–16 GPa, >1900 °C for Sahara 98222 and 17–24 GPa, >2100 °C for Yamato 74445.     

Photoluminescence of baratovite and katayamalite

The photoluminescence (PL) and optical excitation spectra of baratovite in aegirine syenite from Dara-i-Pioz, Tien Shan Mts., Tajikistan and katayamalite in aegirine syenite from Iwagi Islet, Ehime, Japan were obtained at 300 and 80 K. Under short wave (253.7 nm) ultraviolet light, baratovite and katayamalite exhibited bright blue-white luminescence. The PL spectrum of baratovite at 300 K consisted of a wide band with a peak at approximately 406 nm and a full width at half maximum (FWHM) of approximately 6.32k cm⁻¹. The excitation spectrum of the blue-white luminescence from baratovite at 300 K consisted of a prominent band with a peak at approximately 250 nm. The PL and excitation spectra of katayamalite were similar to those of baratovite. The luminescence from these minerals was attributed to the intrinsic luminescence from the TiO₆ center.     

Fragmentation model dependence of collision cascades

Mass depletion of bodies through successive collisional disruptions (i.e., collision cascade) is one of the most important processes in the studies of the asteroids belt, the Edgeworth-Kuiper belt, debris disks, and planetary formation. The collisional disruption is divided into two types, i.e., catastrophic disruption and cratering. Although some studies of the collision cascades neglected the effect of cratering, it is unclear which type of disruption makes a dominant contribution to the collision cascades. In the present study, we construct a simple outcome model describing both catastrophic disruption and cratering, which has some parameters characterizing the total ejecta mass, the mass of the largest fragment, and the power-law exponent of the size distribution of fragments. Using this simple outcome model with parameters, we examine the model dependence of the mass depletion time in collision cascades for neglect of coalescence of colliding bodies due to high collisional velocities. We find the cratering collisions are much more effective in collision cascades than collisions with catastrophic disruption in a wide region of the model parameters. It is also found that the mass depletion time in collision cascades is mainly governed by the total ejecta mass and almost insensitive to the mass of the largest fragment and the power-law exponent of fragments for a realistic parameter region. The total ejecta mass is usually determined by the ratio of the impact energy divided by the target mass (i.e. Q-value) to its threshold value for catastrophic disruption, as well as in our simple model. We derive a mass depletion time in collision cascades, which is determined by of the high-mass end of collision cascades. The mass depletion time derived with our model would be applicable to debris disks and planetary formation.     

Determination of Fluorine and Chlorine by Pyrohydrolysis and Ion Chromatography: Comparison with Alkaline Fusion Digestion and Ion Chromatography

A method to determine F and Cl in silicate materials by employing pyrohydrolysis and ion chromatography (IC) is described. Pyrohydrolysis involved mixing a pulverised sample (∼ 40 mg) with V₂O₅ (∼ 160 mg) and heating to 1100 °C under a wet oxygen flow in a quartz tube. Recovery yields of F and Cl were ∼ 97% using a NaF + NaCl standard solution. Detection limits of the pyrohydrolysis-IC method for silicate samples were 0.36 and 0.69 μg g^-1 for F and Cl, respectively. Fluorine and Cl concentrations were determined in the reference materials JB-2, JB-3 and JA-1 from the GSJ; BCR-2, BHVO-1, BHVO-2, AGV-1 and AGV-2 from the USGS; and NIST SRM 610, 612 and 614 glasses. Precisions (RSD) for determinations of F were 1–13% (except NIST SRM 614) and 2–19% for Cl, and were dependent on the concentration and blank correction. Most results obtained in this study were in good agreement with those of previous studies. In comparison, the Na₂CO₃ + ZnO fusion method at 900 °C showed that the yields of F and Cl by alkaline fusion systematically decreased with fusion duration time. The yields were 84% and 83% for JB-3, inferring that F and Cl were lost in this alkaline fusion.     

Intercomparison and analyses of the climatology of the West African Monsoon in the West African Monsoon Modeling and Evaluation project (WAMME) first model intercomparison experiment

This paper briefly presents the West African Monsoon (WAM) Modeling and Evaluation Project (WAMME) and evaluates WAMME general circulation models’ (GCM) performances in simulating variability of WAM precipitation, surface temperature, and major circulation features at seasonal and intraseasonal scales in the first WAMME experiment. The analyses indicate that models with specified sea surface temperature generally have reasonable simulations of the pattern of spatial distribution of WAM seasonal mean precipitation and surface temperature as well as the averaged zonal wind in latitude-height cross-section and low level circulation. But there are large differences among models in simulating spatial correlation, intensity, and variance of precipitation compared with observations. Furthermore, the majority of models fail to produce proper intensities of the African Easterly Jet (AEJ) and the tropical easterly jet. AMMA Land Surface Model Intercomparison Project (ALMIP) data are used to analyze the association between simulated surface processes and the WAM and to investigate the WAM mechanism. It has been identified that the spatial distributions of surface sensible heat flux, surface temperature, and moisture convergence are closely associated with the simulated spatial distribution of precipitation; while surface latent heat flux is closely associated with the AEJ and contributes to divergence in AEJ simulation. Common empirical orthogonal functions (CEOF) analysis is applied to characterize the WAM precipitation evolution and has identified a major WAM precipitation mode and two temperature modes (Sahara mode and Sahel mode). Results indicate that the WAMME models produce reasonable temporal evolutions of major CEOF modes but have deficiencies/uncertainties in producing variances explained by major modes. Furthermore, the CEOF analysis shows that WAM precipitation evolution is closely related to the enhanced Sahara mode and the weakened Sahel mode, supporting the evidence revealed in the analysis using ALMIP data. An analysis of variability of CEOF modes suggests that the Sahara mode leads the WAM evolution, and divergence in simulating this mode contributes to discrepancies in the precipitation simulation.     

A trial for monitoring temporal variation of seismic velocity using an ACROSS system

Abstract The temporal variation of seismic velocity near the Nojima Fault, which ruptured during the 1995 Kobe earthquake (Hyogo-ken Nanbu earthquake), was detected using an accurately controlled routine-operated seismic source (ACROSS). The source generates elastic waves by a centrifugal force of an eccentric mass rotating around an axis. The mass is driven with an AC servomotor whose angular position is accurately controlled with reference to a very accurate global positioning system (GPS) clock. The error of the mass' position is less than 0.002 radian and does not accumulate. As a result, the source generates sinusoidal waves of very narrow spectral peaks enabling their detection with an excellent signal-to-noise ratio. Although the stability of the rotation is quite excellent, a large daily variation was found, which seems to be caused by changes in atmospheric temperature. The daily variation was 10% in amplitude and 0.1 radian in phase of the signal observed at the 800 m borehole seismometer. A significant variation was found to be due to that of coupling between the rotational source and the foundation made of reinforced concrete in which the source was situated. In order to make a correction on the signal of the 800 m borehole seismometer, the vibration of the foundation was measured and modeled assuming a rigid body movement. The correction successfully reduced the daily variation by approximately 90%, resulting in a variation of 1% in amplitude and 0.01 radian in phase. The phase variation of 0.01 radian corresponds to 100 μs and less than 0.1% in velocity over 1000 m between the source and the receiver.     

Upward longwave radiation from a non-black urban canopy

The longwave upward radiation was calculated for an urban canopy by using a Monte Carlo model. The effects of the urban geometry were examined in terms of the fractional roof area, the height of the buildings and the emissivity. The urban canopy consists of identically sized buildings and the ground surfaces. The model allows for the temperature differences between the buildings and the ground surface and for multiple reflections in the canyon.The Monte Carlo results show that neglect of the geometric effects causes significant errors in calculated upward radiation: calculations with area-weighting of the radiation emitted from flat homogeneous surfaces are not appropriate. The upward flux is a nonlinear function of the fractional roof area, which may be approximated by a function of the square or cube of the fractional roof area. Neglect of the reflections by non-black surfaces (emissivity<1) underestimates the upward flux by a few percent for a canopy of emissivity=0.9. Radiation effects due to multiple reflections in the canyon are parameterized by use of the view factor and the fractional roof area. The parameterization scheme yields accurate results.     

Phase and group velocities and Q of mantle Love and Rayleigh waves of the first two modes and their azimuthal dependences for the 1963 Kurile Islands earthquake

Phase and group velocities and Q of mantle Love and Rayleigh waves from the 1963 Kurile Islands earthquake (M_w = 8.5) were determined over 37 great circle paths by a time variable filtering technique, in a period range 100–500 s for the fundamental modes and 100–275 s for the first higher modes. The preliminary reference Earth model (PREM) explains reasonably well the average dispersion results for the fundamental Love and Rayleigh waves. There exists a small, but significant inconsistency between the observation and the model for the first higher Love and Rayleigh waves. The Q structure of PREM is inconsistent with the observation for the fundamental Love waves, but explains other observations reasonably well. The dispersion of each mode shows a clear azimuthal dependence from which the four azimuthal windows were established. The phase and group velocity measurements for each window were, in general, shown to be mutually consistent. The azimuthal variations are largest for the first higher Rayleigh waves, indicating strong lateral heterogeneity in the structure of the low velocity zone. The first of the four windows is characterized by the largest fraction of Precambrian shields and the second window by the largest fraction of normal oceans. A comparison of these two windows may give some insight into deep lateral heterogeneity between continents and oceans. The observed phase and group velocities of the first window are systematically higher than those of the second window for the fundamental Love and Rayleigh waves at periods up to 400 s, and for the first higher Love and Rayleigh waves up to 175 s. Their differences are greatest for the first higher Rayleigh waves and least for the fundamental Rayleigh waves. Although the fundamental Rayleigh waves show the least velocity differences, their persistence up to a period of longer than 300 s is in striking contrast with some of the pure path phase velocities derived earlier for continents and oceans. A set of models for continents and oceans. PEM-C and PEM-O are not consistent with our observation. The third azimuthal window is characterized by trench-marginal seas and the fourth window by mountainous areas, typically the Asian high plateaus from northern China to the Middle East through Tibet. A comparison of these two windows gives some information about deep structural differences between subduction zones and continental collision zones, both belonging to plate convergence zones. For the fundamental and the first higher Love waves, the phase and group velocities for the third window are markedly low, whereas those for the fourth window are somewhat comparable to those for the second window. Slow Rayleigh waves are evident for two windows, with the fourth window apparently being the slowest for the fundamental Rayleigh above 200 s and for the first higher Rayleigh. For the fundamental Rayleigh waves, the third window is very slow below 200 s, but becomes progressively fast as the period increases and tends to be the fastest window around 400 s, suggesting a deep seated high velocity anomaly beneath trench-marginal seas. The dispersion characteristics of the fourth window indicate a thick high velocity lid with an extensive low velocity zone beneath it. The shield-like lithosphere, coupled with an extensive low velocity zone, may be a characteristic feature of continental collision zones. The particle motion of the fundamental Love waves was found not to be purely transverse to a great-circle connecting the epicenter to a station. The departure from the purely transverse motion is systematic among different periods, different G arrivals (G₂, G₃,…) and different stations, which may be interpreted as being due to lateral refraction.     

Interstadial coral reef terraces and relative sea-level changes during marine oxygen isotope stages 3–4, Kikai Island, central Ryukyus, Japan

Coral reef terraces are one of the best recorders of relative sea-level changes during the last glacial cycle. Thus far, knowledge of relative sea-level record based on coral reefs during the marine Oxygen Isotope Stage (OIS) 3 has been limited to studies of the Huon Peninsula, Papua New Guinea. High-precision a α-spectrometric ²³⁰Th/²³⁴U dating demonstrated an offlapping sequence of five coral reef complexes, ages of which are 66, 64, 62, 55 and 52 ka, in the northern part of Kikai Island, central Ryukyus of Japan. Interstadial reefs, characterized by deepening-upward sequences of coral assemblages, recorded three hemicycles from transgression to highstand at 52, 62, and 66 ka, during which these reefs were drowned. These highstands in the relative sea-level record can be correlated with the eustatic record reconstructed from the Huon reef terraces and with the interstadials 14, 18, and 19 of the GISP 2 oxygen isotope record. This consistency confirms the Huon sea-level record of OIS 3 and implies that the eustatic sea level responded to the millennial-scale climate changes even during the glacial period of OIS 4.