Umbral flashes (UFs) are emissions in the core of chromospheric lines caused by upward propagating waves steepening into shocks. UFs are followed by an expanding blue shifted umbral wave (UW) and red-shifted plasma returning to the initial state. Here we use 5 s cadence images acquired at $\pm$0.04~nm off the Halpha line center by the Visible Imaging Spectrometer (VIS) installed on the Goode Solar Telescope (GST) to detect the origin of UFs and UWs in a sunspot with a uniform umbra free of LBs and clusters of umbral dots. The data showed that UFs do not randomly originate over the umbra. Instead, they appear to be repeatedly triggered at locations with the lowest umbral intensity and the most powerful oscillations of Halpha-0.04 nm intensity. GST magnetic field measurements using Near Infra-Red Imaging Spectropolarimeter (NIRIS) also showed that the dominant location of prevalent UF origin is co-spatial associated with the strongest fields in the umbra. Interface Region Imaging Spectrograph 149.0 nm images showed that no bright UV loops were anchored in the umbra in general and near the UF patches in particular suggesting that UFs and UWs alone can not be responsible for the origin of warm coronal loops. We thus conclude that the existence of locations with prevalent origin of UFs confirms the idea that they may be driven by a sub-surface source located near the axis of a flux rope, while the presence of several UFs trigger centers may indicate the complex structure of a sunspot umbra.

We analyze ground-based chromospheric data acquired at a high temporal cadence of 2 s in wings of the H$\alpha$ spectral line using Goode Solar Telescope (GST) operating at the Big Bear Solar Observatory. We inspected a 30 minute long H$\alpha$-0.08~nm data set to find that rapid blue-shifted H$\alpha$ excursions (RBEs), which are a cool component of type II spicules, experience very rapid morphological changes on the time scales of the order of 1 second. Unlike typical reconnection jets, RBEs very frequently appear \textit{in situ} without any clear evidence of H$\alpha$ material being injected from below. Their evolution includes inverted "Y", "V", "N", and parallel splitting (doubling) patterns as well as sudden formation of a diffuse region followed by branching. We also find that the same feature may undergo several splitting episodes within about 1 min time interval.

Using magnetographic data provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, we analyzed the structure of magnetic fields and vertical electric currents in six active regions (ARs) with different level of flare activity. We found that electric currents are well balanced over the entire AR: for all of them the current imbalance is below 0.1%, which means that any current system is closed within an AR. Decomposition of the transverse magnetic field vector into two components allowed us to reveal the existence of large-scale vortex structures of the azimuthal magnetic field component around main sunspots of ARs. In each AR, we found a large-scale electric current system occupying a vast area of an AR, which we call distributed electric current. For ARs obeying the Hale polarity law and the hemispheric helicity sign rule, the distributed current is directed upward in the leading part of an AR and it appears to be closing back to the photosphere in the following part of an AR through the corona and chromosphere. Our analysis of the time variations of the magnitude of the distributed electric currents showed that low-flaring ARs exhibit small variations of the distributed currents in the range of $\pm 20 \times 10^{12}$ A, whereas the highly flaring ARs exhibited significant slow variations of the distributed currents in the range of $30-95 \times 10^{12}$ A. Intervals of the enhanced flaring appear to be co-temporal with smooth enhancements of the distributed electric current.

Magnetic field data provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory were utilized to explore the changes in the magnetic energy of four active regions (ARs) during their emergence. We found that at the very early stage of the magnetic flux emergence, an abrupt steepening of the magnetic power spectrum takes place leading to rapid increase of the absolute value of the negative spectra power index $\alpha$ in $E(k)$ $\sim$ $k^{\alpha}$. As the emergence proceeds, the energy increases at all scales simultaneously implying that elements of all sizes do appear in the photosphere. Meanwhile, the energy gain at scales larger than $\approx$10 Mm prevails as compared to that at smaller scales. Both direct (i.e., fragmentation of large structures into smaller ones) and inverse (i.e., merging of small magnetic features into larger elements) cascades are readily observed during the emergence. However, in the case of inverse cascade, the total energy gained at large scales exceeds the energy loss at smaller scales assuming simultaneous appearance of large-scale magnetic entities from beneath the photosphere. We conclude that most of the time the energy may grow at all scales. We also cannot support the point of view regarding the dominant role of the inverse cascade in the formation of an AR. Although the coalescence of small magnetic elements into larger pores and sunspots is observed, our analysis shows that the prevailed energy contribution to an AR comes from emergence of large-scale structures.

A process of diffusion of small-scale magnetic elements inside four active regions (ARs) was analyzed. Line-of-sight magnetograms acquired by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO) during a two-day time interval around the AR culmination time were utilized. Small magnetic elements of size of 3-100 squared HMI pixels with the field strength above the detection threshold of 30 Mx cm$^{-2}$ were detected and tracked. The turbulent diffusion coefficient was retrieved using the pair-separation technique. Comparison with the previously reported quiet-sun (QS) diffusivity was performed. It was found that: i) dispersion of small-scale magnetic elements inside the AR area occurs in the regime close to normal diffusion, whereas well-pronounced super-diffusion is observed in QS; ii) the diffusivity regime operating in an AR (the magnitude of the spectral index and the range of the diffusion coefficient) does not seem to depend on the individual properties of an AR, such as total unsigned magnetic flux, state of evolution, and flaring activity. We conclude that small-scale magnetic elements inside an AR do not represent an undisturbed photosphere, but they rather are intrinsic part of the whole coherent magnetic structure forming an active region. Moreover, turbulence of small-scale elements in an AR is not closely related to processes above the photosphere, but it rather carries the footprint of the sub-photospheric dynamics.

