Technical Paper

Korea Pathfinder Lunar Orbiter Magnetometer Instrument and Initial Data Processing

Wooin Jo1, Ho Jin1,, Hyeonhu Park1, Yunho Jang1, Seongwhan Lee2, Khan-Hyuk Kim1, Ian Garrick-Bethell3, Jehyuck Shin4, Seul-Min Baek4, Junhyun Lee1, Derac Son5, Eunhyeuk Kim6
Author Information & Copyright
1School of Space Research, Kyung Hee University, Yongin 17104, Korea
2NARA Space Technology, Seoul 07245, Korea
3Department of Earth and Planetary Sciences, University of California, Santa Cruz, CA 95065, USA
4Korea Astronomy and Space Science Institute, Daejeon 34055, Korea
5Sensorpia Inc., Daejeon 34054, Korea
6Korea Aerospace Research Institute, Daejeon 34133, Korea
Corresponding Author : Tel: +82-31-201-3865, E-mail:

© Copyright 2023 The Korean Space Science Society. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Nov 02, 2023; Revised: Nov 10, 2023; Accepted: Nov 13, 2023

Published Online: Dec 31, 2023


The Korea Pathfinder Lunar Orbiter (KPLO), the first South Korea lunar exploration probe, successfully arrived at the Moon on December, 2022 (UTC), following a 4.5-month ballistic lunar transfer (BLT) trajectory. Since the launch (4 August, 2022), the KPLO magnetometer (KMAG) has carried out various observations during the trans-lunar cruise phase and a 100 km altitude lunar polar orbit. KMAG consists of three fluxgate magnetometers capable of measuring magnetic fields within a ± 1,000 nT range with a resolution of 0.2 nT. The sampling rate is 10 Hz. During the originally planned lifetime of one year, KMAG has been operating successfully while performing observations of lunar crustal magnetic fields, magnetic fields induced in the lunar interior, and various solar wind events. The calibration and offset processes were performed during the TLC phase. In addition, reliabilities of the KMAG lunar magnetic field observations have been verified by comparing them with the surface vector mapping (SVM) data. If the KPLO’s mission orbit during the extended mission phase is close enough to the lunar surface, KMAG will contribute to updating the lunar surface magnetic field map and will provide insights into the lunar interior structure and lunar space environment.

Keywords: Korea Pathfinder Lunar Orbiter Magnetometer (KPLO); lunar exploration; magnetometer; payload; data calibration


The Korean Pathfinder Lunar Orbiter (KPLO) was launched on August 4, 2022, by SpaceX’s Falcon 9. Before the lunar orbit insertion (LOI), KPLO completed the trans-lunar cruise (TLC) journey during 4.5 months on a ballistic lunar transfer (BLT) trajectory. In December 27, 2022, KPLO successfully settled down in a lunar polar orbit at an altitude of 100 ± 20 km after multiple LOI maneuvers (Song et al. 2023). Since January 1, 2023, KPLO has been operating in this polar orbit, passing over various regions of the Moon per its one-year mission plan.

The major scientific objective is to conduct a magnetic investigation near the Moon. Lunar surface magnetic anomalies are one of our scientific targets (Hood et al. 1979; Baek et al. 2017; Garrick-Bethell & Kelley 2019), with lunar magnetic induction another major scientific objective (Dyal & Parkin 1971; Russell et al. 1981; Hood et al. 1999; Haviland et al. 2019). The ARTEMIS P1/P2 (THEMIS B/C) spacecraft operation provides an excellent reference to study the lunar magnetic induction by KPLO magnetometer (KMAG) observation. The technical objective is a demonstration of multi-sensor operation of the KMAG instrument, which is composed of commercial-grade parts.

The expected KPLO lifetime was one year, but it has been extended over two more years because of sufficient fuel and stable health status. During the first and second years, KPLO operates at an altitude of 100 km, which is too high to investigate the lunar surface magnetic field because the magnetic field intensity is less than 3 nT at that altitude (Richmond & Hood 2008). However, the KMAG shows reasonable observation results when compared with surface vector mapping (SVM) data (Tsunakawa et al. 2010). During passage through the magnetic anomaly regions, KMAG data shows comparable magnitude and shape at the corresponding locations. Moreover, several times we have captured the lunar-induced field phenomenon from the lunar interior driven by rapid changes in the interplanetary magnetic field (IMF). Lunar magnetic induction has been suspected since the Apollo era (Dyal et al. 1976). Similar to the Apollo-era method, we use two ARTEMIS spacecraft in a highly elliptical lunar orbit as a reference for the IMF to identify lunar-induced responses in KPLO measurements close to the Moon.

To minimize magnetic interference from the spacecraft, previous lunar missions like Lunar Prospector (LP) and Kaguya placed magnetometers at the end of a boom, with boom lengths of 3.6 m (Binder 1998) and 12 m (Kato et al. 2010), respectively. However, KMAG has a relatively short 1.2 m boom. Consequently, KMAG was equipped with three triaxial fluxgate sensors for multi-sensor measurement and redundancy purposes. Additionally, to monitor the spacecraft-generated magnetic field, magnetoresistive (MR) sensor was installed within the spacecraft body, providing a measurement range of ± 60,000 nT and a resolution of 100 nT.

Utilizing these four magnetometers, KMAG observes the magnetic field around the Moon and processes the raw data into the final data products. The scientific data obtained from the KPLO mission will be released to the public through the Korea Aerospace Research Institute (KARI)’s planetary data system (KPDS) (Kim et al. 2017). Since KPDS complies with the PDS4 standard, used as a standard for planetary exploration scientific data, KMAG scientific data will be open as the PDS4 standard. Accordingly, the scientific data generated is divided into three levels: Raw, Partially Processed (PP), and Calibrated (CAL), accompanied by metadata explaining them as a product.

