Automatic Generation of Artificial Space Weather Forecast Product Based on Sequence-to-sequence Model

Both analyzing a large amount of space weather observed data and alleviating personal experience bias are significant challenges in generating artificial space weather forecast products.With the use of natural language generation methods based on the sequence-to-sequence model,space weather forecast texts can be automatically generated.To conduct our generation tasks at a fine-grained level,a taxonomy of space weather phenomena based on descriptions is presented.Then,our MDH(Multi-Domain Hybrid)model is proposed for generating space weather summaries in two stages.This model is composed of three sequence-to-sequence-based deep neural network sub-models(one Bidirectional Auto-Regressive Transformers pre-trained model and two Transformer models).Then,to evaluate how well MDH performs,quality evaluation metrics based on two prevalent automatic metrics and our innovative human metric are presented.The comprehensive scores of the three summaries generating tasks on testing datasets are 70.87,93.50,and 92.69,respectively.The results suggest that MDH can generate space weather summaries with high accuracy and coherence,as well as suitable length,which can assist forecasters in generating high-quality space weather forecast products,despite the data being starved.

Image Feature Extraction and Matching of Augmented Solar Images in Space Weather

Augmented solar images were used to research the adaptability of four representative image extraction and matching algorithms in space weather domain.These include the scale-invariant feature transform algorithm,speeded-up robust features algorithm,binary robust invariant scalable keypoints algorithm,and oriented fast and rotated brief algorithm.The performance of these algorithms was estimated in terms of matching accuracy,feature point richness,and running time.The experiment result showed that no algorithm achieved high accuracy while keeping low running time,and all algorithms are not suitable for image feature extraction and matching of augmented solar images.To solve this problem,an improved method was proposed by using two-frame matching to utilize the accuracy advantage of the scale-invariant feature transform algorithm and the speed advantage of the oriented fast and rotated brief algorithm.Furthermore,our method and the four representative algorithms were applied to augmented solar images.Our application experiments proved that our method achieved a similar high recognition rate to the scale-invariant feature transform algorithm which is significantly higher than other algorithms.Our method also obtained a similar low running time to the oriented fast and rotated brief algorithm,which is significantly lower than other algorithms.

Two methods for separating the magnetospheric solar wind charge exchange soft X-ray emission from the diffuse X-ray background

Solar wind charge exchange(SWCX)is the process of solar wind high-valence ions exchanging charges with neutral components and generating soft X-rays.Recently,detecting the SWCX emission from the magnetosphere is proposed as a new technique to study the magnetosphere using panoramic soft X-ray imaging.To better prepare for the data analysis of upcoming magnetospheric soft X-ray imaging missions,this paper compares the magnetospheric SWCX emission obtained by two methods in an XMM-Newton observation,during which the solar wind changed dramatically.The two methods differ in the data used to fit the diffuse X-ray background(DXB)parameters in spectral analysis.The method adding data from the ROSAT All-Sky Survey(RASS)is called the RASS method.The method using the quiet observation data is called the Quiet method,where quiet observations usually refer to observations made by the same satellite with the same target but under weaker solar wind conditions.Results show that the spectral compositions of magnetospheric SWCX emission obtained by the two methods are very similar,and the changes in intensity over time are highly consistent,although the intensity obtained by the RASS method is about 2.68±0.56 keV cm^(-2)s^(-1)sr^(-1)higher than that obtained by the Quiet method.Since the DXB intensity obtained by the RASS method is about 2.84±0.74 keV cm^(-2)s^(-1)sr^(-1)lower than that obtained by the Quiet method,and the linear correlation coefficient between the difference of SWCX and DXB obtained by the two methods in diffe rent energy band is close to-1,the diffe rences in magnetospheric SWCX can be fully attributed to the diffe rences in the fitted DXB.The difference between the two methods is most significant when the energy is less than 0.7 keV,which is also the main energy band of SWCX emission.In addition,the difference between the two methods is not related to the SWCX intensity and,to some extent,to solar wind conditions,because SWCX intensity typically va ries with the solar wind.In summary,both methods are robust and reliable,and should be considered based on the best available options.

Assessment of International GNSS Service Global Ionosphere Map products over China region based on measurements from the Crustal Movement Observation Network of China

