Maximal size of chondrules in shock wave heating model: Stripping of liquid surface in a hypersonic rarefied gas flow

Abstract— We examined partially molten dust particles that have a solid core and a surrounding liquid mantle, and estimated the maximal size of chondrules in a framework of the shock wave heating model for chondrule formation. First, we examined the dynamics of the liquid mantle by analytically solving the hydrodynamics equations for a core‐mantle structure via a linear approximation. We obtained the deformation, internal flow, pressure distribution in the liquid mantle, and the force acting on the solid core. Using these results, we estimated conditions in which liquid mantle is stripped off from the solid core. We found that when the particle radius is larger than about 1–2 mm, the stripping is expected to take place before the entire dust particle melts. So chondrules larger than about 1–2 mm are not likely to be formed by the shock wave heating mechanism. Also, we found that the stripping of the liquid mantle is more likely to occur than the fission of totally molten particles. Therefore, the maximal size of chondrules may be determined by the stripping of the liquid mantle from the partially molten dust particles in the shock waves. This maximal size is consistent with the sizes of natural chondrules.     

Ca‐Al‐Fe‐rich inclusion in the Vigarano CV3 chondrite

Abstract– An anomalous Ca‐Al‐Fe‐rich spherical inclusion (CAFI) was found in the Vigarano CV3 chondrite. The CAFI has an igneous texture and contains large amounts of almost pure and coarse‐grained hercynite grains (approximately 56 vol%) as well as refractory phases such as grossite and perovskite. However, melilite and Mg‐spinel, which are common in ordinary Ca‐Al‐rich inclusions, are very rare (<1 vol%). Another unique characteristic of the CAFI is the presence in its core of dmitryivanovite (CaAl₂O₄), which was formed by shock metamorphism of a low‐pressure form of CaAl₂O₄ that was originally crystallized from a molten droplet. The fine‐grained hercynite and unidentified aluminous phase in the rim of the CAFI may have been produced from grossite during aqueous alteration in the Vigarano parent body.     

Sound velocity measurements in dhcp-FeH up to 70 GPa with inelastic X-ray scattering: Implications for the composition of the Earth's core

We have determined the density evolution of the sound velocity of dhcp-FeH_x (x ≈ 1) up to 70 GPa at room temperature, by inelastic X-ray scattering and by X-ray diffraction. We find that the variation of V_P with density is different for the ferromagnetic and nonmagnetic dhcp-FeH_x, and that only nonmagnetic dhcp-FeH_x follows Birch's law. Combining our results with Birch's law for iron and assuming an ideal two-component mixing model, we obtain an upper bound of the hydrogen content in the Earth's inner core, 0.23(6) wt.% H, corresponding to FeH_0.13(3). The iron alloy with 0.23(6) wt.% H can satisfy the density, and compressional and shear sound velocities of the PREM inner core, assuming that there are no other light elements in the inner core.     

Computation of Solar Radiative Fluxes by 1D and 3D Methods Using Cloudy Atmospheres Inferred from A-train Satellite Data

This study used realistic representations of cloudy atmospheres to assess errors in solar flux estimates associated with 1D radiative transfer models. A scene construction algorithm, developed for the EarthCARE mission, was applied to CloudSat, CALIPSO and MODIS satellite data thus producing 3D cloudy atmospheres measuring 61 km wide by 14,000 km long at 1 km grid-spacing. Broadband solar fluxes and radiances were then computed by a Monte Carlo photon transfer model run in both full 3D and 1D independent column approximation modes. Results were averaged into 1,303 (50 km)² domains. For domains with total cloud fractions A _c  < 0.7 top-of-atmosphere (TOA) albedos tend to be largest for 3D transfer with differences increasing with solar zenith angle. Differences are largest for A _c  > 0.7 and characterized by small bias yet large random errors. Regardless of A _c , differences between 3D and 1D transfer rarely exceed ±30 W m^?2 for net TOA and surface fluxes and ±10 W m^?2 for atmospheric absorption. Horizontal fluxes through domain sides depend on A _c with ～20% of cases exceeding ±30 W m^?2; the largest values occur for A _c  > 0.7. Conversely, heating rate differences rarely exceed ±20%. As a cursory test of TOA radiative closure, fluxes produced by the 3D model were averaged up to (20 km)² and compared to values measured by CERES. While relatively little attention was paid to optical properties of ice crystals and surfaces, and aerosols were neglected entirely, ～30% of the differences between 3D model estimates and measurements fall within ±10 W m^?2; this is the target agreement set for EarthCARE. This, coupled with the aforementioned comparison between 3D and 1D transfer, leads to the recommendation that EarthCARE employ a 3D transport model when attempting TOA radiative closure.     