KMAG’s initial operations started 4 hr after launch. Early operation allowed KMAG to acquire sufficient observation data to confirm data-processing pipelines for each data level during the BLT cruise. To verify the accuracy of KMAG’s data processing results, we compared its observational data with magnetic field measurements from the DSCOVR spacecraft operating in the upstream solar wind. The comparison was performed taking into account the time-delayed responses of magnetic field due to the different observation sites (Maggiolo et al. 2017). Afterward, the data-processing pipeline was applied to KMAG lunar orbit. We also verified the lunar orbit data by comparing it with data from the ARTEMIS P1 and P2 spacecraft, that are in elliptical orbits around the Moon.

In this paper, we present the instrument overview in Section 2 and the data processing in Sections 3 and 4. Sections 5 and 6 describe the results and future works.


KMAG is the first magnetometer instrument developed in South Korea for space exploration beyond the Earth, and it was designed to carry out only the most basic functions. The instrument specifications and performance are described in a previous paper (Lee et al. 2021). In the design stage, the magnetic cleanliness requirement was insufficient to obtain clean observation, and the boom length was also relatively short because of the limitations of mass and volume (Park et al. 2022). Therefore, KMAG has three fluxgate sensors installed inside the boom to use a multi-sensing technique. Based on the observation data from the nominal lunar operation, these attempts are inadequate to obtain a clean signal. Nevertheless, data-reduction processes such as using a high-pass filter and calculating moving averages allow the generation of a reliable dataset.

KMAG consists of two units: the MAG unit and fluxgate magnetometer control electronics (FCE) unit, as shown in Fig. 1. The MAG unit has three identical three-axis fluxgate sensors in the 1.2-meter carbon fiber-reinforced polymer (CFRP) boom that includes a hinge structure. The FCE unit has four boards: a low-voltage power supply (LVPS), an on-board computer (OBC), and digital and analog boards for fluxgate sensor operation. The OBC and LVPS boards have the same dimensions as a 1U CubeSat standard size (Fig. 2), and all electronic parts are commercial-grade except the connector and wire harness between the KMAG and the spacecraft. These major electronic parts are listed in Table 1.

Fig. 1. KMAG instrument configuration. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer.
Download Original Figure
Fig. 2. FCE unit boards. (a) OBC, (b) LVPS, (c) digital, and (d) analog for fluxgate magnetometer. The dimensions of the OBC and the LVPS boards are the same as those of a standard 1U CubeSat. FCE, fluxgate magnetometer control electronics; OBC, on-board computer; LVPS, low-voltage power supply.
Download Original Figure
Table 1. Major parts list of FCE
Board Lib ref Part name Manufacturer
OBC PC-104 (52P) connector ESQ-126-38-G-D Samtec
Supervisory circuits TPS3838K33DBVR Texas Instruments
Mixed-signal microcontrollers MSP430FR5992IPMR Texas Instruments
Triple schmitt-trigger buffer SN74LVC3G17DCTR Texas Instruments
Triple schmitt-trigger buffer SN74LVC1G125DBVR Texas Instruments
Bi-directional, level-shifting TXS0108EPWR Texas Instruments
Line driver AM26LV31ESDREP Texas Instruments
Schmitt-trigger buffer SN74LVC2G17DBVR Texas Instruments
Differential line receiver AM26LV32EMDREP Texas Instruments
Single UART SC16IS740IPW/Q900 NXP Semiconductors
RTC DS1391U-33+ Maxim Integrated
MOSFET IRF7509TRPBF Infineon Technologies
OpAmp LMV324IDT ST Microelectronics
MR sensor HMC1053 Honeywell
Quartz crystal ECS-73-20-4X-DN ECS Inc.
Crystal CMR200T32768DZBT/DZFT Citizen
LVPS Gate SN54LS08J Texas Instruments
OpAmp AD822ANZ Analog Devices
Voltage regulator LM117HV Texas Instruments
DC-DC filter_50V FGDS-2A_50V_Filter Gaia
DC-DC ±5V bipolar MGDB-04-H-C Gaia
DC-DC 5V single MGDS-04-H-C Gaia
HSSR-7111 HSSR-7111 Broadcom
Temperature sensor AD590LH Analog Devices
Analog digital Microcontroller ATmega8L-8AU Microchip Technology
OpAmp AD8676 Analog Devices
ADC AD7193BRUZ Analog Devices

The diodes, resistors, and ceramic capacitors were obtained from various manufacturers, KEMET, Vishay, Panasonics, and Murata. All parts are commercial grade.

FCE, fluxgate magnetometer control electronics; OBC, on-board computer; LVPS, low-voltage power supply.

Download Excel Table
2.1 Fluxgate & Electronics

The OBC board has just a few lightweight functions that include RS422 communication, obtaining data at a rate of 10 Hz, and maintaining time sync with the spacecraft. The main processor used a commercial grade mixed-signal microcontroller. In the case of the flight firmware, we tried to reduce the command and data handling items to carry out essential functions only.

The OBC board has one anisotropic magnetoresistive (AMR) sensor (HMC1053) for the ground and in-orbit test purposes. Because of the temperature variation of the AMR sensor, there is a temperature senor on the board. It has proven to be a good reference (beyond expectations) to monitor large magnetic disturbances from the spacecraft for the in-orbit calibration.