The global ionosphere maps(GIM)provided by the International GNSS Service(IGS)are extensively utilized for ionospheric morphology monitoring,scientific research,and practical application.Assessing the credibility of GIM products in data-sparse regions is of paramount importance.In this study,measurements from the Crustal Movement Observation Network of China(CMONOC)are leveraged to evaluate the suitability of IGS-GIM products over China region in 2013-2014.The indices of mean error(ME),root mean square error(RMSE),and normalized RMSE(NRMSE)are then utilized to quantify the accuracy of IGS-GIM products.Results revealed distinct local time and latitudinal dependencies in IGS-GIM errors,with substantially high errors at nighttime(NRMSE:39%)and above 40°latitude(NRMSE:49%).Seasonal differences also emerged,with larger equinoctial deviations(NRMSE:33.5%)compared with summer(20%).A preliminary analysis implied that the irregular assimilation of sparse IGS observations,compounded by China’s distinct geomagnetic topology,may manifest as error variations.These results suggest that modeling based solely on IGS-GIM observations engenders inadequate representations across China and that a thorough examination would proffer the necessary foundation for advancing regional total electron content(TEC)constructions.

An Improved HVQ Algorithm for Compression and Rendering of Space Environment Volume Data with Multi-correlated Variables

Volume visualization can not only illustrate overall distribution but also inner structure and it is an important approach for space environment research.Space environment simulation can produce several correlated variables at the same time.However,existing compressed volume rendering methods only consider reducing the redundant information in a single volume of a specific variable,not dealing with the redundant information among these variables.For space environment volume data with multi-correlated variables,based on the HVQ-1d method we propose a further improved HVQ method by compositing variable-specific levels to reduce the redundant information among these variables.The volume data associated with each variable is divided into disjoint blocks of size 43 initially.The blocks are represented as two levels,a mean level and a detail level.The variable-specific mean levels and detail levels are combined respectively to form a larger global mean level and a larger global detail level.To both global levels,a splitting based on a principal component analysis is applied to compute initial codebooks.Then,LBG algorithm is conducted for codebook refinement and quantization.We further take advantage of progressive rendering based on GPU for real-time interactive visualization.Our method has been tested along with HVQ and HVQ-1d on high-energy proton flux volume data,including>5,>10,>30 and>50 MeV integrated proton flux.The results of our experiments prove that the method proposed in this paper pays the least cost of quality at compression,achieves a higher decompression and rendering speed compared with HVQ and provides satisficed fidelity while ensuring interactive rendering speed.

Analysis of the joint detection capability of the SMILE satellite and EISCAT-3D radar

The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)satellite is a small magnetosphere–ionosphere link explorer developed cooperatively between China and Europe.It pioneers the use of X-ray imaging technology to perform large-scale imaging of the Earth’s magnetosheath and polar cusp regions.It uses a high-precision ultraviolet imager to image the overall configuration of the aurora and monitor changes in the source of solar wind in real time,using in situ detection instruments to improve human understanding of the relationship between solar activity and changes in the Earth’s magnetic field.The SMILE satellite is scheduled to launch in 2025.The European Incoherent Scatter Sciences Association(EISCAT)-3D radar is a new generation of European incoherent scatter radar constructed by EISCAT and is the most advanced ground-based ionospheric experimental device in the high-latitude polar region.It has multibeam and multidirectional quasi-real-time three-dimensional(3D)imaging capabilities,continuous monitoring and operation capabilities,and multiple-baseline interferometry capabilities.Joint detection by the SMILE satellite and the EISCAT-3D radar is of great significance for revealing the coupling process of the solar wind–magnetosphere–ionosphere.Therefore,we performed an analysis of the joint detection capability of the SMILE satellite and EISCAT-3D,analyzed the period during which the two can perform joint detection,and defined the key scientific problems that can be solved by joint detection.In addition,we developed Web-based software to search for and visualize the joint detection period of the SMILE satellite and EISCAT-3D radar,which lays the foundation for subsequent joint detection experiments and scientific research.

The Soft X-ray Imager(SXI)on the SMILE Mission