Uncertainty Estimate of Surface Irradiances Computed with MODIS-, CALIPSO-, and CloudSat-Derived Cloud and Aerosol Properties

Differences of modeled surface upward and downward longwave and shortwave irradiances are calculated using modeled irradiance computed with active sensor-derived and passive sensor-derived cloud and aerosol properties. The irradiance differences are calculated for various temporal and spatial scales, monthly gridded, monthly zonal, monthly global, and annual global. Using the irradiance differences, the uncertainty of surface irradiances is estimated. The uncertainty (1σ) of the annual global surface downward longwave and shortwave is, respectively, 7?W?m^?2 (out of 345?W?m^?2) and 4?W?m^?2 (out of 192?W?m^?2), after known bias errors are removed. Similarly, the uncertainty of the annual global surface upward longwave and shortwave is, respectively, 3?W?m^?2 (out of 398?W?m^?2) and 3?W?m^?2 (out of 23?W?m^?2). The uncertainty is for modeled irradiances computed using cloud properties derived from imagers on a sun-synchronous orbit that covers the globe every day (e.g., moderate-resolution imaging spectrometer) or modeled irradiances computed for nadir view only active sensors on a sun-synchronous orbit such as Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation and CloudSat. If we assume that longwave and shortwave uncertainties are independent of each other, but up- and downward components are correlated with each other, the uncertainty in global annual mean net surface irradiance is 12?W?m^?2. One-sigma uncertainty bounds of the satellite-based net surface irradiance are 106?W?m^?2 and 130?W?m^?2.     

Elemental and isotopic abundances of lithium in chondrule constituents in the Allende meteorite

Elemental and isotopic abundances of lithium in chondrule constituents in the Allende CV3 meteorite were determined using secondary ion mass spectrometry. Olivines and mesostasis dominated by a feldspathic phase are depleted in Li ( and , respectively). In contrast, low-Ca pyroxenes and mesostasis dominated by a Na-rich phase are enriched in Li ( and , respectively) and the interchondrule matrix is generally enriched in Li ( on average). The Li isotopic abundance of olivine ranges from to 21. The spatial distributions of elemental and isotopic abundances of Li in olivines within individual chondrules exhibit no systematic pattern. This suggests that the distribution of Li in olivine was not disturbed during aqueous alteration or thermal metamorphism on the Allende meteorite parent body. Although mesostasis is the last crystallizing phase from a chondrule melt and is expected to be enriched in Li, in the Allende meteorite it is generally depleted in Li. We suggest that during aqueous alteration on the CV asteroid, Li in mesostasis was leached out by aqueous fluids. The Li-enriched Na-rich mesostasis was probably produced later by infiltration of Na-rich fluids. It seems likely that aqueous fluids sequestered alkali elements from the Allende-chondrite region in the CV parent asteroid, although significant amounts of Li are preserved in ferrous olivine in the interchondrule matrix.     

Crustal structure beneath Aso Caldera,Southwest Japan,as derived from receiver function analysis

