Search

We present the first quiet Sun spectropolarimetric observations obtained with the Visible SpectroPolarimeter (ViSP) at the $4-$m Daniel K. Inouye Solar Telescope (DKIST). We recorded observations in a wavelength range that includes the magnetically sensitive Fe I $6301.5/6302.5$ $\AA$ doublet. With an estimated spatial resolution of 0.08'', this represents the highest spatial resolution full-vector spectropolarimetric observations ever obtained of the quiet Sun. We identified $53$ small-scale magnetic elements, including $47$ magnetic loops and $4$ unipolar magnetic patches, with linear and circular polarisation detected in all of them. Of particular interest is a magnetic element in which the polarity of the magnetic vector appears to change three times in only $400$ km and which has linear polarisation signals throughout. We find complex Stokes $V$ profiles at the polarity inversion lines of magnetic loops and discover degenerate solutions, as we are unable to conclusively determine whether these arise due to gradients in the atmospheric parameters or smearing of opposite polarity signals. We analyse a granule which notably has linear and circular polarisation signals throughout, providing an opportunity to explore its magnetic properties. On this small scale we see the magnetic field strength range from $25$ G at the granular boundary to $2$ kG in the intergranular lane (IGL), and sanity check the values with the weak and strong field approximations. A value of $2$ kG in the IGL is among the highest measurements ever recorded for the internetwork., Comment: Accepted for publication in ApJL

The solar physics community is entering a golden era that is ripe with next-generation ground- and space-based facilities. With ever-increasing resolving power stemming from the newest observational telescopes, it becomes more challenging to obtain (near-)simultaneous measurements at high spatial, temporal and spectral resolutions, while operating at the diffraction limit of these new facilities. Hence, in recent years there has been increased interest in the capabilities integral field units (IFUs) offer towards obtaining the trifecta of high spatial, temporal and spectral resolutions contemporaneously. To date, IFUs developed for solar physics research have focused on mid-optical and infrared measurements. Here, we present an IFU prototype that has been designed for operation within the near-ultraviolet to mid-optical wavelength range, hence providing additional spectral coverage to the instrument suites developed to date. The IFU was constructed as a low-budget proof-of-concept for the upcoming 2m class Indian National Large Solar Telescope and employs circular cross-section fibres to guide light into a Czerny-Turner configuration spectrograph, with the resulting spectra captured using a high quantum efficiency scientific CMOS camera. Mapping of each input fibre allows for the reconstruction of two-dimensional spectral images, with frame rates exceeding 20 per second possible while operating in a non-polarimetric configuration. The science verification data presented here highlights the suitability of fibre-fed IFUs operating at near-ultraviolet wavelengths for solar physics research. Importantly, the successful demonstration of this type of instrument paves the way for further technological developments to make a future variant suitable for upcoming ground-based and space-borne telescope facilities., Comment: 53 pages, 16 figures, accepted for publication by Solar Physics

Using high cadence observations from the Hydrogen-alpha Rapid Dynamics camera imaging system on the Dunn Solar Telescope, we present an investigation of the statistical properties of transverse oscillations in spicules captured above the solar limb. At five equally separated atmospheric heights, spanning approximately 4900-7500 km, we have detected a total of 15 959 individual wave events, with a mean displacement amplitude of 151 +/- 124 km, a mean period of 54 +/- 45 s, and a mean projected velocity amplitude of 21 +/- 13 km s^-1. We find that both the displacement and velocity amplitudes increase with height above the solar limb, ranging from 132 +/- 111 km and 17.7 +/- 10.6 km s^-1 at 4900 km, and 168 +/- 125 km and 26.3 +/- 14.1 km s^-1 at 7500 km, respectively. Following the examination of neighboring oscillations in time and space, we find 45% of the waves to be upwardly propagating, 49% to be downwardly propagating, and 6% to be standing, with mean absolute phase velocities for the propagating waves on the order of 75-150 km s^-1. While the energy flux of the waves propagating downwards does not appear to depend on height, we find the energy flux of the upwardly propagating waves decreases with atmospheric height at a rate of -13 200 +/- 6500 W m^-2 /Mm. As a result, this decrease in energy flux as the waves propagate upwards may provide significant thermal input into the local plasma., Comment: Accepted for publication in ApJ with 15 pages and 9 figures

Solar pores are efficient magnetic conduits for propagating magnetohydrodynamic wave energy into the outer regions of the solar atmosphere. Pore observations often contain isolated and/or unconnected structures, preventing the statistical examination of wave activity as a function of atmospheric height. Here, using high resolution observations acquired by the Dunn Solar Telescope, we examine photospheric and chromospheric wave signatures from a unique collection of magnetic pores originating from the same decaying sunspot. Wavelet analysis of high cadence photospheric imaging reveals the ubiquitous presence of slow sausage mode oscillations, coherent across all photospheric pores through comparisons of intensity and area fluctuations, producing statistically significant in-phase relationships. The universal nature of these waves allowed an investigation of whether the wave activity remained coherent as they propagate. Utilizing bi-sector Doppler velocity analysis of the Ca II 8542 {\AA} line, alongside comparisons of the modeled spectral response function, we find fine-scale 5 mHz power amplification as the waves propagate into the chromosphere. Phase angles approaching zero degrees between co-spatial bi-sectors spanning different line depths indicate standing sausage modes following reflection against the transition region boundary. Fourier analysis of chromospheric velocities between neighboring pores reveals the annihilation of the wave coherency observed in the photosphere, with examination of the intensity and velocity signals from individual pores indicating they behave as fractured wave guides, rather than monolithic structures. Importantly, this work highlights that wave morphology with atmospheric height is highly complex, with vast differences observed at chromospheric layers, despite equivalent wave modes being introduced into similar pores in the photosphere., Comment: Accepted for publication in ApJ

Using high cadence observations from the Hydrogen-alpha Rapid Dynamics camera imaging system on the Dunn Solar Telescope, we present an investigation of the statistical properties of transverse oscillations in spicules captured above the solar limb. At five equally separated atmospheric heights, spanning approximately 4900-7500 km, we have detected a total of 15 959 individual wave events, with a mean displacement amplitude of 151 +/- 124 km, a mean period of 54 +/- 45 s, and a mean projected velocity amplitude of 21 +/- 13 km s^-1. We find that both the displacement and velocity amplitudes increase with height above the solar limb, ranging from 132 +/- 111 km and 17.7 +/- 10.6 km s^-1 at 4900 km, and 168 +/- 125 km and 26.3 +/- 14.1 km s^-1 at 7500 km, respectively. Following the examination of neighboring oscillations in time and space, we find 45% of the waves to be upwardly propagating, 49% to be downwardly propagating, and 6% to be standing, with mean absolute phase velocities for the propagating waves on the order of 75-150 km s^-1. While the energy flux of the waves propagating downwards does not appear to depend on height, we find the energy flux of the upwardly propagating waves decreases with atmospheric height at a rate of -13 200 +/- 6500 W m^-2 /Mm. As a result, this decrease in energy flux as the waves propagate upwards may provide significant thermal input into the local plasma., Comment: Accepted for publication in ApJ with 15 pages and 9 figures