The Soft X-ray Imager(SXI)is part of the scientific payload of the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.SMILE is a joint science mission between the European Space Agency(ESA)and the Chinese Academy of Sciences(CAS)and is due for launch in 2025.SXI is a compact X-ray telescope with a wide field-of-view(FOV)capable of encompassing large portions of Earth’s magnetosphere from the vantage point of the SMILE orbit.SXI is sensitive to the soft X-rays produced by the Solar Wind Charge eXchange(SWCX)process produced when heavy ions of solar wind origin interact with neutral particles in Earth’s exosphere.SWCX provides a mechanism for boundary detection within the magnetosphere,such as the position of Earth’s magnetopause,because the solar wind heavy ions have a very low density in regions of closed magnetic field lines.The sensitivity of the SXI is such that it can potentially track movements of the magnetopause on timescales of a few minutes and the orbit of SMILE will enable such movements to be tracked for segments lasting many hours.SXI is led by the University of Leicester in the United Kingdom(UK)with collaborating organisations on hardware,software and science support within the UK,Europe,China and the United States.

Deformations at Earth’s dayside magnetopause during quasi-radial IMF conditions:Global kinetic simulations and Soft X-ray Imaging

The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)is a joint mission of the European Space Agency(ESA)and the Chinese Academy of Sciences(CAS).Primary goals are investigating the dynamic response of the Earth's magnetosphere to the solar wind(SW)impact via simultaneous in situ magnetosheath plasma and magnetic field measurements,X-Ray images of the magnetosheath and magnetic cusps,and UV images of global auroral distributions.Magnetopause deformations associated with magnetosheath high speed jets(HSJs)under a quasi-parallel interplanetary magnetic field condition are studied using a threedimensional(3-D)global hybrid simulation.Soft X-ray intensity calculated based on both physical quantities of solar wind proton and oxygen ions is compared.We obtain key findings concerning deformations at the magnetopause:(1)Magnetopause deformations are highly coherent with the magnetosheath HSJs generated at the quasi-parallel region of the bow shock,(2)X-ray intensities estimated using solar wind h+and self-consistentO7+ions are consistent with each other,(3)Visual spacecraft are employed to check the discrimination ability for capturing magnetopause deformations on Lunar and polar orbits,respectively.The SMILE spacecraft on the polar orbit could be expected to provide opportunities for capturing the global geometry of the magnetopause in the equatorial plane.A striking point is that SMILE has the potential to capture small-scale magnetopause deformations and magnetosheath transients,such as HSJs,at medium altitudes on its orbit.Simulation results also demonstrate that a lunar based imager(e.g.,Lunar Environment heliospheric X-ray Imager,LEXI)is expected to observe a localized brightening of the magnetosheath during HSJ events in the meridian plane.These preliminary results might contribute to the pre-studies for the SMILE and LEXI missions by providing qualitative and quantitative soft X-ray estimates of dayside kinetic processes.

Influence of topography on the fine structures of stratospheric gravity waves:An analysis using COSMIC-2 temperature data

We derive the potential energy of gravity waves(GWs)in the upper troposphere and stratosphere at 45°S-45°N from December 2019 to November 2022 by using temperature profiles retrieved from the Constellation Observing System for Meteorology,Ionosphere,and Climate-2(COSMIC-2)satellite.Owing to the dense sampling of COSMIC-2,in addition to the strong peaks of gravity wave potential energy(GWPE)above the Andes and Tibetan Plateau,we found weak peaks above the Rocky,Atlas,Caucasus,and Tianshan Mountains.The land-sea contrast is responsible for the longitudinal variations of the GWPE in the lower and upper stratosphere.At 40°N/S,the peaks were mainly above the topographic regions during the winter.At 20°N/S,the peaks were a slight distance away from the topographic regions and might be the combined effect of nontopographic GWs and mountain waves.Near the Equator,the peaks were mainly above the regions with the lowest sea level altitude and may have resulted from convection.Our results indicate that even above the local regions with lower sea level altitudes compared with the Andes and Tibetan Plateau,the GWPE also exhibits fine structures in geographic distributions.We found that dissipation layers above the tropopause jet provide the body force to generate secondary waves in the upper stratosphere,especially during the winter months of each hemisphere and at latitudes of greater than 20°N/S.

Ensemble Numerical Simulations of Realistic SEP Events and the Inspiration for Space Weather Awareness

