Mingantu Spectral Radioheliograph for Solar and Space Weather Studies

The Chinese Spectral Radioheliograph (CSRH) covering 400 MHz-15 GHz frequency range was constructed during 2009–2016 in Mingantu Observing Station, National Astronomical Observatories, Chinese Academy of Sciences at Zhengxiangbaiqi, Inner Mongolia of China. The CSRH is renamed as MingantU SpEctral Radioheliograph (MUSER) after its accomplishment. Currently, MUSER consists of two arrays spreading over three spiral-shaped arms. The maximum baseline length is ∼3 km in both east-west and north-south directions. The MUSER array configuration is optimized to meet the needs of observing the full-disk Sun over ultrawide wavebands with images of high temporal, spatial and spectral resolutions and high dynamic range. The low frequency array, called MUSER-I, covers 400 MHz-2.0 GHz with 40 antennas of 4.5-m-diameter each and the high frequency array, called MUSER-II, covers 2–15 GHz with 60 antennas of 2-m-diameter each. The MUSER-I can obtain full-disk solar radio images in 64 frequency channels with a time cadence of 25 ms and a spatial resolution of 51.6″ to 10.3″ (corresponding to the frequency range 400 MHz to 2 GHz), whereas the MUSER-II can obtain full-disk solar images in 520 channels with a time cadence of 206.25 ms and a spatial resolution of 10.3″ to 1.3 (corresponding to the frequency range 2 to 15 GHz). A dynamic range of 25 dB can be obtained with snapshot images produced with the MUSER. An extension of MUSER in the further lower frequency range covering 30–400 MHz with an array of 224 logarithm-periodic dipole antennas (LPDAs) has been approved and will be completed during the next 4 years. The MUSER, as a dedicated solar instrument, has the following advantages providing simultaneous images over a wide frequency range with a unique high temporal-spatial-spectral resolutions; high-performing ultrawide-band dual-polarization feeds for wide-band signal collection; advanced high data-rate, large-scale digital correlation receiver for multiple-frequency and faster snapshot observations; and applications of new technologies such as using optical fiber to obtain remote antenna and wide-band analog signal transmission. The MUSER thus provides a unique opportunity to measure solar magnetic fields and trace dynamic evolution of energetic electrons in several radio frequencies, which, in turn, will help to have better understandings of the origin of various solar activities and the basic drivers of space weather.