This deposit contains the LPJmL model simulation outputs of mean monthly root-zone soil moisture ("LPJmL_rzsm_hist.zip") and the data underlying the plot in Fig. 3 ("Fig_3_plotdata.xlsx") in the article: Wang-Erlandsson, L., Tobian, A., van der Ent, R. J., Fetzer, I., te Wierik, S., Porkka, M., Staal, A., Jaramillo, F., Dahlmann, H., Singh, C., Greve, P., Gerten, D., Keys, P.W., Gleeson, T, Cornell, S. E., Steffen, W., Bai, X., Rockström, J., (2022): A planetary boundary for green water. Nature Reviews Earth & Environment. For method description, please refer to the article. Files in the folder "LPJmL_rzsm_hist.zip" are named as NN_rzsm_hist_M.tif, in which: NN refers to the name of an Earth system model, of which the outputs were used as forcing in the LPJmL runs rzsm refers to "root zone soil moisture" hist refers to historical period 1850-2014 M refers to name of month (jan for January etc.)

With the advent of next generation high resolution telescopes, our understanding of how the magnetic field is organized in the internetwork (IN) photosphere is likely to advance.We aim to evaluate the extent to which we can retrieve information about the magnetic vector in the IN photosphere using inversions. We use snapshots produced from high resolution 3D magnetohydrodynamic (MHD) simulations and employ the Stokes Inversions based on Response functions (SIR) code to produce synthetic observables in the near infrared spectral window observed by the GREGOR Infrared Spectrograph (GRIS), which contains the highly magnetically sensitive photospheric Fe I line pair at 15648.52 A and 15652.87 A. We perform nearly 14 million inversions to test how well the true MHD atmospheric parameters can be constrained. Finally, we degrade the synthetic Stokes vectors spectrally and spatially to GREGOR resolutions and examine how this influences observations, considering the impact of stray light, spatial resolution and signal-to-noise (S-to-N). We find the depth-averaged parameters can be recovered by the inversions of the undegraded profiles, and by adding gradients to magnetic field strength, inclination and line of sight velocity we show an improvement in the chi squared value is achieved. We evaluate the extent to which we can constrain these parameters at various optical depths, with the kinematic and thermodynamic parameters sensitive deeper in the atmosphere than the magnetic parameters. We find the S-to-N and spatial resolution play a significant role in determining how the atmosphere appears and the magnetic and kinematic parameters are invariant upon inclusion of unpolarized stray light. We studied a linear polarization feature which resembles those recently observed by GRIS, appearing as loop-like structures with similar magnetic flux density., Comment: 13 pages, 11 figures, accepted for publication in A&A. Abstract abridged for mailing list

With the advent of next generation high resolution telescopes, our understanding of how the magnetic field is organized in the internetwork (IN) photosphere is likely to advance.We aim to evaluate the extent to which we can retrieve information about the magnetic vector in the IN photosphere using inversions. We use snapshots produced from high resolution 3D magnetohydrodynamic (MHD) simulations and employ the Stokes Inversions based on Response functions (SIR) code to produce synthetic observables in the near infrared spectral window observed by the GREGOR Infrared Spectrograph (GRIS), which contains the highly magnetically sensitive photospheric Fe I line pair at 15648.52 A and 15652.87 A. We perform nearly 14 million inversions to test how well the true MHD atmospheric parameters can be constrained. Finally, we degrade the synthetic Stokes vectors spectrally and spatially to GREGOR resolutions and examine how this influences observations, considering the impact of stray light, spatial resolution and signal-to-noise (S-to-N). We find the depth-averaged parameters can be recovered by the inversions of the undegraded profiles, and by adding gradients to magnetic field strength, inclination and line of sight velocity we show an improvement in the chi squared value is achieved. We evaluate the extent to which we can constrain these parameters at various optical depths, with the kinematic and thermodynamic parameters sensitive deeper in the atmosphere than the magnetic parameters. We find the S-to-N and spatial resolution play a significant role in determining how the atmosphere appears and the magnetic and kinematic parameters are invariant upon inclusion of unpolarized stray light. We studied a linear polarization feature which resembles those recently observed by GRIS, appearing as loop-like structures with similar magnetic flux density., Comment: 13 pages, 11 figures, accepted for publication in A&A. Abstract abridged for mailing list

The umbral regions of sunspots and pores in the solar photosphere are generally dominated by $3$~mHz oscillations, which are due to $p$-modes penetrating the magnetic region. In these locations, wave power is also significantly reduced with respect to the quiet Sun. However, here we study a pore where the power of the oscillations in the umbra is not only comparable, or even larger than that of the quiet Sun, but the main dominant frequency is not $3$~mHz as expected, but instead $5$~mHz. By combining Doppler velocities and spectropolarimetry and analysing the relationship between magnetic field strength and frequency, the resultant $B-\omega$ diagram reveals distinct ridges which are remarkably clear signatures of resonant MHD oscillations confined within the pore umbra. In addition to velocity oscillations, we demonstrate that these modes are also accompanied by magnetic oscillations, as predicted from MHD theory. The novel technique of $B-\omega$ analysis, proposed in this Letter, opens an exciting new avenue for identifying MHD wave modes in the umbral regions of both pores and sunspots., Comment: Accepted in A&A

