We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304\AA\ images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology.

We investigated the relationship between the spectral structures of type II solar radio bursts in the hectometric and kilometric wavelength ranges and solar energetic particles (SEPs). To examine the statistical relationship between type II bursts and SEPs, we selected 26 coronal mass ejection (CME) events with similar characteristics (e.g., initial speed, angular width, and location) observed by the Large Angle and Spectrometric Coronagraph (LASCO), regardless of the characteristics of the corresponding type II bursts and the SEP flux. Then, we compared associated type II bursts observed by the Radio and Plasma Wave Experiment (WAVES) onboard the Wind spacecraft and the SEP flux observed by the Geostationary Operational Environmental Satellite (GOES) orbiting around the Earth. We found that the bandwidth of the hectometric type II bursts and the peak flux of the SEPs has a positive correlation (with a correlation coefficient of 0.64). This result supports the idea that the nonthermal electrons of type II bursts and the nonthermal ions of SEPs are generated by the same shock and suggests that more SEPs may be generated for a wider or stronger CME shock with a longer duration. Our result also suggests that considering the spectral structures of type II bursts can improve the forecasting accuracy for the peak flux of gradual SEPs.

Coronal mass ejections (CMEs) cause disturbances in the environment of the Earth when they arrive at the Earth. However, the prediction of the arrival of CMEs still remains a challenge. We have developed an interplanetary scintillation (IPS) estimation system based on a global magnetohydrodynamic (MHD) simulation of the inner heliosphere to predict the arrival time of CMEs. In this system, the initial speed of a CME is roughly derived from white light coronagraph observations. Then, the propagation of the CME is calculated by a global MHD simulation. The IPS response is estimated by the three-dimensional density distribution of the inner heliosphere derived from the MHD simulation. The simulated IPS response is compared with the actual IPS observations made by the Institute for Space-Earth Environmental Research, Nagoya University, and shows good agreement with that observed. We demonstrated how the simulation system works using a halo CME event generated by a X9.3 flare observed on September 5, 2017. We find that the CME simulation that best estimates the IPS observation can more accurately predict the time of arrival of the CME at the Earth. These results suggest that the accuracy of the CME arrival time can be improved if our current MHD simulations include IPS data.

We comparatively studied the long-term variation (1992-2017) in polar brightening observed with the Nobeyama Radioheliograph, the polar solar wind velocity with interplanetary scintillation observations at the Institute for Space-Earth Environmental Research, and the coronal hole distribution computed by potential field calculations of the solar corona using synoptic magnetogram data obtained at Kitt Peak National Solar Observatory. First, by comparing the solar wind velocity (V) and the brightness temperature (T_b) in the polar region, we found good correlation coefficients (CCs) between V and T_b in the polar regions, CC = 0.91 (0.83) for the northern (southern) polar region, and we obtained the V-T_b relationship as V =12.6 (T_b-10,667)^{1/2}+432. We also confirmed that the CC of V-T_b is higher than those of V-B and V-B/f, where B and f are the polar magnetic field strength and magnetic flux expansion rate, respectively. These results indicate that T_b is a more direct parameter than B or B/f for expressing solar wind velocity. Next, we analyzed the long-term variation of the polar brightening and its relation to the area of the polar coronal hole (A). As a result, we found that the polar brightening matches the probability distribution of the predicted coronal hole and that the CC between T_b and A is remarkably high, CC = 0.97. This result indicates that the polar brightening is strongly coupled to the size of the polar coronal hole. Therefore, the reasonable correlation of V-T_b is explained by V-A. In addition, by considering the anti-correlation between A and f found in a previous study, we suggest that the V-T_b relationship is another expression of the Wang-Sheeley relationship (V-1/f) in the polar regions.

Polar brightening of the Sun at radio frequencies has been studied for almost fifty years and yet a disagreement persists between solar atmospheric models and observations. Some observations reported brightening values much smaller than the expected values obtained from the models, with discrepancies being particularly large at millimeter wavelengths. New clues to calibrate the atmospheric models can be obtained with the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) radio interferometer. In this work, we analyzed the lower limit of the polar brightening observed at 100 and 230 GHz by ALMA, during its Science Verification period, 2015 December 16-20. We find that the average polar intensity is higher than the disk intensity at 100 and 230 GHz, with larger brightness intensities at the South pole in eight of the nine maps analyzed. The observational results were compared with calculations of the millimetric limb brightnening emission for two semi-empirical atmospheric models, FAL- C (Fontenla et al. 1993) and SSC (Selhorst et al. 2005a). Both models presented larger limb intensities than the averaged observed values. The intensities obtained with the SSC model were closer to the observations, with polar brightenings of 10.5% and 17.8% at 100 and 230 GHz, respectively. This discrepancy may be due to the presence of chromospheric features (like spicules) at regions close to the limb.

