In the present work we expand the study of time-dependent ionization previously identified to be of pivotal importance for acoustic waves in solar magnetic flux tube simulations. We focus on longitudinal tube waves (LTW) known to be an important heating agent of solar magnetic regions. Our models also consider new results of wave energy generation as well as an updated determination of the mixing length of convection now identified as 1.8 scale heights in the upper solar convective layers. We present 1-D wave simulations for the solar chromosphere by studying tubes of different spreading as function of height aimed at representing tubes in environments of different magnetic filling factors. Multi-level radiative transfer has been applied to correctly represent the total chromospheric emission function. The effects of time-dependent ionization are significant in all models studied. They are most pronounced behind strong shocks and in low density regions, i.e., the middle and high chromosphere. Concerning our models of different tube spreading, we attained pronounced differences between the various types of models, which were largely initiated by different degrees of dilution of the wave energy flux as well as the density structure partially shaped by strong shocks, if existing. Models showing a quasi-steady rise of temperature with height are obtained via monochromatic waves akin to previous acoustic simulations. However, longitudinal flux tube waves are identified as insufficient to heat the solar transition region and corona in agreement with previous studies.

We investigate the linear correlation coefficient between the intensities at different wavelengths in photospheric and chromospheric spectral lines. Waves that propagate vertically through the stratified solar atmosphere affect different wavelengths at different times. This leads to a characteristic pattern of (non-)coherence of the intensity at various wavelengths. We derived the correlation matrices for several photospheric and chromospheric spectral lines from observations. For comparison with the observations, we calculate correlation matrices for spectra from LTE modeling approaches, 1-D NLTE simulations, and a 3-D MHD simulation run. We apply the correlation method also to temperature maps at different optical depth layers. We find that all photospheric spectral lines show a similar pattern: a pronounced asymmetry of the correlation between line core and red or blue wing. The pattern cannot be reproduced with a simulation of the granulation pattern, but with waves that travel upwards. All chromospheric spectral lines show a more complex pattern. In the case of Ca II H, the 1-D NLTE simulations of monochromatic waves produce a correlation matrix that qualitatively matches to the observations. The photospheric signature is well reproduced in the matrix derived from the 3-D MHD simulation. The correlation matrices of observed photospheric and chromospheric spectral lines are highly structured with characteristic and different patterns in every spectral line. The comparison with matrices derived from simulations and simple modeling suggests that the main driver of the detected patterns are upwards propagating waves. Application of the correlation method to 3-D temperature cubes seems to be a promising tool for a detailed comparison of simulation results and observations in future studies.

We analyze a 1-hour time series of Ca II H intensity spectra and polarimetric spectra around 630 nm to derive the signature of the chromospheric heating and to investigate its relation to magnetic fields. We derived several characteristic quantities of Ca II H to define the chromospheric atmosphere properties. We study the power of the Fourier transform at different wavelengths and their phase relations. We perform local thermodynamic equilibrium inversions of the data to obtain the magnetic field, once including the Ca spectra. We find that the emission in the Ca II H line core at locations without detectable photospheric polarization signal is due to waves that propagate from low forming continuum layers in the line wing up to the line core. The phase differences of intensity oscillations at different wavelengths indicate standing waves for v < 2 mHz and propagating waves for higher frequencies. The waves steepen into shocks in the chromosphere. On average, shocks are both preceded and followed by intensity reductions. In field-free regions, the profiles show emission about half of the time. The correlation between wavelengths and the decorrelation time is significantly higher in the presence of magnetic fields than for field-free areas. The average Ca II H profile in the presence of magnetic fields contains emission features symmetric to the line core and an asymmetric contribution, where mainly the blue H2V emission peak is increased. We find that acoustic waves steepening into shocks are responsible for the emission in the Ca II H line core for locations without photospheric magnetic fields. We suggest using wavelengths in the line wing of Ca II H, where LTE still applies, to compare theoretical heating models with observations.

It is widely believed that the heating of the chromosphere in quiet-Sun internetwork regions is provided by dissipation of acoustic waves that are excited by the convective motions close to the top of the convection zone and in the photospheric overshoot layer. This view lately became challenged by observations suggesting that the acoustic energy flux into the chromosphere is too low, by a factor of at least ten. Based on a comparison of TRACE data with synthetic image sequences for a three-dimensional simulation extending from the top layers of the convection zone to the middle chromosphere, we come to the contradicting conclusion that the acoustic flux in the model provides sufficient energy for heating the solar chromosphere of internetwork regions. The role of a weak magnetic field and associated electric current sheets is also discussed.