INTRODUCTION
Solar radio bursts are rich with information as they are associated with different types solar eruptions, such as solar flares, coronal mass ejections (CMEs), and various thermal and nonthermal processes. The solar eruptions are believed to arise due to the sudden energy release process because of topological reorganization of solar magnetic field or magnetic reconnection (Benz, 2009). Radio bursts are hence prompt indicators of those solar activities. X-ray observations show that non-thermal particles are highly related to the energy release process and demonstrate a significant amount of total flare energy during initial phase of solar flares (Lin (2008)). Normally, the radio bursts are observed in an ultrawide frequency range, starting from decimetric wavelength range, extending down to a few tens of KHz and up through several GHz even up to mm-wavelengths in a short time scale (e.g., Benz, 2009). Especially, the radio observations covering centimeter and decimetric wavelengths are important as they can reveal key informations about the energy release, and particle acceleration and transportation (Bastian et al., 1998;Gary and Keller 2004;Aschwanden 2005;Pick and Vilmer 2008;Chernov et al., 2014).
The radio burst emissions were commonly generated by different mechanisms. If the plasma mechanism were assumed as the generation mechanism, one can infer from statistical study of radio dynamic spectra during different flare events that the electrons are accelerated from a region with electron density of n acc e ∼ 3 × 10 9 to 10 11 cm −3 . The corresponding plasma frequency is about ] p ∼ 500 MHz to 3.4 GHz, where electron beams are accelerated and propagated in either upward (type III bursts, down to ∼10 s KHz) and/or downward (Reverse Slope type III bursts, up to microwave) directions (Aschwanden and Benz 1997;Tan et al., 2016). However, the present solar radio imaging observations are only available at a few discrete frequencies at 70 MHz for Gauribidanur Radioheliograph (Ramesh et al., 1998), in the range 150-450 MHz for Nancay Radioheliograph (NRH, Radioheliograph Pick and Vilmer, 2008;Kerdraon and Delouis 1997), and at 17/34 GHz for Nobeyama Radioheliograph (NoRH, Nakajima et al., 1994). Those closed facilities in the past include Culgoora radioheliograph at three metric-wave frequencies (Sheridan et al., 1973), the Clark Lake Array in the range 1.5-12.5 MHz (Erickson et al., 1984), etc. Some radioheliographs were upgraded in the recent years, eg. Siberian Solar Radio Telescope at 5.7 GHz (SSRT, Grechnev et al., 2003) and the Expanded Owens Valley Solar Array in a frequency range of 1-18 GHz (EOVSA). There are also some other radio telescopes designed to investigate astronomical objects and the sun is also one of their purposes, eg. Very Large Array (Napier et al., 1983), the Giant Metrewave Radio Telescope (GMRT, 150-1450 MHz) in India, the Low-Frequency Array (LOFAR, in Europe, the Murchison Widefield Array (Tingay et al., 2013) in Australia, and the Long Wavelength Array (LWA, 10-88 MHz). Though the succession of these solar-dedicated radio imaging instruments has greatly advanced the science of the Sun as described in the related reviews (e.g., Bastian et al., 1998;Pick and Vilmer 2008;Benz 2009;Chernov 2011), a key parameter has been lacking, i.e., high quality spectral imaging of the Sun over an ultra-wide band frequency range, covering from the primary energy release site and to the beyond.
Many studies on the associations between hard X-ray (HXR) and radio emissions demand this need. For example, Vilmer et al., 2002 studied and compared locations of RHESSI HXR and NRH radio sources observed for a flare event on 20 February 2002. Since the direct link between HXR and radio emissions was not available from observations, it was hard to understand the details of magnetic reconnection and the energy conversion processes, including their temporal and spatial evolutions. Therefore, it is crucial to image radio emissions in wide centimetric and decimetric wavelengths in order to cover the whole acceleration sites and primary propagation regions of nonthermal electrons which are responsible for the HXR emission. Trottet et al. (2006) carried out a detailed analysis of the HXR images observed by the Yohkoh and the radio images observed by the NRH in a flare on 5 November 1998. They found the evidence that the HXR and radio-emitting electrons were produced by the same accelerator. However, the available observations and analysis did not allow them to make a definite connection between the HXR source region and the radio source sources in the middle corona. Therefore, the spectral imaging observations in frequency range higher than 400 MHz are really important in this regard. As for the famous flare event on 13 December 2006, many radio fine structures were registered by the Chinese Spectral Broadband Radio Spectrometer (SBRS) in the microwave range. Deduced from the microwave Zebra Patterns (ZPs), it was found that the scaleheight ratio between plasma density and magnetic field decreased by a factor of about 2 before and after the flare maximum (Yan et al., 2007). A further statistical investigation from 74 ZPs confirmed the above result (Yu et al., 2012). However, due to lack of observations of the relevant coronal magnetic field structures, it is difficult to present an exact interpretation of the above results. The radio imaging-spectroscopy over a wide frequency range provide such observations. The spectral observations of small-scale microwave bursts (SMBs) in solar flares provides the evidence of small-scale elementary energy releasing activities and electron accelerations in the flaring source regions (Tan, 2013), however, in order to understand the nature of SMBs, we need imaging-spectroscopy in the related frequency range, which may provide the real position of the source region, the relevant magnetic field, and the relationship with the physical mechanisms. The microwave type III pairs (Aschwanden and Benz 1997;Tan et al., 2016) indicate that the frequency where magnetic reconnection and particle acceleration taking place occurred around 0.3-3.4 GHz, but we need imaging observations at the corresponding frequencies to deduce the magnetic fields and the relevant topological structures. Thus, the solar radio imaging-spectroscopy can provide the following crucial information for the solar bursts: 1) positions, 2) topological structures, 3) coronal magnetic fields, and 4) spatio-temporal evolutions. With these informations, we can deduce the primary energy release, particle acceleration, and the mass-energy transportation, furthermore, predicting the occurrence of the solar bursts and their corresponding influences on space weather.
Several general purpose radio telescopes can also be used to observe the Sun, such as the VLA (Napier et al., 1983), the Low-Frequency Array (van Haarlem et al., 2013) and the (Swarup, 1991). The upgraded VLA has been used for solar studies (e.g., Chen et al., 2013;Chen et al., 2015). However, as they were primarily built for non-solar radio astronomical objectives, only about 2% time is available for solar observations, and their field of view is normally too small to cover the whole solar disk. In fact, the solar-dedicated instrument of imaging spectroscopy should have high temporal, spatial, and spectral resolutions and at least with field of view more than 32 arc-minute simultaneously (Bastian et al., 1998;Hudson and Vilmer, 2007;Pick and Vilmer, 2008). The Chinese Spectral Radioheliograph (Yan et al., 2004) and the Frequency-Agile Solar Radiotelescope (FASR, Bastian, 2003;Gary, 2003) proposed to realize this goal. The FASR is still not granted for construction yet and the Expanded Owens Valley Solar Array (Gary et al., 2018) has been developed as a pathfinder for the FASR (Nita et al., 2016). The microwave spectral imaging observations of the well known X8.2-class limb flare on 2017 September 10 by EOVSA have been extensively studied to indicate the nonthermal emissions by flareand shock-accelerated electrons (Gary et al., 2018;Fleishman, et al., 2020;Karlický et al., 2020), or the detection of nonthermal emission at conjugate flux rope footpoints showing the solid evidence of particle transport along the erupting magnetic flux rope during the early impulsive phase .
The Chinese solar physics community had planned to build a radioheliograph since 1960s. Some pre-studies were carried out on proposals for radioheliograph in either centimeter-band (Hu et al., 1984) or millimeter-band (Fu et al., 1997), but none of these had been implemented. Following these lines, it was suggested to build a Chinese Spectral radioheliograph in the decimetric to centimeter wavelength range (Yan et al., 2004). It was later recommended as one of the two major ground-based facilities by Chinese solar physics community in 2006. A 2-element interferometer prototype was built and tested for tackling the key technologies in -2005(Yan et al., 2009). The Chinese Spectral Radioheliograph was officially supported in 2009 as a National Major Scientific Research Facility Program of China. The site survey was pursued at Mingantu town in Inner Mongolia of China. The radio quiet zone protection of 10 km radius centered at Mingantu Observing Station has been setup since 2008. The construction was fulfilled during 2009-2016, and the instrument was renamed as MingantU SpEctral Radioheliograph (MUSER) after its accomplishment, as seen in Figure 1. The brief description and progress of MUSER project while it was under construction as well as a few initial results from MUSER were reported by Yan et al. (2013); Yan et al. (2016). The Meridian-II project among the Major National Infrastructure Projects for Science and Technology under "13th 5-year plan" program (2016-2020) has been approved and a Solar and Interplanetary subsystem as a new part in the Meridian-II project will include a logarithm-periodic dipole antenna (LPDA) array at metric and decametric wave range to be built in the Mingantu Observing Station (Blanc et al., 2020).
The key technical issues consist of implementing high cadence imaging with an aperture synthesis system at about two order higher multiple frequencies than the presently-available radioheliographs over an ultrawide frequency band, and the data processing for such a large data volume from MUSER. We introduce the MUSER system description in § 1. The calibration and synthesis imaging are described in § 2. The observational results are presented in § 3 and the future plans are briefly described in § 4. Finally, we summarize our conclusions in § 5.