The LVPS board has two main components: the regulated voltage supply and the boom actuator control. The boom actuator circuit consists of two MOSFET Optocouplers (HSSR 7111) to supply high current to a Frangibolt (FC2) actuator as shown in Fig. 3. This circuit has redundancy.

Fig. 3. Boom Frangibolt power control circuit. This is the only one redundant circuit in all of the KMAG’s electronics. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer.
Download Original Figure

The analog and digital boards incorporate three identical fluxgate sensors operating circuits for the three magnetometers, respectively. This means that one damaged circuit does not affect the others. The data acquisition continues during power-on, and the observation data is transferred by the OBC 10 Hz trigger signal. All circuit voltages (± 5 analog, + 5 digital) are provided by the LVPS through an internal PC-104 stacked connector.

The fluxgate magnetometer was developed by Sensorpia Inc. The shape of the sensor core is a tiny racetrack (Son 2012). Each KMAG fluxgate magnetometer has three sensors that are installed in the magnetometer housing for three-axis observations as shown in Fig. 4.

Fig. 4. KMAG fluxgate magnetometer (left). The magnetometer is attached to the sensor guide bar (right). The sensor guide bar makes it easy to align the three magnetometers. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer.
Download Original Figure
2.2 Mechanical Parts and Wire Harness

The mass of the KMAG instrument is 3.5 kg including wire and harness. Most of the mechanical structures use AL6061 alloy. The FCE box has a stack hold structure, as shown in Fig. 5. The electronics boards are assembled in each structure and then these structures are connected by screws. The FCE box has eleven connectors for all of the KMAG instrument function flows including communication with the spacecraft. The wire harness diagram is shown in Fig. 6.

Fig. 5. FCE box: The size is 150 × 135 × 82 mm, and the stack-type structure is a three-layer configuration. LVPS (bottom), OBC (middle), analog and digital boards for the fluxgate magnetometer (top). FCE, fluxgate magnetometer control electronics; LVPS, low-voltage power supply; OBC, on-board computer.
Download Original Figure
Fig. 6. The KMAG wire harness diagram. The harness length is about 1.4 m between the FCE and MAG unit. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer; FCE, fluxgate magnetometer control electronics.
Download Original Figure

The hinge structure includes the boom rotation mechanism, electric connectors, and position sensors. For the rotation angle information, we used one micro switch (11HM1, Honeywell, Charlotte, NC, USA) and one potentiometer (6538, Bourns, Riverside, CA, USA). The deployment locking mechanism involves a locking rod inserted into the rotation axis structure when the boom rotation ends, as shown in Fig. 7.

Fig. 7. Hinge unit assembly. The deployment angle is 135°, and the mechanical spring is used to deploy the boom.
Download Original Figure

As shown in Fig. 8, the actuator assembly uses a TiNi FC2 Frangibolt and one micro switch. The micro switch assembly has a push spring pin to initiate the deployment force. This deployment mechanism was tested twelve times during the development stage. There were two failures, one caused by a defective epoxy bond and the other, by a mechanical design problem. The contact surface of separate plates has a slope for easy release. There are two connectors in this unit: a main power supply that comes from the LVPS, and an enable-plug to cut off the input power during unstable power supply conditions. The enable-plug connector was installed at the final stage just before the launch.

Fig. 8. The actuator assembly uses a TiNi FC2 Frangibolt and one micro switch.
Download Original Figure

The KMAG did not use an instrument base plate (i.e., an interface structure) between the KMAG and the spacecraft because of the mass limitation as shown in Fig. 9. Therefore, we repeatedly carried out a variety of tests such as thermal and vibration to verify the safety of the direct attachment on the spacecraft top panel.

Fig. 9. KMAG installed on the KPLO top plate. The length from the end of the boom to the end of the hinge plate is 1,318 mm. KMAG, KPLO Magnetometer; KPLO, Korea Pathfinder Lunar Orbiter.
Download Original Figure


The KMAG data for over a year, since August 5, 2022, is available from both the TLC phase and nominal mission phase, except for about 20 days. This data has been archived by the KPLO deep space ground system (KDGS) in KARI, which supports the KPLO system rehearsals and mission operations (Song et al. 2021). The science data management subsystem (SDMS), one of the subsystems of KDGS, transmits science telemetry (TM) data. This data is sent to instrument science operation centers (SOCs) by the KPDS sFTP server, with the removal of channel access data units (CADU) header, tail, and padding. In addition, KPDS will release SOCs’ processed scientific data to the public with the PDS4 standard, established by NASA’s planetary data system, widely recognized as the standard for scientific data in the field of planetary exploration.

Hence, KMAG SOC processes the TM data into higher-levels of scientific data in compliance with the PDS4 standard. Ultimately, KMAG SOC generates data processed at the Raw, PP, and CAL levels according to the data levels defined in Table 2.

Table 2. KMAG data levels
Data level Description
TM Raw TM data as received at the ground receiving station (binary).
Raw (Level 0) Converting TM data (binary) to ASCII format and scientific units.
PP (Level 1) The orthogonality of the magnetometer sensor is corrected. In addition, zero-offset determination and coordinate transform are applied.
CAL (Level 2) Research-grade. End-user level data products. This data is the final data calibrated for spacecraft-generated fields.

There are three levels of science data: Raw, PP, and CAL.

KMAG, Korea Pathfinder Lunar Orbiter Magnetometer; TM, science telemetry; PP, Partially Processed; CAL, Calibrated.

Download Excel Table

