The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

The SunPy Project developed a 13-question survey to understand the software and hardware usage of the solar physics community. 364 members of the solar physics community, across 35 countries, responded to our survey. We found that 99$\pm$0.5% of respondents use software in their research and 66% use the Python scientific software stack. Students are twice as likely as faculty, staff scientists, and researchers to use Python rather than Interactive Data Language (IDL). In this respect, the astrophysics and solar physics communities differ widely: 78% of solar physics faculty, staff scientists, and researchers in our sample uses IDL, compared with 44% of astrophysics faculty and scientists sampled by Momcheva and Tollerud (2015). 63$\pm$4% of respondents have not taken any computer-science courses at an undergraduate or graduate level. We also found that most respondents utilize consumer hardware to run software for solar-physics research. Although 82% of respondents work with data from space-based or ground-based missions, some of which (e.g. the Solar Dynamics Observatory and Daniel K. Inouye Solar Telescope) produce terabytes of data a day, 14% use a regional or national cluster, 5% use a commercial cloud provider, and 29% use exclusively a laptop or desktop. Finally, we found that 73$\pm$4% of respondents cite scientific software in their research, although only 42$\pm$3% do so routinely.

Extreme ultraviolet (EUV) waves are large-scale propagating disturbances observed in the solar corona, frequently associated with coronal mass ejections and flares. They appear as faint, extended structures propagating from a source region across the structured solar corona. Since their discovery, over two hundred papers discussing their properties, causes and physical nature have been published. However, despite this their fundamental properties and the physics of their interactions with other solar phenomena are still not understood. To further the understanding of EUV waves, we have constructed the Automated Wave Analysis and REduction (AWARE) algorithm for the measurement of EUV waves. AWARE is implemented in two stages. In the first stage, we use a new type of running difference image, the running difference persistence image, which enables the efficient isolation of propagating, brightening wavefronts as they propagate across the corona. In the second stage, AWARE detects the presence of a wavefront, and measures the distance, velocity and acceleration of that wavefront across the Sun. The fit of propagation models to the wave progress isolated in the first stage is achieved using the Random Sample and Consensus (RANSAC) algorithm. AWARE is tested against simulations of EUV wave propagation, and is applied to measure EUV waves in observational data from the Atmospheric Imaging Assembly (AIA). We also comment on unavoidable systematic errors that bias the estimation of wavefront velocity and acceleration. In addition, the full AWARE software suite comes with a package that creates simulations of waves propagating across the disk from arbitrary starting points.

Solar flares often display pulsating and oscillatory signatures in the emission, known as quasi-periodic pulsations (QPP). QPP are typically identified during the impulsive phase of flares, yet in some cases, their presence is detected late into the decay phase. Here, we report extensive fine structure QPP that are detected throughout the large X8.2 flare from 2017 September 10. Following the analysis of the thermal pulsations observed in the GOES/XRS and the 131 A channel of SDO/AIA, we find a pulsation period of ~65 s during the impulsive phase followed by lower amplitude QPP with a period of ~150 s in the decay phase, up to three hours after the peak of the flare. We find that during the time of the impulsive QPP, the soft X-ray source observed with RHESSI rapidly rises at a velocity of approximately 17 km/s following the plasmoid/coronal mass ejection (CME) eruption. We interpret these QPP in terms of a manifestation of the reconnection dynamics in the eruptive event. During the long-duration decay phase lasting several hours, extended downward contractions of collapsing loops/plasmoids that reach the top of the flare arcade are observed in EUV. We note that the existence of persistent QPP into the decay phase of this flare are most likely related to these features. The QPP during this phase are discussed in terms of MHD wave modes triggered in the post-flaring loops.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

