Observations of the solar photosphere show spatially-compact large-amplitude Doppler velocity events with short lifetimes. In data from the Imaging Magnetograph eXperiment (IMaX) on the first flight of the Sunrise balloon in 2009, events with velocities in excess of 4{\sigma} from the mean can be identified in both intergranular downflow lanes and granular upflows. We show that the statistics of such events are consistent with the random superposition of strong convective flows and p-mode coherence patches. Such coincident superposition has implications for the interpretation of acoustic wave sources in the solar photosphere, and may be important to the interpretation of spectral line profiles formed in solar photosphere.

Radiative transfer coupled with highly realistic simulations of the solar atmosphere is routinely used to infer the physical properties underlying solar observations. Due to its computational efficiency, the method of short-characteristics is often employed, despite it introducing numerical diffusion as an interpolation artifact. In this paper, we quantify the effect of the numerical diffusion on the spatial resolution of synthesize emergent intensity images, and derive a closed form analytical model of the diffusion error as a function of viewing angle when using linear interpolation. We demonstrate that the image degradation adversely affects the comparison between simulated data and observations, for observations away from disk-center, unless the simulations are computed at much higher intrinsic resolution than the observations. We also show that the diffusion error is readily avoided by interpolating the simulation solution on a viewing-angle aligned grid prior to computing the radiative transfer. Doing this will be critical for comparisons with observations using the upcoming large aperture telescopes --- the Daniel K. Inouye Solar Telescope and the European Solar Telescope.

The motions of small-scale magnetic flux elements in the solar photosphere can provide some measure of the Lagrangian properties of the convective flow. Measurements of these motions have been critical in estimating the turbulent diffusion coefficient in flux-transport dynamo models and in determining the Alfven wave excitation spectrum for coronal heating models. We examine the motions of internetwork flux elements in a 24 hour long Hinode/NFI magnetogram sequence with 90 second cadence, and study both the scaling of their mean squared displacement and the shape of their displacement probability distribution as a function of time. We find that the mean squared displacement scales super-diffusively with a slope of about 1.48. Super-diffusive scaling has been observed in other studies for temporal increments as small as 5 seconds, increments over which ballistic scaling would be expected. Using high-cadence MURaM simulations, we show that the observed super-diffusive scaling at short temporal increments is a consequence of random changes in the barycenter positions due to flux evolution. We also find that for long temporal increments, beyond granular lifetimes, the observed displacement distribution deviates from that expected for a diffusive process, evolving from Rayleigh to Gaussian. This change in the distribution can be modeled analytically by accounting for supergranular advection along with motions due to granulation. These results complicate the interpretation of magnetic element motions as strictly advective or diffusive on short and long timescales and suggest that measurements of magnetic element motions must be used with caution in turbulent diffusion or wave excitation models. We propose that passive trace motions in measured photospheric flows may yield more robust transport statistics.

Supergranulation is characterized by horizontally divergent flows with typical length scales of 32 Mm in the solar photosphere. Unlike granulation, the size of which is comparable to both the thickness of the radiative boundary layer and local scale height of the plasma in the photosphere, supergranulation does not reflect any obvious length scale of the solar convection zone. Early suggestions that the depth of second helium ionization is important are not supported by numerical simulations. Thus the origin of the solar supergranulation remains largely a mystery. Moreover, observations of flows in the photosphere using either Doppler imaging or correlation or feature tracking show a monotonic decrease in power at scales larger than supergranulation. Both local area and global spherical shell simulations of solar convection by contrast show the opposite, a power law increase in horizontal flow amplitudes to low wavenumber. Here we examine this disparity, and investigate how the solar supergranulation may arise as a consequence of strong photospheric driving and non-local heat transport by cool diving plumes. Using three dimensional anelastic simulations with surface driving, we show that the kinetic energy of largest convective scales in the upper layers of a stratified domain reflects the depth of transition from strong buoyant driving to adiabatic stratification below. We interpret the observed monotonic decrease in solar convective power at scales larger than supergranulation to be a consequence of this rapid transition, and show how the supergranular scale can be understood as the largest buoyantly driven mode of convection in the Sun.