We use observations of line-of-sight magnetograms from Helioseismic and Magnetic Imager (HMI) on board of Solar Dynamics Observatory (SDO) to investigate polarity separation, magnetic flux, flux emergence rate, twist and tilt of solar emerging active regions. Functional dependence of polarity separation and maximum magnetic flux of an active region is in agreement with a simple model of flux emergence as the result of buoyancy forces. Our investigation did not reveal any strong dependence of emergence rate on twist properties of active regions.

We analyse the spatial distribution of current helicity in solar active regions. A comparison of current helicity maps derived from three different instruments (Helioseismic and Magnetc Imager on board the Solar Dynamics Observatory, SDO/HMI, Spectro-Polarimeter on board the Hinode, and Solar Magnetic Field Telescope at the Huairou Solar Observing Station, China, HSOS/SMFT) is carried out. The comparison showed an excellent correlation between the maps derived from the spaceborne instruments and moderate correlation between the maps derived from SDO/HMI and HSOS/SMFT vector magnetograms. The results suggest that the obtained maps characterize real spatial distribution of current helicity over an active region. To analyse intermittency of current helicity, we traditionally use the high-order structure functions and flatness function approach. The slope of a flatness function within some range of scales - the flatness exponent - is a measure of the degree of intermittency. SDO/HMI vector magnetograms for 3 ARs (NOAA 11158, 12494, and 12673) were used to calculate the flatness exponent time variations. All three ARs exhibited emergence of a new magnetic flux during the observational interval. The flatness exponent indicated the increase of intermittency 12-20 hours before the emergence of a new flux. We suppose that this behaviour can indicate subphotospheric fragmentation or distortion of the pre-existed electric current system by emerging magnetic flux.

Self-organization is a property of dissipative nonlinear processes that are governed by an internal driver and a positive feedback mechanism, which creates regular geometric and/or temporal patterns and decreases the entropy, in contrast to random processes. Here we investigate for the first time a comprehensive number of 16 self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous {\sl order out of chaos}, during the evolution from an initially disordered system to an ordered stationary system, via quasi-periodic limit-cycle dynamics, harmonic mechanical resonances, or gyromagnetic resonances. The internal driver can be gravity, rotation, thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational instability, the Rayleigh-B\'enard convection instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or loss-cone instability. Physical models of astrophysical self-organization processes involve hydrodynamic, MHD, and N-body formulations of Lotka-Volterra equation systems.

To explore the magnetic flux dispersion in the undisturbed solar photosphere, magnetograms acquired by Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO) were utilized. Two areas, a coronal hole area (CH) and an area of super-granulation pattern, SG, were analyzed. We explored the displacement and separation spectra and the behavior of the turbulent diffusion coefficient, $K$. The displacement and separation spectra are very similar to each other. Small magnetic elements (of size 3-100 squared pixels and the detection threshold of 20 Mx sm$^{-2}$) in both CH and SG areas disperse in the same way and they are more mobile than the large elements (of size 20-400 squared pixels and the detection threshold of 130 Mx sm$^{-2}$). The regime of super-diffusivity is found for small elements ($\gamma \approx 1.3 $ and $K$ growing from $\sim$100 to $\sim$ 300 km$^2$ s$^{-1}$). Large elements in the CH area are scanty and show super-diffusion with $\gamma \approx 1.2$ and $K$ = (62-96) km$^2$ s$^{-1}$ on rather narrow range of 500-2200 km. Large elements in the SG area demonstrate two ranges of linearity and two diffusivity regimes: sub-diffusivity on scales (900-2500) km with $\gamma=0.88$ and $K$ decreasing from $\sim$130 to $\sim$100 km$^2$ s$^{-1}$, and super-diffusivity on scales (2500-4800) km with $\gamma \approx 1.3$ and $K$ growing from $\sim$140 to $\sim$200 km$^2$ s$^{-1}$. Comparison of our results with the previously published shows that there is a tendency of saturation of the diffusion coefficient on large scales, i.e., the turbulent regime of super-diffusivity is gradually replaced by normal diffusion.

Seven-year long seeing-free observations of solar magnetic fields with the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) were used to study the sources of the solar mean magnetic field, SMMF, defined as the net line-of-sight magnetic flux divided over the solar disk area. To evaluate the contribution of different regions to the SMMF, we separated all the pixels of each SDO/HMI magnetogram into three subsets: weak (B_W), intermediate (B_I), and strong (B_S) fields. The B_W component represents areas with magnetic flux densities below the chosen threshold; the B_I component is mainly represented by network fields, remains of decayed active regions (ARs), and ephemeral regions. The B_S component consists of magnetic elements in ARs. To derive the contribution of a subset to the total SMMF, the linear regression coefficients between the corresponding component and the SMMF were calculated. We found that: i) when the threshold level of 30 Mx cm^-2 is applied, the B_I and B_S components together contribute from 65% to 95% of the SMMF, while the fraction of the occupied area varies in a range of 2-6% of the disk area; ii) as the threshold magnitude is lowered to 6 Mx cm^-2, the contribution from B_I+B_S grows to 98%, and the fraction of the occupied area reaches the value of about 40% of the solar disk. In summary, we found that regardless of the threshold level, only a small part of the solar disk area contributes to the SMMF. This means that the photospheric magnetic structure is an intermittent, inherently porous medium, resembling a percolation cluster. These findings suggest that the long-standing concept that continuous vast unipolar areas on the solar surface are the source of the SMMF may need to be reconsidered.