The current fleet of space-based solar observatories offers us a wealth of opportunities to study solar flares over a range of wavelengths. Significant advances in our understanding of flare physics often come from coordinated observations between multiple instruments. Consequently, considerable efforts have been, and continue to be made to coordinate observations among instruments (e.g. through the Max Millennium Program of Solar Flare Research). However, there has been no study to date that quantifies how many flares have been observed by combinations of various instruments. Here we describe a technique that retrospectively searches archival databases for flares jointly observed by RHESSI, SDO/EVE (MEGS-A and -B), Hinode/(EIS, SOT, and XRT), and IRIS. Out of the 6953 flares of GOES magnitude C1 or greater that we consider over the 6.5 years after the launch of SDO, 40 have been observed by six or more instruments simultaneously. Using each instrument's individual rate of success in observing flares, we show that the numbers of flares co-observed by three or more instruments are higher than the number expected under the assumption that the instruments operated independently of one another. In particular, the number of flares observed by larger numbers of instruments is much higher than expected. Our study illustrates that these missions often acted in cooperation, or at least had aligned goals. We also provide details on an interactive widget now available in SSWIDL that allows a user to search for flaring events that have been observed by a chosen set of instruments. This provides access to a broader range of events in order to answer specific science questions. The difficulty in scheduling coordinated observations for solar-flare research is discussed with respect to instruments projected to begin operations during Solar Cycle 25, such as DKIST, Solar Orbiter, and Parker Solar Probe.

Solar flare emission at X-ray and extreme ultraviolet (EUV) energies can cause substantial enhancements in the electron density in the Earth's lower ionosphere. It is now become clear that flares exhibit quasi-periodic pulsations with timescales of minutes at X-ray energies, but to date, it has not been known if the ionosphere is sensitive to this variability. Here, using a combination of Very Low Frequency (24 kHz) measurement together with space-based X-ray and EUV observations, we report pulsations of the ionospheric D-region, which are synchronized with a set of pulsating flare loops. Modeling of the ionosphere show that the D-region electron density varies by up to an order of magnitude over the timescale of the pulsations ($\sim$20 mins). Our results reveal that the Earth's ionosphere is more sensitive to small-scale changes in solar soft X-ray flux than previously thought, and implies that planetary ionospheres are closely coupled to small-scale changes in solar/stellar activity.

In this Letter we report the detection of chromospheric 3-minute oscillations in disk-integrated EUV irradiance observations during a solar flare. A wavelet analysis of detrended Lyman-alpha (from GOES/EUVS) and Lyman continuum (from SDO/EVE) emission from the 2011 February 15 X-class flare (SOL2011-02-15T01:56) revealed a $\sim$3-minute period present during the flare's main phase. The formation temperature of this emission locates this radiation to the flare's chromospheric footpoints, and similar behaviour is found in the SDO/AIA 1600\AA\ and 1700\AA\ channels, which are dominated by chromospheric continuum. The implication is that the chromosphere responds dynamically at its acoustic cutoff frequency to an impulsive injection of energy. Since the 3-minute period was not found at hard X-ray energies (50-100 keV) in RHESSI data we can state that this 3-minute oscillation does not depend on the rate of energization of non-thermal electrons. However, a second period of 120 s found in both hard X-ray and chromospheric emission is consistent with episodic electron energization on 2-minute timescales. Our finding on the 3-minute oscillation suggests that chromospheric mechanical energy should be included in the flare energy budget, and the fluctuations in the Lyman-alpha line may influence the composition and dynamics of planetary atmospheres during periods of high activity.

Quasi-periodic pulsations (QPP) seen in the time derivative of the GOES soft X-ray light curves are analyzed for the near-limb X3.2 event on 14 May 2013. The pulsations are apparent for a total of at least two hours from the impulsive phase to well into the decay phase, with a total of 163 distinct pulses evident to the naked eye. A wavelet analysis shows that the characteristic time scale of these pulsations increases systematically from $\sim$25 s at 01:10 UT, the time of the GOES peak, to $\sim$100 s at 02:00 UT. A second ridge in the wavelet power spectrum, most likely associated with flaring emission from a different active region, shows an increase from $\sim$40 s at 01:40 UT to $\sim$100 s at 03:10 UT. We assume that the QPP that produced the first ridge result from vertical kink-mode oscillations of the newly formed loops following magnetic reconnection in the coronal current sheet. This allows us to estimate the magnetic field strength as a function of altitude given the density, loop length, and QPP time scale as functions of time determined from the GOES light curves and RHESSI images. The calculated magnetic field strength of the newly formed loops ranges from about $\sim$500 G at an altitude of 24 Mm to a low value of $\sim$10 G at 60 Mm, in general agreement with the expected values at these altitudes. Fast sausage mode oscillations are also discussed and cannot be ruled out as an alternate mechanism for producing the QPP.