The formation of shocks within the solar atmosphere remains one of the few observable signatures of energy dissipation arising from the plethora of magnetohydrodynamic waves generated close to the solar surface. Active region observations offer exceptional views of wave behavior and its impact on the surrounding atmosphere. The stratified plasma gradients present in the lower solar atmosphere allow for the potential formation of many theorized shock phenomena. In this study, using chromospheric Ca II 854.2nm spectropolarimetric data of a large sunspot, we examine fluctuations in the plasma parameters in the aftermath of powerful shock events that demonstrate polarimetric reversals during their evolution. Modern inversion techniques are employed to uncover perturbations in the temperatures, line-of-sight velocities, and vector magnetic fields occurring across a range of optical depths synonymous with the shock formation. Classification of these non-linear signatures is carried out by comparing the observationally-derived slow, fast, and Alfv\'en shock solutions to the theoretical Rankine-Hugoniot relations. Employing over 200,000 independent measurements, we reveal that the Alfv\'en (intermediate) shock solution provides the closest match between theory and observations at optical depths of log(tau) = -4, consistent with a geometric height at the boundary between the upper photosphere and lower chromosphere. This work uncovers first-time evidence of the manifestation of chromospheric intermediate shocks in sunspot umbrae, providing a new method for the potential thermalization of wave energy in a range of magnetic structures, including pores, magnetic flux ropes, and magnetic bright points., Comment: 17 pages, 9 figures, accepted by ApJ

Stokes inversion codes are crucial in returning properties of the solar atmosphere, such as temperature and magnetic field strength. However, the success of such algorithms to return reliable values can be hindered by the presence of oscillatory phenomena within magnetic wave guides. Returning accurate parameters is crucial to both magnetohydrodynamics studies and solar physics in general. Here, we employ a simulation featuring propagating MHD waves within a flux tube with a known driver and atmospheric parameters. We invert the Stokes profiles for the 6301 $\unicode{0xc5}$ and 6302 $\unicode{0xc5}$ line pair emergent from the simulations using the well-known Stokes Inversions from Response functions (SIR) code to see if the atmospheric parameters can be returned for typical spatial resolutions at ground-based observatories. The inversions return synthetic spectra comparable to the original input spectra, even with asymmetries introduced in the spectra from wave propagation in the atmosphere. The output models from the inversions match closely to the simulations in temperature, line-of-sight magnetic field and line-of-sight velocity within typical formation heights of the inverted lines. Deviations from the simulations are seen away from these height regions. The inversion results are less accurate during passage of the waves within the line formation region. The original wave period could be recovered from the atmosphere output by the inversions, with empirical mode decomposition performing better than the wavelet approach in this task., Comment: Accepted for publication in Phil. Trans. R. Soc. A, 20 pages, 4 figures

The formation of shocks within the solar atmosphere remains one of the few observable signatures of energy dissipation arising from the plethora of magnetohydrodynamic waves generated close to the solar surface. Active region observations offer exceptional views of wave behavior and its impact on the surrounding atmosphere. The stratified plasma gradients present in the lower solar atmosphere allow for the potential formation of many theorized shock phenomena. In this study, using chromospheric Ca II 854.2nm spectropolarimetric data of a large sunspot, we examine fluctuations in the plasma parameters in the aftermath of powerful shock events that demonstrate polarimetric reversals during their evolution. Modern inversion techniques are employed to uncover perturbations in the temperatures, line-of-sight velocities, and vector magnetic fields occurring across a range of optical depths synonymous with the shock formation. Classification of these non-linear signatures is carried out by comparing the observationally-derived slow, fast, and Alfv\'en shock solutions to the theoretical Rankine-Hugoniot relations. Employing over 200,000 independent measurements, we reveal that the Alfv\'en (intermediate) shock solution provides the closest match between theory and observations at optical depths of log(tau) = -4, consistent with a geometric height at the boundary between the upper photosphere and lower chromosphere. This work uncovers first-time evidence of the manifestation of chromospheric intermediate shocks in sunspot umbrae, providing a new method for the potential thermalization of wave energy in a range of magnetic structures, including pores, magnetic flux ropes, and magnetic bright points., Comment: 17 pages, 9 figures, accepted by ApJ

The planetary boundaries framework proposes quantified guardrails to human modification of global environmental processes that regulate the stability of the planet and has been considered in sustainability science, governance, and corporate management. However, the planetary boundary for human freshwater use has been critiqued as a singular measure that does not reflect all types of human interference with the complex global water cycle and Earth System. We suggest that the water planetary boundary will be more scientifically robust and more useful in decision-making frameworks if it is redesigned to consider more specifically how climate and living ecosystems respond to changes in the different forms of water on Earth: atmospheric water, frozen water, groundwater, soil moisture, and surface water. This paper provides an ambitious scientific road map to define a new water planetary boundary consisting of sub-boundaries that account for a variety of changes to the water cycle.

The planetary boundaries framework proposes quantified guardrails to human modification of global environmental processes that regulate the stability of the planet and has been considered in sustainability science, governance, and corporate management. However, the planetary boundary for human freshwater use has been critiqued as a singular measure that does not reflect all types of human interference with the complex global water cycle and Earth System. We suggest that the water planetary boundary will be more scientifically robust and more useful in decision-making frameworks if it is redesigned to consider more specifically how climate and living ecosystems respond to changes in the different forms of water on Earth: atmospheric water, frozen water, groundwater, soil moisture, and surface water. This paper provides an ambitious scientific road map to define a new water planetary boundary consisting of sub-boundaries that account for a variety of changes to the water cycle.

The planetary boundaries framework proposes quantified guardrails to human modification of global environmental processes that regulate the stability of the planet and has been considered in sustainability science, governance, and corporate management. However, the planetary boundary for human freshwater use has been critiqued as a singular measure that does not reflect all types of human interference with the complex global water cycle and Earth System. We suggest that the water planetary boundary will be more scientifically robust and more useful in decision-making frameworks if it is redesigned to consider more specifically how climate and living ecosystems respond to changes in the different forms of water on Earth: atmospheric water, frozen water, groundwater, soil moisture, and surface water. This paper provides an ambitious scientific road map to define a new water planetary boundary consisting of sub-boundaries that account for a variety of changes to the water cycle.

The planetary boundaries framework proposes quantified guardrails to human modification of global environmental processes that regulate the stability of the planet and has been considered in sustainability science, governance, and corporate management. However, the planetary boundary for human freshwater use has been critiqued as a singular measure that does not reflect all types of human interference with the complex global water cycle and Earth System. We suggest that the water planetary boundary will be more scientifically robust and more useful in decision-making frameworks if it is redesigned to consider more specifically how climate and living ecosystems respond to changes in the different forms of water on Earth: atmospheric water, frozen water, groundwater, soil moisture, and surface water. This paper provides an ambitious scientific road map to define a new water planetary boundary consisting of sub-boundaries that account for a variety of changes to the water cycle.

