In addition to being driven by tidal winds,the sporadic E(Es)layers are modulated by gravity waves(GWs),although the effects are not yet comprehensively understood.In this article,we discuss the effects of mesoscale G...In addition to being driven by tidal winds,the sporadic E(Es)layers are modulated by gravity waves(GWs),although the effects are not yet comprehensively understood.In this article,we discuss the effects of mesoscale GWs on the Es layers determined by using a newly developed model,MISE-1D(one-dimensional Model of Ionospheric Sporadic E),with low numerical dissipation and high resolution.Driven by the wind fields resolved by the high-resolution version of the Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension(WACCM-X),the MISE-1D simulation revealed that GWs significantly influence the evolution of the Es layer above 100 km but have a very limited effect at lower altitudes.The effects of GWs are diverse and complex,generally including the generation of fluctuating wavelike structures on the Es layer with frequencies similar to those of the GWs.The mesoscale GWs can also cause increases in the density of Es layers,or they can disperse or diffuse the Es layers and increase their thickness.In addition,the presence of GWs is a key factor in sustaining the Es layers in some cases.展开更多
Sporadic E(Es)layers in the ionosphere are characterized by intense plasma irregularities in the E region at altitudes of 90-130 km.Because they can significantly influence radio communications and navigation systems,...Sporadic E(Es)layers in the ionosphere are characterized by intense plasma irregularities in the E region at altitudes of 90-130 km.Because they can significantly influence radio communications and navigation systems,accurate forecasting of Es layers is crucial for ensuring the precision and dependability of navigation satellite systems.In this study,we present Es predictions made by an empirical model and by a deep learning model,and analyze their differences comprehensively by comparing the model predictions to satellite RO measurements and ground-based ionosonde observations.The deep learning model exhibited significantly better performance,as indicated by its high coefficient of correlation(r=0.87)with RO observations and predictions,than did the empirical model(r=0.53).This study highlights the importance of integrating artificial intelligence technology into ionosphere modelling generally,and into predicting Es layer occurrences and characteristics,in particular.展开更多
The algebraic reconstruction technique(ART),multiplicative algebraic reconstruction technique(MART),and simultaneous iterative reconstruction technique(SIRT)are computational methodologies extensively utilized within ...The algebraic reconstruction technique(ART),multiplicative algebraic reconstruction technique(MART),and simultaneous iterative reconstruction technique(SIRT)are computational methodologies extensively utilized within the field of computerized ionospheric tomography(CIT)to facilitate three-dimensional reconstruction of the ionospheric morphology.However,reconstruction accuracy elicits recurrent disputes over its practical application,and people usually attribute this issue to incomplete and uneven coverage of the measurements.The Thermosphere Ionosphere Electrodynamics General Circulation Model(TIEGCM)offers a reasonable physics-based ionospheric background and is widely utilized in ionospheric research.We use the TIEGCM simulations as the targeted ionosphere because the current measurements are far from able to realistically reproduce the ionosphere in detail.Optimized designations of satellite measurements are conducted to investigate the limiting performance of CIT methods in ionospheric reconstruction.Similar to common practice,electron density distributions from outputs of the International Reference Ionosphere(IRI)model are used as the iterative initial value in CIT applications.The outcomes suggest that despite data coverage,iterative initial conditions also play an essential role in ionospheric reconstruction.In particular,in the longitudinal sectors where the iterative initial height of the F2-layer peak electron density(hmF2)differs substantially from the background densities,none of the three CIT methods can reproduce the exact background profile.When hmF2 is close but the ionospheric F2-layer peak density(NmF2)is different between the targeted background and initial conditions,the MART performs better than the ART and SIRT,as evidenced by the correlation coefficients of MART being above 0.97 and those of ART and SIRT being below 0.85.In summary,this investigation reveals the potential uncertainties in traditional CIT reconstruction,particularly when realistic hmF2 or NmF2 values differ substantially from the initial CIT conditions.展开更多
The formation of an embedded electron current sheet within the magnetotail plasma sheet has been poorly understood.In this article,we present an electron current layer detected at the edge of the magnetotail plasma sh...The formation of an embedded electron current sheet within the magnetotail plasma sheet has been poorly understood.In this article,we present an electron current layer detected at the edge of the magnetotail plasma sheet.The ions were demagnetized inside the electron current layer,but the electrons were still frozen in with the magnetic field line.Thus,this decoupling of ions and electrons gave rise to a strong Hall electric field,which could be the reason for the formation of the embedded thin current layer.The magnetized electrons,the absence of the nongyrotropic electron distribution,and negligible energy dissipation in the layer indicate that magnetic reconnection had not been triggered within the embedded thin current layer.The highly asymmetric plasma on the two sides of the current layer and low magnetic shear across it could suppress magnetic reconnection.The observations indicate that the embedded electric current layer,probably generated by the Hall electric field,even down to electron scale,is not a sufficient condition for magnetic reconnection.展开更多
As a major fault in the northeastern Qinghai-Xizang Plateau,the Haiyuan fault zone is important for understanding the regional deformation.Aiming at the differences in the slip rate and locking degree obtained from di...As a major fault in the northeastern Qinghai-Xizang Plateau,the Haiyuan fault zone is important for understanding the regional deformation.Aiming at the differences in the slip rate and locking degree obtained from different studies,this study constructs a refined block model(including Qilian,Alxa,Ordos,Xining,Haiyuan,and Lanzhou blocks)and uses the grid search and simulated annealing methods to invert GPS data for slip rate and locking degree of the Haiyuan fault zone.The results are as follows:(1)The sinistral slip rates in the western,middle,and eastern segments are 4.93-5.22 mm/a,1.52-4.94 mm/a,and 0.43-1.18 mm/a,decreasing eastward on the whole,while the compression rates are 0.45-1.26 mm/a,0.58-2.62 mm/a,and3.52-4.48 mm/a,increasing eastward on the whole.(2)The locking depth of the western segment increases from about 5 km to about 20 km eastward;the middle segment decreases and then increases eastward;the eastern segment concentrates at about 20 km(PHI is about 0.86).(3)The slip deficit is relatively higher in the Lenglongling,Jinqianghe,Maomaoshan,and Liupanshan faults(averaging about 3.42 mm/a,4.16 mm/a,4.23 mm/a,and 3.43 mm/a within 20 km).(4)The Qilian,Alxa,Xining,Lanzhou,and Haiyuan blocks rotate clockwise,while the Ordos block rotates counterclockwise.Additionally,by comparing different block models,the Haiyuan block should be considered independently.The Haiyuan fault zone adjusts surrounding block movements and uplifts Liupanshan mountain tectonically.The results can provide important references for understanding the regional earthquake risk and deformation mechanism.展开更多
Electrostatic ion-cyclotron wave(EICW)and ion-Bernstein wave(IBW)are two typical normal modes of plasma propagating vertically to a magnetic field.These modes satisfy the extreme conditionsω《k_(||)v_(te)andω》k_(||...Electrostatic ion-cyclotron wave(EICW)and ion-Bernstein wave(IBW)are two typical normal modes of plasma propagating vertically to a magnetic field.These modes satisfy the extreme conditionsω《k_(||)v_(te)andω》k_(||)v_(te),whereω,k_(||),and v_(te)are the wave angular frequency,wave number parallel to the magnetic field,and the electron thermal speed.Between these extremities,there exist IBW-like waves that have not been studied.In this paper,the co-excitation of the EICW and IBW-like wave in an experiment is presented.It is shown that,in an argon(Ar)discharge,the EICW appears in the fundamental and second harmonic frequency bands of the ioncyclotron frequency.In contrast,the IBW-like wave appears in the higher harmonic frequency bands and in the case of a relatively strong magnetic field.In a helium(He)discharge,the IBW-like wave is excited in the fundamental frequency band,but the wave behaviour becomes complicated in the harmonic frequency bands.The measured dispersion relations of the excited EICW are in agreement with the results of the fluid and kinetic models in some range of frequencies,but those of the excited IBW-like wave are only in qualitative agreement with the kinetic model of the IBW.This study is important for identifying the IBW-like mode that may have potential application in plasma heating via Landau damping.展开更多
The solar wind's interaction with the Moon has traditionally been understood through the Moon's absorption of solar wind particles and the formation of a plasma cavity on its nightside,known as the lunar wake....The solar wind's interaction with the Moon has traditionally been understood through the Moon's absorption of solar wind particles and the formation of a plasma cavity on its nightside,known as the lunar wake.This study reveals unexpected,large-scale perturbations in the solar wind upstream of the Moon,using 11 years of data from the OMNI and ARTEMIS(Acceleration,Reconnection,Turbulence and Electrodynamics of Moon's Interaction with the Sun)missions(2012-2023).We find systematic moonward deviations of~tens of km/s in a direction perpendicular to the solar wind(moonward),at altitudes of up to 1000 km,particularly when the interplanetary magnetic field(IMF)lines are oblique to the solar wind(30°<θ<60°)and connected to the lunar dayside.The longer the duration of the interaction,the greater the moonward deviation.These perturbations can be explained by neither solar wind pickup of the reflected ions,nor lunar wake dynamics.Instead,they appear to correlate with magnetic connectivity between the ARTEMIS probes and the lunar surface,suggesting a more complex solar wind interaction than previously thought.展开更多
Gravity Recovery and Climate Experiment(GRACE) observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 2...Gravity Recovery and Climate Experiment(GRACE) observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 26,2004.This article focuses on investigating some gravita-tional effects caused by this huge earthquake.We computed the geoid height changes,the equivalent water height(EWH) changes,and the gravity changes using the GRACE Level-2 monthly spherical harmonic(SH) solutions released by University of Texas Center for Space Research(UTCSR).The GRACE results agree well with the prediction by a dislocation model and are consistent with the results obtained by some previous scholars.In particular,we calculated the three components of the gravity gradient variations and found that they can recover the seismic-related signature more sensitively due to a certain degree of amplification of the signals.A positive-negative-positive mode predominates in the spatial distribution of the horizontal components of the gravity gradient variations,which is possibly attributed to the anomalies in the crustal density distribution caused by the uplift-subduction effect of the dip-slip earthquake.Moreover,the latitude components of the gravity gradient changes show strong suppression of the north-south stripes,which is due to the along-orbit measurements of the two GRACE satellites.We conclude that the posi-tive-negative-positive mode in latitude gravity gradient changes would be a more sensitive fea-ture to detect the deformations of some major dip-slip earthquakes by GRACE data.展开更多
As one of the seven scientific payloads on board the Tianwen-1 orbiter,the Mars Orbiter Magnetometer(MOMAG)will measure the magnetic fields of and surrounding Mars to study its space environment and the interaction wi...As one of the seven scientific payloads on board the Tianwen-1 orbiter,the Mars Orbiter Magnetometer(MOMAG)will measure the magnetic fields of and surrounding Mars to study its space environment and the interaction with the solar wind.The instrument consists of two identical triaxial fluxgate magnetometer sensors,mounted on a 3.19 meter-long boom with a seperation of about 90 cm.The dual-magnetometers configuration will help eliminate the magnetic field interference generated by the spacecraft platform and payloads.The sensors are controlled by an electric box mounted inside the orbiter.Each magnetometer measures the ambient vector magnetic field over a wide dynamic range(to 10,000 nT per axis)with a resolution of 1.19 pT.Both magnetometers sample the ambient magnetic field at an intrinsic frequency of 128 Hz,but will operate in a model with alternating frequency between 1 and 32 Hz to meet telemetry allocations.展开更多
Based on gravity frequency shift effect predicted by general relativity theory, this study discusses an approach for determining the gravity potential(geopotential) difference between arbitrary two points P and Q by r...Based on gravity frequency shift effect predicted by general relativity theory, this study discusses an approach for determining the gravity potential(geopotential) difference between arbitrary two points P and Q by remote comparison of two precise optical clocks via optical fiber frequency transfer. After synchronization, by measuring the signal's frequency shift based upon the comparison of bidirectional frequency signals from P and Q oscillators connected with two optical atomic clocks via remote optical fiber frequency transfer technique, the geopotential difference between the two points could be determined, and its accuracy depends on the stabilities of the optical clocks and the frequency transfer comparison technique. Due to the fact that the present stability of optical clocks achieves 1.6×10-18 and the present frequency transfer comparison via optical fiber provides stabilities as high as 10-19 level, this approach is prospective to determine geopotential difference with an equivalent accuracy of 1.5 cm. In addition, since points P and Q are quite arbitrary, this approach may provide an alternative way to determine the geopotential over a continent, and prospective potential to unify a regional height datum system.展开更多
Under the paraxial approximation, the analytical propagation expression of an Airy–Gaussian beam(Ai GB) in uniaxial crystals orthogonal to the optical axis is investigated. The propagation dynamics of the Ai GB is ...Under the paraxial approximation, the analytical propagation expression of an Airy–Gaussian beam(Ai GB) in uniaxial crystals orthogonal to the optical axis is investigated. The propagation dynamics of the Ai GB is given for different ratios of the extraordinary index to the ordinary refractive index. It has been found that the continuity and the self-bending effect of Ai GB become weaker when the ratio increases. From the figure of the maximum intensity of Ai GB, one can see that the maximum intensity is not monotone decreasing due to the anisotropic effect of the crystals. The intensity distribution of Ai GB in different distribution factors is shown. The Ai GB converges toward a Gaussian beam as the distribution factor increases.展开更多
The propagation dynamics of the Airy Gaussian vortex beams in uniaxial crystals orthogonal to the optical axis has been investigated analytically and numerically. The propagation expression of the beams has been obtai...The propagation dynamics of the Airy Gaussian vortex beams in uniaxial crystals orthogonal to the optical axis has been investigated analytically and numerically. The propagation expression of the beams has been obtained. The propagation features of the Airy Gaussian vortex beams are shown with changes of the distribution factor and the ratio of the extraordinary refractive index to the ordinary refractive index. The correlations between the ratio and the maximum intensity value during the propagation, and its appearing distance have been investigated.展开更多
Lushan Earthquake (-Mw 6.6) occurred in Sichuan Province of China on 20 April 2013, was the largest earthquake in Longmenshan fault belt since 2008 Wenchuan Earthquake. To better understand its rupture pattern, we f...Lushan Earthquake (-Mw 6.6) occurred in Sichuan Province of China on 20 April 2013, was the largest earthquake in Longmenshan fault belt since 2008 Wenchuan Earthquake. To better understand its rupture pattern, we focused on the influences of fault parameters on fault slips and performed fault slip inversion using Akaike's Bayesian Information Criterion (ABIC) method. Based on GPS coseismic data, our inverted results showed that the fault slip was mainly confined at depths. The maximum slip amplitude is about 0.7 m, and the scalar seismic moment is about 9.47x10TM N.m. Slip pattern reveals that the earthquake occurred on the thrust fault with large dip-slip and small strike-slip, such a simple fault slip represents no second sub-event occurred. The Coulomb stress changes (ACFF) matched the most aftershocks with negative anomalies. The in- verted results demonstrated that the source parameters have significant impacts on fault slip distri- bution, especially on the slip direction and maximum displacement.展开更多
Scientists pay great attention to different-time-scale signals in the lengllh of day (LOD) variations △LOD, which provide signatures of the Earth's interior structure, couplings among different layers, and potenti...Scientists pay great attention to different-time-scale signals in the lengllh of day (LOD) variations △LOD, which provide signatures of the Earth's interior structure, couplings among different layers, and potential excitations of ocean and atmosphere. In this study, based on the ensemble empirical mode decomposition (EEMD), we analyzed the latest time series of △LOD data spanning from January 1962 to March 2015. We observed the signals with periods and amplitudes of about 0.5 month and 0.19 ms, 1.0 month and 0.19 ms, 0.5 yr and 0.22 ms, 1.0 yr and 0.18 ms, 2.28 yr and 0.03 ms, 5.48 yr and 0.05 ms, respectively, in coincidence with the results of predecessors. In addition, some signals that were previously not definitely observed by predecessors were detected in this study, with periods and amplitudes of 9.13 d and 0.12 ms, 13.69 yr and 0.10 ms, respectively. The mechanisms of the LOD fluctuations of these two signals are still open.展开更多
Previous studies have reached consensus that low velocity zones are widespread in the crust beneath Yunnan region. However, the relationships between the low velocity zones and large faults, earthquake distribution ar...Previous studies have reached consensus that low velocity zones are widespread in the crust beneath Yunnan region. However, the relationships between the low velocity zones and large faults, earthquake distribution are less investigated by available studies. By analyzing the seismic ambient noise recorded by Yunnan Seismic Networks and Tengchong volcano array, we construct a 3D crustal shear wave velocity model for the Yunnan region, which provides more details of the distribution of intra-crustal low velocity zones all over Yunnan. The distribution of low velocity zones shows different features at different depths. At shallow depths, the results are well correlated with near surface geological features. With increasing depth, the low velocity zones are gradually concentrated on the northern part of our study area, most likely reflecting variations in crustal thickness beneath the Yunnan region. The low velocity zones are truncated at depth by several large faults in Yunnan. It is interesting that most strong earthquakes(Ms≥5.0) occurred in Yunnan are distributed in low velocity zones or the transition zones between low and high velocity anomalies within the upper-to-middle crust. The crustal structure is composed of a brittle, seismically active upper-to-middle crust and a warm, aseismic lower crust.展开更多
Surface mass anomalies estimated by mass concentration(mascon)approach using Gravity Recovery and Climate Experiment(GRACE)observations with regularization constraints generally present higher spatial resolution than ...Surface mass anomalies estimated by mass concentration(mascon)approach using Gravity Recovery and Climate Experiment(GRACE)observations with regularization constraints generally present higher spatial resolution than the spheric harmonic(SH)solutions.To analyze the influence of different types of constraints on the estimation of mascon solutions,we carried out a closed-loop simulation experiment to estimate surface mass anomalies over South America based on simulated GRACE intersatellite geopotential differences.Tikhonov regularization with spatial constraint(SC),uniform weighting constraint(UWC),and a prior information constraint(APC)were employed to stabilize the mascon solutions,and the corresponding optimal regularization parameters were determined based on the minimum residual root-mean-square(RMS)criterion.The results show that mascon solutions estimated under different types of constraints are consistent and equivalent when the optimal regularization parameters are selected.The spatial distributions and main characteristics of regional surface mass anomalies estimated by the three types of constraints agree well,and the values of residual RMS with different constraints are very close.But due to the smoothing effect of regularization,the signal strength of mascon solutions is a bit weaker than that of original true signal,especially in the regions with strong signals.In addition,due to the ill-conditioned problem is more serious for higher grid resolution,the relative contribution of the three types of constraints to the final mascon solutions would be stronger.The results show that the averages of relative contribution percentages of these constraints for 2°×2° mascon grids are 80%-90%,while the corresponding values for 4°×4° mascon grids are 30%-60%.However,based on the minimum residual RMS criterion,the accuracy of estimation results is not affected by the type of constraints and their relative contribution to the final mascon solutions.展开更多
Readout electronics is developed for a prototype time-of-flight(TOF) ion composition spectrometer for in situ measurement of the mass/charge distributions of major ion species from 200 to 100 ke V/e in space plasma.By...Readout electronics is developed for a prototype time-of-flight(TOF) ion composition spectrometer for in situ measurement of the mass/charge distributions of major ion species from 200 to 100 ke V/e in space plasma.By utilizing a constant fraction discriminator(CFD) and time-to-digital converter(TDC), challenging dynamic range measurements were performed with high time resolution and event rates. CFD was employed to discriminate the TOF signals from the micro-channel plate and channel electron multipliers. TDC based on the combination of counter and OR-gate delay chain was designed in a highreliability flash field programmable gate array. Owing to the non-uniformity of the delay chain, a correction algorithm based on integral nonlinearity compensation was implemented to reduce the time uncertainty. The test results showed that the electronics achieved a low timingerror of < 200 ps in the input range from 35 to 500 m V for the CFD, and a time resolution of ~550 ps with time uncertainty < 180 ps after correction and a time range of6.4 ls for the TDC. The TOF spectrum from an electron beam experiment of the impacting N_2 gas further indicated the good performance of this readout electronic.展开更多
The space-borne fluxgate magnetometer(FGM)requires regular in-flight calibration to obtain its zero offset.Recently,Wang GQ and Pan ZH(2021a)developed a new method for the zero offset calibration based on the properti...The space-borne fluxgate magnetometer(FGM)requires regular in-flight calibration to obtain its zero offset.Recently,Wang GQ and Pan ZH(2021a)developed a new method for the zero offset calibration based on the properties of Alfvén waves.They found that an optimal offset line(OOL)exists in the offset cube for a pure Alfvén wave and that the zero offset can be determined by the intersection of at least two nonparallel OOLs.Because no pure Alfvén waves exist in the interplanetary magnetic field,calculation of the zero offset relies on the selection of highly Alfvénic fluctuation events.Here,we propose an automatic procedure to find highly Alfvénic fluctuations in the solar wind and calculate the zero offset.This procedure includes three parts:(1)selecting potential Alfvénic fluctuation events,(2)obtaining the OOL,and(3)determining the zero offset.We tested our automatic procedure by applying it to the magnetic field data measured by the FGM onboard the Venus Express.The tests revealed that our automatic procedure was able to achieve results as good as those determined by the Davis-Smith method.One advantage of our procedure is that the selection criteria and the process for selecting the highly Alfvénic fluctuation events are simpler.Our automatic procedure could also be applied to find fluctuation events for the Davis-Smith method.展开更多
As critical component of hydrologic cycle, basin discharge is a key issue for understanding the hydrological and climatologic related to water and energy cycles. Combining GRACE gravity field models with ET from GLDAS...As critical component of hydrologic cycle, basin discharge is a key issue for understanding the hydrological and climatologic related to water and energy cycles. Combining GRACE gravity field models with ET from GLDAS models and precipitation from GPCP, discharge of the Yellow River basin are estimated from the water balance equation. While comparing the results with discharge from GLDAS model and in situ measurements, the results reveal that discharge from Mosaic and CLM GLDAS model can partially represent the river discharge and the discharge estimation from water balance equation could reflect the discharge from precipitation over the Yellow River basin.展开更多
Urban faults in Shenzhen are potential threats to city security and sustainable development. In consideration of the importance of the Shenzhen fault zone, the author provide a detailed interpretation on gravity data ...Urban faults in Shenzhen are potential threats to city security and sustainable development. In consideration of the importance of the Shenzhen fault zone, the author provide a detailed interpretation on gravity data model. Bouguer gravity covering the whole Shenzhen City was calculated with a 1-km resolution. Wavelet multi-scale analysis(MSA) was applied to the Bouguer gravity data to obtain the multilayer residual anomalies corresponding to different depths. In addition, 2D gravity models were constructed along three profiles. The Bouguer gravity anomaly shows an NE-striking high-low-high pattern from northwest to southeast, strongly related to the main faults. According to the results of MSA, the correlation between gravity anomaly and faults is particularly significant from 4 to 12 km depth. The residual gravity with small amplitude in each layer indicates weak tectonic activity in the crust. In the upper layers, positive anomalies along most of faults reveal the upwelling of high-density materials during the past tectonic movements. The multilayer residual anomalies also yield important information about the faults, such as the vertical extension and the dip direction. The maximum depth of the faults is about 20 km. In general, NE-striking faults extend deeper than NW-striking faults and have a larger dip angle.展开更多
基金supported by the Project of Stable Support for Youth Teams in Basic Research Field,Chinese Academy of Sciences(CASGrant No.YSBR-018)+2 种基金the B-type Strategic Priority Program of CAS(Grant No.XDB41000000)the National Natural Science Foundation of China(Grant No.42204165)the National Key Research and Development Program(Grant No.2022YFF0504400).
文摘In addition to being driven by tidal winds,the sporadic E(Es)layers are modulated by gravity waves(GWs),although the effects are not yet comprehensively understood.In this article,we discuss the effects of mesoscale GWs on the Es layers determined by using a newly developed model,MISE-1D(one-dimensional Model of Ionospheric Sporadic E),with low numerical dissipation and high resolution.Driven by the wind fields resolved by the high-resolution version of the Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension(WACCM-X),the MISE-1D simulation revealed that GWs significantly influence the evolution of the Es layer above 100 km but have a very limited effect at lower altitudes.The effects of GWs are diverse and complex,generally including the generation of fluctuating wavelike structures on the Es layer with frequencies similar to those of the GWs.The mesoscale GWs can also cause increases in the density of Es layers,or they can disperse or diffuse the Es layers and increase their thickness.In addition,the presence of GWs is a key factor in sustaining the Es layers in some cases.
基金supported by the Project of Stable Support for Youth Team in Basic Research Field,CAS(grant No.YSBR-018)the National Natural Science Foundation of China(grant Nos.42188101,42130204)+4 种基金the B-type Strategic Priority Program of CAS(grant no.XDB41000000)the National Natural Science Foundation of China(NSFC)Distinguished Overseas Young Talents Program,Innovation Program for Quantum Science and Technology(2021ZD0300301)the Open Research Project of Large Research Infrastructures of CAS-“Study on the interaction between low/mid-latitude atmosphere and ionosphere based on the Chinese Meridian Project”.The project was supported also by the National Key Laboratory of Deep Space Exploration(Grant No.NKLDSE2023A002)the Open Fund of Anhui Provincial Key Laboratory of Intelligent Underground Detection(Grant No.APKLIUD23KF01)the China National Space Administration(CNSA)pre-research Project on Civil Aerospace Technologies No.D010305,D010301.
文摘Sporadic E(Es)layers in the ionosphere are characterized by intense plasma irregularities in the E region at altitudes of 90-130 km.Because they can significantly influence radio communications and navigation systems,accurate forecasting of Es layers is crucial for ensuring the precision and dependability of navigation satellite systems.In this study,we present Es predictions made by an empirical model and by a deep learning model,and analyze their differences comprehensively by comparing the model predictions to satellite RO measurements and ground-based ionosonde observations.The deep learning model exhibited significantly better performance,as indicated by its high coefficient of correlation(r=0.87)with RO observations and predictions,than did the empirical model(r=0.53).This study highlights the importance of integrating artificial intelligence technology into ionosphere modelling generally,and into predicting Es layer occurrences and characteristics,in particular.
基金supported by the National Key R&D Program of China (Grant No. 2022YFF0503702)the National Natural Science Foundation of China (Grant No. 42074186)
文摘The algebraic reconstruction technique(ART),multiplicative algebraic reconstruction technique(MART),and simultaneous iterative reconstruction technique(SIRT)are computational methodologies extensively utilized within the field of computerized ionospheric tomography(CIT)to facilitate three-dimensional reconstruction of the ionospheric morphology.However,reconstruction accuracy elicits recurrent disputes over its practical application,and people usually attribute this issue to incomplete and uneven coverage of the measurements.The Thermosphere Ionosphere Electrodynamics General Circulation Model(TIEGCM)offers a reasonable physics-based ionospheric background and is widely utilized in ionospheric research.We use the TIEGCM simulations as the targeted ionosphere because the current measurements are far from able to realistically reproduce the ionosphere in detail.Optimized designations of satellite measurements are conducted to investigate the limiting performance of CIT methods in ionospheric reconstruction.Similar to common practice,electron density distributions from outputs of the International Reference Ionosphere(IRI)model are used as the iterative initial value in CIT applications.The outcomes suggest that despite data coverage,iterative initial conditions also play an essential role in ionospheric reconstruction.In particular,in the longitudinal sectors where the iterative initial height of the F2-layer peak electron density(hmF2)differs substantially from the background densities,none of the three CIT methods can reproduce the exact background profile.When hmF2 is close but the ionospheric F2-layer peak density(NmF2)is different between the targeted background and initial conditions,the MART performs better than the ART and SIRT,as evidenced by the correlation coefficients of MART being above 0.97 and those of ART and SIRT being below 0.85.In summary,this investigation reveals the potential uncertainties in traditional CIT reconstruction,particularly when realistic hmF2 or NmF2 values differ substantially from the initial CIT conditions.
基金the National Natural Science Founda-tion of China(NSFC,Grant No.42174181)and the Key Research Program of Frontier Sciences,CAS(Grant No.QYZDJ-SSW-DQC010).
文摘The formation of an embedded electron current sheet within the magnetotail plasma sheet has been poorly understood.In this article,we present an electron current layer detected at the edge of the magnetotail plasma sheet.The ions were demagnetized inside the electron current layer,but the electrons were still frozen in with the magnetic field line.Thus,this decoupling of ions and electrons gave rise to a strong Hall electric field,which could be the reason for the formation of the embedded thin current layer.The magnetized electrons,the absence of the nongyrotropic electron distribution,and negligible energy dissipation in the layer indicate that magnetic reconnection had not been triggered within the embedded thin current layer.The highly asymmetric plasma on the two sides of the current layer and low magnetic shear across it could suppress magnetic reconnection.The observations indicate that the embedded electric current layer,probably generated by the Hall electric field,even down to electron scale,is not a sufficient condition for magnetic reconnection.
基金supported by the National Natural Science Foundation of China(42474003,42074007)the Fundamental Research Funds for the Central Universities(2042023kfyq01)。
文摘As a major fault in the northeastern Qinghai-Xizang Plateau,the Haiyuan fault zone is important for understanding the regional deformation.Aiming at the differences in the slip rate and locking degree obtained from different studies,this study constructs a refined block model(including Qilian,Alxa,Ordos,Xining,Haiyuan,and Lanzhou blocks)and uses the grid search and simulated annealing methods to invert GPS data for slip rate and locking degree of the Haiyuan fault zone.The results are as follows:(1)The sinistral slip rates in the western,middle,and eastern segments are 4.93-5.22 mm/a,1.52-4.94 mm/a,and 0.43-1.18 mm/a,decreasing eastward on the whole,while the compression rates are 0.45-1.26 mm/a,0.58-2.62 mm/a,and3.52-4.48 mm/a,increasing eastward on the whole.(2)The locking depth of the western segment increases from about 5 km to about 20 km eastward;the middle segment decreases and then increases eastward;the eastern segment concentrates at about 20 km(PHI is about 0.86).(3)The slip deficit is relatively higher in the Lenglongling,Jinqianghe,Maomaoshan,and Liupanshan faults(averaging about 3.42 mm/a,4.16 mm/a,4.23 mm/a,and 3.43 mm/a within 20 km).(4)The Qilian,Alxa,Xining,Lanzhou,and Haiyuan blocks rotate clockwise,while the Ordos block rotates counterclockwise.Additionally,by comparing different block models,the Haiyuan block should be considered independently.The Haiyuan fault zone adjusts surrounding block movements and uplifts Liupanshan mountain tectonically.The results can provide important references for understanding the regional earthquake risk and deformation mechanism.
基金supported by National Natural Science Foundation of China (No.11975229)。
文摘Electrostatic ion-cyclotron wave(EICW)and ion-Bernstein wave(IBW)are two typical normal modes of plasma propagating vertically to a magnetic field.These modes satisfy the extreme conditionsω《k_(||)v_(te)andω》k_(||)v_(te),whereω,k_(||),and v_(te)are the wave angular frequency,wave number parallel to the magnetic field,and the electron thermal speed.Between these extremities,there exist IBW-like waves that have not been studied.In this paper,the co-excitation of the EICW and IBW-like wave in an experiment is presented.It is shown that,in an argon(Ar)discharge,the EICW appears in the fundamental and second harmonic frequency bands of the ioncyclotron frequency.In contrast,the IBW-like wave appears in the higher harmonic frequency bands and in the case of a relatively strong magnetic field.In a helium(He)discharge,the IBW-like wave is excited in the fundamental frequency band,but the wave behaviour becomes complicated in the harmonic frequency bands.The measured dispersion relations of the excited EICW are in agreement with the results of the fluid and kinetic models in some range of frequencies,but those of the excited IBW-like wave are only in qualitative agreement with the kinetic model of the IBW.This study is important for identifying the IBW-like mode that may have potential application in plasma heating via Landau damping.
基金supported by the National Natural Science Foundation of China(Grant No.42474227,42241106,42388101)financial support through the German Ministry for Economy and Technology and the German Center for Aviation and Space(DLR)under contract 50 OC 0302
文摘The solar wind's interaction with the Moon has traditionally been understood through the Moon's absorption of solar wind particles and the formation of a plasma cavity on its nightside,known as the lunar wake.This study reveals unexpected,large-scale perturbations in the solar wind upstream of the Moon,using 11 years of data from the OMNI and ARTEMIS(Acceleration,Reconnection,Turbulence and Electrodynamics of Moon's Interaction with the Sun)missions(2012-2023).We find systematic moonward deviations of~tens of km/s in a direction perpendicular to the solar wind(moonward),at altitudes of up to 1000 km,particularly when the interplanetary magnetic field(IMF)lines are oblique to the solar wind(30°<θ<60°)and connected to the lunar dayside.The longer the duration of the interaction,the greater the moonward deviation.These perturbations can be explained by neither solar wind pickup of the reflected ions,nor lunar wake dynamics.Instead,they appear to correlate with magnetic connectivity between the ARTEMIS probes and the lunar surface,suggesting a more complex solar wind interaction than previously thought.
基金supported by the National Natural Science Foundation of China (Nos. 40974015,40637034)the Fund of Key Laboratory of Geodynamic Geodesy, Chinese Academy of Sciences (No. 09-18)the Fund of Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,China (No. 07-12)
文摘Gravity Recovery and Climate Experiment(GRACE) observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 26,2004.This article focuses on investigating some gravita-tional effects caused by this huge earthquake.We computed the geoid height changes,the equivalent water height(EWH) changes,and the gravity changes using the GRACE Level-2 monthly spherical harmonic(SH) solutions released by University of Texas Center for Space Research(UTCSR).The GRACE results agree well with the prediction by a dislocation model and are consistent with the results obtained by some previous scholars.In particular,we calculated the three components of the gravity gradient variations and found that they can recover the seismic-related signature more sensitively due to a certain degree of amplification of the signals.A positive-negative-positive mode predominates in the spatial distribution of the horizontal components of the gravity gradient variations,which is possibly attributed to the anomalies in the crustal density distribution caused by the uplift-subduction effect of the dip-slip earthquake.Moreover,the latitude components of the gravity gradient changes show strong suppression of the north-south stripes,which is due to the along-orbit measurements of the two GRACE satellites.We conclude that the posi-tive-negative-positive mode in latitude gravity gradient changes would be a more sensitive fea-ture to detect the deformations of some major dip-slip earthquakes by GRACE data.
基金supported by the Chinese National Space Administration,and the science team of MOMAG is supported by the Strategic Priority Program(Grant No.XDB41000000)the key Research Program of Frontier Sciences(No.QYZDB-SSW-DQC015)of the Chinese Academy of Sciences.
文摘As one of the seven scientific payloads on board the Tianwen-1 orbiter,the Mars Orbiter Magnetometer(MOMAG)will measure the magnetic fields of and surrounding Mars to study its space environment and the interaction with the solar wind.The instrument consists of two identical triaxial fluxgate magnetometer sensors,mounted on a 3.19 meter-long boom with a seperation of about 90 cm.The dual-magnetometers configuration will help eliminate the magnetic field interference generated by the spacecraft platform and payloads.The sensors are controlled by an electric box mounted inside the orbiter.Each magnetometer measures the ambient vector magnetic field over a wide dynamic range(to 10,000 nT per axis)with a resolution of 1.19 pT.Both magnetometers sample the ambient magnetic field at an intrinsic frequency of 128 Hz,but will operate in a model with alternating frequency between 1 and 32 Hz to meet telemetry allocations.
基金supported by the National Natural Science Foundation of China (Nos. 41631072, 41721003, 41574007, and 41429401)the Discipline Innovative Engineering Plan of Modern Geodesy and Geodynamics (No. B17033)+1 种基金the DAAD Thematic Network Project (No. 57173947)the International Space Science Institute (ISSI) 2017–2019
文摘Based on gravity frequency shift effect predicted by general relativity theory, this study discusses an approach for determining the gravity potential(geopotential) difference between arbitrary two points P and Q by remote comparison of two precise optical clocks via optical fiber frequency transfer. After synchronization, by measuring the signal's frequency shift based upon the comparison of bidirectional frequency signals from P and Q oscillators connected with two optical atomic clocks via remote optical fiber frequency transfer technique, the geopotential difference between the two points could be determined, and its accuracy depends on the stabilities of the optical clocks and the frequency transfer comparison technique. Due to the fact that the present stability of optical clocks achieves 1.6×10-18 and the present frequency transfer comparison via optical fiber provides stabilities as high as 10-19 level, this approach is prospective to determine geopotential difference with an equivalent accuracy of 1.5 cm. In addition, since points P and Q are quite arbitrary, this approach may provide an alternative way to determine the geopotential over a continent, and prospective potential to unify a regional height datum system.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11374108 and 10904041)the Foundation for the Author of Guangdong Provincial Excellent Doctoral Dissertation+6 种基金China(Grant No.SYBZZXM201227)the Foundation of Cultivating Outstanding Young Scholars("ThousandHundredTen"Program)of Guangdong Province in Chinathe Fund from the Key Laboratory of Geospace EnvironmentUniversity of Science and Technology of ChinaChinese Academy of Sciences
文摘Under the paraxial approximation, the analytical propagation expression of an Airy–Gaussian beam(Ai GB) in uniaxial crystals orthogonal to the optical axis is investigated. The propagation dynamics of the Ai GB is given for different ratios of the extraordinary index to the ordinary refractive index. It has been found that the continuity and the self-bending effect of Ai GB become weaker when the ratio increases. From the figure of the maximum intensity of Ai GB, one can see that the maximum intensity is not monotone decreasing due to the anisotropic effect of the crystals. The intensity distribution of Ai GB in different distribution factors is shown. The Ai GB converges toward a Gaussian beam as the distribution factor increases.
基金supported by the National Natural Science Foundation of China(Grant Nos.11374108,11374107,10904041,and 11547212)the Foundation of Cultivating Outstanding Young Scholars of Guangdong Province,China+2 种基金the CAS Key Laboratory of Geospace Environment,University of Science and Technology of Chinathe National Training Program of Innovation and Entrepreneurship for Undergraduates(Grant No.2015093)the Science and Technology Projects of Guangdong Province,China(Grant No.2013B031800011)
文摘The propagation dynamics of the Airy Gaussian vortex beams in uniaxial crystals orthogonal to the optical axis has been investigated analytically and numerically. The propagation expression of the beams has been obtained. The propagation features of the Airy Gaussian vortex beams are shown with changes of the distribution factor and the ratio of the extraordinary refractive index to the ordinary refractive index. The correlations between the ratio and the maximum intensity value during the propagation, and its appearing distance have been investigated.
基金supported by the 973 Project of China (No.2013CB733303)the National Natural Science Foundation of China (No.41474093)the Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy of Ministry of Education,China (No.12-02-08)
文摘Lushan Earthquake (-Mw 6.6) occurred in Sichuan Province of China on 20 April 2013, was the largest earthquake in Longmenshan fault belt since 2008 Wenchuan Earthquake. To better understand its rupture pattern, we focused on the influences of fault parameters on fault slips and performed fault slip inversion using Akaike's Bayesian Information Criterion (ABIC) method. Based on GPS coseismic data, our inverted results showed that the fault slip was mainly confined at depths. The maximum slip amplitude is about 0.7 m, and the scalar seismic moment is about 9.47x10TM N.m. Slip pattern reveals that the earthquake occurred on the thrust fault with large dip-slip and small strike-slip, such a simple fault slip represents no second sub-event occurred. The Coulomb stress changes (ACFF) matched the most aftershocks with negative anomalies. The in- verted results demonstrated that the source parameters have significant impacts on fault slip distri- bution, especially on the slip direction and maximum displacement.
基金supported by National 973 Project China (2013CB733305)National Natural Science Foundation of China (NSFCs) (41174011,41429401,41210006,41128003,41021061)
文摘Scientists pay great attention to different-time-scale signals in the lengllh of day (LOD) variations △LOD, which provide signatures of the Earth's interior structure, couplings among different layers, and potential excitations of ocean and atmosphere. In this study, based on the ensemble empirical mode decomposition (EEMD), we analyzed the latest time series of △LOD data spanning from January 1962 to March 2015. We observed the signals with periods and amplitudes of about 0.5 month and 0.19 ms, 1.0 month and 0.19 ms, 0.5 yr and 0.22 ms, 1.0 yr and 0.18 ms, 2.28 yr and 0.03 ms, 5.48 yr and 0.05 ms, respectively, in coincidence with the results of predecessors. In addition, some signals that were previously not definitely observed by predecessors were detected in this study, with periods and amplitudes of 9.13 d and 0.12 ms, 13.69 yr and 0.10 ms, respectively. The mechanisms of the LOD fluctuations of these two signals are still open.
基金financially supported by the National 973-Project(No.2013CB733303)the National Natural Science Foundation of China(No.41474093)the Key Natural Science Foundation of Hubei Province(No.2014CFA110)
文摘Previous studies have reached consensus that low velocity zones are widespread in the crust beneath Yunnan region. However, the relationships between the low velocity zones and large faults, earthquake distribution are less investigated by available studies. By analyzing the seismic ambient noise recorded by Yunnan Seismic Networks and Tengchong volcano array, we construct a 3D crustal shear wave velocity model for the Yunnan region, which provides more details of the distribution of intra-crustal low velocity zones all over Yunnan. The distribution of low velocity zones shows different features at different depths. At shallow depths, the results are well correlated with near surface geological features. With increasing depth, the low velocity zones are gradually concentrated on the northern part of our study area, most likely reflecting variations in crustal thickness beneath the Yunnan region. The low velocity zones are truncated at depth by several large faults in Yunnan. It is interesting that most strong earthquakes(Ms≥5.0) occurred in Yunnan are distributed in low velocity zones or the transition zones between low and high velocity anomalies within the upper-to-middle crust. The crustal structure is composed of a brittle, seismically active upper-to-middle crust and a warm, aseismic lower crust.
基金funded by the National Key Research and Development Program of China(Grant No.2018YFC1503503)the National Natural Science Foundation of China(Grant Nos.41974015,42061134007,41474019)。
文摘Surface mass anomalies estimated by mass concentration(mascon)approach using Gravity Recovery and Climate Experiment(GRACE)observations with regularization constraints generally present higher spatial resolution than the spheric harmonic(SH)solutions.To analyze the influence of different types of constraints on the estimation of mascon solutions,we carried out a closed-loop simulation experiment to estimate surface mass anomalies over South America based on simulated GRACE intersatellite geopotential differences.Tikhonov regularization with spatial constraint(SC),uniform weighting constraint(UWC),and a prior information constraint(APC)were employed to stabilize the mascon solutions,and the corresponding optimal regularization parameters were determined based on the minimum residual root-mean-square(RMS)criterion.The results show that mascon solutions estimated under different types of constraints are consistent and equivalent when the optimal regularization parameters are selected.The spatial distributions and main characteristics of regional surface mass anomalies estimated by the three types of constraints agree well,and the values of residual RMS with different constraints are very close.But due to the smoothing effect of regularization,the signal strength of mascon solutions is a bit weaker than that of original true signal,especially in the regions with strong signals.In addition,due to the ill-conditioned problem is more serious for higher grid resolution,the relative contribution of the three types of constraints to the final mascon solutions would be stronger.The results show that the averages of relative contribution percentages of these constraints for 2°×2° mascon grids are 80%-90%,while the corresponding values for 4°×4° mascon grids are 30%-60%.However,based on the minimum residual RMS criterion,the accuracy of estimation results is not affected by the type of constraints and their relative contribution to the final mascon solutions.
基金supported by the National Key Scientific Instrument and Equipment Development Projects of the National Natural Science Foundation of China(No.41327802)China Mars Project
文摘Readout electronics is developed for a prototype time-of-flight(TOF) ion composition spectrometer for in situ measurement of the mass/charge distributions of major ion species from 200 to 100 ke V/e in space plasma.By utilizing a constant fraction discriminator(CFD) and time-to-digital converter(TDC), challenging dynamic range measurements were performed with high time resolution and event rates. CFD was employed to discriminate the TOF signals from the micro-channel plate and channel electron multipliers. TDC based on the combination of counter and OR-gate delay chain was designed in a highreliability flash field programmable gate array. Owing to the non-uniformity of the delay chain, a correction algorithm based on integral nonlinearity compensation was implemented to reduce the time uncertainty. The test results showed that the electronics achieved a low timingerror of < 200 ps in the input range from 35 to 500 m V for the CFD, and a time resolution of ~550 ps with time uncertainty < 180 ps after correction and a time range of6.4 ls for the TDC. The TOF spectrum from an electron beam experiment of the impacting N_2 gas further indicated the good performance of this readout electronic.
基金financially supported by the Shenzhen Science and Technology Research Program(grant no.JCYJ2021032412141-2034)the National Natural Science Foundation of China(grant nos.41804157,41774171,41774167,41974205,41804157,and 41904156)+3 种基金financial support from the pre-research Project on Civil Aerospace Technologies(grant no.D020103)funded by the China National Space Administrationthe 111 Project(grant no.B18017)the Chinese Academy of Sciences Center for Excellence in Comparative Planetologythe Macao Foundation.
文摘The space-borne fluxgate magnetometer(FGM)requires regular in-flight calibration to obtain its zero offset.Recently,Wang GQ and Pan ZH(2021a)developed a new method for the zero offset calibration based on the properties of Alfvén waves.They found that an optimal offset line(OOL)exists in the offset cube for a pure Alfvén wave and that the zero offset can be determined by the intersection of at least two nonparallel OOLs.Because no pure Alfvén waves exist in the interplanetary magnetic field,calculation of the zero offset relies on the selection of highly Alfvénic fluctuation events.Here,we propose an automatic procedure to find highly Alfvénic fluctuations in the solar wind and calculate the zero offset.This procedure includes three parts:(1)selecting potential Alfvénic fluctuation events,(2)obtaining the OOL,and(3)determining the zero offset.We tested our automatic procedure by applying it to the magnetic field data measured by the FGM onboard the Venus Express.The tests revealed that our automatic procedure was able to achieve results as good as those determined by the Davis-Smith method.One advantage of our procedure is that the selection criteria and the process for selecting the highly Alfvénic fluctuation events are simpler.Our automatic procedure could also be applied to find fluctuation events for the Davis-Smith method.
基金funded by the National 973 Project China (2013CB733302)National Natural Science Foundation of China (41504014, 41474019)
文摘As critical component of hydrologic cycle, basin discharge is a key issue for understanding the hydrological and climatologic related to water and energy cycles. Combining GRACE gravity field models with ET from GLDAS models and precipitation from GPCP, discharge of the Yellow River basin are estimated from the water balance equation. While comparing the results with discharge from GLDAS model and in situ measurements, the results reveal that discharge from Mosaic and CLM GLDAS model can partially represent the river discharge and the discharge estimation from water balance equation could reflect the discharge from precipitation over the Yellow River basin.
基金supported by the National Natural Science Foundation of China (Nos.41504015,41429401)the National 973 Project of China (No.2013CB733302)+2 种基金 China Postdoctoral Science Foundation (No.2015M572146)the National High Technology Research and Development Program of China (No.2011AA060503)the Surveying and Mapping Basic Research Program of National Administration of Surveying,Mapping and Geoinformation (No.15-01-08)
文摘Urban faults in Shenzhen are potential threats to city security and sustainable development. In consideration of the importance of the Shenzhen fault zone, the author provide a detailed interpretation on gravity data model. Bouguer gravity covering the whole Shenzhen City was calculated with a 1-km resolution. Wavelet multi-scale analysis(MSA) was applied to the Bouguer gravity data to obtain the multilayer residual anomalies corresponding to different depths. In addition, 2D gravity models were constructed along three profiles. The Bouguer gravity anomaly shows an NE-striking high-low-high pattern from northwest to southeast, strongly related to the main faults. According to the results of MSA, the correlation between gravity anomaly and faults is particularly significant from 4 to 12 km depth. The residual gravity with small amplitude in each layer indicates weak tectonic activity in the crust. In the upper layers, positive anomalies along most of faults reveal the upwelling of high-density materials during the past tectonic movements. The multilayer residual anomalies also yield important information about the faults, such as the vertical extension and the dip direction. The maximum depth of the faults is about 20 km. In general, NE-striking faults extend deeper than NW-striking faults and have a larger dip angle.