Context. Solar observatories are providing the world-wide community with a wealth of data, covering large time ranges, multiple viewpoints, and returning large amounts of data. In particular, the large volume of SDO data presents challenges: it is available only from a few repositories, and full-disk, full-cadence data for reasonable durations of scientific interest are difficult to download practically due to their size and download data rates available to most users. From a scientist's perspective this poses three problems: accessing, browsing and finding interesting data as efficiently as possible. Aims. To address these challenges, we have developed JHelioviewer, a visualisation tool for solar data based on the JPEG2000 compression standard and part of the open source ESA/NASA Helioviewer Project. Since the first release of JHelioviewer, the scientific functionality of the software has been extended significantly, and the objective of this paper is to highlight these improvements. Methods. The JPEG2000 standard offers useful new features that facilitate the dissemination and analysis of high-resolution image data and offers a solution to the challenge of efficiently browsing petabyte-scale image archives. The JHelioviewer software is open source, platform independent and extendable via a plug-in architecture. Results. With JHelioviewer, users can visualise the Sun for any time period between September 1991 and today. They can perform basic image processing in real time, track features on the Sun and interactively overlay magnetic field extrapolations. The software integrates solar event data and a time line display. As a first step towards supporting science planning of the upcoming Solar Orbiter mission, JHelioviewer offers a virtual camera model that enables users to set the vantage point to the location of a spacecraft or celestial body at any given time.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

The nature of quasi-periodic pulsations in solar flares is poorly constrained, and critically the general prevalence of such signals in solar flares is unknown. Therefore, we perform a large-scale search for evidence of signals consistent with quasi-periodic pulsations in solar flares, focusing on the 1 - 300s timescale. We analyse 675 M- and X-class flares observed by GOES in 1-8\AA\ soft X-rays between 2011 February 1 and 2015 December 31. Additionally, over the same era we analyse Fermi/GBM 15-25 keV X-ray data for each of these flares that was associated with a Fermi/GBM solar flare trigger, a total of 261 events. Using a model comparison method, we determine whether there is evidence for a substantial enhancement in the Fourier power spectrum that may be consistent with a QPP signature, based on three tested models; a power-law plus a constant, a broken power-law plus constant, and a power-law-plus-constant with an additional QPP signature component. From this, we determine that ~30% of GOES events and ~8% of Fermi/GBM events show strong signatures consistent with classical interpretations of QPP. For the remaining events either two or more tested models cannot be strongly distinguished from each other, or the events are well-described by single power-law or broken power-law Fourier power spectra. For both instruments, a preferred characteristic timescale of ~5-30 s was found in the QPP-like events, with no dependence on flare magnitude in either GOES or GBM data. We also show that individual events in the sample show similar characteristic timescales in both GBM and GOES datasets. We discuss the implications of these results for our understanding of solar flares and possible QPP mechanisms.

Quasi-periodic pulsations (QPP) are often observed in X-ray emission from solar flares. To date, it is unclear what their physical origins are. Here, we present a multi-instrument investigation of the nature of QPP during the impulsive and decay phases of the X1.0 flare of 28 October 2013. We focus on the character of the fine structure pulsations evident in the soft X-ray time derivatives and compare this variability with structure across multiple wavelengths including hard X-ray and microwave emission. We find that during the impulsive phase of the flare, high correlations between pulsations in the thermal and non-thermal emissions are seen. A characteristic timescale of ~20s is observed in all channels and a second timescale of ~55s is observed in the non-thermal emissions. Soft X-ray pulsations are seen to persist into the decay phase of this flare, up to 20 minutes after the non-thermal emission has ceased. We find that these decay phase thermal pulsations have very small amplitude and show an increase in characteristic timescale from ~40s up to ~70s. We interpret the bursty nature of the co-existing multi-wavelength QPP during the impulsive phase in terms of episodic particle acceleration and plasma heating. The persistent thermal decay phase QPP are most likely connected with compressive MHD processes in the post-flare loops such as the fast sausage mode or the vertical kink mode.

