Electron acceleration in interaction of magnetic islands in large temporal-spatial turbulent magnetic reconnection

A new combined Fermi, betatron, and turbulent electron acceleration mechanism is proposed in interaction of magnetic islands during turbulent magnetic reconnection evolution in explosive astrophysical phenomena at large temporal-spatial scale(LTSTMR), the ratio of observed current sheets thickness to electron characteristic length, electron Larmor radius for low-β and electron inertial length for high-β, is on the order of 10^(10)–10^(11); the ratio of observed evolution time to electron gyroperiod is on the order of 10~7–10~9).The original combined acceleration model is known to be one of greatest importance in the interaction of magnetic islands; it assumes that the continuous kinetic-dynamic temporal-spatial scale evolution occurs as two separate independent processes.In this paper, we reconsider the combined acceleration mechanism by introducing a kinetic-dynamic-hydro full-coupled model instead of the original micro-kinetic or macro-dynamic model.We investigate different acceleration mechanisms in the vicinity of neutral points in magnetic islands evolution, from the stage of shrink and breakup into smaller islands(kinetic scale), to the stage of coalescence and growth into larger islands(dynamic scale), to the stages of constant and quasi-constant(contracting-expanding) islands(hydro scale).As a result, we give for the first time the acceleration efficiencies of different types of acceleration mechanisms in magnetic islands' interactions in solar atmosphere LTSTMR activities(pico-, 10^(–2)–10~5 m; nano-, 10~5–10~6 m; micro-, 10~6–10~7 m; macro-, 10~7–10~8 m; large-,10~8–10~9 m).

基于U形卷积神经网络的震相识别与到时拾取方法研究

精确获取震相到时是地震定位和地震走时成像等研究的重要基础.近年来,随着地震台站的不断加密,地震台网监测到的地震数量成倍增长,发展快速、准确、适用性强的震相到时自动拾取算法是地震行业的迫切需求.本文在前人工作基础上,发展了Pg、Sg震相自动识别与到时拾取的U网络算法(Unet_cea),使用汶川余震和首都圈地震台网记录的89344个不同震级、不同信噪比的样本进行训练和测试.研究表明,U网络能够较好地识别Pg、Sg震相类型和拾取到时,Pg、Sg震相的正确识别率分别为81%和79.1%,与人工标注到时的均方根误差分别为0.41s和0.54s.U网络在命中率、均方根误差等性能指标上均明显优于STA/LTA和峰度分析自动拾取方法.研究获得的最优模型可以为区域地震台网的自动处理提供辅助.

Parallel Numerical Analysis on the Rheology of the Martian Ice-Rock Mixture

The distribution and amount of ground ice on Mars is an important issue to be addressed for the future exploration of the planet.The occurrence of interstitial ice in Martian frozen ground is indicated by landforms,such as fluidized ejecta craters,softened terrain,and fretted channels.How-ever,experimental data on the rheology of ice-rock mixture under Martian physical conditions are sparse,and the amount of ground ice that is required to produce the viscous deformation observed in Martian ice-related landforms is still unknown.In our study,we put forward a three-dimensional non-Newtonian viscous finite element model to investigate the behavior of ice-rock mixtures numeri-cally.The randomly distributed tetrahedral elements are generated in regular domain to represent the natural distribution of ice-rock materials.Numerical simulation results show that when the volume of rock is less than 40%,the rheology of the mixture is dominated by ice,and there is occurrence of a brittle-ductile transition when ice fraction reaches a certain value.Our preliminary results contribute to the knowledge of the determination of the rheology and ice content in Martian ice-rock mixture.The presented model can also be utilized to evaluate the amount of ground ice on Mars.

Melt Evolution above a Spontaneously Retreating Subducting Slab in a Three-Dimensional Model

Dehydration of the subducting slab favors the melting of the surrounding mantle. Water content and melt evolution atop a spontaneously retreating subducting slab are reported in a three- dimensional （3-D） model. We find that fluids, including water and melts in the rocks, vary substantially along the trench, which cannot be found in two-dimensional （2-D） models. Their maxima along the subducting slab are mainly located at about 50 to 70 and 120 to 140 kin. Volumetric melt production rate changes spatially and episodically atop the slab, which may reflect the intensity and variations of volcanoes.

Influences of Lower-Mantle Properties on the Formation of Asthenosphere in Oceanic Upper Mantle