The planetary boundaries framework proposes quantified guardrails to human modification of global environmental processes that regulate the stability of the planet and has been considered in sustainability science, governance, and corporate management. However, the planetary boundary for human freshwater use has been critiqued as a singular measure that does not reflect all types of human interference with the complex global water cycle and Earth System. We suggest that the water planetary boundary will be more scientifically robust and more useful in decision-making frameworks if it is redesigned to consider more specifically how climate and living ecosystems respond to changes in the different forms of water on Earth: atmospheric water, frozen water, groundwater, soil moisture, and surface water. This paper provides an ambitious scientific road map to define a new water planetary boundary consisting of sub-boundaries that account for a variety of changes to the water cycle.

We report on a search for short-period intensity variations in the green-line FeXIV 530.3 nm emission from the solar corona during the 21 August 2017 total eclipse viewed from Idaho in the United States. Our experiment was performed with a much more sensitive detection system, and with better spatial resolution, than on previous occasions (1999 and 2001 eclipses), allowing fine details of quiet coronal loops and an active-region loop system to be seen. A guided 200-mm-aperture Schmidt-Cassegrain telescope was used with a state-of-the-art CCD camera having 16-bit intensity discrimination and a field-of-view 0.43 degree x 0.43 degree that encompassed approximately one third of the visible corona. The camera pixel size was 1.55 arcseconds, while the seeing during the eclipse enabled features of approx. 2 arcseconds (1450 km on the Sun) to be resolved. A total of 429 images were recorded during a 122.9 second portion of the totality at a frame rate of 3.49 images per second. In the analysis, we searched particularly for short-period intensity oscillations and travelling waves, since theory predicts fast-mode magneto-hydrodynamic (MHD) waves with short periods may be important in quiet coronal and active-region heating. Allowing first for various instrumental and photometric effects, we used a wavelet technique to search for periodicities in some 404 000 pixels in the frequency range 0.5-1.6 Hz (periods: 2 second to 0.6 second). We also searched for travelling waves along some 65 coronal structures. However, we found no statistically significant evidence in either. This negative result considerably refines the limit that we obtained from our previous analyses, and it indicates that future searches for short-period coronal waves may be better directed towards Doppler shifts as well as intensity oscillations.

In this paper we describe observations acquired by satellite instruments (Hinode/SOT and IRIS) and ground-based telescopes (ROSA@DST) during two consecutive C7.0 and X1.6 flares occurred in active region NOAA 12205 on 2014 November 7. The analysis of these data show the presence of continuum enhancements during the evolution of the events, observed both in ROSA images and in IRIS spectra. Moreover, we analyze the role played by the evolution of the $\delta$ sunspots of the active region in the flare triggering, indicating the disappearance of a large portion of penumbra around these sunspots., Comment: 3 pages, 2 figures, to be published in "Nuovo Cimento C" as proceeding of the Third Meeting of the Italian Solar and Heliospheric Community

© 2018 IEEE. Quantum computing's value proposition of an exponential speedup in computing power for certain applications has propelled a vast array of research across the globe. While several different physical implementations of device level qubits are being investigated, semiconductor spin qubits have many similarities to scaled transistors. In this article, we discuss the device/integration of full 300mm based spin qubit devices. This includes the development of (i) a 28 Si epitaxial module ecosystem for growing isotopically pure substrates with among the best Hall mobility at these oxide thicknesses, (ii) a custom 300mm qubit testchip and integration/device line, and (iii) a novel dual nested gate integration process for creating quantum dots.

We report on a search for short-period intensity variations in the green-line FeXIV 530.3 nm emission from the solar corona during the 21 August 2017 total eclipse viewed from Idaho in the United States. Our experiment was performed with a much more sensitive detection system, and with better spatial resolution, than on previous occasions (1999 and 2001 eclipses), allowing fine details of quiet coronal loops and an active-region loop system to be seen. A guided 200-mm-aperture Schmidt-Cassegrain telescope was used with a state-of-the-art CCD camera having 16-bit intensity discrimination and a field-of-view 0.43 degree x 0.43 degree that encompassed approximately one third of the visible corona. The camera pixel size was 1.55 arcseconds, while the seeing during the eclipse enabled features of approx. 2 arcseconds (1450 km on the Sun) to be resolved. A total of 429 images were recorded during a 122.9 second portion of the totality at a frame rate of 3.49 images per second. In the analysis, we searched particularly for short-period intensity oscillations and travelling waves, since theory predicts fast-mode magneto-hydrodynamic (MHD) waves with short periods may be important in quiet coronal and active-region heating. Allowing first for various instrumental and photometric effects, we used a wavelet technique to search for periodicities in some 404 000 pixels in the frequency range 0.5-1.6 Hz (periods: 2 second to 0.6 second). We also searched for travelling waves along some 65 coronal structures. However, we found no statistically significant evidence in either. This negative result considerably refines the limit that we obtained from our previous analyses, and it indicates that future searches for short-period coronal waves may be better directed towards Doppler shifts as well as intensity oscillations.

In this paper we describe observations acquired by satellite instruments (Hinode/SOT and IRIS) and ground-based telescopes (ROSA@DST) during two consecutive C7.0 and X1.6 flares occurred in active region NOAA 12205 on 2014 November 7. The analysis of these data show the presence of continuum enhancements during the evolution of the events, observed both in ROSA images and in IRIS spectra. Moreover, we analyze the role played by the evolution of the $\delta$ sunspots of the active region in the flare triggering, indicating the disappearance of a large portion of penumbra around these sunspots., Comment: 3 pages, 2 figures, to be published in "Nuovo Cimento C" as proceeding of the Third Meeting of the Italian Solar and Heliospheric Community

© 2018 IEEE. Quantum computing's value proposition of an exponential speedup in computing power for certain applications has propelled a vast array of research across the globe. While several different physical implementations of device level qubits are being investigated, semiconductor spin qubits have many similarities to scaled transistors. In this article, we discuss the device/integration of full 300mm based spin qubit devices. This includes the development of (i) a 28 Si epitaxial module ecosystem for growing isotopically pure substrates with among the best Hall mobility at these oxide thicknesses, (ii) a custom 300mm qubit testchip and integration/device line, and (iii) a novel dual nested gate integration process for creating quantum dots.