DESCRIPTION OF MUSER
The MUSER is a solar-dedicated radio interferometric array with high temporal, spatial and spectral resolutions, and it can simultaneously perform the spectral and imaging observations of the full Sun in a wide frequency range (Yan et al., 2004;Yan et al., 2009). The main characteristics and performance are listed in Table 1. The MUSER currently covers the centimetric to decimetric wave range, including MUSER-I operating in the frequency range of 400 MHz-2 GHz and MUSER-II in 2-15 GHz. The MUSER-I array contains 40 antennas with 4.5 m diameter each, and the MUSER-II array contains 60 antennas with 2 m diameter each. All the 100 antennas are located on three logspiral arms with the maximum baseline length of about 3 km in both north-south and east-west directions. The MUSER antennas, numbered as IA#, IB#, IC# from 1 to 13 for MUSER-I array and HA#, HB#, HC# from 1 to 20 for MUSER-II array in local coordinates with the central antenna IA0 as the reference point, are schematically shown in black dots in Figure 2. IA0's location is E 115°15'1.8", N 42°12′42.6″, with an altitude of 1365 m. The central part of the MUSER arrays is shown in Figure 1. Figure 3 shows the schematic system block diagram for MUSER.
The signal processing for MUSER-I (400 MHz-2 GHz) and MUSER-II (2-15 GHz) is almost same. Additionally two 20-m antennas (400 MHz-1 GHz) were established for interferometry experiment in 2011 at the same site. They are proposed to work as a part of MUSER-I array for calibrations. Figure 3 shows that the solar radio signal (400 MHz-2 GHz, 2-15 GHz) as received by the MUSER-I and MUSER-II antennas with broadband feeds, front-end LNAs and optic transmitters. The signal is then transmitted through optic fibers to indoor analog receivers with an Intermediate frequency (IF) output of 400 MHz bandwidth (50-450 MHz), which covers the whole bandwidth by sweeping four times for MUSER-I and 33 times for MUSER-II. It takes 25 ms to cover the 400 MHz-2 GHz bandwidth of MUSER-I by sweeping. For MUSER-II, the sweeping time is ∼200 ms. The observing mode implemented in MUSER allows us to observe at a given frequency in full 25 ms without band or polarization switching. So it provides the flexibility to investigate useful information by setting the working mode of the instrument. The 400-MHz IF signal is received by a digital correlation receiver. Firstly, it is sampled with 1 Gsps Analog to Digital converters (ADC), then goes into the polyphase filter bank (PFB). The PFB generates 16 complex baseband signals with ∼2-25 MHz simultaneously. All baseband signals with same frequency are computed to send out the correlation data. The time delay compensation and fringe stopping are implemented in the digital correlation receiver (Liu et al., 2019). Both MUSER-II and MUSER-I has the same architecture in their digital correlation receivers. The whole correlation procedure is controlled by a monitoring subsystem. Among the various issues of aperture synthesis techniques for MUSER (Taylor et al., 1999;Thompson et al., 2017), one of the key problems is to develop a high performance ultra-wide band feed for the reflector antennas. Such feed should have wide impedance bandwidth, low profile, symmetrical radiation patterns, and fixed phase center over the whole frequency band. So far, the eleven feed (Olsson et al., 2006) can be used for observing the slow-evolving astrophysical objects since its original version only has a maximum return loss of about −5 dB (or Voltage Standing Wave Ratio-VSWR is about 3.53). However, the solar radio bursts and accordingly polarizations always vary quickly. The isolation is thus vital for correct observations of the polarizations, and the return loss or VSWR for solar observations should be further reduced in order to meet the requirement for observing fast-changing solar radio burst signals. We have successfully developed the ultra-wide band feeds for both MUSER-I and MUSER-II with the VSWR less than 1.5 over most of the frequency range, wide (∼ 133%) impedance bandwidth and good radiation characteristics (Li et al., 2015a;Li et al., 2015b). Figure 4 shows the polarization degree measurement of one antenna element (IC8) in frequency of 400 MHz-2 GHz as an example. It shows that in almost 90% frequency band the polarization degree, or (P L − P R )/(P L + P R ) (where P L and P R represent the intensity of left and right polarization) is less than 5%, which represents good isolation performance. Similar performance have been measured for other antenna elements. Time delays among the different antennas of MUSER array have been measured for calibration. The measured RMS errors of time delay compensations for MUSER-I are normally <1 ns. The result is very robust as similar results obtained in multiple measurements that were carried out at the interval of 1 year later . The robust satisfactory <1 ns RMS accuracy of time delay compensations for MUSER-II has also been obtained. During the test observations for the signals of either the satellites, the Sun, or Cygnus A, correlation fringes have been obtained successfully for all baselines (Wang et al., 2013a). For every tri-antenna composition among the MUSER-I or MUSER-II arrays, the residuals of the phase closures were measured with a value of around two degrees for both geostationary and GPS satellites. Fringe stopping has been achieved for all baselines when observing the quiet Sun. These experiments validated the system design and demonstrated the system performance. The two 20 m antennas operating in 400 MHz-1 GHz has also been incorporated into MUSER-I for calibrations. By simulating different cases for observing the quiet Sun, solar active regions and radio bursts, it is shown that the shortest baselines for MUSER-I are important for the quiet Sun image recovering (Du, et al., 2015), which are important for MUSER data analysis.