First, Raw data is extracted from TM binary data. This data contains magnetic field dataset and monitoring dataset. Secondly, the PP data undergoes instrument calibrations such as orthogonality correction and the sensor offset determination. Last, we perform a coordinate transformation correcting for spacecraft orientation using the SPICE toolkit and kernels. At this point, the PP data becomes a useful magnetic field dataset. Finally, the CAL level data (the end-user data) is achieved by calibration of spacecraft noise fields. This data is available for scientific research.

To create a single data product, we generate not only the primary scientific data, but also accompanying supplementary metadata for each level. Hence, we produce metadata in XML syntax; this includes details about the product, mission, and observations, as well as descriptions of the scientific data.

The format of the data product and the data transmission process were verified by creating an interface control document (ICD) and interface test procedure (TP) with KARI before the release of KPDS. Consequently, the processed scientific data and metadata are combined into a single product for multiple dates. The entire data-processing pipeline is shown in Fig. 10. Finally, the data-processing pipeline provides processed magnetic field data products at Raw, PP, and CAL levels.

Fig. 10. The KMAG data-processing pipeline. KMAG TM data are processed at the Raw, PP, and CAL levels, and level products generated with metadata are provided to KPDS. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer; TM, science telemetry; PP, Partially Processed; CAL, Calibrated; KPDS, Korea Aerospace Research Institute (KARI)’s planetary data system.
Download Original Figure


The final goal is to process the scientific data into CAL data, which is research-grade. Data calibration involves two main procedures: ground calibration, which is performed before launch, and in-flight calibration, which is carried out during the cruise phase. On the ground, we conducted calibration tests related to the characteristics of the instrument itself. From the ground calibration test utilizing Helmholtz coils, we derived equations that correct sensor offset and non-orthogonality. Additionally, we implemented a thermal test to evaluate the temperature dependency of the sensor. The temperature offset was found to vary by approximately 0.1 nT/10°C. However, the temperature calibration differs slightly depending on the slope of the temperature variation. In the lunar orbit, the temperature of the magnetometer sensor fluctuates by up to 10°C per orbit, resulting in a minor effect on the magnetic field change. Therefore, thermal correction was not applied to the initial public data.

Corrections of some parameters that are difficult to predict during ground testing are performed by in-orbit calibration. The elements of in-orbit calibration performed during the cruise are sensor offset variations, spacecraft orientation changes, and spacecraft noise fields. Since the study on removing spacecraft-generated magnetic fields is in progress, it will be covered in a later paper on noise removal.

4.1 Offset Determination

In the initial measurements of KMAG, the maximum offset among the three magnetometers reached ~70 nT, which is believed to be caused by the magnetization effects from long-term storage prior to launch and the different operational environment. The zero offset determination is one of the most important stages in the entire data processing because the offset can change over a long period and an inaccurate offset can lead to a totally different result. To find the offset, previous missions employed various methods. For example, there are the spin-averaging technique (Farrell et al. 1995; Kepko et al. 1996) and the minimum variance analysis in the Alfvénic characteristics of solar wind, called the Davis-Smith (D-S) method (Davis & Smith 1968; Belcher 1973).

The KPLO is a three-axis stabilized spacecraft without spin. Thus, the KMAG adopts the D-S method, used in previous missions such as Mariner 5 (Belcher & Davis 1971) and Kaguya (Takahashi et al. 2009), which is applicable to KMAG employing three sets of three-axis fluxgate sensors. In the case of KMAG, the D-S method is applied individually to each fluxgate sensor to derive offsets along the three axes. This method can also be utilized as a correlation technique, as suggested by Leinweber et al. (2008).

The D-S method is based on two main assumptions: (1) During the observation period, the IMF intensity is stable compared to the other three vector components, indicating Alfvénic perturbations (Ness et al. 1964) and (2) the offset remains constant throughout that period. Hence, this method is suitable for solar wind with Alfvénic characteristics, where magnetic field variations primarily manifest in their direction rather than in their magnitude. Accordingly, the D-S method proposes to determine the zero-offset as the value that minimizes the variance of the squares of the total magnetic field.

The D-S method can be applied over intervals where the offset may be assumed to be constant. The time interval should be long enough to include averaged background magnetic field variations, while short enough to configure the entire period with a sufficient number of intervals. Therefore, we determine constant zero-offset within the intervals characterized by the Alfvén wave properties over the entire period.

The KMAG magnetic field measurement, M, can be expressed as the sum of the actual background magnetic field, B, spacecraft magnetic field, N, and sensor offset, O. To determine the sensor offset, the spacecraft-generated field has to be considered; but this will be done during higher-level data processing. Instead, we can select the dataset for which the spacecraft magnetic field (N) strength is much smaller than other two field (M and O) strengths. This condition can be confirmed by the same observation of all three magnetometers.

In this case, the measured magnetic field M is represented as the sum of the actual background magnetic field B and the sensor offset O:

M = B + O .

The D-S method calculates a constant offset O that minimizes the variance of B2. Since B = MO, the variance of Bi2 within a single interval is defined by

V i = [ ( M i O i ) ( M i O i ) ( M i O i ) ( M i O i ) ] 2 ,

where the parameter within angle brackets represents average value and the parameter with subscript i represents the time series values over the ith interval. Thus, the condition for minimizing the variance of B2 can be formulated as an equation in terms of Oi:

O i ( C i ) = 1 2 ( M i 2 M i M i 2 M i ) ,

where Ci denotes the covariance matrix given by

C i = ( M x 2 M x 2 M x M y M x M y M x M z M x M z M x M y M x M y M y 2 M y 2 M y M z M y M z M x M z M x M z M y M z M y M z M z 2 M z 2 )

Eq. (3) provides an expression for obtaining a constant offset considering variance only for a single short interval. To account for variance across the entire block comprising multiple intervals, we have