This paper presents SunPy (version 0.5), a community-developed Python package for solar physics. Python, a free, cross-platform, general-purpose, high-level programming language, has seen widespread adoption among the scientific community, resulting in the availability of a large number of software packages, from numerical computation (NumPy, SciPy) and machine learning (scikit-learn) to visualisation and plotting (matplotlib). SunPy is a data-analysis environment specialising in providing the software necessary to analyse solar and heliospheric data in Python. SunPy is open-source software (BSD licence) and has an open and transparent development workflow that anyone can contribute to. SunPy provides access to solar data through integration with the Virtual Solar Observatory (VSO), the Heliophysics Event Knowledgebase (HEK), and the HELiophysics Integrated Observatory (HELIO) webservices. It currently supports image data from major solar missions (e.g., SDO, SOHO, STEREO, and IRIS), time-series data from missions such as GOES, SDO/EVE, and PROBA2/LYRA, and radio spectra from e-Callisto and STEREO/SWAVES. We describe SunPy's functionality, provide examples of solar data analysis in SunPy, and show how Python-based solar data-analysis can leverage the many existing tools already available in Python. We discuss the future goals of the project and encourage interested users to become involved in the planning and development of SunPy.

The nature of quasi-periodic pulsations in solar and stellar flares remains debated. Recent work has shown that power-law-like Fourier power spectra, also referred to as 'red' noise processes, are an intrinsic property of solar and stellar flare signals, a property that many previous studies of this phenomenon have not accounted for. Hence a re-evaluation of the existing interpretations and assumptions regarding QPP is needed. Here we adopt a Bayesian method for investigating this phenomenon, fully considering the Fourier power law properties of flare signals. Using data from the PROBA2/LYRA, Fermi/GBM, Nobeyama Radioheliograph and Yohkoh/HXT instruments, we study a selection of flares from the literature identified as QPP events. Additionally we examine optical data from a recent stellar flare that appears to exhibit oscillatory properties. We find that, for all but one event tested, an explicit oscillation is not required in order to explain the observations. Instead, the flare signals are adequately described as a manifestation of a power law in the Fourier power spectrum, rather than a direct signature of oscillating components or structures. However, for the flare of 1998 May 8, strong evidence for the existence of an explicit oscillation with P ~ 14-16 s is found in the 17 GHz radio data and the 13-23 keV Yohkoh HXT data. We conclude that, most likely, many previously analysed events in the literature may be similarly described in terms of power laws in the flare Fourier power spectrum, without the need to invoke a narrowband, oscillatory component. As a result the prevalence of oscillatory signatures in solar and stellar flares may be less than previously believed. The physical mechanism behind the appearance of the observed power laws is discussed.

The dynamics of regions of the solar corona are investigated using Atmospheric Imaging Assembly (AIA) 171\AA\ and 193\AA\ data. The coronal emission from the quiet Sun, coronal loop footprints, coronal moss, and from above a sunspot is studied. It is shown that the mean Fourier power spectra in these regions can be described by a power law at lower frequencies that tails to flat spectrum at higher frequencies, plus a Gaussian-shaped contribution that varies depending on the region studied. This Fourier spectral shape is in contrast to the commonly-held assumption that coronal time-series are well described by the sum of a long time-scale background trend plus Gaussian-distributed noise, with some specific locations also showing an oscillatory signal. The implications of this discovery to the field of coronal seismology and the automated detections of oscillations are discussed. The power law contribution to the shape of the Fourier power spectrum is interpreted as being due to the summation of a distribution of exponentially decaying emission events along the line of sight. This is consistent with the idea that the solar atmosphere is heated everywhere by small energy deposition events.

This paper presents measurements of the energy radiated by the lower solar atmosphere, at optical, UV, and EUV wavelengths, during an X-class solar flare (SOL2011-02-15T01:56) in response to an injection of energy assumed to be in the form of nonthermal electrons. Hard X-ray observations from RHESSI were used to track the evolution of the parameters of the nonthermal electron distribution to reveal the total power contained in flare accelerated electrons. By integrating over the duration of the impulsive phase, the total energy contained in the nonthermal electrons was found to be $>2\times10^{31}$ erg. The response of the lower solar atmosphere was measured in the free-bound EUV continua of H I (Lyman), He I, and He II, plus the emission lines of He II at 304\AA\ and H I (Ly$\alpha$) at 1216\AA\ by SDO/EVE, the UV continua at 1600\AA\ and 1700\AA\ by SDO/AIA, and the WL continuum at 4504\AA, 5550\AA, and 6684\AA, along with the Ca II H line at 3968\AA\ using Hinode/SOT. The summed energy detected by these instruments amounted to $\sim3\times10^{30}$ erg; about 15% of the total nonthermal energy. The Ly$\alpha$ line was found to dominate the measured radiative losses. Parameters of both the driving electron distribution and the resulting chromospheric response are presented in detail to encourage the numerical modelling of flare heating for this event, to determine the depth of the solar atmosphere at which these line and continuum processes originate, and the mechanism(s) responsible for their generation.

