We evaluate validity of NLFFF extrapolation performed with Optimization class (OPTI) codes. While explaining inevitable for OPTI partial non-solenoidality caused by the gas pressure notable role in pressure balance at photospheric heights and by mathematical aspects related to optimization and BVP (boundary value problem), we justify elimination of the non-solenoidal component (postprocessing) from the OPTI result obtained. In essence, postprocessing converts the entire non-solenoidal part into a solenoidal force part, which possibly reflects factual deviation of magnetic field from its force-free pproximation on the photosphere and in the solar corona. Two forms of postprocessing have been analyzed in this paper.Postprocessing I eliminates the non-solenoidal component without changing transverse field at the measurement level, and Postprocessing II leaves the field normal component unchanged. Extrapolation, postprocessing, and then comparison of metric and energy characteristics are performed over AR 11158 active region for a small fragment of its evolution containing the February X-class flare.Our version of OPTI code showed that free energy decreased by $\sim 10^{32}$ erg within 1 hour, which corresponds to theoretical estimations of the flare-caused magnetic energy loss. This result differs significantly from the one in Sun et al. (2012).Therefore, we also comment on some features of our OPTI code implementation, which may cause significant differences between our results and those obtained using the Wiegelmann (2004) version of OPTI code in study by Sun et al. (2012).

In this study, we investigate magnetic properties of umbra of magneto-conjugate leading and following sunspots, i.e. connected through magnetic field lines. We established dependences between individual sunspot umbra field characteristics, and between these characteristics on the umbra area ($S$) separately for sunspot pairs, for which the minimal angle between the umbra magnetic field line of the leading ($L$) sunspot and the positive normal line to the Sun surface is smaller, than that in the following ($F$) sunspot ($\alpha_{min-L}<\alpha_{min-F}$; such sunspot pairs are the bulk) and, on the contrary, when $\alpha_{min-L}>\alpha_{min-F}$. The $\alpha_{min-L}(S_L)$, $\alpha_{min-F}(S_F)$, $B_{max-L}(S_L)$ and $B_{max-F}(S_F)$ dependences are shown to have similar behavior features, and are quantitatively close for two sunspot groups with a different asymmetry of the sunspot magnetic field connecting them (here, $B_{max-L,F}(S_L)$ is the magnetic induction maximum induction in umbrae of the leading and the following sunspots). The dependence of mean values of angles within umbra $<\alpha_{L,F}>$ on the sunspot umbra area $S_{L,F}$ and on the mean value of magnetic induction in umbra $<B_{L,F}>$ appeared different for two cases. Also, in the bulk of the investigated sunspot pairs, the leading sunspot was shown to appear closer to the polarity inversion line between the sunspots, than the following one. This result and the conclusion that, in the bulk of the investigated pairs of the magnetically conjugate sunspots, $\alpha_{min-L}<\alpha_{min-F}$ are closely coupled.

In this paper, we came to conclusion that there is a significant systematic error in the SDO/HMI vector magnetic data, which reveals itself in a significant deviation of the lines of the knot magnetic fields from the radial direction. The value of this deviation demonstrates a clear dependence on the distance to the disk center. This paper suggests a method for correction of the vector magnetograms that eliminates the detected systematic error.

Calculation of the potential field inside a three-dimensional box with the normal magnetic field component given on all boundaries is needed for estimation of important quantities related to the magnetic field such as free energy and relative helicity. In this work we present an analysis of three methods for calculating potential field inside a three-dimensional box. The accuracy and performance of the methods are tested on artificial models with a priori known solutions.

From vector measurements of the photospheric magnetic field with the SDO/HMI instrument, we studied the field variations within the 2011 June 7 eruptive event related to the filament eruption (FE), flare, and coronal mass ejection (CME). We analyzed the variations in modulus (B), in radial (Br) and transversal (Bt) magnetic induction components, as well as in the field line inclination angle ({\alpha}) to the radial direction from the Sun center were. The conclusion is that the most probable cause for the FE onset being the trigger related to the flares and CME, as well as for Stage 1 of the filament dramatic acceleration, is {\guillemotleft}magnetic flux cancellation{\guillemotright} in several sites inside and around the filament channel. We assumed that the flare and the stage of the filament fastest acceleration are related to a spot emergence near the south-eastern filament footpoint, in which the field polarity and value satisfy a possibility of this field magnetic reconnection with the field surrounding the filament. We studied in detail the time variations in the field line inclination angles within different sites of the eruptive event. During the filament slow emergence, a decrease in the field line inclination angles was revealed to occur in its channel neighboring, whereas the field line inclination angles dramatically grew in the neutral line neighboring within the flare region after the flare onset. The flare ribbons at all the stages of their existence are shown to be located under the photosphere sites with the field modulus local maxima and with minima of the field line inclination angles. We show that, near the photospheric magnetic field polarity inversion line (PIL), the azimuth decreases after the flare onset, which means a decrease in the angle \b{eta} between PIL and the magnetic induction vector projection onto the sky plane.