i C i · 0 = i 1 2 ( M i 2 M i M i 2 M i ) .

To determine offsets that change by a maximum of 10 nT/month, we apply Eq. (5) over every 8 hr block. Additionally, each interval comprising the block is set at a 10 min duration. Consequently, we apply the D-S method to each of the three fluxgate sensors, deriving three offsets per sensor. After the offset subtraction, we obtain zero-level measurements.

4.2 Coordinate Transformation

The subsequent coordinate transformation involves converting magnetic field measurements from KMAG’s sensor coordinate system to a spacecraft coordinate system, and then to Earth-centered and Moon-centered coordinate systems such as geocentric solar ecliptic (GSE), geocentric solar magnetospheric (GSM), selenocentric solar ecliptic (SSE), and selenographic (SEL) coordinates.

The KMAG sensor coordinate systems and the spacecraft coordinate systems are presented in Fig. 11. The origin of the spacecraft bus reference coordinate system is located at the center point of the propulsion rails and the top of the spacecraft adapter. The +X axis extends from the origin in the direction opposite to the propulsion module installed on the spacecraft. The +Z axis, also originating from the origin, runs parallel to the optical cameras and payload devices. The +Y axis completes the right-handed Cartesian coordinate system with the +X and +Z axes. KMAG’s sensor coordinate system is oriented with a 90° rotation (with respect to the spacecraft frame) about the Z-axis, followed by a negative rotation corresponding to the boom deployment angles about the X-axis. The transformation equation from magnetic field measurements in sensor coordinates BMAG to spacecraft coordinates BS/C has been designed and validated by ground testing. This equation includes the boom deployment angle φ.

Fig. 11. Coordinate systems of the spacecraft and KMAG sensors. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer.
Download Original Figure
( B x B y B z ) S / C = ( 0 1 0 1 0 0 0 0 1 ) · ( 1 0 0 0 cos ( φ ) s i n ( φ ) 0 s i n ( φ ) cos ( φ ) ) · ( B x B y B z ) M A G

The boom deployment angle of the KMAG is set at 135°. During the operation, we receive the boom deployment angle as 1 Hz data to check for errors in the boom deployment angle. To obtain the equation of conversion to the spacecraft coordinates, we substitute this angle into Eq. (6).

Subsequently, using SPICE kernels such as the camera-matrix kernel (ck), which contains spacecraft attitude information, and the frame kernel (fk), which contains frame information, we obtain Euler angles to transform the spacecraft coordinates into various meaningful coordinate systems. To finally convert the spacecraft frame to the desired frame, we need fks for the target frame. We either use an existing fk or define a new fk. The newly defined fk is verified with Kaguya data in various coordinates. In conclusion, we obtain the magnetic field measurements in the desired frame by substituting Euler angles into the coordinate transformation equation. For instance, the ZYX Euler angles α, β, and γ that are used in the coordinate transformation equations to convert magnetic field observations from spacecraft coordinates BS/C to GSE coordinates BGSE are applied as follows.

( B x B y B z ) G S E = ( cos ( α ) s i n ( α ) 0 sin ( α ) cos ( α ) 0 0 0 1 ) · ( cos ( β ) 0 s i n ( β ) 0 1 0 s i n ( β ) 0 cos ( β ) ) · ( 1 0 0 0 cos ( γ ) s i n ( γ ) 0 sin ( γ ) cos ( γ ) ) · ( B x B y B z ) S / C

As a result, we generate observation data in GSE, GSM, SSE and SEL coordinate systems. We also confirmed the reliabilities of the coordinate transformation equation by applying it to Kaguya data.


5.1 Interplanetary Magnetic Field (IMF) Observations in Trans-Lunar Cruise (TLC) Phase

The data accumulated during 4.5 months in the TLC phase were suitable to carry out calibration tests because there were relatively few spacecraft maneuvers. In the case of general KMAG data processing, we first convert the binary TM file delivered through KPDS sFTP into ASCII and scientific units to confirm the data status and check errors. The monitoring data are stable and there are a few missing data points (1.52% for about a year). Then, the ground calibration factors like an initial sensor offset correction and orthogonality corrections determined on the ground are applied.

Afterward, the removal of the spacecraft fields becomes the priority. The noise field is removed by high-pass filtering, but a more precise removal method is currently under investigation. The next step is to determine the zero offset for each sensor.

As previously mentioned, we apply the D-S method to 8 hr blocks consisting of 10 min intervals, resulting in three offsets per day for each sensor axis. The daily averaged offsets over the cruise phase are summarized in Table 3. The offsets for the representative dates with consistent spacecraft orientation remain stable within a range of 0.1–0.7 nT. The offsets exhibit an average variation of ± 2 nT for a month, with a maximum variation of ± 5 nT noted in August (shortly after the launch of KPLO). The ± 5 nT variations can be attributed to the initial phase of stabilization operations. Referring to Table 3 and Fig. 12, the offsets have changed by a maximum of approximately 5 nT over 4.5 months. Comparing this to the offset changes of Kaguya at 2 nT/5 month and THEMIS at 0.2 nT/6 month (Auster et al. 2008), the KMAG offsets display relatively larger variations. This suggests a substantial amount of spacecraft magnetic interference on the KMAG sensors, and the main reason is presumed to be the short boom length. Following the offset calibration, a coordinate transformation is performed using SPICE kernels and toolkits.