The solar energetic particle(SEP)event is a kind of hazardous space weather phenomena,so its quantitative forecast is of great importance from the aspect of space environmental situation awareness.We present here a set of SEP forecast tools,which consists of three components:(1)a simple polytropic solar wind model to estimate the background solar wind conditions at the inner boundary of 0.1 AU(about 20 R⊙);(2)an ice-cream-cone model to estimate the erupted coronal mass ejection(CME)parameters;and(3)the improved Particle Acceleration and Transport in the Heliosphere(i PATH)model to calculate particle fluxes and energy spectra.By utilizing the above models,we have simulated six realistic SEP events from 2010 August 14 to 2014 September 10,and compared the simulated results to the Geostationary Operational Environmental Satellites(GOES)spacecraft observations.The results show that the simulated fluxes of>10 Me V particles agree with the observations while the simulated fluxes of>100 Me V particles are higher than the observed data.One of the possible reasons is that we have adopted a simple method in the model to calculate the injection rate of energetic particles.Furthermore,we have conducted the ensemble numerical simulations over these events and investigated the effects of different background solar wind conditions at the inner boundary on SEP events.The results imply that the initial CME density plays an important role in determining the power spectrum,while the effect of varying background solar wind temperature is not significant.Naturally,we have examined the influence of CME initial density on the numerical prediction results for virtual SEP cases with different CME ejection speeds.The result shows that the effect of initial CME density variation is inversely associated with CME speed.

Study of the Modified Gaussian Model on olivine diagnostic spectral features and its applications in space weathering experiments