We present results derived from the analysis of spectropolarimetric measurements of active region AR12546, which represents one of the largest sunspots to have emerged onto the solar surface over the last $20$ years. The region was observed with full-Stokes scans of the Fe I 617.3 nm and Ca II 854.2 nm lines with the Interferometric BIdimensional Spectrometer (IBIS) instrument at the Dunn Solar Telescope over an uncommon, extremely long time interval exceeding three hours. Clear circular polarization (CP) oscillations localized at the umbra-penumbra boundary of the observed region were detected. Furthermore, the multi-height data allowed us to detect the downward propagation of both CP and intensity disturbances at $2.5-3$~mHz, which was identified by a phase delay between these two quantities. These results are interpreted as a propagating magneto-hydrodynamic surface mode in the observed sunspot., Comment: 8 Pages, 4 Figures, ApJ Accepted

Magneto-hydrodynamic (MHD) Alfv\'en waves have been a focus of laboratory plasma physics and astrophysics for over half a century. Their unique nature makes them ideal energy transporters, and while the solar atmosphere provides preferential conditions for their existence, direct detection has proved difficult as a result of their evolving and dynamic observational signatures. The viability of Alfv\'en waves as a heating mechanism relies upon the efficient dissipation and thermalization of the wave energy, with direct evidence remaining elusive until now. Here we provide the first observational evidence of Alfv\'en waves heating chromospheric plasma in a sunspot umbra through the formation of shock fronts. The magnetic field configuration of the shock environment, alongside the tangential velocity signatures, distinguish them from conventional umbral flashes. Observed local temperature enhancements of 5% are consistent with the dissipation of mode-converted Alfv\'en waves driven by upwardly propagating magneto-acoustic oscillations, providing an unprecedented insight into the behaviour of Alfv\'en waves in the solar atmosphere and beyond., Comment: Letter: 7 pages, 4 figures. Supplementary Material: 22 pages, 10 figures

We present results derived from the analysis of spectropolarimetric measurements of active region AR12546, which represents one of the largest sunspots to have emerged onto the solar surface over the last $20$ years. The region was observed with full-Stokes scans of the Fe I 617.3 nm and Ca II 854.2 nm lines with the Interferometric BIdimensional Spectrometer (IBIS) instrument at the Dunn Solar Telescope over an uncommon, extremely long time interval exceeding three hours. Clear circular polarization (CP) oscillations localized at the umbra-penumbra boundary of the observed region were detected. Furthermore, the multi-height data allowed us to detect the downward propagation of both CP and intensity disturbances at $2.5-3$~mHz, which was identified by a phase delay between these two quantities. These results are interpreted as a propagating magneto-hydrodynamic surface mode in the observed sunspot., Comment: 8 Pages, 4 Figures, ApJ Accepted

Magneto-hydrodynamic (MHD) Alfv\'en waves have been a focus of laboratory plasma physics and astrophysics for over half a century. Their unique nature makes them ideal energy transporters, and while the solar atmosphere provides preferential conditions for their existence, direct detection has proved difficult as a result of their evolving and dynamic observational signatures. The viability of Alfv\'en waves as a heating mechanism relies upon the efficient dissipation and thermalization of the wave energy, with direct evidence remaining elusive until now. Here we provide the first observational evidence of Alfv\'en waves heating chromospheric plasma in a sunspot umbra through the formation of shock fronts. The magnetic field configuration of the shock environment, alongside the tangential velocity signatures, distinguish them from conventional umbral flashes. Observed local temperature enhancements of 5% are consistent with the dissipation of mode-converted Alfv\'en waves driven by upwardly propagating magneto-acoustic oscillations, providing an unprecedented insight into the behaviour of Alfv\'en waves in the solar atmosphere and beyond., Comment: Letter: 7 pages, 4 figures. Supplementary Material: 22 pages, 10 figures

Oscillation amplitudes are generally smaller within magnetically active regions like sunspots and plage, when compared to their surroundings. Such magnetic features, when viewed in spatially-resolved powermaps, appear as regions of suppressed power due to reductions in the oscillation amplitudes. Employing high spatial- and temporal-resolution observations from the Dunn Solar Telescope (DST) in New Mexico, we study the power suppression in a region of evolving magnetic fields adjacent to a pore. By utilising wavelet analysis, we study for the first time, how the oscillatory properties in this region change as the magnetic field evolves with time. Image sequences taken in the blue continuum, G-band, Ca~\textsc{ii}~K and H$\alpha$ filters were used in this study. It is observed that the suppression found in the chromosphere occupies a relatively larger area confirming previous findings. Also, the suppression is extended to structures directly connected to the magnetic region and is found to get enhanced as the magnetic field strength increased with time. The dependence of the suppression on the magnetic field strength is greater at longer periods and higher formation heights. Furthermore, the dominant periodicity in the chromosphere was found to be anti-correlated with increases in the magnetic field strength., Comment: accepted for publication in ApJ (in press)

Sunspots on the surface of the Sun are the observational signatures of intense manifestations of tightly packed magnetic field lines, with near-vertical field strengths exceeding 6,000 G in extreme cases. It is well accepted that both the plasma density and the magnitude of the magnetic field strength decrease rapidly away from the solar surface, making high-cadence coronal measurements through traditional Zeeman and Hanle effects difficult since the observational signatures are fraught with low-amplitude signals that can become swamped with instrumental noise. Magneto-hydrodynamic (MHD) techniques have previously been applied to coronal structures, with single and spatially isolated magnetic field strengths estimated as 9-55 G. A drawback with previous MHD approaches is that they rely on particular wave modes alongside the detectability of harmonic overtones. Here we show, for the first time, how omnipresent magneto-acoustic waves, originating from within the underlying sunspot and propagating radially outwards, allow the spatial variation of the local coronal magnetic field to be mapped with high precision. We find coronal magnetic field strengths of 32 +/- 5 G above the sunspot, which decrease rapidly to values of approximately 1 G over a lateral distance of 7000 km, consistent with previous isolated and unresolved estimations. Our results demonstrate a new, powerful technique that harnesses the omnipresent nature of sunspot oscillations to provide magnetic field mapping capabilities close to a magnetic source in the solar corona., Comment: Published in Nature Physics, 12, 179 (2016; DOI:10.1038/nphys3544). 22 pages (Letter + Supplementary Material), including 7 figures

