Advanced Search



ISSN  2096-3955

CN  10-1502/P

Browse Articles


Statistical study of “trunk-like” heavy ion structures in the inner magnetosphere
Hao Zhang, YaBing Wang, JianYong Lu
Recently Published, doi: 10.26464/epp2022032
[Abstract](165) [FullText HTML](7) [PDF 695KB](10)
We present a statistical study of “trunk-like” structures observed in He+ and O+ in the inner magnetosphere. The main characteristic of these structures is that the energy of the peak flux decreases earthward. Using observations from the Helium Oxygen Proton Electron (HOPE) instrument onboard Van Allen Probe A, we identify the trunks observed from November 2012 to June 2019 and obtain the universal time, L shell, magnetic local time (MLT), and energy information of each trunk’s root and tip. We then investigate the behavior of trunks in terms of their frequency of occurrence, temporal evolution, spatial and energy distribution, and trunk dependence on different geomagnetic indices. We find that (1) the trunks are always located at L = 1.5−4.0 and have a preferential location mainly concentrated at MLT = 18−24, (2) the sector within MLT = 14−16 is a forbidden zone without trunk roots, and (3) the energy of He+ trunks is the largest near dusk and gradually decreases in the counterclockwise direction, whereas the energy of O+ trunks is relatively evenly distributed with MLT and L. The differences between He+ and O+ trunks are probably due to the different charge exchange and Coulomb collision lifetime. The dependence on different geomagnetic indices indicates that the trunk structures occur more frequently during relatively quiet periods.
No seamount subduction, no magmatic arc?
XiaoBo He, JunTao Tao, Yi Cao, FaBin Pan, EnYuan He, LingMin Cao, YiXian Zheng
Recently Published, doi: 10.26464/epp2022031
[Abstract](216) [FullText HTML](21) [PDF 2641KB](20)
Water is essential for the formation of a magmatic arc by lowering the melting temperature of materials in the mantle wedge. As such, it is logical to attribute the absence of a magmatic arc to insufficient water released from the subducting plate, although a number of other factors may cause volcanic arc quiescence as well, such as a slab window or flat slab subduction. In this contribution, we present a possible but testable correlation between the occurrence of a magmatic arc and seamount subduction in light of bathymetric data obtained near trenches. This correlation, if it holds true, in turn means that a magmatic arc is unlikely to occur when the subducting slabs have not been severely fractured and that one of the main reasons for excluding effects such as the slab window or flat slab subduction may be that the plate is not accompanied by seamounts. Therefore, the role that seamount subduction plays in recycling water back into the mantle deserves more attention from the earth sciences community.