We compare four different methods of calculating uncertainty estimates in fitting parameterized models to RHESSI X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the GOES X1.3 class flare of 19 January 2005, and the other from the X4.8 flare of 23 July 2002. The four methods give approximately the same uncertainty estimates for the 19 January 2005 spectral fit parameters, but lead to very different uncertainty estimates for the 23 July 2002 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that differ greatly depending on the shape of the hypersurface. For the 23 July 2002 flare data, there is a 95% probability that the low energy cutoff lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. The low-energy cutoff for the 19 January 2005 flare is more tightly constrained to 107+/-4 keV (68% probability). Using the Bayesian/MCMC approach, we also estimate for the first time probability density functions (PDFs) for the total number of flare accelerated electrons and the energy they carry. For the 23 July 2002 event, these PDFs are asymmetric with long tails orders of magnitude higher than the most probable value, caused by the poorly constrained value of the low-energy cutoff. For the 19 January 2005 flare spectrum, the PDFs for the total number of flare accelerated electrons and their energy are much more symmetric and narrow.

Recently observed oscillations in the solar atmosphere have been interpreted and modeled as magnetohydrodynamic wave modes. This has allowed the estimation of parameters that are otherwise hard to derive, such as the coronal magnetic-field strength. This work crucially relies on the initial detection of the oscillations, which is commonly done manually. The volume of Solar Dynamics Observatory (SDO) data will make manual detection inefficient for detecting all of the oscillating regions. An algorithm is presented which automates the detection of areas of the solar atmosphere that support spatially extended oscillations. The algorithm identifies areas in the solar atmosphere whose oscillation content is described by a single, dominant oscillation within a user-defined frequency range. The method is based on Bayesian spectral analysis of time-series and image filtering. A Bayesian approach sidesteps the need for an a-priori noise estimate to calculate rejection criteria for the observed signal, and it also provides estimates of oscillation frequency, amplitude and noise, and the error in all these quantities, in a self-consistent way. The algorithm also introduces the notion of quality measures to those regions for which a positive detection is claimed, allowing simple post-detection discrimination by the user. The algorithm is demonstrated on two Transition Region and Coronal Explorer (TRACE) datasets, and comments regarding its suitability for oscillation detection in SDO are made.

Across all disciplines that work with image data - from astrophysics to medical research and historic preservation - there is a growing need for efficient ways to browse and inspect large sets of high-resolution images. We present the development of a visualization software for solar physics data based on the JPEG 2000 image compression standard. Our implementation consists of the JHelioviewer client application that enables users to browse petabyte-scale image archives and the JHelioviewer server, which integrates a JPIP server, metadata catalog and an event server. JPEG 2000 offers many useful new features and has the potential to revolutionize the way high-resolution image data are disseminated and analyzed. This is especially relevant for solar physics, a research field in which upcoming space missions will provide more than a terabyte of image data per day. Providing efficient access to such large data volumes at both high spatial and high time resolution is of paramount importance to support scientific discovery.

Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.

A Bayesian probability based approach is applied to the problem of detecting and parameterizing oscillations in the upper solar atmosphere for the first time. Due to its statistical origin, this method provides a mechanism for determining the number of oscillations present, gives precise estimates of the oscillation parameters with a self-consistent statistical error analysis, and allows the oscillatory model signals to be reconstructed within these errors. A highly desirable feature of the Bayesian approach is the ability to resolve oscillations with extremely small frequency separations. The code is applied to SOHO/CDS (Solar and Heliospheric Observatory/Coronal Diagnostic Spectrometer) O V 629A observations and resolves four distinct P4, P5, P6 and P7 p-modes within the same sunspot transition region. This suggests that a spectrum of photospheric p-modes is able to propagate into the upper atmosphere of the Sun and Sun-like stars, and places precise observational constraints on models of umbral eigen modes.