This study presents an “Earth observation-based” method for estimating root zone storage capacity – a critical, yet uncertain parameter in hydrological and land surface modelling. By assuming that vegetation optimises its root zone storage capacity to bridge critical dry periods, we were able to use state-of-the-art satellite-based evaporation data computed with independent energy balance equations to derive gridded root zone storage capacity at global scale. This approach does not require soil or vegetation information, is model independent, and is in principle scale independent. In contrast to a traditional look-up table approach, our method captures the variability in root zone storage capacity within land cover types, including in rainforests where direct measurements ofroot depths otherwise are scarce. Implementing the estimated root zone storage capacity in the global hydrological model STEAM (Simple Terrestrial Evaporation to Atmosphere Model) improved evaporation simulation overall, and in particular during the least evaporating months in sub-humid to humid regions with moderate to high seasonality. Our results suggest that several forest types are able to create a large storage to buffer for severe droughts (with a very long return period), in contrast to, for example, savannahs and woody savannahs (medium length return period), as well as grasslands, shrublands, and croplands (very short return period). The presented method to estimate root zone storage capacity eliminates the need for poor resolution soil and rooting depth data that form a limitation for achieving progress in the global land surface modelling community.

This study presents an “Earth observation-based” method for estimating root zone storage capacity – a critical, yet uncertain parameter in hydrological and land surface modelling. By assuming that vegetation optimises its root zone storage capacity to bridge critical dry periods, we were able to use state-of-the-art satellite-based evaporation data computed with independent energy balance equations to derive gridded root zone storage capacity at global scale. This approach does not require soil or vegetation information, is model independent, and is in principle scale independent. In contrast to a traditional look-up table approach, our method captures the variability in root zone storage capacity within land cover types, including in rainforests where direct measurements ofroot depths otherwise are scarce. Implementing the estimated root zone storage capacity in the global hydrological model STEAM (Simple Terrestrial Evaporation to Atmosphere Model) improved evaporation simulation overall, and in particular during the least evaporating months in sub-humid to humid regions with moderate to high seasonality. Our results suggest that several forest types are able to create a large storage to buffer for severe droughts (with a very long return period), in contrast to, for example, savannahs and woody savannahs (medium length return period), as well as grasslands, shrublands, and croplands (very short return period). The presented method to estimate root zone storage capacity eliminates the need for poor resolution soil and rooting depth data that form a limitation for achieving progress in the global land surface modelling community.

Sunspots on the surface of the Sun are the observational signatures of intense manifestations of tightly packed magnetic field lines, with near-vertical field strengths exceeding 6,000 G in extreme cases. It is well accepted that both the plasma density and the magnitude of the magnetic field strength decrease rapidly away from the solar surface, making high-cadence coronal measurements through traditional Zeeman and Hanle effects difficult since the observational signatures are fraught with low-amplitude signals that can become swamped with instrumental noise. Magneto-hydrodynamic (MHD) techniques have previously been applied to coronal structures, with single and spatially isolated magnetic field strengths estimated as 9-55 G. A drawback with previous MHD approaches is that they rely on particular wave modes alongside the detectability of harmonic overtones. Here we show, for the first time, how omnipresent magneto-acoustic waves, originating from within the underlying sunspot and propagating radially outwards, allow the spatial variation of the local coronal magnetic field to be mapped with high precision. We find coronal magnetic field strengths of 32 +/- 5 G above the sunspot, which decrease rapidly to values of approximately 1 G over a lateral distance of 7000 km, consistent with previous isolated and unresolved estimations. Our results demonstrate a new, powerful technique that harnesses the omnipresent nature of sunspot oscillations to provide magnetic field mapping capabilities close to a magnetic source in the solar corona., Comment: Published in Nature Physics, 12, 179 (2016; DOI:10.1038/nphys3544). 22 pages (Letter + Supplementary Material), including 7 figures

Oscillation amplitudes are generally smaller within magnetically active regions like sunspots and plage, when compared to their surroundings. Such magnetic features, when viewed in spatially-resolved powermaps, appear as regions of suppressed power due to reductions in the oscillation amplitudes. Employing high spatial- and temporal-resolution observations from the Dunn Solar Telescope (DST) in New Mexico, we study the power suppression in a region of evolving magnetic fields adjacent to a pore. By utilising wavelet analysis, we study for the first time, how the oscillatory properties in this region change as the magnetic field evolves with time. Image sequences taken in the blue continuum, G-band, Ca~\textsc{ii}~K and H$\alpha$ filters were used in this study. It is observed that the suppression found in the chromosphere occupies a relatively larger area confirming previous findings. Also, the suppression is extended to structures directly connected to the magnetic region and is found to get enhanced as the magnetic field strength increased with time. The dependence of the suppression on the magnetic field strength is greater at longer periods and higher formation heights. Furthermore, the dominant periodicity in the chromosphere was found to be anti-correlated with increases in the magnetic field strength., Comment: accepted for publication in ApJ (in press)

We present observational evidence of compressible magnetohydrodynamic wave modes propagating from the solar photosphere through to the base of the transition region in a solar magnetic pore. High cadence images were obtained simultaneously across four wavelength bands using the Dunn Solar Telescope. Employing Fourier and wavelet techniques, sausage-mode oscillations displaying significant power were detected in both intensity and area fluctuations. The intensity and area fluctuations exhibit a range of periods from 181-412s, with an average period ~290s, consistent with the global p-mode spectrum. Intensity and area oscillations present in adjacent bandpasses were found to be out-of-phase with one another, displaying phase angles of 6.12 degrees, 5.82 degrees and 15.97 degrees between 4170 Angstrom continuum - G-band, G-band - Na I D1 and Na I D1 - Ca II K heights, respectively, reiterating the presence of upwardly-propagating sausage-mode waves. A phase relationship of ~0 degrees between same-bandpass emission and area perturbations of the pore best categorises the waves as belonging to the `slow' regime of a dispersion diagram. Theoretical calculations reveal that the waves are surface modes, with initial photospheric energies in excess of 35000 W/m^2. The wave energetics indicate a substantial decrease in energy with atmospheric height, confirming that magnetic pores are able to transport waves that exhibit appreciable energy damping, which may release considerable energy into the local chromospheric plasma., Comment: 14 pages, 5 figures, accepted into ApJ

We analyse high temporal and spatial resolution time-series of spectral scans of the Halpha line obtained with the CRisp Imaging SpectroPolarimeter (CRISP) instrument mounted on the Swedish Solar Telescope. The data reveal highly dynamic, dark, short-lived structures known as Rapid Redshifted and Blueshifted Excursions (RREs, RBEs) that are on-disk absorption features observed in the red and blue wings of spectral lines formed in the chromosphere. We study the dynamics of RREs and RBEs by tracking their evolution in space and time, measuring the speed of the apparent motion, line-of-sight Doppler velocity, and transverse velocity of individual structures. A statistical study of their measured properties shows that RREs and RBEs have similar occurrence rates, lifetimes, lengths, and widths. They also display non-periodic, non-linear transverse motions perpendicular to their axes at speeds of 4 - 31 km/s. Furthermore, both types of structures either appear as high speed jets and blobs that are directed outwardly from a magnetic bright point with speeds of 50 - 150 km/s, or emerge within a few seconds. A study of the different velocity components suggests that the transverse motions along the line-of-sight of the chromospheric flux tubes are responsible for the formation and appearance of these redshifted/blueshifted structures. The short lifetime and fast disappearance of the RREs/RBEs suggests that, similar to type II spicules, they are rapidly heated to transition region or even coronal temperatures. We speculate that the Kelvin-Helmholtz instability triggered by observed transverse motions of these structures may be a viable mechanism for their heating., Comment: 8 pages, 6 figures, accepted in ApJ

Water Management, Civil Engineering and Geosciences

We present observational evidence of compressible magnetohydrodynamic wave modes propagating from the solar photosphere through to the base of the transition region in a solar magnetic pore. High cadence images were obtained simultaneously across four wavelength bands using the Dunn Solar Telescope. Employing Fourier and wavelet techniques, sausage-mode oscillations displaying significant power were detected in both intensity and area fluctuations. The intensity and area fluctuations exhibit a range of periods from 181-412s, with an average period ~290s, consistent with the global p-mode spectrum. Intensity and area oscillations present in adjacent bandpasses were found to be out-of-phase with one another, displaying phase angles of 6.12 degrees, 5.82 degrees and 15.97 degrees between 4170 Angstrom continuum - G-band, G-band - Na I D1 and Na I D1 - Ca II K heights, respectively, reiterating the presence of upwardly-propagating sausage-mode waves. A phase relationship of ~0 degrees between same-bandpass emission and area perturbations of the pore best categorises the waves as belonging to the `slow' regime of a dispersion diagram. Theoretical calculations reveal that the waves are surface modes, with initial photospheric energies in excess of 35000 W/m^2. The wave energetics indicate a substantial decrease in energy with atmospheric height, confirming that magnetic pores are able to transport waves that exhibit appreciable energy damping, which may release considerable energy into the local chromospheric plasma., Comment: 14 pages, 5 figures, accepted into ApJ

We analyse high temporal and spatial resolution time-series of spectral scans of the Halpha line obtained with the CRisp Imaging SpectroPolarimeter (CRISP) instrument mounted on the Swedish Solar Telescope. The data reveal highly dynamic, dark, short-lived structures known as Rapid Redshifted and Blueshifted Excursions (RREs, RBEs) that are on-disk absorption features observed in the red and blue wings of spectral lines formed in the chromosphere. We study the dynamics of RREs and RBEs by tracking their evolution in space and time, measuring the speed of the apparent motion, line-of-sight Doppler velocity, and transverse velocity of individual structures. A statistical study of their measured properties shows that RREs and RBEs have similar occurrence rates, lifetimes, lengths, and widths. They also display non-periodic, non-linear transverse motions perpendicular to their axes at speeds of 4 - 31 km/s. Furthermore, both types of structures either appear as high speed jets and blobs that are directed outwardly from a magnetic bright point with speeds of 50 - 150 km/s, or emerge within a few seconds. A study of the different velocity components suggests that the transverse motions along the line-of-sight of the chromospheric flux tubes are responsible for the formation and appearance of these redshifted/blueshifted structures. The short lifetime and fast disappearance of the RREs/RBEs suggests that, similar to type II spicules, they are rapidly heated to transition region or even coronal temperatures. We speculate that the Kelvin-Helmholtz instability triggered by observed transverse motions of these structures may be a viable mechanism for their heating., Comment: 8 pages, 6 figures, accepted in ApJ

Water Management, Civil Engineering and Geosciences

We use ground-based images of high spatial and temporal resolution to search for evidence of nanoflare activity in the solar chromosphere. Through close examination of more than 10^9 pixels in the immediate vicinity of an active region, we show that the distributions of observed intensity fluctuations have subtle asymmetries. A negative excess in the intensity fluctuations indicates that more pixels have fainter-than-average intensities compared with those that appear brighter than average. By employing Monte Carlo simulations, we reveal how the negative excess can be explained by a series of impulsive events, coupled with exponential decays, that are fractionally below the current resolving limits of low-noise equipment on high-resolution ground-based observatories. Importantly, our Monte Carlo simulations provide clear evidence that the intensity asymmetries cannot be explained by photon-counting statistics alone. A comparison to the coronal work of Terzo et al. (2011) suggests that nanoflare activity in the chromosphere is more readily occurring, with an impulsive event occurring every ~360s in a 10,000 km^2 area of the chromosphere, some 50 times more events than a comparably sized region of the corona. As a result, nanoflare activity in the chromosphere is likely to play an important role in providing heat energy to this layer of the solar atmosphere., Comment: 7 pages, 3 figures, accepted into ApJ

We use ground-based images of high spatial and temporal resolution to search for evidence of nanoflare activity in the solar chromosphere. Through close examination of more than 10^9 pixels in the immediate vicinity of an active region, we show that the distributions of observed intensity fluctuations have subtle asymmetries. A negative excess in the intensity fluctuations indicates that more pixels have fainter-than-average intensities compared with those that appear brighter than average. By employing Monte Carlo simulations, we reveal how the negative excess can be explained by a series of impulsive events, coupled with exponential decays, that are fractionally below the current resolving limits of low-noise equipment on high-resolution ground-based observatories. Importantly, our Monte Carlo simulations provide clear evidence that the intensity asymmetries cannot be explained by photon-counting statistics alone. A comparison to the coronal work of Terzo et al. (2011) suggests that nanoflare activity in the chromosphere is more readily occurring, with an impulsive event occurring every ~360s in a 10,000 km^2 area of the chromosphere, some 50 times more events than a comparably sized region of the corona. As a result, nanoflare activity in the chromosphere is likely to play an important role in providing heat energy to this layer of the solar atmosphere., Comment: 7 pages, 3 figures, accepted into ApJ