MUSER OBSERVATIONS
The MUSER array began observations in 2014, routinely observed the Sun from 2016, roughly in a time range from 01:00 to 08:00 UT. In total, the MUSER has accumulated more than 390 TB observational data. Considering a large amount of raw data, it has not been available in an online archive so far. But the data are public and interested users can get in touch with us if they are looking for some specific data. There were fewer spectral events observed by MUSER-II in the high frequency band. Wang et al. (2019) proposed a new restoration method for the data processing of MUSER-II observations at 4.2 GHz, which shows a better quality than using the deconvolution algorithm for the production of radio images. Here, more attentions were paid to a lower frequency band, which show more fruitful spectral structures. During 2014-2019, a total number of 85 solar radio burst events have been registered by the MUSER-I, as shown in Table 2. More than 60 radio burst events contain the fine structures. The MUSER data are processed using the radio astronomy software, Common Astronomy Software Applications (CASA) and the own developed Fourier routines. The data processing code has been made available at https://github.com/astroitlab/museros (Mei et al., 2018). A brief introduction of the pipeline for data processing was discussed in Mei et al. (2018) and Chen et al. (2019). Next, we will present three typical solar radio burst events observed by MUSER.

A Radio Burst Event on 11 November 2014
During 04:22-04:24UT On 11 Nov 2014, a radio burst event was recorded by MUSER-I array at 400 MHz-2 GHz. Figure 5 presents a comparison of time profiles of GOES SXR flux, the NoRP radio flux at 2 GHz and MUSER-I radio fluxes at several frequencies in the range of 400 MHz-2 GHz. It was wrongly attributed to a C-class flare (started from 04:22 UT and peaked at 04:49UT) near the disk center according to SGD event list (http://www.swpc.noaa.gov/products/goes-x-ray-flux). The Nobeyama Radio Polarimeters (NoRP) also demonstrated strong burst signals at 1, 2, and 3.75 GHz during 04:22-04:24 UT but weaker signals at higher frequencies. While NoRH (Nakajima  The radio images obtained from MUSER-I at 1.7 GHz show that the radio source is located at the east limb of the Sun, not in the solar disk area, as shown on top panels in Figure 6 at different times before, during and after the radio burst. In order to identify the correctness of the radio sources, we checked the EUV images observed by AIA/SDO (Lemen et al., 2012) during the same time period and at the same location.
The EUV images at wavelength of 94 Å and 131 Å are presented in the middle and bottom panels of Figure 6. It is found that there were indeed loop opening at 94 Å and 131 Å corresponding to hot plasmas whereas responses at other EUV wavelengths were not so obvious. A CME event with a signature starting around 04:36 UT was observed by SOHO/LASCO (Brueckner et al., 1995), which indicates the CME at that location ( Figure 7). As there was no radio burst at higher frequencies at the solar east limb, we propose that the decimetric radio burst starting at 04:22 UT in the east limb was due to the solar eruptions taking place in the backside of the Sun but it was very close to the east limb. As radio bursts occurred at higher altitude which can be observed by MUSER-I in the Earth direction. Its occurrence during the impulsive phase of the C-class flare in the solar disk was just a coincidence and there should be no physical connection to the C-class flare on the solar disk. This FIGURE 7 | The CME process recored by the SOHO/LASCO C2 at 04:29, 04:59 and 05:29UT after the radio burst. The CME and radio burst locations are indicated by the arrows, respectively.
Frontiers in Astronomy and Space Sciences | www.frontiersin.org March 2021 | Volume 8 | Article 584043 demonstrate the importance of the image spectroscopy observation of the solar radio burst (Yan et al., 2016).

A Radio Burst Event on 17 December 2014
The spectral imaging of quasi-periodic pulsations (QPPs) overlaid on a type IV microwave continuum was recorded by MUSER at frequencies of 1.2-2.0 GHz for the first time. It took place during 04:25-04:45 UT in an M8.7 flare (start at 04:25 UT, peak at 04:51 UT) in active region AR12242 on 2014 December 17.
The flare region has circular ribbons over multiple-scale loop structures as revealed by the EUV images of AIA/SDO. The flarerelated loops can be classified into three groups: small-scale lowlying loops, intermediate dome-like structure, and a group of large-scale loops.
The radio observations (shown in Figure 8) indicates that 1) the temporal profile at 2 GHz is similar to that of NoRP, but with more spikes superimposed on an increasing-to-decreaing intensity profile; 2) the intensity profile integrated from the radio images of the source region at 1.7 GHz matches well with the temporal profile at 1.7 GHz obtained from the spectrum; 3) the period of the radio QPP at 2.0 GHz is about FIGURE 8 | (A) The dynamic spectrum of quasi-periodic pulsations overlaid on type IV solar radio continuum was recorded by MUSER-I in a time range of 04: 25-04:45 UT and a frequency range of 1.2-2.0 GHz. (B) and (C) Images of the radio source with a 50% intensity level observed at multi frequencies from 1.2 GHz(red) to 2.0 GHz (black), overlaid on the AIA images of three combined wavelengths at the start (04:25 UT) and peak (04:32 UT) time of the radio burst. The radio source centroids were marked by the plus signs. (D) The radio flux curves from NoRP, the spectrum of MUSER and the intensity of radio source from MUSER imaging (green).
Frontiers in Astronomy and Space Sciences | www.frontiersin.org March 2021 | Volume 8 | Article 584043 8 121s from the wavelet analysis; 4) the size of the radio source varies from small to large and then small at a given intensity and frequency, which corresponds to the evolutions of radio emission processes; 5) the radio sources locate in the middle of the active region, over the positive magnetic field; 6) the radio sources at six frequencies line up nicely and move along the similar expanding direction of flaring loop. The detailed data processing of spectral and imaging observations of the radio source, the main features and complicated interactions among different loops, and the physical relationships from various multi-wavelength observations to the possible mechanisms of QPPs have been presented in Chen et al. (2019).

