Ground magnetic survey on Mars from the Zhurong rover

Mars’ magnetic field has been measured at large scale by orbiting spacecraft and at very small scale via Martian meteorites. Here we report on a ground magnetic survey on metre to kilometre scales. The Zhurong rover made vector measurements at 16 sites along a 1,089 m track in the Utopia Basin on Mars. It recorded an extremely weak magnetic field, with an order of the average intensity less than that inferred from orbit, in contrast to the large magnetic field in Elysium Planitia measured by InSight. A spacecraft measurement samples an area with radius comparable to its altitude, while a ground measurement samples an area with radius comparable to the depth of the magnetized body. The weak magnetic field measured by Zhurong indicates no magnetization anomalies for a depth of many kilometres around and below the rover’s traverse. We suggest two possible explanations for the weak magnetic field: the entire Utopia Basin may have remained unmagnetized since its formation about 4 billion years ago or that the 5-km-radius ghost crater where Zhurong landed may have been been demagnetized by impact. The Zhurong rover measured in situ the Martian magnetic field vector at 16 sites. An average horizontal intensity of 11.2 ± 10.9 nT is recorded, substantially lower than the ~80 nT measured from orbit. Demagnetization over the entire Utopia Basin or at small scales around the Zhurong site might explain this difference.

of the spectrum, palaeomagnetic analyses of samples from a Martian meteorite (ALH84001) revealed evidence of the ancient magnetic field at very small length scales (~1 mm), supporting the operation of an early Martian dynamo before 4.1 Ga or possibly 3.9 Ga (ref.7).However, the origin of Martian crustal magnetization and the history and attributes of the ancient dynamo remain topics of active interest.Many of the outstanding key questions require new high-resolution regional data sets such as local magnetic field measurements 8 .The NASA InSight fluxgate magnetometer (IFG) measured the ground-level magnetic field at a single landing site 9 .The surface magnetic intensity (~2,000 nT) measured by IFG is very high by the standard of a terrestrial crustal field  47 .The rover landed in the southern Utopia Basin (marked) at 25.066° N, 109.926°E, northwest of the InSight landing site.b, The modelled total magnetic field strength at an altitude of 200 km (ref.14).c, The corresponding geologic map 48 whose key for the colours and geologic units is provided as Supplementary text.Note that resurfacing events in localized areas, such as around the Zhurong landing site, were not indicated in this standard geologic map.Major impact basin boundaries of Isidis (Is) and Utopia are shown with solid black (topography 49,50 ) and dashed black (demagnetization 21 ) lines.

Article
https://doi.org/10.1038/s41550-023-02008-7onboard the Zhurong rover consists of two identical tri-axial fluxgate magnetometers (Fig. 2a) with a dynamic range of ±65,000 nT, a noise level at 0.01 nT/√Hz @1 Hz.One is mounted on the mast, and the other at the intersection of the mast and the rover deck.The upper fluxgate can rotate about a vertical axis.The horizontal magnetic field associated with the rover is removed by rotating the upper fluxgate while keeping the lower one fixed and by repeating the process with the rover pointing in different directions.It should be noted, however, that it is not possible to calibrate the vertical component in this way, so only the horizontal components are discussed in this paper.
The RoMAG measures a combination of the internal fields, rover-generated interference and the external fields.The convention in geomagnetism is to refer to the crustal magnetic field (and core dynamo-generated field where one exists) as 'internal' and those originating in the upper atmosphere and beyond (ionosphere, magnetosphere, Sun, etc.) as 'external'.In this study, rover-generated fields were removed by a two-step calibration method, and the external contribution of the measurements is estimated using observations from the InSight lander and an ionospheric dynamo model (Methods).