The absorption features of olivine in visible and near-infrared(VNIR)reflectance spectra are the key spectral parameters in its mineralogical studies.Generally,these spectral parameters can be obtained by exploiting the Modified Gaussian Model(MGM)with a proper continuum removal.However,different continua may change the deconvolution results of these parameters.This paper investigates the diagnostic spectral features of olivine with diverse chemical compositions.Four different continuum removal methods with MGM for getting the deconvolution results are presented and the regression equations for predicting the Mg-number(Fo#)are introduced.The results show that different continua superimposed on the mineral absorption features will make the absorption center shift,as well as the obvious alterations in shape,width and strength of the absorption band.Additionally,it is also found that the logarithm of a second-order polynomial continuum can match the overall shape of the spectrum in logarithmic space,and the improved regression equations applied to estimate the chemical composition of olivine-dominated spectra also have a better performance.As an application example,the improved approach is applied to pulse laser irradiated olivine grains to simulate and study the space weathering effects on olivine diagnostic spectral features.The experiments confirm that space weathering can make the absorption band center shift toward longer wavelength.Therefore,the Fo#estimated from remote sensing spectra may be less than its actual chemical composition.These results may provide valuable information for revealing the difference between the spectra of olivine grains and olivine-dominated asteroids.

Statistical properties of kinetic-scale magnetic holes in terrestrial space

Kinetic-scale magnetic holes(KSMHs)are structures characterized by a significant magnetic depression with a length scale on the order of the proton gyroradius.These structures have been investigated in recent studies in near-Earth space,and found to be closely related to energy conversion and particle acceleration,wave-particle interactions,magnetic reconnection,and turbulence at the kineticscale.However,there are still several major issues of the KSMHs that need further study—including(a)the source of these structures(locally generated in near-Earth space,or carried by the solar wind),(b)the environmental conditions leading to their generation,and(c)their spatio-temporal characteristics.In this study,KSMHs in near-Earth space are investigated statistically using data from the Magnetospheric Multiscale mission.Approximately 200,000 events were observed from September 2015 to March 2020.Occurrence rates of such structures in the solar wind,magnetosheath,and magnetotail were obtained.We find that KSMHs occur in the magnetosheath at rates far above their occurrence in the solar wind.This indicates that most of the structures are generated locally in the magnetosheath,rather than advected with the solar wind.Moreover,KSMHs occur in the downstream region of the quasi-parallel shock at rates significantly higher than in the downstream region of the quasi-perpendicular shock,indicating a relationship with the turbulent plasma environment.Close to the magnetopause,we find that the depths of KSMHs decrease as their temporal-scale increases.We also find that the spatial-scales of the KSMHs near the subsolar magnetosheath are smaller than those in the flanks.Furthermore,their global distribution shows a significant dawn-dusk asymmetry(duskside dominating)in the magnetotail.

Determination of the Kelvin-Helmholtz Wave Parameters on the Magnetopause in Single-spacecraft Observations

Kelvin-Helmholtz(K-H) waves are formed from the triggering of the K-H instability on the magnetopause,which is a candidate mechanism for solar wind entry into the magnetosphere,especially under northward interplanetary magnetic field conditions.In this study,a K-H wave event was identified from the observation of probe B of the Time History of Events and Macroscale Interactions during Substorms(THEMIS) mission on 15 May 2008.A new method to determine the wave parameters of the K-H waves in single-spacecraft observations is proposed.The dominant wave period is determined by three kinds of spectrograms for three key parameters,namely the ion density,the ion temperature,and the z component of magnetic field.The phase velocity is estimated by calculating the center-of-mass velocity of the detected K-H vortex region.This approximation is validated by comparison with other alternative methods.The method to determine the wave parameters is a first step to further study K-H wave properties and their relationship with interplanetary conditions.

Theoretical study of the ionospheric dynamo region inside the South Atlantic Anomaly

The South Atlantic Anomaly(SAA)is a region where the geomagnetic field is significantly lower than that of the surrounding area.On the basis of the models of CHAOS-7.8,Mass Spectrometer Incoherent Scatter Model(NRLMSISE-00),and International Reference Ionosphere 2016(IRI-2016),we theoretically investigated the lower and upper boundaries of the ionospheric dynamo region inside the SAA.In the ionospheric dynamo region,electrons are coupled with magnetic field lines,whereas ions are decoupled from magnetic field lines.Our results showed that the ionospheric dynamo region inside the SAA is higher and larger than that outside the SAA.We also studied the boundary variations of the dynamo region inside the SAA depending on the seasons and solar activities.We found that the dynamo region inside the SAA is the highest and largest in the summer of the southern hemisphere at solar maximum.The larger and higher altitude range of the ionospheric dynamo region in the SAA can contribute to the stronger ionospheric currents in this region.

Clarification on Polarity of Bipolar Electric Field Solitary Structures in Space Plasmas with Satellite Observation

The bipolar electric field solitary(EFS)structures observed frequently in space plasmas by satellites have two different polarities,first positive electric field peak then negative(i.e.,positive/negative)and first negative then positive peak(i.e.,negative/positive).We provide the physical explanation on the polarity of observed bipolar EFS structures with an electrostatic ion fluid model.The results show that if initial electric field𝐸E_(0)>0,the polarity of the bipolar EFS structure will be positive/negative;and if𝐸E_(0)<0,the polarity of the bipolar EFS structure will be negative/positive.However,for a fixed polarity of the EFS,either positive/negative or negative/positive,if the satellite is located at the positive side of the EFS,the observed polarity should be positive/negative,if the satellite is located at the negative side of the EFS,the observed polarity should be negative/positive.Therefore,we provide a method to clarify the natural polarity of the EFS with observed polarity by satellites.Our results are significant to understand the physical process in space plasma with the satellite observation.

Spectral deconvolution analysis on Olivine-Orthopyroxene mixtures with simulated space weathering modifications

Olivine and pyroxene are important mineral end-members for studying the surface material compositions of mafic bodies.The profiles of visible and near-infrared spectra of olivine-orthopyroxene mixtures systematically vary with their composition ratios.In our experiments,we combine the RELAB spectral database with new spectral data obtained from some assembled olivine-orthopyroxene mixtures.We found that the commonly-used band area ratio(BAR,Cloutis et al.)does not work well on our newly obtained spectral data.To investigate this issue,an empirical procedure based on fitted results by a modified Gaussian model is proposed to analyze the spectral curves.Following the new empirical procedure,the endmember abundances can be estimated with a 15%accuracy with some prior mineral absorption features.In addition,the mixture samples configured in our experiments are also irradiated by pulsed lasers to simulate and investigate the space weathering effects.Spectral deconvolution results confirm that low-content olivine on celestial bodies is difficult to measure and estimate.Therefore,the olivine abundance of space weathered materials may be underestimated from remote sensing data.This study may be utilized to quantify the spectral relationship of olivine-orthopyroxene mixtures and further reveal their correlation between the spectra of ordinary chondrites and silicate asteroids.

Space Environment Measurements and Applications on the Tianhe Core Module of China Space Station

Tianhe Core Module of China Space Station(CSS)equips a set of instruments consisting of a Particle Fluxes and Solar Activity Detector(PFSAD)and two Atmospheric Density Multi-directional Detectors(ADMDs).The PFSAD is to measure X-rays from the Sun and energetic particles in the low-latitude and low altitude regions,including electrons,protons,and helium ions.The ADMDs are to measure thermospheric atmospheric density.The instruments provide real-time data of the orbital space environment,including solar flares,energetic particle variation and thermospheric density enhancement.All the data contribute to the CSS space weather service for mission control and astronaut’s safety.The paper gives preliminary analyses of the space environment measurements from the PFSAD and the ADMDs.By further analysis,the 1024-channel fine spectra of the solar X-ray can be used to study the mechanism of solar flares and their impacts on the Earth’s atmosphere.Data accumulation will be helpful for analyzing mid-term and long-term variations of the South Atlantic Anomaly and atmosphere density.Furthermore,the data are useful to calibrate previous empirical models and establish new models to study the space environment.

Estimation of surface water content at the Tianwen-1 Zhurong landing site

China’s first Mars rover,Zhurong,successfully landed in the south of Utopia Planitia.The surface water content at the landing area can provide constraints on mineral formation conditions and help us better understand the evolution of the Martian aqueous and geological environment.In this work,the surface kinetic temperature of the Zhurong landing area was derived by analyzing data from the Mars Express Observatoire pour la Minéralogie,l’Eau,les Glaces et l’Activité(OMEGA)spectrometer.Using the Discrete Ordinate Radiative Transfer(DISORT)model,we performed atmospheric correction and thermal correction for the OMEGA data to obtain the surface effective single-particle absorption thickness(ESPAT)parameter to evaluate the surface water content.The surface water content distribution at the landing area was relatively uniform at a lateral scale of~10 km.At the Zhurong landing site,the surface water content in the topmost layer(a few hundred micrometers)of the regolith was 5−8 weight percent water(wt%H_(2)O),assuming surface particle sizes of<45μm,or 1.6−2.5 wt%H_(2)O,assuming surface particle sizes in the range of 125−250μm.The Mars Surface Composition Detector(MarSCoDe)onboard Zhurong also observed significant H_(2)O/OH signals in the landing area.Our results provide an important regional context for the hydration state of the area and can be further verified by the H content derived from the Laser-Induced Breakdown Spectrometer(LIBS)data of MarSCoDe.

Change Ratios of Magnetic Helicity and Magnetic Free Energy During Major Solar Flares

Magnetic helicity is an important concept in solar physics,with a number of theoretical statements pointing out the important role of magnetic helicity in solar flares and coronal mass ejections(CMEs).Here we construct a sample of 47 solar flares,which contains 18 no-CME-associated confined flares and 29 CME-associated eruptive flares.We calculate the change ratios of magnetic helicity and magnetic free energy before and after these 47 flares.Our calculations show that the change ratios of magnetic helicity and magnetic free energy show distinct different distributions in confined flares and eruptive flares.The median value of the change ratios of magnetic helicity in confined flares is-0.8%,while this number is-14.5%for eruptive flares.For the magnetic free energy,the median value of the change ratios is-4.3%for confined flares,whereas this number is-14.6%for eruptive flares.This statistical result,using observational data,is well consistent with the theoretical understandings that magnetic helicity is approximately conserved in the magnetic reconnection,as shown by confined flares,and the CMEs take away magnetic helicity from the corona,as shown by eruptive flares.

Study of the geomagnetic field's regional gradients in Chinese continent using three-dimensional surface Spline model

We combined domestic ground-based and satellite magnetic measurements to create a regional three-dimensional surface Spline(3DSS)gradient model of the main geomagnetic field over the Chinese continent.To improve the precision of the model,we considered the data gap between the ground and satellite data.We compared and analyzed the results of the Taylor polynomial,surface Spline,and CHAOS-6(the CHAMP,?rsted and SAC-C model of Earth’s magnetic field)gradient models.Results showed that the gradients in the south-north and east-west directions of the four models were consistent.The 3DSS model was able to express not only gradients at different altitudes,but also average gradients inside the research area.The two Spline models were able to capture more information on gradient anomalies than were the fitted models.Strong local anomalies were observed in northern Xinjiang,Beijing,and the junction area between Jiangsu and Zhejiang,and the total intensity F decreased whereas the altitude increased.The gradient decreased by 21.69%in the south-north direction and increased by 11.78%in the east-west direction.In addition,the altitude gradient turned from negative to positive while the altitude increased.The Spline model and the two fitted models differed mainly in the field sources they expressed and the modeling theory.