Asthenosphere is a venerable concept based on geological intuition of Reginald Daly nearly 100 years ago. There have been various explanations for the existence of the asthenosphere. The concept of a plume-fed asthenosphere has been around for a few years due to the ideas put forth by Yamamoto et al.. Using a two-dimensional Cartesian code based on finite-volume method, we have investigated the influences of lower-mantle physical properties on the formation of a low-viscosity zone in the oceanic upper mantle in regions close to a large mantle upwelUng. The rheological law is Newtonian and depends on both temperature and depth. An extended-Boussinesq model is assumed for the energetics and the olivine to spinel, the spinel to perovskite and perovskite to post-perovskite （ppv） phase transitions are considered. We have compared the differences in the behavior of hot upweilings passing through the transition zone in the mid-mantle for a variety of models, starting with constant physical properties in the lower-mantle and culminating with complex models which have the post-perovskite phase transition and depth-dependent coefficient of thermal expansion and thermal conductivity. We found that the formation of the asthenosphere in the upper mantle in the vicinity of large upwellings is facilitated in models where both depth-dependent thermal expansivity and conductivity are included. Models with constant thermal expansivity and thermal conductivity do not produce a hot low-viscosity zone, resembling the asthenosphere. We have also studied the influences of a cylindrical model and found similar results as the Cartesian model with the important difference that upper-mantle temperatures were much cooler than the Cartesian model by about 600 to 700 K. Our findings argue for the potentially important role played by lower-mantle material properties on the development of a plume-fed asthenosphere in the oceanic upper mantle.

Regional Stress Fields under Tibet from 3D Global Flow Simulation

Tibetan area is the most active continental collision zone on earth. Several major earthquakes occurred around the boundaries of Tibetan plateau and caused massive damages and casualties. The dynamics of this area is not well understood due to the complex structure of Tibet and its surrounding area. In this study, a 3D global flow simulation with only viscous rheology is applied to studying the stress distribution in this area, and the interaction between Tibet and its surrounding areas is investigated. Finally, the possibility of combining regional modeling with global models is also discussed.

Submicron size-scale mapping of carbonate effective elastic properties from FIB-SEM images and finite element method

In this paper, an automatic unstructured focused ion beam （FIB） and scanning electron microscopy （SEM） images induced representative volume element （RVE） finite element （FE） method is developed to predict submicron scale carbonate rock effective Young＇s and bulk moduli and Poisson＇s ratio on parallel CPU-GPU platform. Based on high resolution-contrast surface morphology and internal fabric-texture structure images from carbonate rock specimen （covered 0.12-64 μm2 area and 8000 μm3 domain）, the cubic RVE FE models are constructed from different sites through Avizo with user-defined parameters Matlab coding. The effective Young＇s and bulk moduli and Poisson＇s ratio of the different RVEs and porosity and pore size are computed by using periodic boundary condition in the well-known FE software Abaqus. FE mesh sensitivity analysis has been conducted where all moduli converge to a certain constant value at larger FE mesh density. The effect of fabric-texture （pore size, shape, and distribution） on the elastic properties is discussed. The correlations between the computed effective elastic properties and pore size, porosity, RVE size have been established. The simulation results show that the effective Young＇s and bulk moduli and Poisson＇s ratio have strong anisotropic behavior and depend on RVE size, porosity and pore size. The RVE size, porosity and pore size are three independent factors in affecting of the effective elastic moduli, the effect mechanism of porosity and pore size is same while the effect mechanism of RVE size is difference.

Small-Scale Mineralogical Heterogeneity from Variations in Phase Assemblages in the Transition Zone and D" Layer Predicted by Convection Modelling

Small-scale heterogeneity in the deep mantle is concentrated in the upper-mantle transition zone （TZ）, in the depth range 410-660 km and also at the bottom 250 km D＂ region. This encourages a more detailed investigation of the potential for seismic reflectivity imaging by modelling heterogeneous structures in mantle convection models including phase transitions of the TZ and D＂ regions. We applied finite elements with variable spacing near the boundary layers in 2-D cylindrical geometry that allow for sufficient spatial resolution. We investigated several models including an extended Boussinesq （EBA） model, focused on the D＂ region, and a compressible （ALA） model for the TZ region. The results for the D＂ region show typical lens-shaped structures of post-perovskite （PPV） embedded in the perovskite （PV） background mantle, where the thickness of the lenses, at 200-400 km, strongly depends on the Clapeyron slope of the PV-PPV transition. A second phase transition （double crossing） occurs in case the core temperature is higher than the intercept temperature Ti. Our phase-dependent rheology results in contrasting effective viscosity between PV and PPV. Our model results reveal a distinctly clear mechanical weakening of the PPV lenses with about an order of magnitude lower viscosity. The shear wave-speed distributions computed from our convection results are strongly correlated with the heterogeneous distribution of the mineral phase. Gradients in the seismic wave-speed that are the target of seismological reflectivity imaging are clearly revealed. The wave-speed results show a clear resolution of the top and bottom interfaces of the PPV lenses. Our ALA model for the TZ is based on a thermodynamical model for the magnesium end- member of an olivine-pyroxene mantle. The model predicts a much more complex distribution of mineral phases, compared to our D＂ results, in agreement with the greater number of mineral phases involved in the olivine-pyroxene phase diagram for the P, T conditions of the transition zone. Near cold downwelling flows representing subducting lithospheric slabs, where the local geotherm can differ by up to 1 000 K from the horizontal average, and small-scale lateral variations in the mineral phases can occur.