Results
Magnetic measurements in this survey were made at 16 locations along a 1,089 m track during 90 sols from 4 June to 3 September 2021 (Extended Data Fig. 1) in a small ghost crater 12,13 of about 5 km in radius (Fig. 2b) whose centre is located at 25.12° N, 109.97°E. The calibration process, which involves the mast rotations and rover rotations, was conducted to derive the Martian field along its traverse (Extended Data Figs. 2  and 3).The measured magnetic fields of crustal origin were then acquired by removing the external contribution from the Martian fields at the 16 locations (Extended Data Fig. 4).The track of the magnetic survey is depicted in Fig. 3a, and the results of the crustal fields are presented in Fig. 3b-d.The directions of the estimated crustal magnetic field are predominantly southeast.The amplitudes of the north-south horizontal component typically range from 3.4 to 24.0 nT after removing the external field and contributions from the rover itself (Methods).The total intensity, which excludes the uncalibrated vertical component, ranges from 5.2 to 39.8 nT, while the average of horizontal intensities is 11.2 nT with a standard deviation of 10.9 nT.The measurements show a significant variation of the surface field in magnitude and direction along the Zhurong traverse.In the first segment (about 0-200 m, sols 21-50) of the traverse, the magnetic field is extremely weak (<10 nT) and close to zero.In the second segment (about 200-600 m, sols 58-79), the magnetic field is moderately weak at about 30 nT.In the final segment (about 600-1,000 m, sols 87-110), the magnetic field becomes extremely weak again (<10 nT) and close to zero.This suggests that the Zhurong rover landed in an area of exceptionally low magnetic field and that the scale size of the local magnetization is on the order of several hundred metres, which can only be revealed via a ground magnetic survey.One should bear in mind that some, or even all of, the variations may be due to ionospheric contributions that deviate from those estimated in this paper (Methods).The possibility of zero magnetization beneath Zhurong's traverse exists in this sense.
The intensities of the magnetic field measured by RoMAG are much smaller than estimates from the latest orbital data-based model 14 (we shall sometimes refer to ref. 14 as L19 herein) downwards continued to the planet's surface (intensity of 81 nT and horizontal intensity of 55 nT).By contrast, NASA InSight's single-site measurement on Elysium Planitia was about ten times stronger than the orbital prediction (about 300 nT) 9 .The very low values at the Zhurong magnetic survey region are also unlikely to be found in the crustal field of many places on Earth, particularly above volcanics containing magnetic minerals.However, on Mars, there is no induced magnetization because there is no dynamo field at present.Remanence may be absent because of demagnetization or because the rocks were formed after dynamo action ceased.
Zhurong's landing site lies in a geologically simple terrain within the giant impact basin Utopia Planitia in the Northern Hemisphere, where the overall magnetic field is weak 15,16 .The surface of the basin is exceptionally smooth with low crater density, which indicates a surface age of Hesperian to Amazonian 17 .The inferred stratigraphic structure (Fig. 2b) shows that the ancient Noachian basement at a depth of several kilometres has been flooded and filled with Hesperian-aged volcanic lavas with a 1-2 km thickness.The late Hesperian-aged Vastitas Borealis Formation (VBF) sediments, interpreted to be the residue of outflow channel effluents derived from the southern highlands or as a sublimation residue of the putative ocean, have been emplaced on the volcanic Hesperian plain.It is also evident that resurfacing has occurred in the landing area in the middle Amazonian (that is, ~1.6 Ga), which may Zhurong's landing site (https://www.nasa.gov/mission_pages/MRO/spacecraft/sc-instru-ctx.html) and the inferred stratigraphy with estimated layer thickness and a ghost crater beneath Zhurong's landing site (red dot), with the black dashed circle marking the ghost crater previously identified 12,13 . https://doi.org/10.1038/s41550-023-02008-7 represent a volcanic flooding event.Crater counting analysis indicates that the filling of VBF and volcanic materials could have transformed Hesperian craters into ghost craters by 3.45 Ga (ref.18).Zhurong rover radar has revealed Amazonian material tens of metres thick, covered with a regolith layer several metres thick 19 .

Discussion
A straightforward way of explaining an extremely weak magnetic field measured by the Zhurong rover is that the entire Utopia Basin, and the crust beneath it, may have remained unmagnetized since the formation of the Utopia impact basin about 4 billion years ago 20 and that the landing site is representative of surrounding area.The demagnetization radius of the Utopia impact is estimated to be 1,600 km, and the distance of the Zhurong landing site to the impact basin centre is about 1,200 km (ref.21).In other words, the Zhurong measurements were made within a large, demagnetized region in Utopia Basin.The lack of strong magnetic signals above the large impact basins from orbital altitude is obvious.However, up to now, it could not unambiguously be argued that the field is close to zero on the surface.The Martian surface fields measured by Zhurong impose a strong constraint on the remagnetization processes of the crust underlying Utopia Basin after the region was pervasively demagnetized by the giant basin-forming impact 21 .Since negligible subsequent demagnetization is expected from burial by post-dynamo lavas 22 or younger sediments, the extremely weak surface fields from RoMAG measurements suggest that the Hesperian lava flows were never magnetized and that the dynamo may have ceased at or before the original Utopia basin-forming impact in the Noachian.Alternatively, demagnetization associated with the impact in the presence of a dynamo field has been suggested 23 .Further, the possibility remains that the region near RoMAG was demagnetized by a later sizeable impact whose signature is buried.This possibility is discussed below.
Another way of explaining an extremely weak magnetic field measured by the Zhurong rover is that the 5-km-radius ghost crater 12,13 where Zhurong landed was demagnetized by the impact that formed the original crater.Ghost craters are impact craters that have been buried by younger materials such as lavas or sediments.They are often identified by circular structures in the overlying materials, typically in the form of concentric double grabens thought to be formed from compaction around the inferred rim of the buried crater 12 .The double-ringed structure in the Zhurong landing region (Extended Data Fig. 5a) suggests that the landing site is 1-2 km from the rim of the ghost crater 13 , which may have been formed from an impact that penetrated the Hesperian lava units, subsequently buried by VBF (Fig. 2b).How can one explain the difference between the Zhurong measurements and orbital predictions at ground level?Satellite measurements resolve length scales that are 100 times more than the Zhurong's traverse.Either Zhurong landed in an area of exceptionally low magnetization, or the orbital predictions represent the limit of accuracy of the downwards continued magnetic field.If the former, the low-magnetization region could be the ghost crater demagnetized by impact.If the latter, the 50-100 nT prediction could be the result of amplification of errors by the downwards continuation process, which is particularly large for the harmonics of highest degree.The problem of assessing the detailed interpretation of orbital fields at ground level for Utopia and other basins will form a separate future study.
Next, we investigate possible sources of crustal magnetization that would be detectable by the Zhurong rover.The satellite field is often downwards continued to ground level to give an average of the surface field when in situ surface observations are unavailable.This method has the advantage of generality, being independent of the depth of any magnetic source, but it is well known that a magnetic field model constructed from orbital data cannot reliably predict the surface magnetic field strength at a point.Downwards continuation can be viewed by accumulating contributions from an area of interest and with different weights determined by the relevant Green's function for the Laplacian 24,25 .The area of interest is where the weighting is significant, a circle with a radius of order the altitude of the satellite (150 km). Figure 4a-d shows the horizontal magnetic field (see Extended Data Fig. 6 for the corresponding vertical magnetic field) from L19 downwards continued to different altitudes, with a circle denoting the area of significant weighting.It should be highlighted that magnetic fields at orbital altitudes are dominated by wavelengths of crustal magnetism that are comparable to the distance from the source 2 .Downwards continuation of a model from this altitude will thus necessarily miss the shorter-wavelength components that dominate the surface field.Given that the L19 model includes data between 120 and 170 km, these wavelengths (clearly visible in the patterns in Fig. 4) dominate the predicted field at the surface.The crustal magnetic field intensity predicted by the L19 model at 150 km is about 5 nT, which is less than the MAVEN r.m.s.misfit (~9.42 nT).Thus, the crustal field in the L19 model may be primarily noise and consistent with either zero or extremely small crustal magnetic fields at orbital altitude.If there is strong but small-scale magnetization (say, <10 km) around Zhurong's traverse, this will not be clearly visible at orbital altitudes but could cause very strong surface magnetic fields.This is likely the case for the InSight landing site.However, if there is weak or zero crustal magnetization in the region, the surface magnetic field will of course be weak, and the orbital field will also be weak and dominated by non-crustal noise which will be amplified when downwards continued to the surface.A helpful way of understanding the RoMAG measurements is to construct a specific magnetic forwards model of the ghost crater, which has been buried by younger materials such as lavas or sediments.The ghost crater model has the advantage of specificity.It provides an upper-limit estimation of the crust magnetization beneath the landing site constrained by the Zhurong surface measurements.We constructed a forwards model with a weakly magnetized uniform layer (with dimensions of 1,000 km × 1,000 km × 1 km; north, east and vertical components of magnetization M N = 0 A m −1 , M E = 0.1 A m −1 and M Z = 0.1 A m −1 , respectively) representing the layer of Hesperian lava flows.A cylindrical hole was placed in the centre to represent the ghost crater that penetrated the lava layer (Extended Data Fig. 5b, lower panel).The horizontal magnetic field is very weak (~6 nT) at 200 m above the crater centre, as predicted by the forwards model.The strongest horizontal strength (up to 41 nT) along the rim of the hole (Extended Data Fig. 5b, middle panel) is comparable to the maximum strength of the value observed by Zhurong.The weak magnetic field shown by this forwards model is consistent with the crust beneath the Zhurong landing region, which is weakly magnetized at about 0.14 A m −1 (magnitude of magnetization of the uniform layer).The magnitude is likely to represent an upper limit of the crust magnetization beneath the Zhurong rover that is far less than the magnetization of the Southern Hemisphere (~10 A m −1 ) 26 .
The crustal magnetic fields of Mars have been critical in dating the duration and timing of its ancient dynamo [26][27][28][29] .Nonetheless, challenges and unsolved problems remain.The interpretation so far has necessarily been simplistic, with absence of magnetic anomalies indicating structures that were either never magnetized because they were formed when the dynamo was inactive, were subsequently demagnetized by impact or other alteration, or buried by deep, younger lavas erupted when the dynamo was inactive.Virtually all the magnetic interpretation up to now has come from satellites some 150 km or more above the surface of Mars, limiting the spatial resolution to worse than 100 km.It is not surprising, therefore, that satellite predictions are different from point surface measurements because spatial variability on small length scales and/or weak magnetization cannot be captured.
Zhurong's weak values can similarly be explained as natural variations.However, the model derived from orbital magnetic field data predicts anomalies of about 73 nT, which are much larger than the ~10 nT recorded by Zhurong.They could arise from deeper layers that were magnetized at formation.It is possible that an ancient dynamo was operating when the dichotomy formed 27 and some of the weakly magnetized basement was preserved after the Utopia Basin impact.Any new lava flows could be magnetized by the underlying anomalies and would otherwise mask the deeper anomalies.Our crater model requires weak magnetization of the sort that could result from the basin-forming impact, owing to either the absence of a magnetic field during formation 1 or processes associated with the impact such as excavation 27 or cooling in a polarity-changing dynamo field 23 .Indeed, weak or zero magnetization that could occupy most of Utopia Planitia could easily come from the impact that formed the basin, which occurred under the hypothesis that there was no global dynamo magnetic field at that time.Indeed, the magnetization of the highly magnetizable early Hesperian lavas had to already have been weak before the impact that formed the ghost crater, which may provide further constraints on the evolution of the Martian dynamo.
We have started ground magnetic surveys on Mars and entered a new era of more detailed interpretation, using scientific methods similar to those successfully employed for ground magnetic surveys on Earth over the last 100 years.The landing site in a ghost crater is indeed fortunate because, unlike an exposed crater with a steep rim that is potentially hazardous to cross, Zhurong may be able to drive over the outer rim and observe magnetic anomalies that are larger than seen so far.The original 90-sol mission has been completed, but Zhurong is capable of more exploration.In our preliminary extended mission plan, Zhurong will be deployed by further southwards travel of about 20 km, approaching and crossing the outer crater rim and visiting other small volcanic features that may have stronger magnetic anomalies.These would fill the gap in our knowledge on the magnetization of smaller craters and, more importantly, the processes by which all craters are magnetized and demagnetized.
Together with Tianwen-1 orbital observations of space magnetic fields and plasma environments, the Zhurong RoMAG further surface measurements will offer an extraordinary opportunity to study the history of the Martian environment.Concurrent observations between Zhurong RoMAG and InSight IFG may be essential to elucidate the variations of ionospheric currents at different local times and to identify sources of time-varying magnetic fields.

Instrument and operation modes
Full details of the RoMAG instrumentation can be found in a previous study 11 .The RoMAG is a set of magnetometers with two identical tri-axial fluxgate sensors, one mounted near the top of the mast of the Zhurong rover, and the other fixed on the mast base.The fluxgate magnetometers measure magnetic fields with a dynamic range of ±65,000 nT with a noise level of 0.01 nT −1/2 at 1 Hz, and a normal sample rate of 1 Hz and capacity of 32 Hz in burst mode.The first RoMAG measurements were implemented in the rover static mode on 4 June 2021, which occurred at the first observation point, 5.14 m away from the lander.During the following 90 sols, RoMAG acquired raw data of magnetic vector fields at another 15 observation points over a total distance of 1,089 m.
After stopping at the first observation point, the rover started science operations.The starting time and operation sequence on June 2021 are presented in Supplementary Table 1.During the period of mast rotation, RoMAG was the only instrument operating.The same operation sequence was conducted at all 16 observation points, except for slightly different angles and durations for each yaw rotation of the mast (Extended Data Fig. 1).In addition to the 16 regular explorations, we conducted an along-track calibration (Supplementary Table 2) on January 2022 for one time, during which both the mast and rover rotations were implemented at a fixed location to make a cross calibration.The whole along-track calibration process lasted about 1 h, and the RoMAG was also the only instrument operating.

Removal of rover fields and determination of Martian fields
The measured magnetic fields (B o ) were the sum of three contributions: the ambient field (B Mars ), the field originating from the rover mast system (B mast ) and the field generated by the rover body (B body ) including fields from solar array currents.To determine B Mars , we selected data during periods of mast rotations to remove rover fields from the RoMAG measurements that had been processed through a calibration procedure.We only deal with the observed magnetic fields from the top sensor as the bottom one did not rotate with the mast.In addition, owing to the closer distance to the rover body and the mast rotation motor, the dynamic interferences in the bottom sensor were apparently stronger than for the top one, especially during the mast rotations (not shown).Therefore, below we determine the Martian fields based only on the measurements from the top sensor.
It was shown that the measured magnetic field (B o ) is on the order of ~10,000 nT and varies by several hundred nanotesla for each detection.The Multispectral Camera (MSCam), which is also mounted on the mast, is quite close to the top sensor.Although the MSCam is switched off during the mast rotation, magnetic fields from the permanent magnet in the MSCam primarily contribute to the measured magnetic field.There is a steering mechanism that controls the pitch angle of the camera.The pitch angle varies at different detection locations to get a better view, but it does not vary during the mast rotation.Therefore, the Article https://doi.org/10.1038/s41550-023-02008-7magnetic fields generated by the permanent magnet remain steady during the mast yaw rotation, but vary for each measurement at different sites owing to the change of the pitch angle of the camera.For each mast rotation, B Mars (Mars geographic frame), B body (rover frame) and B mast (sensor frame) are constants and can be sufficiently separated by using a simple analytical model with relative sensor rotations, as is common practice in the in-flight calibration of fluxgate magnetometer offsets 5,30,31 .We construct a two-step procedure to separate the three fields using manoeuvres of the mast and the rover.When the mast rotates around its z axis, the B mast in the sensor frame and the combined field (B Mars and B body ) in the rover frame are constants (Extended Data Fig. 2a).The field vector rotations generate smooth variations in their horizontal components (B j , B k ) that are recorded by the RoMAG sensor (Extended Data Fig. 2c), which is presented as a circular rotation in the hodograph of B j − B k (Extended Data Fig. 2d).Using a similar method, we can further separate B Mars and B body using horizontal rotations of the whole rover, in which only B Mars changes its direction in the rover frame (Extended Data Fig. 2b).Considering that the rover rotation was only done once on 5 January 2022, it has to be replaced by comparing the relative changes of the rover's heading directions at observation points in the regular detection, as the RoMAG was conducting measurements only when the rover was still.
To be specific, by assuming that the combined field (B c ) of B Mars and B body near the sensor does not vary during the mast rotation, the observed magnetic field B o in the frame of the rotating sensor can be expressed with a linear equation as where B mast is in the sensor frame (i, j, k), B c is in the rover frame (x, y, z) and C is the rotation matrix defining the transformation from the rover frame to the sensor frame, given by where the Euler angles α, β and γ describe this transformation based on a z-y-x rotation sequence and γ is the mast rotation angle.As the sensor rotates with the z axis, C can be simplified to then equation ( 1) can be rewritten as where the subscripts i, j, k and x, y, z denote corresponding components of the sensor frame and the rover frame, respectively.Moving B mast,j and B mast,k in the second and third equations in equation (4) to the left-hand side respectively, squaring both sides and adding the two equations together, we get This equation represents a circle, where (B mast,j , B mast,k ) are the coordinates of the centre of the circle in the B j − B k plane and √ B 2 c,x + B 2 c,y is the radius.
Expanded equation ( 1), we can reduce the problem to a linear multivariable equation below and solve it by the least-squares method as follows: where N is the number of observations during each mast rotation, C i is the rotation matrix for the ith observation derived by using equation ( 2) and B o i is a 2 × 1 vector of the measured horizontal component for the ith observation.We only consider the horizontal components because the rotations take place in the horizontal plane and thus the vertical component of B c and B mast cannot be decoupled in this case.The four variables that can be calculated by the minimization are B c,x , B c,y , B mast,j and B mast,k , which give the horizontal components of the combined field (the vector sum of B Mars and B body ) and the mast field.
We applied this method to calculate B c (the vector sum of B Mars and B body ) at each of the 16 Zhurong locations in addition to the along-track calibration on 5 January.That is, the measurements from 16 mast rotations give us 16 values of B c (Extended Data Fig. 1).The field component variations of 16 measurements are consistent with the modelled rotation vectors in the horizontal plane (B j -B k plane).There is good agreement between hodographs of combined horizontal fields at the 16 observation points and the model-predicted circles of rotating B c (Extended Data Fig. 2d).Differences in the roundness of above 0.82 (where the roundness of a perfect circle is 1) are probably due to the spatial gradient of B body within the rotation plane of the sensor, corresponding to a value up to 7.37 nT in the rotation plane.
In the second step, when the mast conducts at least another two rotations with different rover yaw angles at a fixed location, B Mars and B body can be decoupled according to where is the transformation matrix from the geographic coordinates of Mars to the rover coordinates, and ω corresponds to the yaw angle of the rover.Under these circumstances, both B Mars and B body can be considered as constants.B body can be acquired by solving equation (7) with different C M2R and the corresponding B c at the fixed location using the similar method as for equation (1).Thus, B body = (65.6 ± 0.8 nT, −82.6 ± 1.1 nT) in the rover frame and B Mars = (−11.2± 0.5 nT, −6.2 ± 0.6 nT) in the northeast-centre (NEC) frame were determined (Supplementary Table 3).
As an independent check, B Mars can also be determined directly from this same rover rotation (and no mast rotation, shown in Supplementary Table 4) because the vector sum of B mast and B body stays constant in the sensor frame.The B Mars determined by these two procedures (B Mars,x = −11.2± 0.5 nT, B Mars,y = −6.2± 0.6 nT for the two-step mast rotation method and B Mars,x = −10.9± 0.5 nT, B Mars,y = −6.1 ± 0.6 nT for the first two rover rotation) were quite consistent in both direction and strength, indicating that the two-step method was valid.The substantial difference of the Martian field (B Mars,H ≈ 7.2 nT) determined by the third Article https://doi.org/10.1038/s41550-023-02008-7 and the first two rover yaw rotations was possibly due to the different external field contribution, which was not reflected in the errors associated with the individual measurements.In addition, the difference in B Mars between the two methods was likely from the variations of the external field during the along-track calibrations (about 1 h duration).Another magnetic source might come from the spatial uniformity of the rover-generated magnetic fields during rotation of the mast, though the influence may not be large due to the small displacement of the sensor 2 centre.The body field acquired during the along-track calibration was used to calculate B body at the 16 detection locations.Although the RoMAG measurements at the 16 locations were performed following the same sequence, the temperatures and the solar array currents had considerable effects on the measurements because they were made at different local times and the batteries were in different charging stages.Such influences of the temperatures and solar array currents on measurements were also found in magnetic fields of InSight IFG data 32 .We developed a similar algorithm to evaluate the interferences from different solar array currents at different locations to account for the variations of B body .To be specific, we quantitatively analysed the correlation between the measurements and the solar array currents (I A , I B and I C ), as well as the temperatures at the sensor (T S ) and in the electronics unit (T E , the temperature sensor was pasted on the printed circuit board of the electronic unit).Supplementary Fig. 1 shows a 90 min time series of three components of the top-sensor measurements in the rover coordinate and the corresponding solar array currents and temperatures.During this time interval, the magnetometers were the only instruments that were switched on.The magnetic field components changed by more than 100 nT.I A decreased by about 2 A, while I B changed by less than 0.2 A and I C did not change during this period.We found that the coefficient between fitting results and the observation was smaller when utilizing both I A and I B than when using I A alone.Therefore, we only used I A to calibrate the magnetic field measurements in this study.Although the temperature dependence might not be linear over a wide range of temperatures, the linearity is justified over such a small temperature range (where T S went from −38 °C to −29 °C and T E increased from 1 °C to 15 °C).Thus, it is reasonable to assume that the magnetic field measurements of each component were linearly dependent on I A , T S and T E .Each component of the measurements can be expressed as where the C's are the coefficients to be determined by the multiple linear regression and i corresponds to each component.The fitting coefficients for each component are listed in Supplementary Table 5. Supplementary Fig. 2 shows scatter plots between the observed magnetic fields and the fitting results based on the fitting coefficients.It was shown that the observed magnetic fields linearly depended strongly on the combination of temperatures and solar array currents.The correlation coefficients reached more than 0.99, with an r.m.s. of 2.2 and 2.4 nT for the B x and B y component, respectively.It was reasonable to assume that the only interferences resulting the solar array currents were contained in B body .The influences of the temperatures on the measurements were mainly contained in B mast and removed in step 1.This is because the effect of the sensor and electronic temperatures on the magnetometer recording is independent of the orientation of the sensor axis.When the mast rotates, the horizontal magnetic field components due to the un-calibrated temperature effect do not change in the sensor frame as long as the temperature remains constant during the rotation.If interference due to temperature is contained in the body field, the two horizontal components will remain unchanged during the mast rotation.On the other hand, the solar array, which is fixed to the rover body, will form a spatial magnetic field distribution and the horizontal components will exhibit a time variation with a sine/cosine-like form in the sensor frame as the sensor rotates with the mast.Therefore, the linear correlation coefficient of I S was used to calibrate B body for the 16 locations on the basis of the solar array current differences between the period of the 16-mast rotation and the along-track calibration.
We made an error estimation for the Martian field according to error propagation analysis for each step.According to equation ( 7), the Martian field can be expressed as where B C and B body are the combined field and body field in the rover coordinates, respectively.C R2M denotes the transformation from the rover to the Martian NEC coordinates.According to the law of uncertainty propagation, the uncertainties in the final Martian fields can be written as where σ B C , σ B ′ body and σ B C B ′ body are the root mean square error (r.m.s.e.) of B C , B ′ body and the covariance of the two variables, respectively.The r.m.s.e. of B C was the fitting error in step 1 (Extended Data Fig. 2e), being calculated as where x i and x o i denote the observation and fitting results (circles).This indicates the degree of deviation from the circle.
A similar method was also used to estimate the errors in step 2, in which B body was derived during the along-track calibration.However, before applying the B body in the 16 explorations, we calibrated B body according to the solar array currents differences as follows: where B ′ body and B body are the body field of the 16 explorations and the along-track calibration, respectively.C 4,i are the fitting coefficients in equation ( 9) of the corresponding components.ΔI are the solar array currents differences between the 16 explorations and the along-track calibration.We first determined the uncertainties of C 4,i ΔI as where σ C 4,i and σ ΔI are the uncertainties of C 4,i and ΔI, respectively.The uncertainties of B ′ body then can be expressed as Substituting σ B ′ body into equation (11), we can obtain the final uncertainties of B Mars .Extended Data Fig. 3

Estimation of the ionospheric contribution at the Zhurong landing site
The diurnal variations of the surface magnetic field were reported based on the InSight magnetic measurements, and the typical Article https://doi.org/10.1038/s41550-023-02008-7peak-to-peak amplitude of the variation in a sol on each component is usually less than 30 nT (refs.9,33,34).
To assess the magnetic field contributed by the diurnal variations mainly produced by Mars' ionospheric currents above the Zhurong rover, we made a comparison by modelling the ionospheric currents over the Zhurong and InSight landing sites and their resulting surface magnetic fields based on the model developed in ref. 35.The average intensity ratio of the modelled surface magnetic field between Zhurong and InSight, along with the measurements at InSight with the same solar longitudes (Ls), was used to estimate the ionospheric current distribution at the 16 observation locations.Although the ionospheric current pattern and its sol-to-sol variability at the Zhurong landing site may be different from those at the InSight landing site owing to their different plasma and neutral environments, we can make an overall comparison of their average intensities.The magnetic field inputs for the model were from MAVEN spacecraft measurements.The background magnetic fields used in the model are from MAVEN spacecraft measurements between 2014 and 2020.Since the magnetic fields above the InSight and Zhurong landing sites vary with local times and altitudes (ref.33 and Supplementary Fig. 3), we calculate mean values of a 10° area in the N-S and E-W directions at each local time (24 local times in total) as the input background magnetic fields.The magnetic fields at each altitude are the averaged values in 2 km intervals (from 130 km to 280 km).
We divided the MAVEN magnetic field measurements from 2014 to 2020 above the Zhurong and InSight landing sites by local time (24 bins) and altitude (from 130 km to 280 km with a 2 km interval).We used the mean values of 10° area in each direction as the model input in the dynamo region (Supplementary Fig. 3).The neutral winds, neutral densities, and electron density above the Zhurong landing site were from the Mars Climate Database (MCD 5.3) 36 .The electron temperature was constructed by fitting the hyperbolic tangent function with a lower and upper boundary of 510 and 3,140 K, respectively 37 .The input parameters above the InSight landing site were from the same data sets as used in ref. 35.We predicted the surface field by assuming the ionospheric current distribution as a steady-state current sheet with large extent in horizontal direction and finite thickness.The surface field was independent from the distance from the sheet and could be given by B surface = μ 0 ∑ j(z) ⋅ , where j(z) is the current density.This assumption will result in an upper limit of the surface field estimation.On the other hand, a line current assumption will lead to a lower limit of the field estimation.The real geometry, of course, lies in an intermediate state.Both the line current and current sheet geometry assumptions were also conducted to estimate the surface field at the InSight landing site 33 .Supplementary Fig. 4 shows a comparison of the modelled mean horizontal (H) components as a function of Ls.It was shown that the average intensity of the H component at the Zhurong landing site was about 83% of that at the InSight landing site during Ls 54-95° (shadow area).This intensity ratio, along with the diurnal variations from InSight IFG observations, was used to evaluate the external fields at the 16 exploration locations.
Next, we investigated the diurnal variations at InSight with respect to local true solar time (LTST) and Ls to evaluate the diurnal fluctuations as well as the sol-to-sol variations within Ls from 54° to 95° (Supplementary Fig. 5).It was shown that the median amplitude during the LTST of interest was less than 15 nT (Supplementary Fig. 5a).The peak amplitudes of H component show variations from 10 to 25 nT during Ls from 54° to 95° (Supplementary Fig. 5b).The peak amplitude, to some extent, represented the upper limit of the amplitude of a sol in the LTST of interest.The standard deviation of the peak amplitudes of B H is about 3.5 nT, which can be treated as the sol-to-sol variability in B H .This could be due to the sol-to-sol variations of the ionospheric currents since both the Thermospheric Neutral Wind 38 and/or the electron densities 39,40 in the Martian ionosphere appear to have sol-to-sol variations.According to the ionospheric current modelling at both Zhurong and InSight, the sol-to-sol variations at Zhurong would also be less than 3.5 nT during the Ls of interest.
We then subtracted 83% of the mean value (north and east component of ionospheric origin, B iono,N = 3.69 nT, B iono , E = −3.8nT) between 10:00 and 16:00 LTST (the interval during which most of Zhurong's measurements were made) observed by InSight at the corresponding Ls from the B Mars observed by Zhurong to obtain the pure crustal fields from location to location.We used the mean value because the differences in the magnetic field variations in a sol as well as sol-to-sol between the Zhurong and InSight landing sites may be considerable owing to their different ionospheric plasma environments.In addition, they were not measured at the same time.We therefore used the mean values observed between 10:00 and 16:00 LTST as a proxy for the external fields at Zhurong.The standard deviation in our estimate of the external field was used, via standard error propagation, to calculate the uncertainties in our derived values of B crust .Note that the uncertainties in B crust due to ionospheric contributions are not completely accounted for in the individual measurements.The B crust and the corresponding uncertainties are shown in the two right columns of Extended Data Fig. 4.
By assuming that the uncertainties in each location obey a normal distribution, we calculate the fraction of the area under each normal curve of the 16 measurements that falls between −5 and 5 nT to indicate whether the possibility of no detectable crustal field exists along Zhurong's traverse.We chose the range of −5 to +5 nT because the average uncertainty of the calculated crustal fields was close to 5 nT.The average probabilities of the 16 curves, which lie in the range between −5 and +5 nT (shadow area) are 21% and 34% for the N-S and E-W component, respectively (Supplementary Fig. 6).

Defining the stratigraphy beneath the landing site
The surface regolith layer of the landing area is composed of pulverized loose materials.Some craters would penetrate through the regolith layer, excavate and emplace larger boulders on the surface.As a result, depending on whether or not rocky materials are present in the crater ejecta, the crater excavation depth (roughly 0.084 times the crater diameter) can constrain the thickness of the regolith layer 41 .We identified ten rocky and ten non-rocky craters with mean diameter between 7 and 82 m in HiRISE images near the landing site and estimated the average regolith thickness to be ~3 m.When catastrophic flooding carrying sediments or younger lava flows resurfaced older units, smaller craters with smaller rim heights would be buried and thus become invisible.This leaves a kink in the measured crater size-frequency distribution in the region, and the diameter of the crater at the kink should correspond to the largest crater whose rim height is equal to the lava flow thickness.In this study, we found the crater diameters that produced the kink in the crater size-frequency distribution to be between ~0.41 and 1.2 km, corresponding to a thickness of the infilling Amazonian materials of around 30-50 m based on the equations in ref. 42.There is a 24.9-km-sized ghost crater near the landing site 18 .The estimated thickness of the filling materials is around ~297 m.It has been hypothesized that the ghost craters have been mostly filled with VBF materials, thus the thickness of the VBF at the landing site is at least 297 − 50 − 3 = 244 m.Estimating the thickness of the Hesperian lava flows at the landing site is challenging as there are no obvious exposed basaltic materials exposed by impact detected by orbital remote sensing observations.A similar challenge is presented by estimating the thickness of the hydrated materials overlying the Noachian basement.Here, we adopted an average thickness of the Hesperian lavas of 1-2 km and a depth of the hydrated mineral layer overlying the Noachian basement of around 4-5.5 km, on the basis of a previous study 43 .
A ghost crater, centred at 25° 6′ N, 109° 54′ E, in Zhurong's landing region was previously identified 12,13 .It was also recognized from the combined images from the context camera (CTX) (Extended Data Fig. 5a).We estimated a diameter of ~10 km for the ghost crater with a depth of ~1 km based on the depth-diameter ratio criterion 44 .

Magnetized layer forwards modelling
We construct a forwards model from a homogeneously magnetized layer with a cylindrical hole placed in the layer representing the Hesperian lava penetrated by a crater where the Zhurong rover is located (Extended Data Fig. 5b).The hole has a diameter and a thickness of 10 and 1 km, respectively.The magnetization vector of the layer [M N , M E , M Z ] = [0, 0.1, 0.1] A m −1 is assumed to be uniform, being less than the maximum magnetization intensity derived from the Mars Global Surveyor (MGS)-based crustal model in the Utopia Basin 26 .The magnetized layer is placed 200 m beneath the Martian surface.We calculate the magnetic field at the Martian surface from the layer divided into 364 triangular prisms.The magnetic field from each prism is calculated based on the polygon method 2,45 .The Martian surface field can be obtained by integrating each prism field.Similar models are also used to investigate the magnetic variations above large craters 46 .

Fig. 1 |
Fig. 1 | Maps of Mars near the landing site in Mercator projection.a, Mars Orbiter Laser Altimeter topography 47 .The rover landed in the southern Utopia Basin (marked) at 25.066° N, 109.926°E, northwest of the InSight landing site.b, The modelled total magnetic field strength at an altitude of 200 km (ref.14).c, The corresponding geologic map 48 whose key for the colours and

Fig. 2 |
Fig. 2 | The Zhurong rover and the relevant geologic context.a, The Zhurong rover (with dimensions of approximately 3.3 m × 3.2 m × 1.85 m) and its two three-component fluxgates 11 , one at the top of the mast, which can rotate, and a fixed one at the intersection of the mast and the rover deck.b, The CTX image of

Fig. 3 |
Fig. 3 | RoMAG measurements along the Zhurong traverse.a, Measured horizontal magnetic fields components of crustal origin at 16 sites along the Zhurong traverse overlaid on an image of surface features from Tianwen-1's high-resolution imaging camera.Coloured dots and numbers show the locations and the sols of magnetic observations, respectively.Black lines indicate horizontal components of the magnetic vector.The values of the measurements along the Zhurong traverse at 16 sites are shown in Extended Data Fig. 4. b,c, The northwards (b) and eastwards (c) magnetic component at 16 observation points in the geographic frame of Mars.The error bar overplotted on each of the 16 points is the uncertainty calculated from the error propagation (Extended Data Fig. 4).d, The distribution of horizontal components in the geographic frame of Mars.The black star marks the vector mean of the horizontal fields at 16 locations.Colour codes in a-d are the same.

Fig. 4 |
Fig. 4 | The satellite-based magnetic field model (L19) and areas influenced by measurements.a-d, Maps of the horizontal magnetic field strength in the Zhurong landing region predicted by the L19 model 14 downwards continued to four different heights of h = 150 km (a), 100 km (b), 50 km (c) and 0 km (d).The corresponding field strength at the Zhurong landing site (25.066°N, 109.926°E; star) is approximately 2.6 nT at h = 150 km, 5.6 nT at h = 100 km, 16.2 nT at summarizes the flowchart that describes each calibration step along with the uncertainty propagation.The Martian surface magnetic fields B Mars at the 16 observations were easily determined by subtracting B body from the combined field B c at the corresponding observation site.The calibrated values of B body and hence B Mars are shown in tabular form in Extended Data Fig. 4. The final B crust along with the uncertainties was acquired by subtracting the external fields estimated from InSight observations from B Mars .

Fig. 1 |
List of Zhurong rover operation parameters during RoMAG measurements.The universal time (UT) dates, days on Mars (Sols), solar longitudes (Ls), local mean solar times (LMST) at the beginning of mast rotation, locations of the Zhurong rover relative to the lander level (LL) frame in which X points north, Y east, and Z down, yaw angles of the rover (Yaw), rotation angles of the mast (θ), and the calculated horizontal strength of the combined field (B C,H ) are provided.Extended Data Fig. 2 | Combined horizontal magnetic fields B c measured during the mast rotations at 16 observation points.a, Sketch of the rotating combined vector B c = B mars + B body , and the fixed field B mast in the sensor frame.b, Sketch of the rotating vector B mars , and the body field B body in the rover frame.c, Variations of field components along with the modeled fields (black lines) in the sensor frame as a function of mast yaw angles during the mast rotations.d, Hodograph of the horizontal components of B c on the B j -B k plane in the sensor frame, where the beginnings and ends of the interval are marked by squares and triangles, respectively.e, The mean fitting errors of the horizontal component of B c (sample size n = 153) with the 95% confidence interval (error bars).Extended Data Fig. 3 | Data processing flow chart.The combined fields (B C ) in the rover frame are determined by the mast rotation.C R2M and C M2R are the transformation matrices between the geographic coordinates of Mars and rover coordinate, in which ω is the yaw angle of rover.B mars is the ambient field and B body is the field generated by the rover body.I A , T S and T E are the solar array current, sensor temperature and electronic temperature, respectively.C are the coefficients to be determined by the multiple linear regression and i corresponds to each magnetic field component.∆I are the solar array currents differences between the 16 explorations and the along-track calibration.B ′ body are the magnetic fields generated by the rover body and calibrated by ∆I and C 4 .Extended Data Fig. 4 | The calculated rover body field, Martian magnetic fields, and crustal fields along with the uncertainties, respectively.The universal time (UT) dates, days on Mars (Sols) are listed in the first two columns.Extended Data Fig. 5 | The ghost crater in the Zhurong's landing region.a, Combined images from Context camera (CTX: D22_035786_2060_ XN_26N250W and CTX: F04_037553_2068_XN_26N250W) enhanced to better identify the ghost crater.The white and black dashed lines indicate the inner and outer rings of the double-ringed structure of a ghost crater, respectively.The Zhurong's landing site and route are denoted by the red dot and the blue curve, respectively.The horizontal white bar indicates the scale bar of the map.b, Modeling results of the ghost crater showing the distributions of vertical magnetic component and horizontal intensity in the top and middle panels, respectively.The cylindrical hole, which has a diameter of 10 km, corresponds to the inner ring of the ghost crater.Extended Data Fig. 6 | The satellite-based magnetic field model (L19) and areas influenced by measurements.Panels (a-d) show maps of the horizontal intensity in the Zhurong landing region predicted by the L19 model 6 downward continued to four different heights h = 150 km, 100 km, 50 km, 0 km.The corresponding vertical field strength at the Zhurong landing site (25.066°N, 109.926°E, the stars) is approximately 4.8 nT at h = 150 km, 9.2 nT at h = 100 km, 21.3 nT at 50 km, and 57.1nT at the Martian surface.The black circle at h = 0 km denotes the approximate ground area averaged by the orbital model, radius 150 km.It includes some of the stronger anomalies.The corresponding circle for Zhurong's ground measurement, assuming a magnetised layer of depth 1 km, is too small to be shown on this scale.The gray curve indicates the topography boundary of Utopia basin at the surface altitude.