Through their lifetime sunspots undergo a change in their area and shape and, as they decay, they fragment into smaller structures. Here, for the first time we analyze the spatial structure of magnetohydrodynamic (MHD) slow body and fast surface modes in observed umbrae as their cross-sectional shape changes. The Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) techniques were used to analyze 3 and 6 hours SDO/HMI time series of Doppler velocities at the photospheric level of approximately circular and elliptically shaped sunspots. Each time series were divided equally into time intervals, to evidence the change of the sunspots' shape. To identify physical wave modes, the POD/DMD modes were cross-correlated with a slow body mode model using the exact shape of the umbra, whereas the shape obtained by applying a threshold level of the mean intensity for every time interval. Our results show that the spatial structure of MHD modes are affected, even by apparently small changes of the umbral shape, especially in the case of the higher-order modes. For the datasets used in our study, the optimal time intervals to consider the influence of the change in the shape on the observed MHD modes is 37 - 60 minutes. The choice of these intervals is crucial to properly quantify the energy contribution of each wave mode to the power spectrum.

Through their lifetime, sunspots undergo a change in their area and shape and, as they decay, they fragment into smaller structures. Here, for the first time we analyze the spatial structure of the magnetohydrodynamic (MHD) slow-body and fast-surface modes in the observed umbrae as their cross-sectional shape changes. The proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) techniques were used to analyze 3 and 6 hr Solar Dynamics Observatory/Helioseismic and Magnetic Imager time series of Doppler velocities at the photospheric level of approximately circular and elliptically shaped sunspots. Each time series was divided into equal time intervals to evidence the change in the shape of the sunspots. To identify the physical wave modes, the POD/DMD modes were cross-correlated with a slow-body mode model using the exact shape of the umbra, whereas the shape obtained by applying a threshold level of the mean intensity for every time interval. Our results show that the spatial structure of MHD modes are affected, even by apparently small changes in the umbral shape, especially in the case of the higher-order modes. For the data sets used in our study, the optimal time intervals to consider the influence of the change in the shape on the observed MHD modes is 37–60 minutes. The choice of these intervals is crucial to properly quantify the energy contribution of each wave mode to the power spectrum.

The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates on the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.

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 the 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 bisector Doppler velocity analysis of the Ca ii 8542 Å 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 line depths 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 waveguides, 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.

Although theoretically predicted, the simultaneous excitation of several resonant modes in sunspots has not been observed. Like any harmonic oscillator, a solar magnetic flux tube can support a variety of resonances, which constitute the natural response of the system to external forcing. Apart from a few single low order eigenmodes in small scale magnetic structures, several simultaneous resonant modes were not found in extremely large sunspots. Here we report the detection of the largest-scale coherent oscillations observed in a sunspot, with a spectrum significantly different from the Sun's global acoustic oscillations, incorporating a superposition of many resonant wave modes. Magnetohydrodynamic numerical modeling agrees with the observations. Our findings not only demonstrate the possible excitation of coherent oscillations over spatial scales as large as 30-40 Mm in extreme magnetic flux regions in the solar atmosphere, but also paves the way for their diagnostic applications in other astrophysical contexts.

Magnetohydrodynamic (MHD) simulations of the solar atmosphere are often performed under the assumption that the plasma is fully ionized. However, in the lower solar atmosphere a reduced temperature often results in only the partial ionization of the plasma. The interaction between the decoupled neutral and ionized components of such a partially ionized plasma produces ambipolar diffusion. To investigate the role of ambipolar diffusion in propagating wave characteristics in the photosphere and chromosphere, we employ the Mancha3D numerical code to model magnetoacoustic waves propagating through the atmosphere immediately above the umbra of a sunspot. We solve the nonideal MHD equations for data-driven perturbations to the magnetostatic equilibrium and the effect of ambipolar diffusion is investigated by varying the simulation to include additional terms in the MHD equations that account for this process. Analyzing the energy spectral densities for simulations with/without ambipolar diffusion, we find evidence to suggest that ambipolar diffusion plays a pivotal role in wave characteristics in the weakly ionized low density regions, hence maximizing the local ambipolar diffusion coefficient. As a result, we propose that ambipolar diffusion is an important mechanism that requires careful consideration into whether it should be included in simulations, and whether it should be utilized in the analysis and interpretation of particular observations of the lower solar atmosphere.

Determining accurate velocity measurements from observations of the Sun is of vital importance to solar physicists who are studying the wave dynamics in the solar atmosphere. Weak chromospheric absorption lines, due to dynamic events in the solar atmosphere, often consist of multiple spectral components. Isolating these components allows for the velocity field of the dynamic and quiescent regimes to be studied independently. However, isolating such components is particularly challenging due to the wide variety of spectral shapes present in the same dataset. MCALF provides a novel method and infrastructure to determine Doppler velocities in a large dataset. Each spectrum is fitted with a model adapted to its specific spectral shape.

Determining accurate plasma Doppler (line-of-sight) velocities from spectroscopic measurements is a challenging endeavour, especially when weak chromospheric absorption lines are often rapidly evolving and, hence, contain multiple spectral components in their constituent line profiles. Here, we present a novel method that employs machine learning techniques to identify the underlying components present within observed spectral lines, before subsequently constraining the constituent profiles through single or multiple Voigt fits. Our method allows active and quiescent components present in spectra to be identified and isolated for subsequent study. Lastly, we employ a Ca ɪɪ 8542 Å spectral imaging dataset as a proof-of-concept study to benchmark the suitability of our code for extracting two-component atmospheric profiles that are commonly present in sunspot chromospheres. Minimization tests are employed to validate the reliability of the results, achieving median reduced χ2-values equal to 1.03 between the observed and synthesized umbral line profiles. This article is part of the Theo Murphy meeting issue ‘High-resolution wave dynamics in the lower solar atmosphere’.

The instrumental advances made in this new era of 4 m class solar telescopes with unmatched spectropolarimetric accuracy and sensitivity will enable the study of chromospheric magnetic fields and their dynamics with unprecedented detail. In this regard, spectropolarimetric diagnostics can provide invaluable insight into magneto-hydrodynamic (MHD) wave processes. MHD waves and, in particular, Alfvénic fluctuations associated with particular wave modes were recently recognized as important mechanisms not only for the heating of the outer layers of the Sun's atmosphere and the acceleration of the solar wind, but also for the elemental abundance anomaly observed in the corona of the Sun and other Sun-like stars (also known as first ionization potential) effect. Here, we take advantage of state-of-the-art and unique spectropolarimetric Interferometric BIdimensional Spectrometer observations to investigate the relation between intensity and circular polarization (CP) fluctuations in a sunspot chromosphere. Our results show a clear link between the intensity and CP fluctuations in a patch which corresponds to a narrow range of magnetic field inclinations. This suggests the presence of Alfvénic perturbations in the sunspot. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.