Table 3. KMAG daily offset (μ ± σ ) in the TLC phase
Date Ox (nT) Oy (nT) Oz (nT)
2022.08.08 77.23 ± 0.32 –8.33 ± 0.66 –3.83 ± 0.06
2022.09.06 78.57 ± 0.06 –11.23 ± 0.12 –0.73 ± 0.12
2022.10.16 77.40 ± 0.17 –10.73 ± 0.38 –1.67 ± 0.06
2022.11.11 73.13 ± 0.06 –10.03 ± 0.21 –0.97 ± 0.06
2022.12.08 75.37 ± 0.10 –7.23 ± 0.19 –0.53 ± 0.29

KMAG, Korea Pathfinder Lunar Orbiter Magnetometer; TLC, trans-lunar cruise.

Download Excel Table
Fig. 12. KMAG monthly averaged offsets in the TLC phase. Error bars denote standard deviation. The offsets have a standard deviation of ±2 nT on average over one month and vary by about 5 nT over 4.5 months. KMAG, Korea Pathfinder Lunar Orbiter Magnetometer; TLC, trans-lunar cruise.
Download Original Figure

We compare IMF observations from the KMAG instrument onboard KPLO with the magnetic field data from DSCOVR, monitoring the solar wind, when KPLO was in the solar wind. This comparison is conducted under an assumption of a constant time delay between both satellites. Fig. 13 shows the KPLO (blue) and DSCOVR (red) trajectories in the GSE equatorial plane for the interval of 00:00-24:00 UT on September 6, 2022. During this 24 hr interval, DSCOVR moved from GSE (x, y, z) = (225.0, –23.5, 6.4) to (225.1, 29.1, 1.0) RE, and KPLO moved from GSE (x, y, z) = (220.6, –19.7, 5.7) to (222.0, –23.5, 6.4) RE. The KPLO’s overall BLT trajectory is indicated by the dashed-dotted line and the orbit of the Moon located at a radial distance of 57.2 RE is presented by the dashed circle. The solid curves represent the model bow shock by Fairfield (1971) and model magnetopause by Shue et al. (1998), respectively. Fig. 14 displays IMF variations over that period measured by KPLO (black) and DSCOVR (gray). The IMF vector quantities are in GSE coordinates. In timeseries plots, the DSCOVR magnetic field data are added to 15 nT and shifted to the right by 5 min for a comparison with KPLO data. The constant time delay of 5 min corresponds to a convection time of the solar wind between DSCOVR and KPLO. The 15 nT offset is marked by a horizontal dashed line in each component and total field intensity. The IMF variations observed at KPLO and DSCOVR exhibit very similar signatures with a mean difference less than 0.2 nT and a standard deviation difference less than 0.1 nT. After applying the D-S method to remove KMAG offset values, it has been confirmed that the KMAG instrument onboard KPLO accurately measures the solar wind magnetic field.

Fig. 13. Trajectories of KPLO (blue) and DSCOVR (red) in GSE coordinate system on September 6, 2022, 00:00–24:00 UT. The last position of each object is indicated by a filled circle. KPLO, Korea Pathfinder Lunar Orbiter; GSE, geocentric solar ecliptic.
Download Original Figure
Fig. 14. Time series of magnetic fields observed by KPLO (black) and DSCOVR (gray) in GSE coordinate system on September 6, 2022. We use 1 sec data for both KPLO and DSCOVR. To compare with KPLO data, a 15 nT offset and a 5 min time delay are applied to the DSCOVR data. (a–c) Each magnetic field component exhibits significant variances, while (d) the total magnetic field remains stable compared to the other three components, indicating Alfvénic perturbations. KPLO, Korea Pathfinder Lunar Orbiter; GSE, geocentric solar ecliptic.
Download Original Figure
5.2 Magnetic Field Observations in Nominal Mission Phase

Fortunately, the KPLO is operating accompanied by the ARTEMIS P1 and P2 spacecraft, which have elliptical trajectories around the Moon and observe the magnetic field. The configuration of these three spacecraft provides a new insight into the lunar magnetic field investigation; moreover, P1 and P2 are the best references for the KMAG observation.

Fig. 15 shows the KPLO (blue) and ARTEMIS P1 (red) trajectories in the GSE equatorial plane and SSE equatorial plane during the 12:00–24:00 UT period on May 4, 2022, respectively, indicating that they are located within the magnetotail. Fig. 16(a)16(d) represents magnetic field data from KMAG observed on May 4, 2023, along with ARTEMIS P1 data including a 15 nT offset. The positions of KPLO, which changes as it follows the Moon’s dawn-dusk orbital plane, are represented in the SEL coordinate system in Fig. 16(e) and 16(f). The spacecraft orientation changes constantly because of nadir pointing. Despite the complicated changes in the attitude of the spacecraft, KMAG at 100 km altitude exhibits similar features within 1–2 nT with respect to ARTEMIS P1. These results show the reliabilities of the KMAG in-flight calibration and KMAG data.

Fig. 15. Position of KPLO (blue) and ARTEMIS P1 (red) in (a) GSE coordinate system and (b) SSE coordinate coordinate system on May 4, 2023, 12:00–24:00 UT. The last position of each object is indicated by a filled circle. KPLO, Korea Pathfinder Lunar Orbiter; GSE, geocentric solar ecliptic; SSE, selenocentric solar ecliptic.
Download Original Figure
Fig. 16. Time series data of KPLO and ARTEMIS P1 in SEL coordinate system on May 4, 2023, 12:00–24:00 UT. We use 1 sec KPLO data and 4 sec ARTEMIS P1 data. To compare with KPLO data, a 15 nT offset is applied to the ARTEMIS P1 data. (a–d) Comparison of magnetic field observations by KPLO (black) and ARTEMIS P1 (gray). (e) The longitude (red) and latitude (blue) of KPLO in SEL coordinate system. (f) The altitude of KPLO on the Moon. KPLO, Korea Pathfinder Lunar Orbiter; SEL, selenographic.
Download Original Figure