Solar irradiance variations over solar rotational time-scales are largely determined by the passage of magnetic structures across the visible solar disk. Variations on solar cycle time scales are thought to be similarly due to changes in surface magnetism with activity. Understanding the contribution of magnetic structures to total solar irradiance and solar spectral irradiance requires assessing their contributions as a function of disk position. Since only relative photometry is possible from the ground, the contrasts of image pixels are measured with respect to a center-to-limb intensity profile. Using nine years of full-disk red and blue continuum images from the Precision Solar Photometric Telescope at the Mauna Loa Solar Observatory (PSPT/MLSO), we examine the sensitivity of continuum contrast measurements to the center-to-limb profile definition. Profiles which differ only by the amount of magnetic activity allowed in the pixels used to determine them yield oppositely signed solar cycle length continuum contrast trends; either agreeing with the result of Preminger et al. (2011) showing negative correlation with solar cycle or disagreeing and showing positive correlation with solar cycle. Changes in the center-to-limb profile shape over the solar cycle are responsible for the contradictory contrast results, and we demonstrate that the lowest contrast structures, internetwork and network, are most sensitive to these. Thus the strengths of the full-disk, internetwork, and network photometric trends depend critically on the magnetic flux density used in the quiet-sun definition. We conclude that the contributions of low contrast magnetic structures to variations in the solar continuum output, particularly to long-term variations, are difficult, if not impossible, to determine without the use of radiometric imaging.

We model the solar horizontal velocity power spectrum at scales larger than granulation using a two-component approximation to the mass continuity equation. The model takes four times the density scale height as the integral (driving) scale of the vertical motions at each depth. Scales larger than this decay with height from the deeper layers. Those smaller are assumed to follow a Kolomogorov turbulent cascade, with the total power in the vertical convective motions matching that required to transport the solar luminosity in a mixing length formulation. These model components are validated using large scale radiative hydrodynamic simulations. We reach two primary conclusions: 1. The model predicts significantly more power at low wavenumbers than is observed in the solar photospheric horizontal velocity spectrum. 2. Ionization plays a minor role in shaping the observed solar velocity spectrum by reducing convective amplitudes in the regions of partial helium ionization. The excess low wavenumber power is also seen in the fully nonlinear three-dimensional radiative hydrodynamic simulations employing a realistic equation of state. This adds to other recent evidence suggesting that the amplitudes of large scale convective motions in the Sun are significantly lower than expected. Employing the same feature tracking algorithm used with observational data on the simulation output, we show that the observed low wavenumber power can be reproduced in hydrodynamic models if the amplitudes of large scale modes in the deep layers are artificially reduced. Since the large scale modes have reduced amplitudes, modes on the scale of supergranulation and smaller remain important to convective heat flux even in the deep layers, suggesting that small scale convective correlations are maintained through the bulk of the solar convection zone.

We present the initial results of an observational study into the variation of the dominant length-scale of quiet solar emission: supergranulation. This length-scale reflects the radiative energy in the plasma of the upper solar chromosphere and transition region at the magnetic network boundaries forming as a result of the relentless interaction of magnetic fields and convective motions of the Sun's interior. We demonstrate that a net difference of ~0.5Mm in the supergranular emission length-scale occurs when comparing observations cycle 22/23 and cycle 23/24 minima. This variation in scale is reproduced in the datasets of multiple space- and ground-based instruments and using different diagnostic measures. By means of extension, we consider the variation of the supergranular length-scale over multiple solar minima by analyzing a subset of the Mt Wilson Solar Observatory (MWO) Ca II K image record. The observations and analysis presented provide a tantalizing look at solar activity in the absence of large-scale flux emergence, offering insight into times of "extreme" solar minimum and general behavior such as the phasing and cross-dependence of different components of the spectral irradiance. Given that the modulation of the supergranular scale imprints itself in variations of the Sun's spectral irradiance, as well as in the mass and energy transport into the entire outer atmosphere, this preliminary investigation is an important step in understanding the impact of the quiet sun on the heliospheric system.

As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance of the flows established in rotating spherical shell segments that capture much of the NSSL with the Curved Spherical Segment (CSS) code using two different diffusion schemes. The CSS code is a new massively parallel modeling tool capable of simulating 3-D compressible MHD convection with a realistic solar stratification in rotating spherical shell segments.

We study the radiative emission of various types of solar features, such as quiet Sun, enhanced network, plage, and bright plage regions, identified on filtergrams taken in the Ca II K line. We analysed fulldisk images obtained with the PSPT, by using three interference filters that sample the Ca II K line with different bandpasses. We studied the dependence of the radiative emission of disk features on the filter bandpass. We also performed a NLTE spectral synthesis of the Ca II K line integrated over the bandpass of PSPT filters. The synthesis was carried out by utilizing both the PRD and CRD with the most recent set of semi empirical atmosphere models in the literature and some earlier atmosphere models. We measured the CLV of intensity values for various solar features identified on PSPT images and compared the results obtained with those derived from the synthesis. We find that CRD calculations derived using the most recent quiet Sun model, on average, reproduce the measured values of the quiet Sun regions slightly more accurately than PRD computations with the same model. This may reflect that the utilized atmospheric model was computed assuming CRD. Calculations with PRD on earlier quiet Sun model atmospheres reproduce measured quantities with a similar accuracy as to that achieved here by applying CRD to the recent model. We also find that the median contrast values measured for most of the identified bright features, disk positions, and filter widths are, on average, a factor 1.9 lower than those derived from PRD simulations performed using the recent bright feature models. The discrepancy between measured and modeled values decreases by 12% after taking into account straylight effects on PSPT images. PRD computations on either the most recent or the earlier atmosphere models of bright features reproduce measurements from plage and bright plage regions with a similar accuracy.