We present the first high-resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of a sunspot at wavelengths of 1.3 mm and 3 mm, obtained during the solar ALMA Science Verification campaign in 2015, and compare them with the predictions of semi-empirical sunspot umbral/penumbral atmosphere models. For the first time millimeter observations of sunspots have resolved umbral/penumbral brightness structure at the chromospheric heights, where the emission at these wavelengths is formed. We find that the sunspot umbra exhibits a radically different appearance at 1.3 mm and 3 mm, whereas the penumbral brightness structure is similar at the two wavelengths. The inner part of the umbra is ~600 K brighter than the surrounding quiet Sun (QS) at 3 mm and is ~700 K cooler than the QS at 1.3 mm, being the coolest part of sunspot at this wavelength. On average, the brightness of the penumbra at 3 mm is comparable to the QS brightness, while at 1.3 mm it is ~1000 K brighter than the QS. Penumbral brightness increases towards the outer boundary in both ALMA bands. Among the tested umbral models, that of Severino et al. (1994) provides the best fit to the observational data, including both the ALMA data analyzed in this study and data from earlier works. No penumbral model amongst those considered here gives a satisfactory fit to the currently available measurements. ALMA observations at multiple mm wavelengths can be used for testing existing sunspot models, and serve as an important input to constrain new empirical models.

The total solar fluxes at 1, 2, 3.75, and 9.4 GHz were observed continuously from 1957 to 1994 at Toyokawa, and from 1994 until now at Nobeyama, Japan with the current Nobeyama Radio Polarimeters. We examined the multi-frequency and long-term datasets, and found that not only the microwave solar flux but also its monthly standard deviation well indicates the long-term variation of solar activity. Furthermore, we found that the microwave spectra at the solar minima of Cycle 20~24 agree with each other. These results show that the average atmospheric structure above the upper chromosphere in the quiet Sun has not varied for half a century, and suggest that the energy input for atmospheric heating from the sub-photosphere to the corona has not changed in the quiet Sun despite significantly differing strengths of magnetic activity in the last five solar cycles.

The polarization characteristics of zebra patterns (ZPs) in type IV solar bursts were studied. We analyzed 21 ZP events observed by the Assembly of Metric-band Aperture Telescope and Real-time Analysis System between 2010 and 2015 and identified the following characteristics: a degree of circular polarization (DCP) in the range of 0%-70%, a temporal delay of 0-70 ms between the two circularly polarized components (i.e., the right- and left-handed components), and dominant ordinary-mode emission in about 81% of the events. For most events, the relation between the dominant and delayed components could be interpreted in the framework of fundamental plasma emission and depolarization during propagation, though the values of DCP and delay were distributed across wide ranges. Furthermore, it was found that the DCP and delay were positively correlated (rank correlation coefficient R = 0.62). As a possible interpretation of this relationship, we considered a model based on depolarization due to reflections at sharp density boundaries assuming fundamental plasma emission. The model calculations of depolarization including multiple reflections and group delay during propagation in the inhomogeneous corona showed that the DCP and delay decreased as the number of reflections increased, which is consistent with the observational results. The dispersive polarization characteristics could be explained by the different numbers of reflections causing depolarization.

We have developed a digital fast Fourier transform (FFT) spectrometer made of an analog-to-digital converter (ADC) and a field-programmable gate array (FPGA). The base instrument has independent ADC and FPGA modules, which allow us to implement different spectrometers in a relatively easy manner. Two types of spectrometers have been instrumented, one with 4.096 GS/s sampling speed and 2048 frequency channels and the other with 2.048 GS/s sampling speed and 32768 frequency channels. The signal processing in these spectrometers has no dead time and the accumulated spectra are recorded in external media every 8 ms. A direct sampling spectroscopy up to 8 GHz is achieved by a microwave track-and-hold circuit, which can reduce the analog receiver in front of the spectrometer. Highly stable spectroscopy with a wide dynamic range was demonstrated in a series of laboratory experiments and test observations of solar radio bursts.