A Radio Burst Event on 22 November 2015
The solar fine radio burst of November 22, 2015 event is shown in the Figure 9. It was occurred just at the beginning of a small SXR flare (indicated by an arrow on the top panel of Figure 9).

EXTENSION OF MUSER TO 30-400 MHZ BAND
The project to extend MUSER to 30-400 MHz frequency regime with 224 LPDAs has been approved under the Meridian-II Project which is a National Science Infrastructure Project under "13th 5-years plan" program (2016-2020). Solar and Interplanetary Subsystem as a new part in Meridian-II project will include metric and decametric wave range with a LPDA array to be built in Mingantu Observing Station. The array configuration for extension of MUSER to lower frequency range of 30-400 MHz with 100 LPDAs has been optimized to achieve minimum RMS number derivation of antenna distributions in axial direction and minimum RMS number derivation of antenna distributions with respect to a Gaussian distribution in the radial direction. Some LPDAs' locations are adjusted to avoid blocking to each other during observations. As the original central part of the MUSER-I and MUSER-II arrays is populated with dense antennas, there is no sufficient space to allocate new LPDAs. So the compact central part of the new array is chosen to be close to the arm B and 31 LPDAs can be allocated there, as shown in the region outlined by blue area in Figure 2. The rest 69 LPDAs are spread over three arms with 20 LPDAs in arm A, 25 LPDAs in arm B and 24 LPDAs in arm C, with red dots indicating the locations of the newly designed 100 LPDAs in Figure 2. Within the compact area located in the red square in Figure 2, a total of 124 LPDAs will be grouped into 16 sub arrays with 1 central subarray   containing 19 LPDAs whereas rest subarray each containing 7 LPDAs. The 16 subarray will form a phased array for calibration as shown in Figure 10. The snapshot uv distribution and the synthesized beam at 200 MHz are shown in Figure 11. The performance of the MUSER array in metric and decametric wave range is listed in Table 3.
The extension of MUSER to lower frequency range of 30-400 MHz will be built in the next a few years. These radio facilities will be very important tools for monitoring solar disturbances from the Sun to the Earth environment and they will play fundamental role for space weather studies and monitoring.