We have analyzed multi-wavelength observations and magnetic-field data for the solar flare of May 10, 2012 (04:18 UT) and have detected a sign inversion of the signal in the line-of-sight magnetic measurements in the umbra of a small spot. This effect is associated, at least partly, with the emergence of a new magnetic field. Almost at the same time, a burst of hard X-rays was recorded, and a wave in the vacuum ultraviolet (EUV) range (a "sunquake") was generated due to the impact of the disturbance in the energy release range on the photosphere. At the beginning of the event, a sigmoid flare was recorded, but it did not spread, as it usually does, along the polarity inversion (neutral) line. SDO/HMI full-vector measurements were used to extrapolate the AR 11476 magnetic field to the corona, and the distribution of vertical currents $j_z$ in the photosphere was obtained. The distribution of currents in the active region shows that the relationship between them and the occurrence of flares is very intricate. We have corroborated that the expected "ideal" behavior of the current system before and after the flare (e.g., see (Sharykin and Kosovichev, 2015)) is observed only in the sigmoid region. The results obtained were compared with the observations of two other flares recorded in this AR on the same day, one of which was similar to the flare under discussion and the other was of different type. Our results confirm that the formation and eruption of large-scale magnetic flux ropes in sigmoid flares are associated with the shear motions in the photosphere and the emergence of twisted magnetic tubes, as well as with the subsequent development of the torus instability.

A method of construction of solution for acoustic-gravity waves (AGW) above a wave source, taking dissipation throughout the atmosphere into account (Dissipative Solution above Source, DSAS), is proposed. The method is to combine three solutions for three parts of the atmosphere: an analytical solution for the upper isothermal part and numerical solutions for the real non-isothermal dissipative atmosphere in the middle part and for the real non-isothermal small dissipation atmosphere in the lower one. In this paper the method has been carried out for the atmosphere with thermal conductivity but without viscosity. The heights of strong dissipation and the total absorption index in the regions of weak and average dissipation are found. For internal gravity waves the results of test calculations for an isothermal atmosphere and calculations for a real non-isothermal atmosphere are shown in graphical form. An algorithm and appropriate code to calculate DSAS, taking dissipation due to finite thermal conductivity into account throughout the atmosphere, are developed. The results of test DSAS calculations for an everywhere isothermal atmosphere are given. The calculation results for DSAS for the real non-isothermal atmosphere are also presented. A method for constructing of the 2x2 Green's matrix fully taking dissipation into account and allowing to find disturbance from some source of AGW in the atmosphere is proposed.

Trapped atmosphere waves, such as IGW waveguide modes and Lamb modes, are described using dissipative solution above source (DSAS) (Dmitrienko and Rudenko, 2016). Accordingly this description, the modes are disturbances penetrating without limit in the upper atmosphere and dissipating their energy throughout the atmosphere; leakage from a trapping region to the upper atmosphere is taken in consideration. The DSAS results are compared to those based on both accurate and WKB approximated dissipationless equations. It is shown that the spatial and frequency characteristics of modes in the upper atmosphere calculated by any of the methods are close to each other and are in good agreement with the observed characteristics of traveling ionospheric disturbances.