The Atacama Large Millimeter-submillimeter Array (ALMA) radio telescope has commenced science observations of the Sun starting in late 2016. Since the Sun is much larger than the field of view of individual ALMA dishes, the ALMA interferometer is unable to measure the background level of solar emission when observing the solar disk. The absolute temperature scale is a critical measurement for much of ALMA solar science, including the understanding of energy transfer through the solar atmosphere, the properties of prominences, and the study of shock heating in the chromosphere. In order to provide an absolute temperature scale, ALMA solar observing will take advantage of the remarkable fast-scanning capabilities of the ALMA 12m dishes to make single-dish maps of the full Sun. This article reports on the results of an extensive commissioning effort to optimize the mapping procedure, and it describes the nature of the resulting data. Amplitude calibration is discussed in detail: a path that utilizes the two loads in the ALMA calibration system as well as sky measurements is described and applied to commissioning data. Inspection of a large number of single-dish datasets shows significant variation in the resulting temperatures, and based on the temperature distributions we derive quiet-Sun values at disk center of 7300 K at lambda=3 mm and 5900 K at lambda=1.3 mm. These values have statistical uncertainties of order 100 K, but systematic uncertainties in the temperature scale that may be significantly larger. Example images are presented from two periods with very different levels of solar activity. At a resolution of order 25 arcsec, the 1.3 mm wavelength images show temperatures on the disk that vary over about a 2000 K range.

We report the discovery of a brightness enhancement in the center of a large sunspot umbra at a wavelength of 3 mm using the Atacama Large Millimeter/sub-millimeter Array (ALMA). Sunspots are amongst the most prominent features on the solar surface, but many of their aspects are surprisingly poorly understood. We analyzed a {\lambda}=3 mm (100 GHz) mosaic image obtained by ALMA, which includes a large sunspot within the active region AR12470 on December 16, 2015. The 3 mm map has a field-of-view and spatial resolution, which is the highest spatial-resolution map of an entire sunspot in this frequency range. We find a gradient of 3 mm brightness from a high value in the outer penumbra to a low value in the inner penumbra/outer umbra. Within the inner umbra, there is a marked increase in 3mm brightness temperature, which we call an umbral brightness enhancement. This enhanced emission corresponds to a temperature excess of 800 K relative to the surrounding inner penumbral region and coincides with excess brightness in the 1330 and 1400 {\AA} slitjaw images of the Interface Region Imaging Spectrograph (IRIS), adjacent to a partial lightbridge. This {\lambda}=3 mm brightness enhancement may be an intrinsic feature of the sunspot umbra at chromospheric heights, such as a manifestation of umbral flashes, or it could be related to a coronal plume since the brightness enhancement was coincident with the footpoint of a coronal loop observed at 171 {\AA}.

Eruptive phenomena such as plasmoid ejections or jets are an important feature of solar activity with the potential for improving our understanding of the dynamics of the solar atmosphere. Such ejections are often thought to be signatures of the outflows expected in regions of fast magnetic reconnection. The 304 A EUV line of Helium, formed at around 10^5 K, is found to be a reliable tracer of such phenomena, but the determination of physical parameters from such observations is not straightforward. We have observed a plasmoid ejection from an X-ray bright point simultaneously at millimeter wavelengths with ALMA, at EUV wavelengths with AIA, in soft X-rays with Hinode/XRT. This paper reports the physical parameters of the plasmoid obtained by combining the radio, EUV and X-ray data. As a result, we conclude that the plasmoid can consist either of (approximately) isothermal 10^5 K plasma that is optically thin at 100 GHz, or else a 10^4 K core with a hot envelope. The analysis demonstrates the value of the additional temperature and density constraints that ALMA provides, and future science observations with ALMA will be able to match the spatial resolution of space-borne and other high-resolution telescopes.

Observations of the Sun at millimeter and submillimeter wavelengths offer a unique probe into the structure, dynamics, and heating of the chromosphere; the structure of sunspots; the formation and eruption of prominences and filaments; and energetic phenomena such as jets and flares. High-resolution observations of the Sun at millimeter and submillimeter wavelengths are challenging due to the intense, extended, low- contrast, and dynamic nature of emission from the quiet Sun, and the extremely intense and variable nature of emissions associated with energetic phenomena. The Atacama Large Millimeter/submillimeter Array (ALMA) was designed with solar observations in mind. The requirements for solar observations are significantly different from observations of sidereal sources and special measures are necessary to successfully carry out this type of observations. We describe the commissioning efforts that enable the use of two frequency bands, the 3 mm band (Band 3) and the 1.25 mm band (Band 6), for continuum interferometric-imaging observations of the Sun with ALMA. Examples of high-resolution synthesized images obtained using the newly commissioned modes during the solar commissioning campaign held in December 2015 are presented. Although only 30 of the eventual 66 ALMA antennas were used for the campaign, the solar images synthesized from the ALMA commissioning data reveal new features of the solar atmosphere that demonstrate the potential power of ALMA solar observations. The ongoing expansion of ALMA and solar-commissioning efforts will continue to enable new and unique solar observing capabilities.