During lunar orbits, KPLO covers all orbital planes, from the dawn-dusk plane to the noon-midnight plane, every six months, maintaining a polar orbit. In particular, near the noon-midnight (LT 12–24 hr) orbital plane, the KPLO spends considerable time in the lunar wake region where the KPLO is within the solar wind. This is an area where the plasma effect of the solar wind is relatively small, and it helps in studying the magnetic field of the Moon. KPLO is located within a solar wind for much of the time, but when it is within the magnetotail (a short period of time per month periodically), a quiet background magnetic field environment is also good for measuring the lunar magnetic field.

With diverse observation orbits, KMAG successfully observed lunar magnetic anomalies as well as solar wind events. Fig. 17 is a trajectory of an orbit in which KPLO passed over lunar magnetic anomalies, including Abel, which is located at 34.6 °S, 85.8 °E with a radius of 137.4 km, on May 5, 2023. KPLO is positioned in the magnetotail, one day apart from Fig. 15. The results obtained by KMAG between 12:23 UT and 13:00 UT are presented in Fig. 18, organized by latitude. The observations are compared with SVM’s 100 km observations based on previous observations by LP and Kaguya. KMAG observed about 0.5 nT of magnetic strength change in magnetic anomalies, similar to SVM.

Fig. 17. Orbit trajectory of KPLO when passing through lunar magnetic anomalies including Abel in SEL coordinates on May 5, 2023, 12:23–13:00 UT. The background magnetic field map shows the SVM’s total magnetic field intensity at an altitude of 100 km. KPLO, Korea Pathfinder Lunar Orbiter; SEL, selenographic; SVM, surface vector mapping.
Download Original Figure
Fig. 18. Comparison of KPLO observed magnetic field and SVM in SEL coordinates when passing over lunar magnetic anomalies including Abel on May 5, 2023, 12:23–13:00 UT. The first and second rows of x-ticks represent latitude and longitude, respectively. (a–c) Each magnetic field component is presented on a 3 nT scale, while (d) the total magnetic field is displayed on a 2 nT scale to provide a more detailed view of magnetic field intensity variations. KPLO, Korea Pathfinder Lunar Orbiter; SVM, surface vector mapping; SEL, selenographic.
Download Original Figure


KMAG operation has been stable for over one year. Using accumulated observations during the cruise orbit and lunar mission orbits, validation of KMAG data reliability has been achieved. Moreover, the data-processing pipeline and data calibration methods are established and updated more precisely. Furthermore, KMAG has successfully observed lunar magnetic fields such as lunar magnetic anomalies and lunar-induced magnetic fields, and solar wind events such as magnetic storms and magnetic clouds, in various environments, which we have been studying further.

However, the spacecraft’s magnetic interferences have not been fully eliminated because of the varied and complicated sources of spacecraft magnetic fields. In particular, the continuous maneuvering of the spacecraft during lunar orbits results in observations that include a significant amount of spacecraft-generated noise. Removing this is crucial for generating final research-grade data. Currently, two causes of the spacecraft noise field that have been investigated include the turning on and off of the valve-drive electronics (VDE) and the battery-charging current, with the possibility of more sources. The step-like and spike-like noise fields caused by the VDE can be removed using the PiCoG method (Magnes et al. 2020). However, the irregular and diverse noise field caused by battery discharging current is being researched, as is the removal of other factors. When these spacecraft field removal techniques are completed, more accurate offset determination and higher-quality data will be obtained.

Furthermore, complex thermal correction coefficients must be studied further on a case-by-case basis. Following a comprehensive temperature calibration, a revised version of the public data will be generated.

The KPLO mission has been extended by an additional three years beyond the initially planned mission duration, now reaching until 2025. During this extended mission period, if lunar magnetic field observations are possible in low orbits, it is expected that more meaningful and detailed studies of the lunar magnetic field, such as low-altitude lunar magnetic field mapping, can be conducted.


This work was supported by (NRF-2020M1A3B7109194) and by the KPLO payload instrument operation program (KASI). The authors wish to thank all of the KPLO team members for their help in the observation. Magnetic field data from NOAA’s DSCOVR spacecraft used as a space weather reference were produced by CDAWeb/NASA. ARTEMIS P1 FGM data was supported by NASA contract NAS5-0299. We acknowledge K. H. Glassmeier, U. Auster, and W. Baumjohann for the use of FGM data provided under the lead of the Technical University of Braunschweig and with financial support through the German Ministry for Economy and Technology and the German Center for Aviation and Space (DLR) under contract 50 OC 0302. IGB acknowledges support from the NASA KPLO Participating Scientist Program. The magnetic field of surface vector mapping (SVM) data are available at



Auster HU, Glassmeier KH, Magnes W, Aydogar O, Baumjohann W, et al., The THEMIS fluxgate magnetometer, Space Sci. Rev. 141, 235-264 (2008).


Baek SM, Kim KH, Garrick-Bethell I, Jin H, Lee HJ, et al., Detailed study of the Mare Crisium northern magnetic anomaly, J. Geophys. Res. Planets 122, 411-430 (2017).


Belcher JW, A variation of the Davis-Smith method for in-flight determination of spacecraft magnetic fields, J. Geophys. Res. 78, 6480-6490 (1973).


Belcher JW, Davis L Jr, Large-amplitude Alfvén waves in the interplanetary medium, 2, J. Geophys. Res. 76, 3534-3563 (1971).