This is the first of four companion papers, which analyze a complex eruptive event of 18 November 2003 in AR 10501 and the causes of the largest Solar Cycle 23 geomagnetic storm on 20 November 2003. Analysis of a complete data set, not considered before, reveals a chain of eruptions to which hard X-ray and microwave bursts responded. A filament in AR 10501 was not a passive part of a larger flux rope, as usually considered. The filament erupted and gave origin to a CME. The chain of events was as follows: i) an eruption at 07:29 accompanied by a not reported M1.2 class flare associated with the onset of a first southeastern CME1, which is not responsible for the superstorm; ii) a confined eruption at 07:41 (M3.2 flare) that destabilized the filament; iii) the filament acceleration (07:56); iv) the bifurcation of the eruptive filament that transformed into a large cloud; v) an M3.9 flare in AR 10501 associated to this transformation. The transformation of the filament could be due to its interaction with the magnetic field in the neighborhood of a null point, located at a height of about 100 Mm above the complex formed by ARs 10501, 10503, and their environment. The CORONAS-F/SPIRIT telescope observed the cloud in 304 A as a large Y-shaped darkening, which moved from the bifurcation region to the limb. The masses and kinematics of the cloud and the filament were similar. Remnants of the filament were not observed in the second southwestern CME2, previously regarded as a source of the 20 November superstorm. These facts do not support a simple scenario, in which the interplanetary magnetic cloud is considered as a flux rope formed from a structure initially associated with the pre-eruption filament in AR 10501. Observations suggest a possible additional eruption above the bifurcation region close to solar disk center between 08:07 and 08:17 that could be the source of the superstorm.

Our analysis in Papers I and II (Grechnev et al., 2014, Solar Phys. 289, 289 and 1279) of the 18 November 2003 solar event responsible for the 20 November geomagnetic superstorm has revealed a complex chain of eruptions. In particular, the eruptive filament encountered a topological discontinuity located near the solar disk center at a height of about 100 Mm, bifurcated, and transformed into a large cloud, which did not leave the Sun. Concurrently, an additional CME presumably erupted close to the bifurcation region. The conjectures about the responsibility of this compact CME for the superstorm and its disconnection from the Sun are confirmed in Paper IV (Grechnev et al., Solar Phys., submitted), which concludes about its probable spheromak-like structure. The present paper confirms the presence of a magnetic null point near the bifurcation region and addresses the origin of the magnetic helicity of the interplanetary magnetic clouds and their connection to the Sun. We find that the orientation of a magnetic dipole constituted by dimmed regions with the opposite magnetic polarities away from the parent active region corresponded to the direction of the axial field in the magnetic cloud, while the pre-eruptive filament mismatched it. To combine all of the listed findings, we come to an intrinsically three-dimensional scheme, in which a spheromak-like eruption originates via the interaction of the initially unconnected magnetic fluxes of the eruptive filament and pre-existing ones in the corona. Through a chain of magnetic reconnections their positive mutual helicity was transformed into the self-helicity of the spheromak-like magnetic cloud.