The absolute brightness temperature of the Sun at millimeter wavelengths is an important diagnostic of the solar chromosphere. Because the Sun is so bright, measurement of this property usually involves the operation of telescopes under extreme conditions and requires a rigorous performance assessment of the telescope. In this study, we establish solar observation and calibration techniques at 2.6-mm wavelength for the Nobeyama 45-m telescope and derive the absolute solar brightness temperature accurately. We tune the superconductor-insulator-superconductor (SIS) receiver by inducing different bias voltages onto the SIS mixer to prevent saturation. Then, we examine the linearity of the receiver system by comparing outputs derived from different tuning conditions. Further, we measure the lunar filled beam efficiency of the telescope using the New Moon, and then derive the absolute brightness temperature of the Sun. The derived solar brightness temperature is 7700+-310 K at 115 GHz. The telescope beam pattern is modeled as a summation of three Gaussian functions and derived using the solar limb. The real shape of the Sun is determined via deconvolution of the beam pattern from the observed map. Such well-calibrated single-dish observations are important for high-resolution chromospheric studies because they provide the absolute temperature scale missing from interferometer observations.

We investigate the upper chromosphere and the transition region of the sunspot umbra using the radio brightness temperature at 34 GHz (corresponding to 8.8-mm observations) as observed by the Nobeyama Radioheliograph (NoRH). Radio free-free emission in the longer millimeter range is generated around the transition region, and its brightness temperature yields the region's temperature and density distribution. We use the NoRH data at 34 GHz by applying the Steer-CLEAN image synthesis. These data and the analysis method enable us to investigate the chromospheric structures in the longer millimeter range with high spatial resolution and sufficient visibilities. We also perform simultaneous observations of one sunspot using the NoRH and the Nobeyama 45-m telescope operating at 115 GHz. We determine that 115-GHz emission mainly originates from the lower chromosphere while 34-GHz emission mainly originates from the upper chromosphere and transition region. These observational results are consistent with the radio emission characteristics estimated from the current atmospheric models of the chromosphere. On the other hand, the observed brightness temperature of the umbral region is almost the same as that of the quiet region. This result is inconsistent with the current sunspot models, which predict a considerably higher brightness temperature of the sunspot umbra at 34 GHz. This inconsistency suggests that the temperature of the region at which the 34 GHz radio emission becomes optically thick should be lower than that predicted by the models.

We estimated the accuracy of coronal magnetic fields derived from radio observations by comparing them to potential field calculations and the DEM measurements using EUV observations. We derived line of sight component of the coronal magnetic field from polarization observations of the thermal bremsstrahlung in the NOAA active region 11150, observed around 3:00 UT on February 3, 2011 using the Nobeyama Radioheliograph at 17 GHz. Because the thermal bremsstrahlung intensity at 17 GHz includes both chromospheric and coronal components, we extracted only the coronal component by measuring the coronal emission measure in EUV observations. In addition, we derived only the radio polarization component of the corona by selecting the region of coronal loops and weak magnetic field strength in the chromosphere along the line of sight. The upper limit of the coronal longitudinal magnetic fields were determined as 100 - 210 G. We also calculated the coronal longitudinal magnetic fields from the potential field extrapolation using the photospheric magnetic field obtained from the Helioseismic and Magnetic Imager (HMI). However, the calculated potential fields were certainly smaller than the observed coronal longitudinal magnetic field. This discrepancy between the potential and the observed magnetic field strengths can be explained consistently by two reasons; (a) the underestimation of the coronal emission measure resulting from the limitation of the temperature range of the EUV observations, (b) the underestimation of the coronal magnetic field resulting from the potential field assumption.

We investigated the polarization characteristics of a zebra pattern (ZP) in a type-IV solar radio burst observed with AMATERAS on 2011 June 21 for the purpose of evaluating the generation processes of ZP. Analyzing highly resolved spectral and polarization data revealed the frequency dependence of the degree of circular polarization and the delay between two polarized components for the first time. The degree of circular polarization was 50-70 percent right-handed and it varied little as a function of frequency. Cross-correlation analysis determined that the left-handed circularly polarized component was delayed by 50-70 ms relative to the right-handed component over the entire frequency range of the ZP and this delay increased with the frequency. We examined the obtained polarization characteristics by using pre-existing ZP models and concluded that the ZP was generated by the double plasma resonance process. Our results suggest that the ZP emission was originally generated in a completely polarized state in the O-mode and was partly converted into the X-mode near the source. Subsequently, the difference between the group velocities of the O-mode and X-mode caused the temporal delay.

The brightness temperature of the radio free-free emission at millimeter range is an effective tool for characterizing the vertical structure of the solar chromosphere. In this paper, we report on the first single-dish observation of a sunspot at 85 and 115 GHz with sufficient spatial resolution for resolving the sunspot umbra using the Nobeyama 45 m telescope. We used radio attenuation material, i.e. a solar filter, to prevent the saturation of the receivers. Considering the contamination from the plage by the side-lobes, we found that the brightness temperature of the umbra should be lower than that of the quiet region. This result is inconsistent with the preexisting atmospheric models. We also found that the brightness temperature distribution at millimeter range strongly corresponds to the ultraviolet (UV) continuum emission at 1700 {\AA}, especially at the quiet region.

We derived the coronal magnetic field, plasma density, and temperature from the observation of polarization and intensity of radio thermal free-free emission using the Nobeyama Radioheliograph (NoRH) and extreme ultraviolet (EUV) observations. We observed a post-flare loop on the west limb 11 April 2013. The line-of-sight magnetic field was derived from the circularly polarized free-free emission observed by NoRH. The emission measure and temperature were derived from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The derived temperature was used to estimate the emission measure from the NoRH radio free-free emission observations. The derived density from NoRH was larger than that determined using AIA, which can be explained by the fact that the low temperature plasma is not within the temperature coverage of the AIA filters used in this study. We also discuss the other observation of the post-flare loops by the EUV Imager onboard the Solar Terrestrial Relations Observatory (STEREO), which can be used in future studies to reconstruct the coronal magnetic field strength. The derived plasma parameters and magnetic field were used to derive the plasma beta, which is a ratio between the magnetic pressure and the plasma pressure. The derived plasma beta is about 5.7*10^(-4) to 7.6*10^(-4) at the loop top region.

The fine spectral structures of solar radio type-I bursts were observed by the solar radio telescope AMATERAS. The spectral characteristics, such as the peak flux, duration, and bandwidth, of the individual burst elements were satisfactorily detected by the highly resolved spectral data of AMATEAS with the burst detection algorithm that is improved in this study. The peak flux of the type-I bursts followed a power-law distribution with a spectral index of 2.9-3.3, whereas their duration and bandwidth were distributed more exponentially. There were almost no correlations between the peak flux, duration, and bandwidth. That means there were no similarity shapes in the burst spectral structures. We defined the growth rate of a burst as the ratio between its peak flux and duration. There was a strong correlation between the growth rate and peak flux. These results suggest that the free energy of type-I bursts that is originally generated by non-thermal electrons is modulated in the subsequent stages of the generation of non-thermal electrons, such as plasma wave generation, radio wave emissions, and propagation. The variation of the time scale of the growth rate is significantly larger than that of the coronal environments. These results can be explained by the situation that the source region may have the inhomogeneity of an ambient plasma environment, such as the boundary of open and closed field lines, and the superposition of entire emitted bursts was observed by the spectrometer.

Coronal and chromospheric magnetic fields are derived from polarization and spectral observations of the thermal free-free emission using the Nobeyama Radioheliograph (NoRH). In magnetized plasma, the ordinary and extraordinary modes of free-free emission have different optical depths. This creates a circularly polarized component in an atmosphere with a temperature gradient. We observed an active region on April 13, 2012 to derive its coronal and chromospheric magnetic fields. The observed degree of circular polarization was between 0.5 % and 1.7 %. The radio circular polarization images were compared with ultraviolet images observed by the Atmospheric Imaging Assembly and the photospheric magnetic field observed by the Helioseismic and Magnetic Imager, both on board the Solar Dynamic Observatory. At the edge of the active region, the radio circular polarization was emitted mainly from coronal loops, and the coronal magnetic field was derived to be about 70 G. At the center of the active region, the chromospheric and coronal components cannot be separated. The derived magnetic field is about 20 % to 50 % of the corresponding photospheric magnetic field, which is an emission-measure-weighted average of the coronal and chromospheric magnetic fields.