Vertical Crustal Deformation Velocity and its Influencing Factors over the Qinghai-Tibet Plateau based on the Satellite Gravity Data
HuRong Duan, JunGang Guo, LingKang Chen, JiaShuang Jiao, and HeTing Jian
Accepted Articles, doi: 10.26464/epp2022034
[Abstract](186) [PDF](7)
The uplift of the Qinghai-Tibet Plateau (TP) strongly influences climate change, both regionally and globally. Surface observation data from this region have limited coverage and are difficult to obtain. Consequently, the vertical crustal deformation velocity (VCDV) distribution of the TP is poorly constrained. In this study, VCDV from the TP was inverted using data from the gravity recovery and climate experiment (GRACE). Based on the assumption that the gravity signal detected by GRACE is mainly composed of hydrological factors and vertical crustal movement, by deducting hydrological factors, the vertical crustal movement could be obtained. From this, the distribution of the vertical crustal deformation velocity across the TP was inverted. The results show that the VCDV of the southern, eastern, and northern TP is ~1.1 mm/a, ~0.5 mm/a, and −0.1 mm/a, respectively, while that of the region between the Qilian Haiyuan Fault and Kunlun Fault is ~0.0 mm/a. These results are consistent with the distributions of crustal deformation, thrust earthquakes/faults, and regional lithospheric activity. Hydrology, crustal thickness, and topographic factors do not change the overall distribution of the VCDV across the TP. The influence of hydrological factors is marked, with the maximum differences being approximately −0.4 mm/a in the northwest and 1.0 mm/a in the central area. The results of this study are significant for understanding the kinematics of the TP.
The anisotropy of hcp iron under inner core conditions: the effect of light elements
Wei Zhang, Mei Tang, ZhenWei Niu
Accepted Articles, doi: 10.26464/epp2022035
[Abstract](206) [PDF](8)
In the last decades, global seismic observations have identified increasingly complex anisotropy of the Earth’s inner core. Plenty of seismic studies have confirmed the anisotropy presents hemispherical variations in inner core. Here we report the effect of light elements on the anisotropy of hcp iron under inner core conditions based on ab initio molecular dynamics calculations. We found that light elements have significant effects on the density, sound velocities, and the anisotropy of hcp-Fe-X (X = C, O, Si, and S) binary alloy. For these binary alloys, the anisotropy depends on combined effects of temperature and the type of light element. Furthermore, there is a certain increase of the anisotropy strength with the increasing temperature. Alloying iron with some light elements such as C or O actually does not improve but reduces the strength of anisotropy of pure iron at high temperatures. Oppositely, light element S can significantly enhance the elastic anisotropy strength of hcp-Fe.
Interpretation of the locally high gravity anomalies using terrestrial gravity data in Bagodo, North Cameroon
Apollinaire Bouba, Kasi Njeudjang, Loudi Yap, Bouba Saido4, Joseph Kamguia, and Charles Tabod Tabod
Accepted Articles, doi: 10.26464/epp2022033
[Abstract](181) [PDF](22)
In this work, we interpreted the gravity data in order to determine the structural characteristics responsible for high gravity anomalies in Bagodo. These anomalies have not been characterized through a local study. Thus, regional-residual separation of the gravity anomalies using the polynomial method is undertaken. Geophysical signatures of near-surface small-extent geological structures are therefore revealed. In order to conduct the quantitative interpretation of the gravity anomalies, one profile was drawn on the residual Bouguer anomaly map and therefore was interpreted using spectral analysis, ideal body solution and 2.5 modeling methods. Our results showed that the intrusive body in Bagodo area consists of two trapezoids blocks. The first and the second block have roofs of about 7.5 and 14 km depths respectively while their bases are about 17km depths. The values are in accordance with those obtained by ideal body which shows two cells with density contrast of 0.3 gcm-3 in comparison to the surrounding rocks. Its density is estimated to about 3gcm-3. The topography of these rocks shows that they are basaltic rocks that would have cooled in fracture zones as an intrusion.
A machine learning-based electron density (MLED) model in the inner magnetosphere
QingHua Zhou, YunXiang Chen, FuLiang Xiao, Sai Zhang, Si Liu, Chang Yang, YiHua He, ZhongLei Gao
Accepted Articles, doi: 10.26464/epp2022036
[Abstract](156) [PDF](5)
The plasma density is an important factor in determining wave-particle interaction in the magnetosphere. We develop a machine learning-based electron density (MLED) model in the inner magnetosphere using a data set of electron density from Van Allen Probes between September 25, 2012 and August 30, 2019. This MLED model is a physics-based nonlinear network that employs fundamental physical principles to represent the variation of electron density. It predicts the plasmapause location under different geomagnetic activities and models the electron density of the plasmasphere and trough separately. We train the model using gradient descent and backpropagation algorithms, which can be widely used to deal effectively with nonlinear relationships among physical quantities in space plasma environments. The model gives explicit expressions with few parameters and describes the association of electron density with geomagnetic activity, solar cycle, and seasonal effects. Under various geomagnetic conditions, the electron densities calculated by this model agree well with the observation and provide a good description of the plasmapause movement. This MLED model, which can be easily incorporated into the previously developed radiation belt model, would be very helpful in modeling and forecasting the radiation belt electron dynamics.
Moonquake hazard distribution and the rolling trajectories of the lunar boulder: I. Inverse algorithm and validation of PGA
Sha Tao, Zhengfeng Zhang, Bei Zhang, Yaolin Shi, David A Yuen, Bojing Zhu
Accepted Articles, doi: 10.26464/epp2022030
[Abstract](192) [PDF](6)
The role of Moonquake hazard distribution (peak ground acceleration, PGA; instrumental intensity Moonquakes scale, IIM) on the co-seismic moment is critical for understanding the lunar quake (including the shallow and deep lunar quakes). Here, we use Lagrangian analysis methods and lunar surface image extraction-calibration technologies to investigate the relationship between the co-seismic induced ground motion (PGA and IIM) and the rolling trajectories of the lunar boulder on the lunar slope. First, the inverse algorithm, the validation of IIM (00.18g_m), the certain critical rolling condition between the PGA and the IIM, and the falling boulder are obtained. Then, the correlation between IIM scale (scale I, ~ 0 m⁄s^2 ; scale II, ~ 0.0017g_m; scale III, 0.0017g_m0.014g_m; scale IV, 0.014 g_m0.039g_m; scale V, 0.039 g_m0.092g_m; scale VI, 0.092g_m0.18g_m) and different parameters (e.g., lunar slope 1°44°; lasting time, t, 0.12s ; fixed PGA coefficient, the ratio of the vertical component to the horizontal component of PGA amplitude, A_y0⁄A_x0 =2) is simulated and discussed. All these results play an important role in understanding lunar quakes' mechanism and the internal structure evolution. The elevation extraction-calibration technologies with the particular software for automatically recognizing lunar-terrain observational images and data (e.g., high-resolution satellite lunar imagery from NASA's lunar reconnaissance orbiter narrow-angle camera) driven IIM with observational trajectories of the lunar boulder are presented in Part II and Part III of this three-paper series studies, respectively. We anticipate our studies to be a key point, give new insights for understanding the role of PGA and IIM, and help assess the stabilities of the rovers (e.g., Zhurong rover, Chang'e series rover) in China's Moon and Mars exploration program.

Display:          |     

Joint observation of the concentric gravity wave event on the Tibetan Plateau
Chang Lai, PengWei Li, JiYao Xu, Wei Yuan, Jia Yue, Xiao Liu, Kogure Masaru, LiLi Qian
2022, 6(3): 219 -227. doi: 10.26464/epp2022029
A concentric gravity wave event was captured by a photographer in Nagarzê County (90.28°N, 28.33°E) between 02:00 and 04:00 (local time) on May 11, 2019. This concentric gravity wave event was also observed by the Suomi National Polar-orbiting Partnership satellite and the all-sky airglow imager at Yangbajing station (90.5°E, 30.1°N). The temporal and spatial information on gravity waves from the photographs provided a rare opportunity to study the propagation of gravity waves over the Tibetan Plateau. According to wind and temperature data from the MERRA-2 reanalysis (Modern-Era Retrospective analysis for Research and Applications, Version 2) and empirical models (NRLMSISE-00 [Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Exosphere] and HWM [horizontal wind model]), we inversely derived the propagation trajectory from the observed wave pattern to the source region by using the ray-tracing method. The source of the concentric gravity wave was identified as deep convection in Bangladesh (90.6°E, 25.0°N). The maximum background wind speed in the propagation direction (31.05 m/s) was less than the phase speed of 53 m/s, which is consistent with the wind-filtering theory.
Runoff variations in the Yangtze River Basin and sub-basins based on GRACE, hydrological models, and in-situ data
WeiLong Rao, WenKe Sun
2022, 6(3): 228 -240. doi: 10.26464/epp2022021
Water budget closure is a method used to study the balance of basin water storage and the dynamics of relevant hydrological components (e.g., precipitation, evapotranspiration, and runoff). When water budget closure is connected with terrestrial water storage change (TWSC) estimated from Gravity Recovery and Climate Experiment (GRACE) data, variations in basin runoff can be understood comprehensively. In this study, total runoff variations in the Yangtze River Basin (YRB) and its sub-basins are examined in detail based on the water budget closure equation. We compare and combine mainstream precipitation and evapotranspiration models to determine the best estimate of precipitation minus evapotranspiration. In addition, we consider human water consumption, which has been neglected in earlier studies, and discuss its impact. To evaluate the effectiveness and accuracy of the combined hydrological models in estimating subsurface runoff, we collect discharge variations derived from in situ observations in the YRB and its sub-basins and compare these data with the models’ final estimated runoff variations. The estimated runoff variations suggest that runoff over the YRB has been increasing, especially in the lower sub-basins and in the post-monsoon season, and is accompanied by apparent terrestrial water loss.
Theoretical calculation of tidal Love numbers of the Moon with a new spectral element method
BinBin Liao, XiaoDong Chen, JianQiao Xu, JiangCun Zhou, HePing Sun
2022, 6(3): 241 -247. doi: 10.26464/epp2022025
The tidal Love numbers of the Moon are a set of nondimensional parameters that describe the deformation responses of the Moon to the tidal forces of external celestial bodies. They play an important role in the theoretical calculation of the Moon’s tidal deformation and the inversion of its internal structure. In this study, we introduce the basic theory for the theoretical calculation of the tidal Love numbers and propose a new method of solving the tidal Love numbers: the spectral element method. Moreover, we explain the mathematical theory and advantages of this method. On the basis of this new method, using 10 published lunar internal structure reference models, the lunar surface and lunar internal tidal Love numbers were calculated, and the influence of different lunar models on the calculated Love numbers was analyzed. Results of the calculation showed that the difference in the second-degree lunar surface Love numbers among different lunar models was within 8.5%, the influence on the maximum vertical displacement on the lunar surface could reach ±8.5 mm, and the influence on the maximum gravity change could reach ±6 μGal. Regarding the influence on the Love numbers inside the Moon, different lunar models had a greater impact on the Love numbers h2 and l2 than on k2 in the lower lunar mantle and core.
Neogene faulting and volcanism in the Victoria Land Basin of the Ross Sea, Antarctica
Mei Yue, JinYao Gao, ChunFeng Li, Chao Zhu, XinZhi Fan, Guochao Wu, ZhongYan Shen, Han Shi, XiaoXian Cai, YiDong Guo
2022, 6(3): 248 -258. doi: 10.26464/epp2022023
The Neogene Terror Rift in the Antarctic Victoria Land Basin (VLB) of the Ross Sea, Antarctica, is composed of the Discovery Graben and the Lee Arch. Many Neogene volcanoes are aligned in the north-south direction in the southern VLB, belonging to the McMurdo Volcanic Group. However, due to multiple glaciations and limited seismic data, the volcanic processes are still unclear in the northern VLB, especially in the Terror Rift. Multichannel seismic profiles were collected at the VLB from the 32nd Chinese National Antarctic Research Expedition (CHINARE). We utilized four seismic profiles from the CHINARE and additional historical profiles, along with gravity and magnetic anomalies, to analyze faults and stratigraphic characteristics in the northern Terror Rift and volcanism in the VLB. Negative flower structures found in the northern Terror Rift suggest that the Terror Rift was affected by dextral strike-slip faults extending from the northern Victoria Land (NVL). After the initial orthogonal tension, the rift transited into an oblique extension, forming a set of downward concaving normal faults and accommodation zones in the Terror Rift. On the Lee Arch, several imbricated normal faults formed and converged into a detachment fault. Under gravitational forces, the strata bent upward and formed a rollover anticline. Many deep faults and thin strata subjected to erosion facilitated volcanic activity. A brittle volcanic region in the VLB was affected by dextral strike-slip movements and east-west extension, resulting in two Neogene volcanic chains that connect three igneous provinces in the VLB: the Hallett, Melbourne, and Erebus Provinces. These two chains contain mud volcanoes with magnetic nuclei, volcanic intrusions, and late-stage volcanic eruptions. Volcanisms have brought about opposite polarities of magnetic anomalies in Antarctica, indicating the occurrence of multiple volcanic activities.
Seismic attenuation compensation with spectral-shaping regularization
QiZhen Du, WanYu Wang, WenHan Sun, Li-Yun Fu
2022, 6(3): 259 -274. doi: 10.26464/epp2022024
Because of the viscoelasticity of the subsurface medium, seismic waves will inherently attenuate during propagation, which lowers the resolution of the acquired seismic records. Inverse-Q filtering, as a typical approach to compensating for seismic attenuation, can efficiently recover high-resolution seismic data from attenuation. Whereas most efforts are focused on compensating for high-frequency energy and improving the stability of amplitude compensation by inverse-Q filtering, low-frequency leakage may occur as the high-frequency component is boosted. In this article, we propose a compensation scheme that promotes the preservation of low-frequency energy in the seismic data. We constructed an adaptive shaping operator based on spectral-shaping regularization by tailoring the frequency spectra of the seismic data. We then performed inverse-Q filtering in an inversion scheme. This data-driven shaping operator can regularize and balance the spectral-energy distribution for the compensated records and can maintain the low-frequency ratio by constraining the overcompensation for high-frequency energy. Synthetic tests and applications on prestack common-reflection-point gathers indicated that the proposed method can preserve the relative energy of low-frequency components while fulfilling stable high-frequency compensation.
On the source of the quasi-Carrington Rotation periodic magnetic variations on the Martian surface: InSight observations and modeling
Hao Luo, AiMin Du, ShaoHua Zhang, YaSong Ge, Ying Zhang, ShuQuan Sun, Lin Zhao, Lin Tian, SongYan Li
2022, 6(3): 275 -283. doi: 10.26464/epp2022022
In a recent paper (Luo H et al., 2022), we found that the peak amplitudes of diurnal magnetic variations, measured during martian days (sols) at the InSight landing site, exhibited quasi Carrington-Rotation (qCR) periods at higher eigenmodes of the natural orthogonal components (NOC); these results were based on ~664 sols of magnetic field measurements. However, the source of these periodic variations is still unknown. In this paper we introduce the neutral-wind driven ionospheric dynamo current model (e.g., Lillis et al., 2019) to investigate the source. Four candidates — the draped IMF, electron density/plasma density, the neutral densities, and the electron temperature in the ionosphere with artificial qCR periodicity, are applied in the modeling to find the main factor likely to be causing the observed surface magnetic field variations that exhibit the same qCR periods. Results show that the electron density/plasma density, which controls the total conductivity in the dynamo region, appears to account for the greatest part of the surface qCR variations; its contribution reaches about 67.6%. The draped IMF, the neutral densities, and the electron temperature account, respectively, for only about 12.9%, 10.3%, and 9.2% of the variations. Our study implies that the qCR magnetic variations on the Martian surface are due primarily to variations of the dynamo currents caused by the electron density variations. We suggest also that the time-varying fields with the qCR period could be used to probe the Martian interior's electrical conductivity structure to a depth of at least 700 km.
A two-dimensional energy balance climate model on Mars
YaoKun Li, JiPing Chao
2022, 6(3): 284 -293. doi: 10.26464/epp2022026
A two-dimensional energy balance climate model has been built to investigate the climate on Mars. The model takes into account the balance among solar radiation, longwave radiation, and energy transmission and can be solved analytically by Legendre polynomials. With the parameters for thermal diffusion and radiation processes being properly specified, the model can simulate a reasonable surface atmospheric temperature distribution but not a very perfect vertical atmospheric temperature distribution compared with numerical results, such as those from the Mars Climate Database. With varying solar radiation in a Martian year, the model can simulate the seasonal variation of the air temperature on Mars. With increasing dust content, the Martian atmosphere gradually warms. However, the warming is insignificant in the cold and warm scenarios, in which the dust mixing ratio varies moderately, whereas the warming is significant in the storm scenario, in which the dust mixing ratio increases dramatically. With an increasing albedo value of either the polar cap or the non-ice region, Mars gradually cools. The mean surface atmospheric temperature decreases moderately with an increasing polar ice albedo, whereas it increases dramatically with an increasing non-ice albedo. This increase occurs because the planetary albedo of the ice regions is smaller than that of the non-ice region.
Groove formation on Phobos from reimpacting orbital ejecta of the Stickney crater
XiangYu Xi, Min Ding, MengHua Zhu
2022, 6(3): 294 -303. doi: 10.26464/epp2022027
Numerous linear grooves have long been recognized as covering the surface of Phobos, but the mechanisms of their formation are still unclear. One possible mechanism is related to the largest crater on Phobos, the Stickney crater, whose impact ejecta may slide, roll, bounce, and engrave groove-like features on Phobos. When the launch velocity is higher than the escape velocity, the impact ejecta can escape Phobos. A portion of these high-velocity ejecta are dragged by the gravitational force of Mars, fall back, and reimpact Phobos. In this research, we numerically test the hypothesis that the orbital ejecta of the Stickney crater that reimpact Phobos could be responsible for a particular subset of the observed grooves on Phobos. We adopt impact hydrocode iSALE-2D (impact-Simplified Arbitrary Lagrangian Eulerian, two-dimensional) to simulate the formation of the Stickney crater and track its impact ejecta, with a focus on orbital ejecta with launch velocities greater than the escape velocity of Phobos. The launch velocity distribution of the ejecta particles is then used to calculate their trajectories in space and determine their fates. For orbital ejecta reimpacting Phobos, we then apply the sliding boulder model to calculate the ejecta paths, which are compared with the observed groove distribution and length to search for causal relationships. Our ejecta trajectory calculations suggest that only ~1% of the orbital ejecta from the Stickney crater can reimpact Phobos. Applying the sliding boulder model, we predict ejecta sliding paths of 9−20 km in a westward direction to the east of the zone of avoidance, closely matching the observed grooves in that region. The best-fit model assumes an ejecta radius of ~150 m and a speed restitution coefficient of 0.3, consistent with the expected ejecta and regolith properties. Our calculations thus suggest the groove class located to the east of the zone of avoidance may have been caused by reimpact orbital ejecta from the Stickney crater.
show more results
A brief review of equatorial ionization anomaly and ionospheric irregularities
Nanan Balan, LiBo Liu, HuiJun Le
2018, 2(4): 257-275doi: 10.26464/epp2018025
Corotating drift-bounce resonance of plasmaspheric electron with poloidal ULF waves
Qiu-Gang Zong, YongFu Wang, Jie Ren, XuZhi Zhou, SuiYan Fu, Robert Rankin, Hui Zhang
2017, 1(1): 2-12doi: 10.26464/epp2017002
Ambient noise surface wave tomography of marginal seas in east Asia
Qing Wang, XiaoDong Song, JianYe Ren
2017, 1(1): 13-25doi: 10.26464/epp2017003
Thermal structures of the Pacific lithosphere from magnetic anomaly inversion
Chun-Feng Li, Jian Wang
2018, 2(1): 52-66doi: 10.26464/epp2018005
A seismic model for crustal structure in North China Craton
TianYu Zheng, YongHong Duan, WeiWei Xu, YinShuang Ai
2017, 1(1): 26-34doi: 10.26464/epp2017004
The first joint experimental results between SURA and CSES
XueMin Zhang, Vladimir Frolov, ShuFan Zhao, Chen Zhou, YaLu Wang, Alexander Ryabov, DuLin Zhai
2018, 2(6): 527-537doi: 10.26464/epp2018051
Reviewing subduction initiation and the origin of plate tectonics: What do we learn from present-day Earth?
Gang Lu, Liang Zhao, Ling Chen, Bo Wan, FuYuan Wu
2021, 5(2): 123-140doi: 10.26464/epp2021014
Different earthquake patterns for two neighboring fault segments within the Haiyuan Fault zone
ZhiKun Ren, ZhuQi Zhang, PeiZhen Zhang
2018, 2(1): 67-73doi: 10.26464/epp2018006
A simulation study of 630 nm and 557.7 nm airglow variations due to dissociative recombination and thermal electrons by high-power HF heating
Tong Dang, JiuHou Lei, XianKang Dou, WeiXing Wan
2017, 1(1): 44-52doi: 10.26464/epp2017006
Probabilistic seismic hazard assessment of Nepal using multiple seismic source models
Md Moklesur Rahman, Ling Bai
2018, 2(4): 327-341doi: 10.26464/epp2018030

Submission Log In

Forgot your password?

Enter your e-mail address to

receive your account information.

Current Issue

Year 2022

Volume 6

Issue 3




All Issues

Supported by Beijing Renhe Information Technology Co. LtdE-mail: