Source author record

Joseph M. Davila

Joseph M. Davila appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2015arXiv

Evidence of thermal conduction suppression in a solar flaring loop by coronal seismology of slow-mode waves

Analysis of a longitudinal wave event observed by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) is presented. A time sequence of 131 A images reveals that a C-class flare occurred at one footpoint of a large loop and triggered an intensity disturbance (enhancement) propagating along it. The spatial features and temporal evolution suggest that a fundamental standing slow-mode wave could be set up quickly after meeting of two initial disturbances from the opposite footpoints. The oscillations have a period of ~12 min and a decay time of ~9 min. The measured phase speed of 500$\pm$50 km/s matches the sound speed in the heated loop of ~10 MK, confirming that the observed waves are of slow mode. We derive the time-dependent temperature and electron density wave signals from six AIA extreme-ultraviolet (EUV) channels, and find that they are nearly in phase.The measured polytropic index from the temperature and density perturbations is 1.64$\pm$0.08 close to the adiabatic index of 5/3 for an ideal monatomic gas. The interpretation based on a 1D linear MHD model suggests that the thermal conductivity is suppressed by at least a factor of 3 in the hot flare loop at 9 MK and above. The viscosity coefficient is determined by coronal seismology from the observed wave when only considering the compressive viscosity dissipation. We find that to interpret the rapid wave damping, the classical compressive viscosity coefficient needs to be enhanced by a factor of 15 as the upper limit.

preprint2015arXiv

Optimizing Global Coronal Magnetic Field Models Using Image-Based Constraints

The coronal magnetic field directly or indirectly affects a majority of the phenomena studied in space physics. It provides energy for coronal heating, controls the release of coronal mass ejections (CMEs), and drives heliospheric and magnetospheric activity, yet the coronal magnetic field itself has proven difficult to measure. This difficulty has prompted a decades-long effort to develop accurate, timely, models of the field - an effort that continues today. We have developed a method for improving global coronal magnetic field models by incorporating the type of morphological constraints which could be derived from coronal images. Here we report promising initial tests of this approach on two theoretical problems, and discuss opportunities for application.

preprint2015arXiv

Spectroscopic Diagnosis of Propagating disturbances in coronal loops: Waves or flows?

The analysis of multiwavelength properties of propagating disturbances (PDs) using Hinode/EIS observations is presented. Quasi-periodic PDs were mostly interpreted as slow magnetoacoustic waves in early studies, but recently suggested to be intermittent upflows of the order of 50-150 km/s based on the Red-Blue (RB) asymmetry analysis of spectral line profiles. Using the forward models, velocities of the secondary component derived from the RB analysis are found significantly overestimated due to the saturation effect when its offset velocities are smaller than the Gaussian width. We developed a different method to examine spectral features of the PDs. This method is assuming that the excessive emission of the PD profile against the background (taken as that prior to the PD) is caused by a hypothetic upflow. The derived LOS velocities of the flow are on the order of 10-30 km/s from the warm (1-1.5 MK) coronal lines, much smaller than those inferred from the RB analysis. This result does not support the flow interpretation but favors of the early wave interpretation.

preprint2014arXiv

Spatiotemporal organization of energy release events in the quiet solar corona

Using data from STEREO and SOHO spacecraft, we show that temporal organization of energy release events in the quiet solar corona is close to random, in contrast to the clustered behavior of flaring times in solar active regions. The locations of the quiet-Sun events follow the meso- and supergranulation pattern of the underling photosphere. Together with earlier reports of the scale-free event size statistics, our findings suggest that quiet solar regions responsible for bulk coronal heating operate in a driven self-organized critical state, possibly involving long-range Alfvénic interactions.

preprint2014arXiv

Validation of Spherically Symmetric Inversion by Use of a Tomographic Reconstructed Three-Dimensional Electron Density of the Solar Corona

Determination of the coronal electron density by the inversion of white-light polarized brightness (pB) measurements by coronagraphs is a classic problem in solar physics. An inversion technique based on the spherically symmetric geometry (Spherically Symmetric Inversion, SSI) was developed in the 1950s, and has been widely applied to interpret various observations. However, to date there is no study about uncertainty estimation of this method. In this study we present the detailed assessment of this method using a three-dimensional (3D) electron density in the corona from 1.5 to 4 Rsun as a model, which is reconstructed by tomography method from STEREO/COR1 observations during solar minimum in February 2008. We first show in theory and observation that the spherically symmetric polynomial approximation (SSPA) method and the Van de Hulst inversion technique are equivalent. Then we assess the SSPA method using synthesized pB images from the 3D density model, and find that the SSPA density values are close to the model inputs for the streamer core near the plane of the sky (POS) with differences generally less than a factor of two or so; the former has the lower peak but more spread in both longitudinal and latitudinal directions than the latter. We estimate that the SSPA method may resolve the coronal density structure near the POS with angular resolution in longitude of about 50 degrees. Our results confirm the suggestion that the SSI method is applicable to the solar minimum streamer (belt) as stated in some previous studies. In addition, we demonstrate that the SSPA method can be used to reconstruct the 3D coronal density, roughly in agreement with that by tomography for a period of low solar activity. We suggest that the SSI method is complementary to the 3D tomographic technique in some cases, given that the development of the latter is still an ongoing research effort.

preprint2013arXiv

Measuring temperature - dependent propagating disturbances in coronal fan loops using multiple SDO/AIA channels and surfing transform technique

A set of co-aligned high resolution images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) is used to investigate propagating disturbances (PDs) in warm fan loops at the periphery of a non-flaring active region NOAA AR 11082. To measure PD speeds at multiple coronal temperatures, a new data analysis methodology is proposed enabling quantitative description of subvisual coronal motions with low signal-to-noise ratios of the order of 0.1 %. The technique operates with a set of one-dimensional "surfing" signals extracted from position-time plots of several AIA channels through a modified version of Radon transform. The signals are used to evaluate a two-dimensional power spectral density distribution in the frequency - velocity space which exhibits a resonance in the presence of quasi-periodic PDs. By applying this analysis to the same fan loop structures observed in several AIA channels, we found that the traveling velocity of PDs increases with the temperature of the coronal plasma following the square root dependence predicted for the slow mode magneto-acoustic wave which seems to be the dominating wave mode in the studied loop structures. This result extends recent observations by Kiddie et al. (Solar Phys., 2012) to a more general class of fan loop systems not associated with sunspots and demonstrating consistent slow mode activity in up to four AIA channels.

preprint2013arXiv

Stochastic coupling of solar photosphere and corona

The observed solar activity is believed to be driven by the dissipation of nonpotential magnetic energy injected into the corona by dynamic processes in the photosphere. The enormous range of scales involved in the interaction makes it difficult to track down the photospheric origin of each coronal dissipation event, especially in the presence of complex magnetic topologies. In this paper, we propose an ensemble-based approach for testing the photosphere - corona coupling in a quiet solar region as represented by intermittent activity in SOHO MDI and STEREO EUVI image sets. For properly adjusted detection thresholds corresponding to the same degree of intermittency in the photosphere and corona, the dynamics of the two solar regions is described by the same occurrence probability distributions of energy release events but significantly different geometric properties. We derive a set of scaling relations reconciling the two groups of results and enabling statistical description of coronal dynamics based on photospheric observations. Our analysis suggests that multiscale intermittent dissipation in the corona at spatial scales > 3 Mm is controlled by turbulent photospheric convection. Complex topology of the photospheric network makes this coupling essentially nonlocal and non-deterministic. Our results are in an agreement with the Parker's coupling scenario in which random photospheric shuffling generates marginally stable magnetic discontinuities at the coronal level, but they are also consistent with an impulsive wave heating involving multiscale Alfvenic wave packets and/or MHD turbulent cascade. A back reaction on the photosphere due to coronal magnetic reconfiguration can be a contributing factor.

preprint2013arXiv

Three-dimensional MHD modeling of propagating disturbances in fan-like coronal loops

Quasi-periodic propagating intensity disturbances (PDs) have been observed in large coronal loops in EUV images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. Motivated by the scenario that the coronal upflows could be cumulative result of numerous individual flow pulses generated by sporadic heating events (nanoflares) at the loop base, we construct a velocity driver with repetitive tiny pulses, whose energy frequency distribution follows the flare power-law scaling. We then perform 3D MHD modeling of an idealized bipolar active region by applying this broadband velocity driver at the footpoints of large coronal loops which appear open in the computational domain. Our model successfully reproduces the PDs with similar features as the observed, and shows that any upflow pulses inevitably excite slow magnetosonic wave disturbances propagating along the loop. We find that the generated PDs are dominated by the wave signature as their propagation speeds are consistent with the wave speed in the presence of flows, and the injected flows rapidly decelerate with height. Our simulation results suggest that the observed PDs and associated persistent upflows may be produced by small-scale impulsive heating events (nanoflares) at the loop base, and that the flows and waves may both contribute to the PDs at lower heights.

preprint2012arXiv

Growing transverse oscillations of a multistranded loop observed by SDO/AIA

The first evidence of transverse oscillations of a multistranded loop with growing amplitudes and internal coupling observed by the Atomspheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) is presented. The loop oscillation event occurred on 2011 March 8, triggered by a CME. The multiwavelength analysis reveals the presence of multithermal strands in the oscillating loop, whose dynamic behaviors are temperature-dependent, showing differences in their oscillation amplitudes, phases and emission evolution. The physical parameters of growing oscillations of two strands in 171 A are measured and the 3-D loop geometry is determined using STEREO-A/EUVI data. These strands have very similar frequencies, and between two 193 A strands a quarter-period phase delay sets up. These features suggest the coupling between kink oscillations of neighboring strands and the interpretation by the collective kink mode as predicted by some models. However, the temperature dependence of the multistarnded loop oscillations was not studied previously and needs further investigation. The transverse loop oscillations are associated with intensity and loop width variations. We suggest that the amplitude-growing kink oscillations may be a result of continuous non-periodic driving by magnetic deformation of the CME, which deposits energy into the loop system at a rate faster than its loss.

preprint2012arXiv

Multiscale dynamics of solar magnetic structures

Multiscale topological complexity of solar magnetic field is among the primary factors controlling energy release in the corona, including associated processes in the photospheric and chromospheric boundaries. We present a new approach for analyzing multiscale behavior of the photospheric magnetic flux underlying this dynamics as depicted by a sequence of high-resolution solar magnetograms. The approach involves two basic processing steps: (1) identification of timing and location of magnetic flux origin and demise events (as defined by DeForest et al., 2007) by tracking spatiotemporal evolution of unipolar and bipolar photospheric regions, and (2) analysis of collective behavior of the detected magnetic events using a generalized version of Grassberger - Procaccia correlation integral algorithm. The scale-free nature of the developed algorithms makes it possible to characterize the dynamics of the photospheric network across a wide range of distances and relaxation times. Three types of photospheric conditions are considered to test the method: a quiet photosphere, a solar active region (NOAA 10365) in a quiescent non-flaring state, and the same active region during a period of M-class flares. The results obtained show (1) the presence of a topologically complex asymmetrically fragmented magnetic network in the quiet photosphere driven by meso- and supergranulation, (2) the formation of non-potential magnetic structures with complex polarity separation lines inside the active region, and (3) statistical signatures of canceling bipolar magnetic structures coinciding with flaring activity in the active region. Each of these effects can represent an unstable magnetic configuration acting as an energy source for coronal dissipation and heating.

preprint2012arXiv

Propagating Intensity Disturbances in Fan-like Coronal Loops: Flows or Waves?

Quasi-periodic intensity disturbances propagating upward along the coronal structure have been extensively studied using EUV imaging observations from SOHO/EIT and TRACE. They were interpreted as either slow mode magnetoacoustic waves or intermittent upflows. In this study we aim at demonstrating that time series of spectroscopic observations are critical to solve this puzzle. Propagating intensity and Doppler shift disturbances in fanlike coronal loops are analyzed in multiple wavelengths using the sit-and-stare observations from Hinode/EIS. We find that the disturbances did not cause the blue-wing asymmetry of spectral profiles in the warm (~1.5 MK) coronal lines. The estimated small line-of-sight velocities also did not support the intermittent upflow interpretation. In the hot (~2 MK) coronal lines the disturbances did cause the blue-wing asymmetry, but the double fits revealed that a high-velocity minor component is steady and persistent, while the propagating intensity and Doppler shift disturbances are mainly due to variations of the core component, therefore, supporting the slow wave interpretation. However, the cause for blueward line asymmetries remains unclear.

preprint2012arXiv

Slow magnetosonic waves and fast flows in active region loops

Recent EUV spectroscopic observations indicate that slow magnetosonic waves are present in active region (AR) loops. Some of the spectral data were also interpreted as evidence of fast (~100-300 km/s) quasi-periodic flows. We have performed three-dimensional magnetohydrodynamic (3D MHD) modeling of a bipolar AR that contains impulsively generated waves and flows in coronal loops. The model AR is initiated with a dipole magnetic field and gravitationally stratified density, with an upflow driven steadily or periodically in localized regions at the footpoints of magnetic loops. The resulting flows along the magnetic field lines of the AR produce higher density loops compared to the surrounding plasma by injection of material into the flux-tubes and the establishment of siphon flow. We find that the impulsive onset of flows with subsonic speeds result in the excitation of damped slow magnetosonic waves that propagate along the loops and coupled nonlinearly driven fast mode waves. The phase speed of the slow magnetosonic waves is close to the coronal sound speed. When the amplitude of the driving pulses is increased we find that slow shock-like wave trains are produced. When the upflows are driven periodically, undamped oscillations are produced with periods determined by the periodicity of the upflows. Based on the results of the 3D MHD model we suggest that the observed slow magnetosonic waves and persistent upflows may be produced by the same impulsive events at the bases of ARs.

preprint2012arXiv

The Relationship Between the Expansion Speed and Radial Speed of CMEs Confirmed Using Quadrature Observations of the 2011 February 15 CME

It is difficult to measure the true speed of Earth-directed CMEs from a coronagraph located along the Sun-Earth line because of the occulting disk. However, the expansion speed (the speed with which the CME appears to spread in the sky plane) can be measured by such a coronagraph. In order to convert the expansion speed to radial speed (which is important for space weather applications) one can use an empirical relationship between the two that assumes an average width for all CMEs. If we have the width information from quadrature observations, we can confirm the relationship between expansion and radial speeds derived by Gopalswamy et al. (2009a). The STEREO spacecraft were in qudrature with SOHO (STEREO-A ahead of Earth by 87o and STEREO-B 94o behind Earth) on 2011 February 15, when a fast Earth-directed CME occurred. The CME was observed as a halo by the Large-Angle and Spectrometric Coronagraph (LASCO) on board SOHO. The sky-plane speed was measured by SOHO/LASCO as the expansion speed, while the radial speed was measured by STEREO-A and STEREO-B. In addition, STEREO-A and STEREO-B images provided the width of the CME, which is unknown from Earth view. From the SOHO and STEREO measurements, we confirm the relationship between the expansion speed (Vexp) and radial speed (Vrad) derived previously from geometrical considerations (Gopalswamy et al. 2009a): Vrad = 1/2 (1 + cot w)Vexp, where w is the half width of the CME. STEREO-B images of the CME, we found that CME had a full width of 76o, so w = 38o. This gives the relation as Vrad = 1.14 Vexp. From LASCO observations, we measured Vexp = 897 km/s, so we get the radial speed as 1023 km/s. Direct measurement of radial speed yields 945 km/s (STEREO-A) and 1058 km/s (STEREO-B). These numbers are different only by 7.6% and 3.4% (for STEREO-A and STEREO-B, respectively) from the computed value.

preprint2011arXiv

Earth-Affecting Solar Causes Observatory (EASCO): A mission at the Sun-Earth L5

Coronal mass ejections (CMEs) and corotating interaction regions (CIRs) as well as their source regions are important because of their space weather consequences. The current understanding of CMEs primarily comes from the Solar and Heliospheric Observatory (SOHO) and the Solar Terrestrial Relations Observatory (STEREO) missions, but these missions lacked some key measurements: STEREO did not have a magnetograph; SOHO did not have in-situ magnetometer. SOHO and other imagers such as the Solar Mass Ejection Imager (SMEI) located on the Sun-Earth line are also not well-suited to measure Earth-directed CMEs. The Earth-Affecting Solar Causes Observatory (EASCO) is a proposed mission to be located at the Sun-Earth L5 that overcomes these deficiencies. The mission concept was recently studied at the Mission Design Laboratory (MDL), NASA Goddard Space Flight Center, to see how the mission can be implemented. The study found that the scientific payload (seven remote-sensing and three in-situ instruments) can be readily accommodated and can be launched using an intermediate size vehicle; a hybrid propulsion system consisting of a Xenon ion thruster and hydrazine has been found to be adequate to place the payload at L5. Following a 2-year transfer time, a 4-year operation is considered around the next solar maximum in 2025.

preprint2011arXiv

LEMUR: Large European Module for solar Ultraviolet Research. European contribution to JAXA's Solar-C mission

Understanding the solar outer atmosphere requires concerted, simultaneous solar observations from the visible to the vacuum ultraviolet (VUV) and soft X-rays, at high spatial resolution (between 0.1" and 0.3"), at high temporal resolution (on the order of 10 s, i.e., the time scale of chromospheric dynamics), with a wide temperature coverage (0.01 MK to 20 MK, from the chromosphere to the flaring corona), and the capability of measuring magnetic fields through spectropolarimetry at visible and near-infrared wavelengths. Simultaneous spectroscopic measurements sampling the entire temperature range are particularly important. These requirements are fulfilled by the Japanese Solar-C mission (Plan B), composed of a spacecraft in a geosynchronous orbit with a payload providing a significant improvement of imaging and spectropolarimetric capabilities in the UV, visible, and near-infrared with respect to what is available today and foreseen in the near future. The Large European Module for solar Ultraviolet Research (LEMUR), described in this paper, is a large VUV telescope feeding a scientific payload of high-resolution imaging spectrographs and cameras. LEMUR consists of two major components: a VUV solar telescope with a 30 cm diameter mirror and a focal length of 3.6 m, and a focal-plane package composed of VUV spectrometers covering six carefully chosen wavelength ranges between 17 and 127 nm. The LEMUR slit covers 280" on the Sun with 0.14" per pixel sampling. In addition, LEMUR is capable of measuring mass flows velocities (line shifts) down to 2 km/s or better. LEMUR has been proposed to ESA as the European contribution to the Solar C mission.

preprint2011arXiv

Underflight calibration of SOHO/CDS and Hinode/EIS with EUNIS-07

Flights of Goddard Space Flight Center's Extreme-Ultraviolet Normal-Incidence Spectrograph (EUNIS) sounding rocket in 2006 and 2007 provided updated radiometric calibrations for SOHO/CDS and Hinode/EIS. EUNIS carried two independent imaging spectrographs covering wavebands of 300-370 A in first order and 170-205 A in second order. After each flight, end-to-end radiometric calibrations of the rocket payload were carried out in the same facility used for pre-launch calibrations of CDS and EIS. During the 2007 flight, EUNIS, SOHO CDS and Hinode EIS observed the same solar locations, allowing the EUNIS calibrations to be directly applied to both CDS and EIS. The measured CDS NIS 1 line intensities calibrated with the standard (version 4) responsivities with the standard long-term corrections are found to be too low by a factor of 1.5 due to the decrease in responsivity. The EIS calibration update is performed in two ways. One is using the direct calibration transfer of the calibrated EUNIS-07 short wavelength (SW) channel. The other is using the insensitive line pairs, in which one member was observed by EUNIS-07 long wavelength (LW) channel and the other by EIS in either LW or SW waveband. Measurements from both methods are in good agreement, and confirm (within the measurement uncertainties) the EIS responsivity measured directly before the instrument's launch. The measurements also suggest that the EIS responsivity decreased by a factor of about 1.2 after the first year of operation. The shape of the EIS SW response curve obtained by EUNIS-07 is consistent with the one measured in laboratory prior to launch. The absolute value of the quiet-Sun He II 304 A intensity measured by EUNIS-07 is consistent with the radiance measured by CDS NIS in quiet regions near the disk center and the solar minimum irradiance obtained by CDS NIS and SDO/EVE recently.

preprint2009arXiv

Absolute radiometric calibration of the EUNIS-06 170-205 A channel and calibration update for CDS/NIS

The Extreme-Ultraviolet Normal-Incidence Spectrograph sounding-rocket payload was flown on 2006 April 12 (EUNIS-06), carrying two independent imaging spectrographs covering wave bands of 300-370 A in first order and 170-205 A in second order, respectively. The absolute radiometric response of the EUNIS-06 long-wavelength (LW) channel was directly measured in the same facility used to calibrate CDS prior to the SOHO launch. Because the absolute calibration of the short-wavelength (SW) channel could not be obtained from the same lab configuration, we here present a technique to derive it using a combination of solar LW spectra and density- and temperature-insensitive line intensity ratios. The first step in this procedure is to use the coordinated, cospatial EUNIS and SOHO/CDS spectra to carry out an intensity calibration update for the CDS NIS-1 waveband, which shows that its efficiency has decreased by a factor about 1.7 compared to that of the previously implemented calibration. Then, theoretical insensitive line ratios obtained from CHIANTI allow us to determine absolute intensities of emission lines within the EUNIS SW bandpass from those of cospatial CDS/NIS-1 spectra after the EUNIS LW calibration correction. A total of 12 ratios derived from intensities of 5 CDS and 12 SW emission lines from Fe Fe X - Fe XIII yield an instrumental response curve for the EUNIS-06 SW channel that matches well to a relative calibration which relied on combining measurements of individual optical components. Taking into account all potential sources of error, we estimate that the EUNIS-06 SW absolute calibration is accurate to about 20%.