We studied the emergence process of 42 active region (ARs) by analyzing the time derivative, R(t), of the total unsigned flux. Line-of-sight magnetograms acquired by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) were used. A continuous piecewise linear fitting to the R(t)-profile was applied to detect an interval, dt_2, of nearly-constant R(t) covering one or several local maxima. The averaged over dt_2 magnitude of R(t) was accepted as an estimate of the maximal value of the flux growth rate, R_MAX, which varies in a range of (0.5-5)x10^20 Mx hour^-1 for active regions with the maximal total unsigned flux of (0.5-3)x10^22 Mx. The normalized flux growth rate, R_N, was defined under an assumption that the saturated total unsigned flux, F_MAX, equals unity. Out of 42 ARs in our initial list, 36 event were successfully fitted and they form two subsets (with a small overlap of 8 events): the ARs with a short (<13 hours) interval dt_2 and a high (>0.024 hour^-1) normalized flux emergence rate, R_N, form the "rapid" emergence event subset. The second subset consists of "gradual" emergence events and it is characterized by a long (>13 hours) interval dt_2 and a low R_N (<0.024 hour^-1). In diagrams of R_MAX plotted versus F_MAX, the events from different subsets are not overlapped and each subset displays an individual power law. The power law index derived from the entire ensemble of 36 events is 0.69+-0.10. The "rapid" emergence is consistent with a "two-step" emergence process of a single twisted flux tube. The "gradual" emergence is possibly related to a consecutive rising of several flux tubes emerging at nearly the same location in the photosphere.

The solar mean magnetic field (SMMF) provided by the Wilcox Solar Observatory (WSO) is compared with the SMMF acquired by the \textit{Helioseismic and Magnetic Imager} (HMI) onboard the \textit{Solar Dynamic Observatory} (SDO). We found that despite the different spectral lines and measurement techniques used in both instruments the Pearson correlation coefficient between these two datasets equals 0.86 while the conversion factor is very close to unity: B(HMI) = 0.99(2)B(WSO). We also discuss artifacts of the SDO/HMI magnetic field measurements, namely the 12 and 24-hour oscillations in SMMF and in sunspots magnetic fields that might be caused by orbital motions of the spacecraft. The artificial harmonics of SMMF reveal significant changes in amplitude and the nearly stable phase. The connection between the 24-hour harmonic amplitude of SMMF and the presence of sunspots is examined. We also found that opposite phase artificial 12 and/or 24-hour oscillations exist in sunspots of opposite polarities.

In order to clarify a possible role of small-scale dynamo in formation of solar magnetic field, we suggest an observational test for small-scale dynamo action based on statistics of anti-Hale sunspot groups. As we have shown, according to theoretical expectations the small-scale dynamo action has to provide a population of sunspot groups which do not follow the Hale polarity law, and the density of such groups on the time-latitude diagram is expected to be independent on the phase of the solar cycle. Correspondingly, a percentage of the anti-Hale groups is expected to reach its maximum values during solar minima. For several solar cycles, we considered statistics of anti-Hale groups obtained by several scientific teams, including ours, to find that the percentage of anti-Hale groups becomes indeed maximal during a solar minimum. Our interpretation is that this fact may be explained by the small-scale dynamo action inside the solar convective zone.

We analyse sunspot oscillations using Interface Region Imaging Spectrograph (IRIS) slit-jaw and spectral data and narrow-band chromospheric images from the New Solar Telescope (NST) for the main sunspot in NOAA AR 11836. We report that the difference between the shock arrival times as measured the Mg II k 2796.35\AA\ and Si IV 1393.76\AA\ line formation levels changes during the observed period and peak-to-peak delays may range from 40~s to zero. The intensity of chromospheric shocks also displays a long term (about 20~min) variations. NST's high spatial resolution \ha\ data allowed us to conclude that in this sunspot umbral flashes (UFs) appeared in the form of narrow bright lanes stretched along the light bridges and around clusters of umbral bright points. Time series also suggested that UFs preferred to appear on the sunspot-center side of light bridges, which may indicate the existence of a compact sub-photospheric driver of sunspot oscillations. The sunspot's umbra as seen in the IRIS chromospheric and transition region data appears bright above the locations of light bridges and the areas where the dark umbra is dotted with clusters of umbral dots. Co-spatial and co-temporal data from the Atmospheric Imaging Assembly on board Solar Dynamics Observatory showed that the same locations were associated with bright footpoints of coronal loops suggesting that the light bridges may play an important role in heating the coronal sunspot loops. Finally, the power spectra analysis showed that the intensity of chromospheric and transition region oscillations significantly vary across the umbra and with height, suggesting that umbral non-uniformities and the structure of sunspot magnetic fields may play a role in wave propagation and heating of umbral loops.

Recent observations of sunspot's umbra suggested that it may be finely structured at a sub-arcsecond scale representing a mix of hot and cool plasma elements. In this study we report the first detailed observations of the umbral spikes, which are cool jet-like structures seen in the chromosphere of an umbra. The spikes are cone-shaped features with a typical height of 0.5-1.0 Mm and a width of about 0.1 Mm. Their life time ranges from 2 to 3 min and they tend to re-appear at the same location. The spikes are not associated with photospheric umbral dots and they rather tend to occur above darkest parts of the umbra, where magnetic fields are strongest. The spikes exhibit up and down oscillatory motions and their spectral evolution suggests that they might be driven by upward propagating shocks generated by photospheric oscillations. It is worth noting that triggering of the running penumbral waves seems to occur during the interval when the spikes reach their maximum height.

An essential parameter for models of coronal heating and fast solar wind acceleration that rely on the dissipation of MHD turbulence is the characteristic energy-containing length $\lambda_{\bot}$ of the squared velocity and magnetic field fluctuations ($u^2$ and $b^2$) transverse to the mean magnetic field inside a coronal hole (CH) at the base of the corona. The characteristic length scale defines directly the heating rate. We use a time series analysis of solar granulation and magnetic field measurements inside two CHs obtained with the New Solar Telescope (NST) at Big Bear Solar Observatory. A data set for transverse magnetic fields obtained with the Solar Optical Telescope/Spectro-Polarimeter (SOT/SP) aboard {\it Hinode} spacecraft was utilized to analyze the squared transverse magnetic field fluctuations $b_t^2$. Local correlation tracking (LCT) was applied to derive the squared transverse velocity fluctuations $u^2$. We find that for $u^2$-structures, Batchelor integral scale $\lambda$ varies in a range of 1800 - 2100 km, whereas the correlation length $\varsigma$ and the $e$-folding length $L$ vary between 660 and 1460 km. Structures for $b_t^2$ yield $\lambda \approx 1600$ km, $\varsigma \approx 640$ km, and $L \approx 620$ km. An averaged (over $\lambda, \varsigma$, and $L$) value of the characteristic length of $u^2$-fluctuations is 1260$\pm$500 km, and that of $b_t^2$ is 950$\pm$560 km. The characteristic length scale in the photosphere is approximately 1.5-50 times smaller than that adopted in previous models (3-30$\times10^3$ km). Our results provide a critical input parameter for current models of coronal heating and should yield an improved understanding of fast solar wind acceleration.

An abstract mathematical concept of fractal organization of certain complex objects received significant attention in astrophysics during last decades. The concept evolved into a broad field including multi-fractality and intermittency, percolation theory, self-organized criticality, theory of catastrophes, etc. Such a strong mathematical and physical approach provide new possibilities for exploring various aspects of astrophysics. In particular, in the solar and stellar magnetism, multi-fractal properties of magnetized plasma turned to be useful for understanding burst-like dynamics of energy release events, conditions for turbulent dynamo action, nature of turbulent magnetic diffusivity, and even the dual nature of solar dynamo. In this review, I will briefly outline how the ideas of multi-fractality are used to explore the above mentioned aspects of solar magnetism.

We used \ha-0.1 nm and magnetic field (at 1.56$\mu$) data obtained with the New Solar Telescope to study the origin of the disk counterparts to type II spicules, so-called rapid blueshifted excursions (RBEs). The high time cadence of our chromospheric (10 s) and magnetic field (45 s) data allowed us to generate x-t plots using slits parallel to the spines of the RBEs. These plots, along with potential field extrapolation, led us to suggest that the occurrence of RBEs is generally correlated with the appearance of new, mixed or unipolar fields in close proximity to network fields. RBEs show a tendency to occur at the interface between large-scale fields and small-scale dynamic magnetic loops and thus are likely to be associated with existence of a magnetic canopy. Detection of kinked and/or inverse "Y" shaped RBEs further confirm this conclusion.

Results of a statistical analysis of solar granulation are presented. A data set of 36 images of a quiet Sun area on the solar disk center was used. The data were obtained with the 1.6 m clear aperture New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO) and with a broad-band filter centered at the TiO (705.7 nm) spectral line. The very high spatial resolution of the data (diffraction limit of 77 km and pixel scale of 0.$"$0375) augmented by the very high image contrast (15.5$\pm$0.6%) allowed us to detect for the first time a distinct subpopulation of mini-granular structures. These structures are dominant on spatial scales below 600 km. Their size is distributed as a power law with an index of -1.8 (which is close to the Kolmogorov's -5/3 law) and no predominant scale. The regular granules display a Gaussian (normal) size distribution with a mean diameter of 1050 km. Mini-granular structures contribute significantly to the total granular area. They are predominantly confined to the wide dark lanes between regular granules and often form chains and clusters, but different from magnetic bright points. A multi-fractality test reveals that the structures smaller than 600 km represent a multi-fractal, whereas on larger scales the granulation pattern shows no multi-fractality and can be considered as a Gaussian random field. The origin, properties and role of the newly discovered population of mini-granular structures in the solar magneto-convection are yet to be explored.

We used photospheric intensity images and magnetic field measurements from the New Solar Telescope in Big Bear and Helioseismic Magnetic Imager on board Solar Dynamics Observatory (SDO) to study the the effect that the new small-scale emerging flux induces on solar granulation. We report that emerging flux appears to leave different types of footprint on solar granulation: i) diffuse irregular patches of increased brightness, ii) well defined filament-like structures and accompanied bright points, and iii) bright point-like features that appear inside granules. We suggest that the type of the footprint depends on the intensity of emerging fields. Stronger fields, emerging as a part of large magnetic structure, create on the solar surface a well defined filamentary pattern with bright points at the ends of the filaments, while weak turbulent fields are associated with bright patches inside the host granule.

Using high resolution off-band \ha\ data from the New Solar Telescope and Morlet wavelet analysis technique, we analyzed transverse motions of type II spicules observed near the North Pole of the Sun. Our new findings are that i) some of the observed type II spicules display kink or an inverse "Y" features, suggesting that their origin may be due to magnetic reconnection, and ii) type II spicules tend to display coherent transverse motions/oscillations. Also, the wavelet analysis detected significant presence of high frequency oscillations in type II spicules, ranging from 30 to 180 s with the the average period of 90 s. We conclude that at least some of type II spicules and their coherent transverse motions may be caused by reconnection between large scale fields rooted in the intergranular lanes and and small-scale emerging dipoles, a process that is know to generate high frequency kink mode MHD waves propagating along the magnetic field lines.

Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent properties we compare kinetic energy spectra on granular and sub-granular scales obtained from infrared TiO observations with the New Solar Telescope (Big Bear Solar Observatory) and from 3D radiative MHD numerical simulations ('SolarBox' code). We find that the numerical simulations require a high spatial resolution with 10 - 25 km grid-step in order to reproduce the inertial (Kolmogorov) turbulence range. The observational data require an averaging procedure to remove noise and potential instrumental artifacts. The resulting kinetic energy spectra show a good agreement between the simulations and observations, opening new perspectives for detailed joint analysis of more complex turbulent phenomena on the Sun, and possibly on other stars. In addition, using the simulations and observations we investigate effects of background magnetic field, which is concentrated in self-organized complicated structures in intergranular lanes, and find an increase of the small-scale turbulence energy and its decrease at larger scales due to magnetic field effects.

Generation and diffusion of the magnetic field on the Sun is a key mechanism responsible for solar activity on all spatial and temporal scales - from the solar cycle down to the evolution of small-scale magnetic elements in the quiet Sun. The solar dynamo operates as a non-linear dynamical process and is thought to be manifest in two types: as a global dynamo responsible for the solar cycle periodicity, and as a small-scale turbulent dynamo responsible for the formation of magnetic carpet in the quiet Sun. Numerous MHD simulations of the solar turbulence did not yet reach a consensus as to the existence of a turbulent dynamo on the Sun. At the same time, high-resolution observations of the quiet Sun from Hinode instruments suggest possibilities for the turbulent dynamo. Analysis of characteristics of turbulence derived from observations would be beneficial in tackling the problem. We analyse magnetic and velocity energy spectra as derived from Hinode/SOT, SOHO/MDI, SDO/HMI and the New Solar Telescope (NST) of Big Bear Solar Observatory (BBSO) to explore the possibilities for the small-scale turbulent dynamo in the quiet Sun.

On the basis of observations of solar granulation obtained with the New Solar Telescope (NST) of Big Bear Solar Observatory, we explored proper motion of bright points (BPs) in a quiet sun area, a coronal hole, and an active region plage. We automatically detected and traced bright points (BPs) and derived their mean-squared displacements as a function of time (starting from the appearance of each BP) for all available time intervals. In all three magnetic environments, we found the presence of a super-diffusion regime, which is the most pronounced inside the time interval of 10-300 seconds. Super-diffusion, measured via the spectral index, $\gamma$, which is the slope of the mean-squared displacement spectrum, increases from the plage area ($\gamma=1.48$) to the quiet sun area ($\gamma=1.53$) to the coronal hole ($\gamma=1.67$). We also found that the coefficient of turbulent diffusion changes in direct proportion to both temporal and spatial scales. For the minimum spatial scale (22 km) and minimum time scale (10 sec), it is 22 and 19 km$^{2}$ s$^{-1}$ for the coronal hole and the quiet sun area, respectively, whereas for the plage area it is about 12 km$^{2}$ s$^{-1}$ for the minimum time scale of 15 seconds. We applied our BP tracking code to 3D MHD model data of solar convection (Stein et al. 2007) and found the super-diffusion with $\gamma=1.45$. An expression for the turbulent diffusion coefficient as a function of scales and $\gamma$ is obtained.

Here we analyze solar activity by focusing on time variations of the number of sunspot groups (SGs) as a function of their modified Zurich class. We analyzed data for solar cycles 2023 by using Rome (cycles 2021) and Learmonth Solar Observatory (cycles 2223) SG numbers. All SGs recorded during these time intervals were separated into two groups. The first group includes small SGs (A, B, C, H, and J classes by Zurich classification) and the second group consists of large SGs (D, E, F, and G classes). We then calculated small and large SG numbers from their daily mean numbers as observed on the solar disk during a given month. We report that the time variations of small and large SG numbers are asymmetric except for the solar cycle 22. In general large SG numbers appear to reach their maximum in the middle of the solar cycle (phase 0.450.5), while the international sunspot numbers and the small SG numbers generally peak much earlier (solar cycle phase 0.290.35). Moreover, the 10.7 cm solar radio flux, the facular area, and the maximum CME speed show better agreement with the large SG numbers than they do with the small SG numbers. Our results suggest that the large SG numbers are more likely to shed light on solar activity and its geophysical implications. Our findings may also influence our understanding of long term variations of the total solar irradiance, which is thought to be an important factor in the Sun - Earth climate relationship.

We investigate the relationship between the monthly averaged maximal speeds of coronal mass ejections (CMEs), international sunspot number (ISSN), and the geomagnetic Dst and Ap indices covering the 1996-2008 time interval (solar cycle 23). Our new findings are as follows. (1) There is a noteworthy relationship between monthly averaged maximum CME speeds and sunspot numbers, Ap and Dst indices. Various peculiarities in the monthly Dst index are correlated better with the fine structures in the CME speed profile than that in the ISSN data. (2) Unlike the sunspot numbers, the CME speed index does not exhibit a double peak maximum. Instead, the CME speed profile peaks during the declining phase of solar cycle 23. Similar to the Ap index, both CME speed and the Dst indices lag behind the sunspot numbers by several months. (3) The CME number shows a double peak similar to that seen in the sunspot numbers. The CME occurrence rate remained very high even near the minimum of the solar cycle 23, when both the sunspot number and the CME average maximum speed were reaching their minimum values. (4) A well-defined peak of the Ap index between 2002 May and 2004 August was co-temporal with the excess of the mid-latitude coronal holes during solar cycle 23. The above findings suggest that the CME speed index may be a useful indicator of both solar and geomagnetic activities. It may have advantages over the sunspot numbers, because it better reflects the intensity of Earth-directed solar eruptions.

We studied bright umbral dots (UDs) detected in a moderate size sunspot and compared their statistical properties to recent MHD models. The study is based on high resolution data recorded by the New Solar Telescope at the Big Bear Solar Observatory and 3D MHD simulations of sunspots. Observed UDs, living longer than 150 s, were detected and tracked in a 46 min long data set, using an automatic detection code. Total 1553 (620) UDs were detected in the photospheric (low chromospheric) data. Our main findings are: i) none of the analyzed UDs is precisely circular, ii) the diameter-intensity relationship only holds in bright umbral areas, and iii) UD velocities are inversely related to their lifetime. While nearly all photospheric UDs can be identified in the low chromospheric images, some small closely spaced UDs appear in the low chromosphere as a single cluster. Slow moving and long living UDs seem to exist in both the low chromosphere and photosphere, while fast moving and short living UDs are mainly detected in the photospheric images. Comparison to the 3D MHD simulations showed that both types of UDs display, on average, very similar statistical characteristics. However, i) the average number of observed UDs per unit area is smaller than that of the model UDs, and ii) on average, the diameter of model UDs is slightly larger than that of observed ones.

We studied flux emergence events of sub-granular scale in a solar active region. New Solar Telescope (NST) of Big Bear Solar Observatory made it possible to clearly observe the photospheric signature of flux emergence with very high spatial (0".11 at 7057{\AA}) and temporal (15 s) resolution. From TiO observations with the pixel scale of 0".0375, we found several elongated granule-like features (GLFs) stretching from the penumbral filaments of a sunspot at a relatively high speed of over 4 km s-1. After a slender arched darkening appeared at a tip of a penumbral filament, a bright point (BP) developed and quickly moved away from the filament forming and stretching a GLF. The size of a GLF was approximately 0.5" wide and 3" long. The moving BP encountered nearby structures after several minutes of stretching, and a well-defined elongated shape of a GLF faded away. Magnetograms from SDO/HMI and NST/IRIM revealed that those GLFs are photospheric indicators of small-scale flux emergence, and their disappearance is related to magnetic cancellation. From two well-observed events, we describe detailed development of the sub-structures of GLFs, and different cancellation processes that each of the two GLFs underwent.

We observe that intergranular jets, originating in the intergranular space surrounding individual granules, tend to be associated with granular fragmentation, in particular, with the formation and evolution of a bright granular lane (BGL) within individual granules. The BGLs have recently been identified as vortex tubes by Steiner et al. We further discover the development of a well-defined bright grain located between the BGL and the dark intergranular lane to which it is connected. Signatures of a BGL may reach the lower chromosphere and can be detected in off-band \ha images. Simulations also indicate that vortex tubes are frequently associated with small-scale magnetic fields. We speculate that the intergranular jets detected in the NST data may result from the interaction between the turbulent small-scale fields associated with the vortex tube and the larger-scale fields existing in the intergranular lanes. The intergranular jets are much smaller and weaker than all previously known jet-like events. At the same time, they appear much more numerous than the larger events, leading us to the speculation that the total energy release and mass transport by these tiny events may not be negligible in the energy and mass-flux balance near the temperature minimum atop the photosphere. The study is based on the photospheric TiO broadband (1.0 nm) filter data acquired with the 1.6 m New Solar Telescope (NST) operating at the Big Bear Solar Observatory. The data set also includes NST off-band \ha images collected through a Zeiss Lyot filter with a passband of 0.025 nm.

We detected 2.8 bright points (BPs) per Mm$^2$ in the Quiet Sun (QS) with the New Solar Telescope (NST) at Big Bear Solar Observatory; using the TiO 705.68 nm spectral line, at an angular resolution ~ 0.1'' to obtain 30 min data sequence. Some BPs formed knots that were stable in time and influenced the properties of the granulation pattern around them. The observed granulation pattern within ~ 3'' of knots presents smaller granules than those observed in a normal granulation pattern; i.e., around the knots a suppressed convection is detected. Observed BPs covered ~ 5% of the solar surface and were not homogeneously distributed. BPs had an average size of 0.22'', they were detectable for 4.28 min in average, and had an averaged contrast of 0.1% in the deep red TiO spectral line.

With the ever increasing influx of high resolution images of the solar surface obtained at a multitude of wavelengths, various processes occurring at small spatial scales have become a greater focus of our attention. Complex small-scale magnetic fields have been reported that appear to have enough stored to heat the chromosphere. While significant progress has been made in understanding small-scale phenomena, many specifics remain elusive. We present here a detailed study of a single event of disappearance of a magnetic dipole and associated chromospheric activity. Based on New Solar Telescope H$\alpha$ data and {\it Hinode} photospheric line-of-sight magnetograms and Ca II H images we report the following. 1) Our analysis indicates that even very small dipoles (elements separated by about 0\arcsec.5 or less) may reach the chromosphere and trigger non-negligible chromospheric activity. 2) Careful consideration of the magnetic environment where the new flux is deposited may shed light on the details of magnetic flux removal from the solar surface. We argue that the apparent collision and disappearance of two opposite polarity elements may not necessarily indicate their cancellation (i.e., reconnection, emergence of a "U" tube or submergence of $ \Omega $ loops). In our case, the magnetic dipole disappeared by reconnecting with overlying large-scale inclined plage fields. 3) Bright points seen in off-band H$\alpha$ images are very well-correlated with the Ca II H bright points, which in turn are co-spatial with G-band bright points. We further speculate that, in general, H$\alpha$ bright points are expected be co-spatial with photospheric BPs, however, a direct comparison is needed to refine their relationship.

We present results of a study of intermittency and multifractality of magnetic structures in solar active regions (ARs). Line-of-sight magnetograms for 214 ARs of different flare productivity observed at the center of the solar disk from January 1997 until December 2006 are utilized. Data from the Michelson Doppler Imager (MDI) instrument on-board the {\it Solar and Heliospheric Observatory} (SOHO) operating in the high resolution mode, the Big Bear Solar Observatory digital magnetograph and {\it Hinode} SOT/SP instrument were used. Intermittency spectra were derived via high-order structure functions and flatness functions. The flatness function exponent is a measure of the degree of intermittency. We found that the flatness function exponent at scales below approximately 10 Mm is correlated to the flare productivity (the correlation coefficient is - 0.63). {\it Hinode} data show that the intermittency regime is extended toward the small scales (below 2 Mm) as compared to the MDI data. The spectra of multifractality, derived from the structure functions and flatness functions, are found to be more broad for ARs of highest flare productivity as compared to that of low flare productivity. The magnetic structure of high-flaring ARs consists of a voluminous set of monofractals, and this set is much richer than that for low-flaring ARs. The results indicate relevance of the multifractal organization of the photospheric magnetic fields to the flaring activity. Strong intermittency observed in complex and high-flaring ARs is a hint that we observe a photospheric imprint of enhanced sub-photospheric dynamics.

We present results of two-hour non-interrupted observations of solar granulation obtained under excellent seeing conditions with the largest aperture ground-based solar telescope - the New Solar Telescope (NST) - of Big Bear Solar Observatory. Observations were performed with adaptive optics correction using a broad-band TiO filter in the 705.7 nm spectral line with a time cadence of 10 s and a pixel size of 0.0375". Photospheric bright points (BPs) were detected and tracked. We find that the BPs detected in NST images are co-spatial with those visible in Hinode/SOT G-band images. In cases where Hinode/SOT detects one large BP, NST detects several separated BPs. Extended filigree features are clearly fragmented into separate BPs in NST images. The distribution function of BP sizes extends to the diffraction limit of NST (77 km) without saturation and corresponds to a log-normal distribution. The lifetime distribution function follows a log-normal approximation for all BPs with lifetime exceeding 100 s. A majority of BPs are transient events reflecting the strong dynamics of the quiet sun: 98.6% of BPs live less than 120 s. The longest registered life time was 44 minutes. The size and maximum intensity of BPs were found to be proportional to their life times.

Line-of-sight magnetograms for 217 active regions (ARs) of different flare rate observed at the solar disk center from January 1997 until December 2006 are utilized to study the turbulence regime and its relationship to the flare productivity. Data from {\it SOHO}/MDI instrument recorded in the high resolution mode and data from the BBSO magnetograph were used. The turbulence regime was probed via magnetic energy spectra and magnetic dissipation spectra. We found steeper energy spectra for ARs of higher flare productivity. We also report that both the power index, $\alpha$, of the energy spectrum, $E(k) \sim k^{-\alpha}$, and the total spectral energy $W=\int E(k)dk$ are comparably correlated with the flare index, $A$, of an active region. The correlations are found to be stronger than that found between the flare index and total unsigned flux. The flare index for an AR can be estimated based on measurements of $\alpha$ and $W$ as $A=10^b (\alpha W)^c$, with $b=-7.92 \pm 0.58$ and $c=1.85 \pm 0.13$. We found that the regime of the fully-developed turbulence occurs in decaying ARs and in emerging ARs (at the very early stage of emergence). Well-developed ARs display under-developed turbulence with strong magnetic dissipation at all scales.

Surface photometry of the quiet Sun has achieved an angular resolution of $0".1$ with the New Solar Telescope at Big Bear Solar Observatory revealing that a disproportionate fraction of the oscillatory events appear above observed bright point-like structures. During the tracking of these structures, we noted that the more powerful oscillatory events are cospatial with them, indicating that observed flux tubes may be the source of many observed oscillatory events.

Low and mid-latitude coronal holes (CHs) observed on the Sun during the current solar activity minimum (from September 21, 2006, Carrington rotation (CR) 2048, until June 26, 2009 (CR 2084)) were analyzed using {\it SOHO}/EIT and STEREO-A SECCHI EUVI data. From both the observations and Potential Field Source Surface (PFSS) modeling, we find that the area occupied by CHs inside a belt of $\pm 40^\circ$ around the solar equator is larger in the current 2007 solar minimum relative to the similar phase of the previous 1996 solar minimum. The enhanced CH area is related to a recurrent appearance of five persistent CHs, which survived during 7-27 solar rotations. Three of the CHs are of positive magnetic polarity and two are negative. The most long-lived CH was being formed during 2 days and existed for 27 rotations. This CH was associated with fast solar wind at 1 AU of approximately 620$\pm 40$ km s$^{-1}$. The 3D MHD modeling for this time period shows an open field structure above this CH. We conclude that the global magnetic field of the Sun possessed a multi-pole structure during this time period. Calculation of the harmonic power spectrum of the solar magnetic field demonstrates a greater prevalence of multi-pole components over the dipole component in the 2007 solar minimum compared to the 1996 solar minimum. The unusual large separation between the dipole and multi-pole components is due to the very low magnitude of the dipole component, which is three times lower than that in the previous 1996 solar minimum.