The Siberian Solar Radio Telescope (SSRT) is a solar-dedicated directly-imaging interferometer observing the Sun at 5.7 GHz. The SSRT operates in the two-dimensional mode since 1996. The imaging principle of the SSRT restricts its opportunities in observations of very bright flare sources, while it is possible to use `dirty' images in studies of low brightness features, which do not overlap with side lobes from bright sources. The interactive CLEAN technique routinely used for the SSRT data provides imaging of active regions but consumes much time and efforts and does not reveal low-brightness features below the CLEAN threshold. The newly developed technique combines the CLEAN routine with the directly imaging capability of the SSRT and provides clean images with an enhanced dynamic range automatically. These elaborations considerably extend the range of tasks, which can be solved with the SSRT. We show some examples of the present opportunities of the SSRT and compare its data with the images produced by the Nobeyama Radioheliograph at 17 GHz as well as observations in different spectral ranges.

We use daily full-disk vector magnetograms from Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) system to synthesize the first Carrington maps of the photospheric vector magnetic field. We describe these maps and make a comparison of observed radial field with the radial field estimate from LOS magnetograms. Further, we employ these maps to study the hemispheric pattern of current helicity density, Hc, during the rising phase of the solar cycle 24. Longitudinal average over the 23 consecutive solar rotations shows a clear signature of the hemispheric helicity rule, i.e. Hc is predominantly negative in the North and positive in South. Although our data include the early phase of cycle 24, there appears no evidence for a possible (systematic) reversal of the hemispheric helicity rule at the beginning of cycle as predicted by some dynamo models. Further, we compute the hemispheric pattern in active region latitudes (-30 deg \le \theta \le 30 deg) separately for weak (100< |B_r| <500 G)and strong (|B_r|>1000 G) radial magnetic fields. We find that while the current helicity of strong fields follows the well-known hemispheric rule (i.e., \theta . Hc < 0), H_c of weak fields exhibits an inverse hemispheric behavior (i.e., \theta . Hc > 0) albeit with large statistical scatter. We discuss two plausible scenarios to explain the opposite hemispheric trend of helicity in weak and strong field region.

The paper presents the possibility of fast and quality azimuth disambiguation of vector magnetogram data regardless of location on the solar disc. The new Super Fast and Quality (SFQ) code of disambiguation is tried out on well-known models of Metcalf et al. (2006), Leka et al. (2009) and artificial model of fixed configuration AR 10930 (Rudenko et al., 2010). We make comparison of Hinode SOT SP vector magnetograms of AR 10930 disambiguated with three codes: SFQ, NPFC (Georgoulis, 2005), and SME (Rudenko et al., 2010). We exemplify the SFQ disambiguation of SDO/HMI measurements of the full disc. The preliminary examination indicates that the SFQ algorithm provides better quality than NPFC and is comparable to SME. In contrast to other codes, SFQ supports relatively high quality of results regardless of the magnetogram proximity to the limb (when being very close to the limb, it is efficient unlike all other algorithms).

On 17 January 2010, STEREO-B observed in extreme ultraviolet (EUV) and white light a large-scale dome-shaped expanding coronal transient with perfectly connected off-limb and on-disk signatures. Veronig et al. (2010, ApJL 716, 57) concluded that the dome was formed by a weak shock wave. We have revealed two EUV components, one of which corresponded to this transient. All of its properties found from EUV, white light, and a metric type II burst match expectations for a freely expanding coronal shock wave including correspondence to the fast-mode speed distribution, while the transient sweeping over the solar surface had a speed typical of EUV waves. The shock wave was presumably excited by an abrupt filament eruption. Both a weak shock approximation and a power-law fit match kinematics of the transient near the Sun. Moreover, the power-law fit matches expansion of the CME leading edge up to 24 solar radii. The second, quasi-stationary EUV component near the dimming was presumably associated with a stretched CME structure; no indications of opening magnetic fields have been detected far from the eruption region.

An algorithm for calculating three gauge-invariant helicities (self-, mutual- and Berger relative helicity) for a magnetic field specified in a rectangular box is described. The algorithm is tested on a well-known force-free model (Low and Lou, 1990) presented in vector-potential form.

Possibilities in principle for satisfactory removal of the 180-azimuthal ambiguity in the transverse field of vector magnetograms and the extrapolation of magnetic fields independently of their position on the solar disk are shown. Revealed here is an exact correspondence between the estimated field and the nonpotential loop structure on the limb. The Metropolis's algorithm modified to work in spherical geometry is used to resolve the azimuthal ambiguity. Based on a version of the optimization method from Rudenko and Myshyakov (2009), we use corrected magnetograms as boundary conditions for magnetic field extrapolation in the nonlinear force-free approximation.

A technique is proposed for constructing three-dimensional solutions comlplying with the self-consistent magnetohydrostatic (MHS) equations and with observations along the line of sight of the magnetic field at the photosphere. The technique is a generalization of a paper Rudenko, {\tt arXiv:astro-ph/0004222} to a potential-field approximation. The solution of the problem under consideration leans a representation of the magnetic field in terms of a scalar function, with its subsequent harmonic expansion in terms of a spherical functional basis that satisfies specified boundary conditions. A numerical realization of the proposed method is expected to permit a real-time modelling of three-dimensional magnetic field, temperature, pressure and density distributions.

This paper is concerned with the Laplace boundary-value problem with the directional derivative, corresponding to the specific nature of measurements of the longitudinal component of the photospheric magnetic field. Boundary conditions are specified by a distribution on the sphere of projection of the magnetic field vector unto a given direction. It is shown that the solution of this problem exists in the form of a spherical harmonic expansion, and uniqueness of this solution is proved. A conceptual sketch of numerical determination of the harmonic series coefficients is given. The field of application of the method is analyzed having regard to the peculiarities of actual data. Finally, we present differences in results derived from extrapolating the magnetic field from a synoptic map and a full-disk magnetogram.