Binder AB, Lunar Prospector: overview, Science 281, 1475-1476 (1998).


Davis L Jr, Smith EJ, The in-flight determination of spacecraft magnetic field zeros, Eos Trans. Am Geophys. Union 49, 257 (1968).


Dyal P, Parkin CW, The Apollo 12 magnetometer experiment: internal lunar properties from transient and steady magnetic field measurements, Proceedings of the 2nd Lunar Science Conference, Houston, TX, USA, 11-14 Jan 1971.


Dyal P, Parkin CW, Daily WD, Structure of the lunar interior from magnetic field measurements, Proceedings of the 7th Lunar Science Conference, Houston, TX, USA 15-19 Mar 1976.


Fairfield DH, Average and unusual locations of the Earth’s magnetopause and bow shock, J. Geophys. Res. 76, 6700-6716 (1971).


Farrell WM, Thompson RF, Lepping RP, Byrnes JB, A method of calibrating magnetometers on a spinning spacecraft, IEEE Trans. Magn. 31, 966-972 (1995).


Garrick-Bethell I, Kelley MR, Reiner gamma: a magnetized elliptical disk on the Moon, Geophys. Res. Lett. 46, 5065-5074 (2019).


Haviland HF, Poppe AR, Fatemi S, Delory GT, de Pater I, Time-dependent hybrid plasma simulations of lunar electromagnetic induction in the solar wind, Geophys. Res. Lett. 46, 4151-4160 (2019).


Hood LL, Coleman PJ Jr, Wilhelms DE, The Moon: sources of the crustal magnetic anomalies, Science 204, 53-57 (1979).


Hood LL, Mitchell DL, Lin RP, Acuna MH, Binder AB, Initial measurements of the lunar induced magnetic dipole moment using Lunar Prospector magnetometer data, Geophys. Res. Lett. 26, 2327-2330 (1999).


Kato M, Sasaki S, Takizawa Y, The Kaguya Project Team, The Kaguya mission overview, Space Sci. Rev. 154, 3-19 (2010).


Kepko EL, Khurana KK, Kivelson MG, Elphic RC, Russell CT, Accurate determination of magnetic field gradients from four point vector measurements. I. Use of natural constraints on vector data obtained from a single spinning spacecraft, IEEE Trans. Magn. 32, 377-385 (1996).


Kim JH, Choi YJ, Kim BY, Development of KPLO science data archive for public release, Proceedings of the 3rd Planetary Data Workshop and the Planetary Geologic Mappers Annual Meeting, Flagstaff, AZ, USA, 12-15 Jun 2017.


Lee H, Jin H, Jeong B, Lee S, Lee S, et al., KMAG: KPLO magnetometer payload, Publ. Astron. Soc. Pac. 133, 034506 (2021).


Leinweber HK, Russell CT, Torkar K, Zhang TL, Angelopoulos V, An advanced approach to finding magnetometer zero levels in the interplanetary magnetic field, Meas. Sci. Technol. 19, 055104 (2008).


Maggiolo R, Hamrin M, De Keyser J, Pitkänen T, Cessateur G, et al., The delayed time response of geomagnetic activity to the solar wind, J. Geophys. Res. Space Phys. 122, 11109-11127 (2017).


Magnes W, Hillenmaier O, Auster HU, Brown P, Kraft S, et al., Space weather magnetometer aboard GEO-KOMPSAT-2A, Space Sci. Rev. 216, 119 (2020).


Ness NF, Scearce CS, Seek JB, Initial results of the Imp 1 magnetic field experiment, J. Geophys. Res. 69, 3531-3569 (1964).


Park HH, Jin H, Kim TY, Kim KH, Lee HJ, et al., Analysis of the KPLO magnetic cleanliness for the KMAG instrument, Adv. Space Res. 69, 1198-1204 (2022).


Richmond NC, Hood LL, A preliminary global map of the vector lunar crustal magnetic field based on Lunar Prospector magnetometer data, J. Geophys. Res. Planets 113, E02010 (2008).


Russell CT, Coleman PJ Jr, Goldstein BE, Measurements of the lunar induced magnetic moment in the geomagnetic tail: evidence for a lunar core?, Proceedings of the 31st Lunar and Planetary Science Conference, Houston, TX, USA, 16-20 Mar 1981.


Shue JH, Song P, Russell CT, Steinberg JT, Chao JK, et al., Magnetopause location under extreme solar wind conditions, J. Geophys. Res. Space Phys. 103, 17691-17700 (1998).


Son DR, Construction of feed-back type flux-gate magnetometer, J. Korean Magn. Soc. 22, 45-48 (2012).


Song YJ, Bae J, Hong S, Bang J, Pohlkamp KM, et al., KARI and NASA JSC collaborative endeavors for joint Korea Pathfinder Lunar Orbiter flight dynamics operations: architecture, challenges, successes, and lessons learned, Aerospace 10, 664 (2023).


Song YJ, Kim YR, Bae J, Park J, Hong S, et al., Overview of the flight dynamics subsystem for Korea Pathfinder Lunar Orbiter mission, Aerospace 8, 222 (2021).


Takahashi F, Shimizu H, Matsushima M, Shibuya H, Matsuoka A, et al., In-orbit calibration of the lunar magnetometer onboard SELENE (KAGUYA), Earth Planets Space 61, 1269-1274 (2009).


Tsunakawa H, Shibuya H, Takahashi F, Shimizu H, Matsushima M, et al., Lunar magnetic field observation and initial global mapping of lunar magnetic anomalies by MAP-LMAG onboard SELENE (Kaguya), Space Sci. Rev. 154, 219-251 (2010).