Parameters of magnetic flux distribution inside low-latitude coronal holes (CHs) were analyzed. A statistical study of 44 CHs based on Solar and Heliospheric Observatory (SOHO)/MDI full disk magnetograms and SOHO/EIT 284\AA images showed that the density of the net magnetic flux, $B_{{\rm net}}$, does not correlate with the associated solar wind speeds, $V_x$. Both the area and net flux of CHs correlate with the solar wind speed and the corresponding spatial Pearson correlation coefficients are 0.75 and 0.71, respectively. A possible explanation for the low correlation between $B_{{\rm net}}$ and $V_x$ is proposed. The observed non-correlation might be rooted in the structural complexity of the magnetic field. As a measure of complexity of the magnetic field, the filling factor, $ f(r)$, was calculated as a function of spatial scales. In CHs, $f(r)$ was found to be nearly constant at scales above 2 Mm, which indicates a monofractal structural organization and smooth temporal evolution. The magnitude of the filling factor is 0.04 from the Hinode SOT/SP data and 0.07 from the MDI/HR data. The Hinode data show that at scales smaller than 2 Mm, the filling factor decreases rapidly, which means a mutlifractal structure and highly intermittent, burst-like energy release regime. The absence of necessary complexity in CH magnetic fields at scales above 2 Mm seems to be the most plausible reason why the net magnetic flux density does not seem to be related to the solar wind speed: the energy release dynamics, needed for solar wind acceleration, appears to occur at small scales below 1 Mm.

Recent studies undoubtedly demonstrate that the magnetic field in the photosphere and corona is an intermittent structure, which offers new views on the underlying physics. In particular, such problems as the existence in the corona of localized areas with extremely strong resistivity (required to explain magnetic reconnection of all scales) and the interchange between small and large scales (required in study of the photosphere/corona coupling), to name a few, can be easily captured by the concept of intermittency. This study is focused on simultaneous time variations of intermittency properties derived in the photosphere, chromosphere and corona. We analyzed data for NOAA AR 10930 acquired between Dec 08, 2006 12:00 UT and Dec 13, 2006 18:45 UT. Photospheric intermittency was inferred from Hinode magnetic field measurements, while intermittency in the transition region and corona was derived from Nobeyama 9 GHz radio polarization measurements, high cadence Hinode/XRT/Be-thin data as well as GOES 1-8\AA flux. Photospheric dynamics and its possible relationship with the intermittency variations were also analyzed by calculating the kinetic vorticity. For this case study we found the following chain of events. Intermittency of the photospheric magnetic field peaked after the specific kinetic vorticity of plasma flows in the AR reached its maximum level (4 hour time delay). In turn, gradual increase of coronal intermittency occurred after the peak of the photospheric intermittency. The time delay between the peak of photospheric intermittency and the occurrence of the first strong (X3.4) flare was approximately 1.3 days. Our analysis seems to suggest that the enhancement of intermittency/complexity first occurs in the photosphere and is later transported toward the corona.

A proxy for horizontal electric currents in the solar photosphere was suggested. For a set of evolving active regions (ARs) observed with {\it Solar and Heliospheric Observatory (SOHO)} Michelson Doppler Imager (MDI) in the high resolution mode, the dissipation spectrum, $k^2E(k)$, and the spatial structure of dissipation, i.e., the Stokes dissipation function $\epsilon(x,y)$, were calculated from the observed $B_z$ component of the magnetic field. These functions allowed us to calculate (a part of) the horizontal electric current density in the photosphere. It was shown that as an active region emerges, large-scale horizontal electric currents are gradually generated and determine a bulk of dissipation. When an active region decays, the large-scale horizontal currents decay faster than the small-scale ones. The density of horizontal currents in active regions is in the range of $<j_h > \sim (0.008 - 0.028)$ A/m$^2$, that is compatible with the density of vertical currents in active regions. We suggest two possible mechanisms for generation of such horizontal currents in the photosphere. One of them is the drift motions of charged particles in the medium of varying plasma pressure gradient in a horizontal plane at the periphery of a sunspot. Such a drift can produce quasi-circular closed horizontal currents around sunspots. Another possibility could be an existence of horizontal axial current inside a highly twisted horizontal magnetic structure laying in the photosphere along the magnetic neutral line. The horizontal currents may contribute significantly to the dynamics of the photosphere/corona coupling, as well as the estimation of non-potentiality of ARs.

The observations of magnetic field variations as a signature of flaring activity is one of the main goal in solar physics. Some efforts in the past give apparently no unambiguous observations of changes. We observed that the scaling laws of the current helicity inside a given flaring active region change clearly and abruptly in correspondence with the eruption of big flares at the top of that active region. Comparison with numerical simulations of MHD equations, indicates that the change of scaling behavior in the current helicity, seems to be associated to a topological reorganization of the footpoint of the magnetic field loop, namely to dissipation of small scales structures in turbulence. It is evident that the possibility of forecasting in real time high energy flares, even if partially, has a wide practical interest to prevent the effects of big flares on Earth and its environment.