Aso Volcano experienced a huge pyroclastic eruption 90 thousand years ago, and formed a large caldera (18 km × 25 km). In order to test the hypothesis of a magma body in the mid and lower crust that has been suggested geophysically and geochemically, we investigated seismic velocity discontinuities and velocity structure beneath Aso Caldera using receiver functions and a genetic algorithm inversion. We confirm the existence of the Moho at depths between 30 km and 35 km and a large velocity anomaly should exist in the deep portion of the crust beneath Aso Caldera, from imaging of receiver functions observed only at stations outside the caldera. As a result of a more detailed examination with GA inversion, a low velocity layer is detected at depths between 10 km and 24 km beneath the western part of the caldera. S-wave velocity of the layer is estimated to be 2.0–2.4 km/s. We estimate that the low velocity layer contains at most 15% melt or 30% aqueous fluid. The layer exists near the Conrad and at the same depths as the swarm of the low frequency earthquakes and a compressional and dilatational deformation source which are expected to be caused by fluid movement beneath the middle-eastern part of the caldera. Fluid contained in the layer might be related with huge pyroclastic eruptions of Aso Volcano.     

基于全球降水数据估计值的地表径流模拟--以长江上游地区为例

为了评价在相对较短时间内针对长江上游地区地表水径流所建立模型的模拟效果,以及检验以GCM模型和其他卫星数据所估算的降水数据作为输入数据的可行性,选择分布式水文模型HSPF以及1987年和1988年的ISLSCP降水数据作为输入数据。模型模拟结果表明:从整个长江上游地区看,在校正期内,5天平均流量的Nash–Sutcliffe相关系数(R2) 为0.94;在验证期内,Nash-Sutcliffe相关系数(R2) 为0.95。此外,该模型对长江上游主要支流的5天平均流量的模拟效果也很好,R2的值在0.46到0.96之间。例外的情况主要发生在沱江和嘉陵江,模型对2年洪水期的峰值流速的估计值偏低,沱江只有实际值的71%,嘉陵江只有实际值的61%。ISLSCP估计的降水比实际测量的降水频繁且程度要弱,这可能是HSPF不能在所有时间和所有区域都具有较好模拟效果的一个主要原因。     

Recognition of a Locked State in Plate Subduction from Microearthquake Seismicity

—A tectonic state of a locked subduction is considered to be a possible source of a future interplate earthquake. Discriminating an actually locked state to verify its extent is therefore essential in constructing an accurate prospect against the forthcoming earthquake. Micorearthquake seismicity is an effective tool for such an analysis because it is considered to be a faithful indicator of the stress state, and is expected to exhibit a characteristic pattern in the area where the locked state in the subduction appears with a certain stress concentration. Focusing on the microearthquake seismicity around the Tokai district in central Japan, where a large interplate earthquake is feared to occur, we tried to identify such an area of locked subduction on the Philippine Sea plate, possibly related to the future earthquake. We investigated the microearthquake seismicity from various perspectives. First, the hypocenter distribution was analyzed to identify the extent of the locked area. The characteristic profile of the distribution was presumed to represent a stress concentrated area induced from the mechanical contact between both plates. The second approach is to interpret stress patterns reflected in focal mechanisms. The locked state was recognized and verified by a comparison of the P-axis distribution pattern with that expected from a model imaging a partially locked subduction. The third approach is to monitor the temporal change of the seismic wave spectrum. Analyzing predominant frequencies of P and S waves and monitoring their changes for a period of 10 years, we found a trend of gradual increase common to both waves. This means an increase of stress drop in microfracturings, and in its turn implies accumulation of stress around the focus area. The rate of the stress change converted from the frequency change was compared with the result derived from a numerical simulation. The simulation, performed on the basis of a constitutive friction law for a stick sliding on the plate interface, computed a changing rate of the maximum shear stress around the locked zone and showed its spatial variation along the subduction axis. Thus the simulated result indicated a certain compatibility with the observed one. Although ambiguities and uncertainties still exist in the study, all the results derived here seem to indicate an identical conclusion that the plate subduction is actually locked in this region at present.