SUMMARY
A dedicated solar radio interferometer MUSER was built in China, covering the frequency range 0.4-15 GHz, targeting high time, space and frequency resolutions with an aim to produce simultaneous solar radio images over a wide frequency band. The images obtained with MUSER can provide a unique opportunity to study different types of solar eruptive activities, such as flares and CMEs processes and their evolutions in decimetric and centimetric radio wavelengths. The few initial results as discussed here provides the evidence for the successful use of MUSER observations yielding the measurements and imaging of solar magnetic fields covering from the solar chromosphere to the higher corona, helping to understand the physics of various solar activities involved and the basic drivers of space weather activities. The MUSER and its extension to metric and decametric wavelengths will further play the role of the newgeneration radioheliographs, and will be the major solardedicated leading radio facility in the world for carrying out solar physics and space weather studies.

DATA AVAILABILITY STATEMENT
Further inquiries can be directed to the corresponding author. Our observed events are listed from the link: http://mos.bao.ac.cn

AUTHOR CONTRIBUTIONS
YY led the design, construction and operation of the MUSER and developed the observational strategy, methodology, and calibration for MUSER system; He also guided the analyses of MUSER data, and wrote the draft manuscript; ZC proposed the methodology for achieving high-performance feeds of MUSER and coordinated the construction of MUSER antenna and feed sub-system; WW designed the configurations for MUSER arrays and developed software for imaging and calibration; FL developed the MUSER digital correlator prototype and coordinated the construction of monitoring and digital correlating sub-system of MUSER; LG coordinated the construction of MUSER analogous receiving subsystem; LC implemented the software for imaging and calibration; CT developed pipeline for spectra analysis and visualization under the guidance of YY and WW; XC conducted the event analysis under the guidance of YY, WW, LC and BT; CS contributed to the MUSER feed performance evaluations and the configuration for future MUSER metric and decametric wave array under the guidance of YY and WW; and BT contributed to the interpretation of MUSER observations. All authors discussed the interpretation of the data, contributed scientific and technical results, and helped prepare the paper.