We investigated the radiative emission of different types of solar features in the spectral range of the Ca II K line. We analyzed full-disk 2k x 2k observations from the PSPT Precision Solar Photometric Telescope. The data were obtained by using three narrow-band interference filters that sample the Ca II K line with different pass bands. Two filters are centered in the line core, the other in the red wing of the line. We measured the intensity and contrast of various solar features, specifically quiet Sun (inter-network), network, enhanced network, plage, and bright plage (facula) regions. Moreover, we compared the results obtained with those derived from the numerical synthesis performed for the three PSPT filters with a widely used radiative code on a set of reference semi-empirical atmosphere models.

We have measured the average radial (cell center to network boundary) profile of the continuum intensity contrast associated with supergranular flows using data from the Precision Solar Photometric Telescope (PSPT) at the Mauna Loa Solar Observatory (MLSO). After removing the contribution of the network flux elements by the application of masks based on Ca II K intensity and averaging over more than 10^5 supergranular cells, we find a ~ 0.1% decrease in red and blue continuum intensity from the supergranular cell centers outward, corresponding to a ~ 1.0 K decrease in brightness temperature across the cells. The radial intensity profile may be caused either by the thermal signal associated with the supergranular flows or a variation in the packing density of unresolved magnetic flux elements. These are not unambiguously distinguished by the observations, and we raise the possibility that the network magnetic fields play an active role in supergranular scale selection by enhancing the radiative cooling of the deep photosphere at the cell boundaries.

We investigate the photometric signature of magnetic flux tubes in the solar photosphere. We developed two dimensional, static numerical models of isolated and clustered magnetic flux tubes. We investigated the emergent intensity profiles at different lines-of-sight for various spatial resolutions and opacity models. We found that both geometric and photometric properties of bright magnetic features are determined not only by the physical properties of the tube and its surroundings, but also by the particularities of the observations, including the line/continuum formation height, the spatial resolution and the image analyses techniques applied. We show that some observational results presented in the literature can be interpreted by considering bright magnetic features to be clusters of smaller elements, rather than a monolithic flux tube.

We have examined images from the Precision Solar Photometric Telescope (PSPT) at the Mauna Loa Solar Observatory (MLSO) in search of latitudinal variation in the solar photospheric intensity. Along with the expected brightening of the solar activity belts, we have found a weak enhancement of the mean continuum intensity at polar latitudes (continuum intensity enhancement $\sim0.1 - 0.2%$ corresponding to a brightness temperature enhancement of $\sim2.5{\rm K}$). This appears to be thermal in origin and not due to a polar accumulation of weak magnetic elements, with both the continuum and CaIIK intensity distributions shifted towards higher values with little change in shape from their mid-latitude distributions. Since the enhancement is of low spatial frequency and of very small amplitude it is difficult to separate from systematic instrumental and processing errors. We provide a thorough discussion of these and conclude that the measurement captures real solar latitudinal intensity variations.

Several studies have investigated the fractal and multifractal nature of magnetic features in the solar photosphere and its variation with the solar magnetic activity cycle. Here we extend those studies by examining the fractal geometry of bright magnetic features at higher atmospheric levels, specifically in the solar chromosphere. We analyze structures identified in CaIIK images obtained with the Precision Solar Photometric Telescopes (PSPTs) at Osservatorio Astronomico di Roma (OAR) and Mauna Loa Solar Observatory (MLSO). Fractal dimension estimates depend on the estimator employed, the quality of the images, and the structure identification techniques used. We examine both real and simulated data and employ two different perimeter-area estimators in order to understand the sensitivity of the deduced fractal properties to pixelization and image quality. The fractal dimension of bright 'magnetic' features in CaIIK images ranges between values of 1.2 and 1.7 for small and large structures respectively. This size dependency largely reflects the importance of image pixelization in the measurement of small objects. The fractal dimension of chromospheric features does not show any clear systematic variation with time over the period examined, the descending phase of solar cycle 23. These conclusions, and the analysis of both real and synthetic images on which they are based, are important in the interpretation of previously reported results.