Aims. We study the formation and evolution of a failed filament eruption observed in NOAA active region 11121 near the southeast limb on November 6, 2010. Methods. We used a time series of SDO/AIA 304, 171, 131, 193, 335, and 94 images, SDO/HMI magnetograms, as well as ROSA and ISOON H_ images to study the erupting active region. Results. We identify coronal loop arcades associated with a quadrupolar magnetic configuration, and show that the expansion and cancellation of the central loop arcade system over the filament is followed by the eruption of the filament. The erupting filament reveals a clear helical twist and develops the same sign of writhe in the form of inverse 298-shape. Conclusions. The observations support the magnetic breakout process in which the eruption is triggered by quadrupolar reconnection in the corona. We propose that the formation mechanism of the inverse -shape flux rope is the magnetohydrodynamic helical kink instability. The eruption has failed because of the large-scale, closed, overlying magnetic loop arcade that encloses the active region., The original document contains color images.

We use images of high spatial and temporal resolution, obtained using both ground- and space-based instrumentation, to investigate the role magnetic field inclination angles play in the propagation characteristics of running penumbral waves in the solar chromosphere. Analysis of a near-circular sunspot, close to the center of the solar disk, reveals a smooth rise in oscillatory period as a function of distance from the umbral barycenter. However, in one directional quadrant, corresponding to the north direction, a pronounced kink in the period-distance diagram is found. Utilizing a combination of the inversion of magnetic Stokes vectors and force-free field extrapolations, we attribute this behaviour to the cut-off frequency imposed by the magnetic field geometry in this location. A rapid, localised inclination of the magnetic field lines in the north direction results in a faster increase in the dominant periodicity due to an accelerated reduction in the cut-off frequency. For the first time we reveal how the spatial distribution of dominant wave periods, obtained with one of the highest resolution solar instruments currently available, directly reflects the magnetic geometry of the underlying sunspot, thus opening up a wealth of possibilities in future magneto-hydrodynamic seismology studies. In addition, the intrinsic relationships we find between the underlying magnetic field geometries connecting the photosphere to the chromosphere, and the characteristics of running penumbral waves observed in the upper chromosphere, directly supports the interpretation that running penumbral wave phenomena are the chromospheric signature of upwardly-propagating magneto-acoustic waves generated in the photosphere., Comment: Accepted for publication in ApJ (11 pages, 5 figures)

We study the formation and evolution of a failed filament eruption observed in NOAA active region 11121 near the southeast limb on November 6, 2010. We use a time series of SDO/AIA 304, 171, 131, 193, 335, 94 {\AA} images, SDO/HMI magnetograms, plus ROSA and ISOON H\alpha images, to study the erupting active region. We identify coronal loop arcades associated with a quadrupolar magnetic configuration, and show that the expansion and cancelation of the central loop arcade system over the filament is followed by the eruption of the filament. The erupting filament reveals a clear helical twist and develops a same sign of writhe in the form of inverse \gamma-shape. The observations support the "magnetic breakout" process with the eruption been triggered by quadrupolar reconnection in the corona. We suggest that the formation mechanism of the inverse \gamma-shape flux rope may be the MHD helical kink instability. The eruption has failed due to the large-scale, closed, overlying magnetic loop arcade that encloses the active region., Comment: 7 Pages and 7 Figures, Accepted into A&A

We use images of high spatial and temporal resolution, obtained using both ground- and space-based instrumentation, to investigate the role magnetic field inclination angles play in the propagation characteristics of running penumbral waves in the solar chromosphere. Analysis of a near-circular sunspot, close to the center of the solar disk, reveals a smooth rise in oscillatory period as a function of distance from the umbral barycenter. However, in one directional quadrant, corresponding to the north direction, a pronounced kink in the period-distance diagram is found. Utilizing a combination of the inversion of magnetic Stokes vectors and force-free field extrapolations, we attribute this behaviour to the cut-off frequency imposed by the magnetic field geometry in this location. A rapid, localised inclination of the magnetic field lines in the north direction results in a faster increase in the dominant periodicity due to an accelerated reduction in the cut-off frequency. For the first time we reveal how the spatial distribution of dominant wave periods, obtained with one of the highest resolution solar instruments currently available, directly reflects the magnetic geometry of the underlying sunspot, thus opening up a wealth of possibilities in future magneto-hydrodynamic seismology studies. In addition, the intrinsic relationships we find between the underlying magnetic field geometries connecting the photosphere to the chromosphere, and the characteristics of running penumbral waves observed in the upper chromosphere, directly supports the interpretation that running penumbral wave phenomena are the chromospheric signature of upwardly-propagating magneto-acoustic waves generated in the photosphere., Comment: Accepted for publication in ApJ (11 pages, 5 figures)

We study the formation and evolution of a failed filament eruption observed in NOAA active region 11121 near the southeast limb on November 6, 2010. We use a time series of SDO/AIA 304, 171, 131, 193, 335, 94 {\AA} images, SDO/HMI magnetograms, plus ROSA and ISOON H\alpha images, to study the erupting active region. We identify coronal loop arcades associated with a quadrupolar magnetic configuration, and show that the expansion and cancelation of the central loop arcade system over the filament is followed by the eruption of the filament. The erupting filament reveals a clear helical twist and develops a same sign of writhe in the form of inverse \gamma-shape. The observations support the "magnetic breakout" process with the eruption been triggered by quadrupolar reconnection in the corona. We suggest that the formation mechanism of the inverse \gamma-shape flux rope may be the MHD helical kink instability. The eruption has failed due to the large-scale, closed, overlying magnetic loop arcade that encloses the active region., Comment: 7 Pages and 7 Figures, Accepted into A&A

High cadence, multi-wavelength observations and simulations are employed for the analysis of solar photospheric magnetic bright points (MBPs) in the quiet Sun. The observations were obtained with the Rapid Oscillations in the Solar Atmosphere (ROSA) imager and the Interferometric BIdimensional Spectrometer (IBIS) at the Dunn Solar Telescope. Our analysis reveals that photospheric MBPs have an average transverse velocity of approximately 1 km/s, whereas their chromospheric counterparts have a slightly higher average velocity of 1.4 km/s. Additionally, chromospheric MBPs were found to be around 63% larger than the equivalent photospheric MBPs. These velocity values were compared with the output of numerical simulations generated using the MURaM code. The simulated results were similar, but slightly elevated, when compared to the observed data. An average velocity of 1.3 km/s was found in the simulated G-band images and an average of 1.8 km/s seen in the velocity domain at a height of 500 km above the continuum formation layer. Delays in the change of velocities were also analysed. Average delays of ~4 s between layers of the simulated data set were established and values of ~29 s observed between G-band and Ca II K ROSA observations. The delays in the simulations are likely to be the result of oblique granular shock waves, whereas those found in the observations are possibly the result of a semi-rigid flux tube., Comment: Accepted for publication in MNRAS Main Journal, 7 pages, 3 figures