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.

In the standard model of solar flares, a large-scale reconnection current sheet is postulated as the central engine for powering the flare energy release and accelerating particles. However, where and how the energy release and particle acceleration occur remain unclear due to the lack of measurements for the magnetic properties of the current sheet. Here we report the measurement of spatially-resolved magnetic field and flare-accelerated relativistic electrons along a current-sheet feature in a solar flare. The measured magnetic field profile shows a local maximum where the reconnecting field lines of opposite polarities closely approach each other, known as the reconnection X point. The measurements also reveal a local minimum near the bottom of the current sheet above the flare loop-top, referred to as a "magnetic bottle". This spatial structure agrees with theoretical predictions and numerical modeling results. A strong reconnection electric field of ~4000 V/m is inferred near the X point. This location, however, shows a local depletion of microwave-emitting relativistic electrons. These electrons concentrate instead at or near the magnetic bottle structure, where more than 99% of them reside at each instant. Our observations suggest that the loop-top magnetic bottle is likely the primary site for accelerating and/or confining the relativistic electrons.

Nonthermal electrons accelerated in solar flares produce electromagnetic emission in two distinct, highly complementary domains - hard X-rays (HXRs) and microwaves (MWs). This paper reports MW imaging spectroscopy observations from the Expanded Owens Valley Solar Array of an M1.2 flare that occurred on 2017 September 9, from which we deduce evolving coronal parameter maps. We analyze these data jointly with the complementary Reuven Ramaty High-Energy Solar Spectroscopic Imager HXR data to reveal the spatially-resolved evolution of the nonthermal electrons in the flaring volume. We find that the high-energy portion of the nonthermal electron distribution, responsible for the MW emission, displays a much more prominent evolution (in the form of strong spectral hardening) than the low-energy portion, responsible for the HXR emission. We show that the revealed trends are consistent with a single electron population evolving according to a simplified trap-plus-precipitation model with sustained injection/acceleration of nonthermal electrons, which produces a double-powerlaw with steadily increasing break energy.

We investigate the performance of the generalized Spectral Kurtosis (SK) estimator in detecting and discriminating natural and artificial very short duration transients in the 2-bit sampling time domain Very-Long-Baseline Interferometry (VLBI) data. We demonstrate that, while both types of transients may be efficiently detected, their natural or artificial nature cannot be distinguished if only a time domain SK analysis is performed. However, these two types of transients become distinguishable from each other in the spectral domain, after a 32-bit FFT operation is performed on the 2-bit time domain voltages. We discuss the implication of these findings on the ability of the Spectral Kurtosis estimator to automatically detect bright astronomical transient signals of interests -- such as pulsar or fast radio bursts (FRB) -- in VLBI data streams that have been severely contaminated by unwanted radio frequency interference.

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.

Energy stored in the magnetic field in the solar atmosphere above active regions is a key driver of all solar activity (e.g., solar flares and coronal mass ejections), some of which can affect life on Earth. Radio observations provide a unique diagnostic of the coronal magnetic fields that make them a critical tool for the study of these phenomena, using the technique of broadband radio imaging spectropolarimetry. Observations with the ngVLA will provide unique observations of coronal magnetic fields and their evolution, key inputs and constraints for MHD numerical models of the solar atmosphere and eruptive processes, and a key link between lower layers of the solar atmosphere and the heliosphere. In doing so they will also provide practical "research to operations" guidance for space weather forecasting.

The spectrum of gyrosynchrotron emission from solar flares generally peaks in the microwave range. Its optically-thin, high-frequency component, above the spectral peak, is often used for diagnostics of the nonthermal electrons and the magnetic field in the radio source. Under favorable conditions, its low-frequency counterpart brings additional, complementary information about these parameters as well as thermal plasma diagnostics, either through gyrosynchrotron self-absorption, free-free absorption by the thermal plasma, or the suppression of emission through the so-called Razin effect. However, their effects on the low-frequency spectrum are often masked by spatial nonuniformity. To disentangle the various contributions to low-frequency gyrosynchrotron emission, a combination of spectral and imaging data is needed. To this end, we have investigated Owens Valley Solar Array (OVSA) multi-frequency images for 26 solar bursts observed jointly with Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) during the first half of 2002. For each, we examined dynamic spectra, time- and frequency-synthesis maps, RHESSI images with overlaid OVSA contours, and a few representative single-frequency snapshot OVSA images. We focus on the frequency dependence of microwave source sizes derived from the OVSA images and their effect on the low-frequency microwave spectral slope. We succeed in categorizing 18 analyzed events into several groups. Four events demonstrate clear evidence of being dominated by gyrosynchrotron self-absorption, with an inferred brightness temperature of $\geq10^8$~K. The low-frequency spectra in the remaining events are affected to varying degree by Razin suppression. We find that many radio sources are rather large at low frequencies, which can have important implications for solar energetic particle production and escape.

We report the first science results from the newly completed Expanded Owens Valley Solar Array (EOVSA), which obtained excellent microwave imaging spectroscopy observations of SOL2017-09-10, a classic partially-occulted solar limb flare associated with an erupting flux rope. This event is also well-covered by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in hard X-rays (HXRs). We present an overview of this event focusing on microwave and HXR data, both associated with high-energy nonthermal electrons, and discuss them within the context of the flare geometry and evolution revealed by extreme ultraviolet (EUV) observations from the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory (SDO/AIA). The EOVSA and RHESSI data reveal the evolving spatial and energy distribution of high-energy electrons throughout the entire flaring region. The results suggest that the microwave and HXR sources largely arise from a common nonthermal electron population, although the microwave imaging spectroscopy provides information over a much larger volume of the corona.

Understanding nonthermal particle generation, transport, and escape in solar flares requires detailed quantification of the particle evolution in the realistic 3D domain where the flare takes place. Rather surprisingly, apart of standard flare scenario and integral characteristics of the nonthermal electrons, not much is known about actual evolution of nonthermal electrons in the 3D spatial domain. This paper attempts to begin to remedy this situation by creating sets of evolving 3D models, the synthesized emission from which matches the evolving observed emission. Here we investigate two contrasting flares: a dense, "coronal-thick-target" flare SOL2002-04-12T17:42, that contained a single flare loop observed in both microwave and X-ray, and a more complex flare, SOL2015-06-22T17:50, that contained at least four distinct flaring loops needed to consistently reproduce the microwave and X-ray emission. Our analysis reveals differing evolution pattern of the nonthermal electrons in the dense and tenuous loops; however, both of which imply the central role of resonant wave-particle interaction with turbulence. These results offer new constraints for theory and models of the particle acceleration and transport in solar flares.

The well-established notion of a "common population" of the accelerated electrons simultaneously producing the hard X-ray (HXR) and the microwave (MW) emission during the flare impulsive phase has been challenged by some studies reporting the discrepancies between the HXR-inferred and the MW-inferred electron energy spectra. The traditional methods of their spectral inversion have some problems that can be mainly attributed to the unrealistic and the oversimplified treatment of the flare emission. To properly address this problem, we use a Non-linear Force Free Field (NLFFF) model extrapolated from an observed photospheric magnetogram as input to the three-dimensional, multi-wavelength modeling platform GX Simulator, and create a unified electron population model that can simultaneously reproduce the observed HXR and MW observations. We model the end of the impulsive phase of the 2015-06-22 M6.5 flare, and constrain the modeled electron spatial and energy parameters using observations made by the highest-resolving instruments currently available in two wavelengths, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) for HXR and the Expanded Owens Valley Solar Array (EOVSA) for MW. Our results suggest that the HXR-emitting electron population model fits the standard flare model with a broken power-law spectrum (E_break ~ 200 keV) that simultaneously produces the HXR footpoint emission and the MW high frequency emission. The model also includes an "HXR invisible" population of nonthermal electrons that are trapped in a large volume of magnetic field above the HXR-emitting loops, which is observable by its gyrosynchrotron (GS) radiation emitting mainly in MW low frequency range.

Ever-increasing multi-frequency imaging of solar observations suggests that solar flares often involve more than one magnetic fluxtube. Some of the fluxtubes are closed, while others can contain open field. The relative proportion of nonthermal electrons among those distinct loops is highly important for understanding the energy release, particle acceleration, and transport. The access of nonthermal electrons to the open field is further important as the open field facilitates the solar energetic particle (SEP) escape from the flaring site, and thus controls the SEP fluxes in the solar system, both directly and as seed particles for further acceleration. The large-scale fluxtubes are often filled with a tenuous plasma, which is difficult to detect in either EUV or X-ray wavelengths; however, they can dominate at low radio frequencies, where a modest component of nonthermal electrons can render the source optically thick and, thus, bright enough to be observed. Here we report detection of a large-scale `plume' at the impulsive phase of an X-class solar flare, SOL2001-08-25T16:23, using multi-frequency radio data from Owens Valley Solar Array. To quantify the flare spatial structure, we employ 3D modeling utilizing force-free-field extrapolations from the line-of-sight SOHO/MDI magnetograms with our modeling tool GX Simulator. We found that a significant fraction of the nonthermal electrons accelerated at the flare site low in the corona escapes to the plume, which contains both closed and open field. We propose that the proportion between the closed and open field at the plume is what determines the SEP population escaping into interplanetary space.

Solar flares are complicated physical phenomena that are observable in a broad range of the electromagnetic spectrum, from radiowaves to $\gamma$-rays. For a more comprehensive understanding of flares, it is necessary to perform a combined multi-wavelength analysis using observations from many satellites and ground-based observatories. For efficient data search, integration of different flare lists and representation of observational data, we have developed an Interactive Multi-Instrument Database of Solar Flares (https://solarflare.njit.edu/). The web accessible database is fully functional and allows the user to search for uniquely-identified flare events based on their physical descriptors and availability of observations by a particular set of instruments. Currently, the data from three primary flare lists (GOES, RHESSI and HEK) and a variety of other event catalogs (Hinode, Fermi GBM, Konus-Wind, OVSA flare catalogs, CACTus CME catalog, Filament eruption catalog) and observing logs (IRIS and Nobeyama coverage) are integrated, and an additional set of physical descriptors (temperature and emission measure) is provided along with an observing summary, data links, and multi-wavelength light curves for each flare event since January, 2002. We envision that this new tool will allow researchers to significantly speed up the search of events of interest for statistical and case studies.

We describe in general terms the practical use in astronomy of a higher-order statistical quantity called Spectral Kurtosis (SK), and describe the first implementation of SK-enabled firmware in the F-engine (Fourier transform-engine) of a digital FX correlator for Expanded Owens Valley Solar Array (EOVSA). The development of the theory for SK is summarized, leading to an expression for generalized SK that is applicable to both SK spectrometers and those not specifically designed for SK. We also give the means for computing both the SK estimator and thresholds for its application as a discriminator of RFI contamination. Tests of the performance of EOVSA as an SK spectrometer are shown to agree precisely with theoretical expectations, and the methods for configuring the correlator for correct SK operation are described.

Recently, in a few case studies we demonstrated that gyrosynchrotron microwave emission can be detected directly from the acceleration region when the trapped electron component is insignificant. For the statistical study reported here, we have identified events with steep (narrowband) microwave spectra that do not show a significant trapped component and at the same time show evidence of source uniformity, which simplifies the data analysis greatly. Initially, we identified a subset of more than 20 radio bursts with such narrow spectra, having low- and high-frequency spectral indices larger than 3 in absolute value. A steep low-frequency spectrum implies that the emission is nonthermal (for optically-thick thermal emission, the spectral index cannot be steeper than 2), and the source is reasonably dense and uniform. A steep high-frequency spectrum implies that no significant electron trapping occurs; otherwise a progressive spectral flattening would be observed. Roughly half of these radio bursts have RHESSI data, which allows for detailed, joint diagnostics of the source parameters and evolution. Based on an analysis of radio-to-X-ray spatial relationships, timing, and spectral fits, we conclude that the microwave emission in these narrowband bursts originates directly from the acceleration regions, which have relatively strong magnetic field, high density, and low temperature. In contrast, the thermal X-ray emission comes from a distinct loop with smaller magnetic field, lower density, but higher temperature. Therefore, these flares occurred likely due to interaction between two (or more) magnetic loops.

We obtain analytical approximations for the expectation and variance of the Spectral Kurtosis estimator in the case of Gaussian and coherent transient time domain signals mixed with a quasi-stationary Gaussian background, which are suitable for practical estimations of their signal-to-noise ratio and duty-cycle relative to the instrumental integration time. We validate these analytical approximations by means of numerical simulations and demonstrate that such estimates are affected by statistical uncertainties that, for a suitable choice of the integration time, may not exceed a few percent. Based on these analytical results, we suggest a multiscale Spectral Kurtosis spectrometer design optimized for real-time detection of transient signals, automatic discrimination based on their statistical signature, and measurement of their properties.

We present detailed 3D modeling of a dense, coronal thick target X-ray flare using the GX Simulator tool, photospheric magnetic measurements, and microwave imaging and spectroscopy data. The developed model offers a remarkable agreement between the synthesized and observed spectra and images in both X-ray and microwave domains, which validates the entire model. The flaring loop parameters are chosen to reproduce the emission measure, temperature, and the nonthermal electron distribution at low energies derived from the X-ray spectral fit, while the remaining parameters, unconstrained by the X-ray data, are selected such as to match the microwave images and total power spectra. The modeling suggests that the accelerated electrons are trapped in the coronal part of the flaring loop, but away from where the magnetic field is minimal, and, thus, demonstrates that the data are clearly inconsistent with electron magnetic trapping in the weak diffusion regime mediated by the Coulomb collisions. Thus, the modeling supports the interpretation of the coronal thick-target sources as sites of electron acceleration in flares and supplies us with a realistic 3D model with physical parameters of the acceleration region and flaring loop.

We developed an efficient algorithm integrated in our 3D modeling tool, GX Simulator (Nita et al. 2015), allowing quick computation of the synthetic intensity and polarization maps of solar active regions (AR) in the ALMA spectral range. The algorithm analyzes the photospheric input (white light and magnetogram) to classify a given photospheric pixel to belong to a given photospheric structure. Then, a 1D chromospheric model (Fontenla et al. 2009) is added on top of each pixel, which forms a chromospheric model of the AR. Next step is computation of the mm and sub-mm emission produced from this chromosphere model. A huge advantage of this approach is that emission from any given AR can be synthesized very fast, on the order of a few minutes after the AR selection. Using the GX Simulator tool it is also possible to produce synthetic maps of the microwave (gyroresonance) and EUV emission from the same AR model and compare them with the ALMA synthetic maps and with the corresponding observed microwave and/or EUV data.

This paper presents a solely thermal flare, which we detected in the microwave range from the thermal gyro- and free-free emission it produced. An advantage of analyzing thermal gyro emission is its unique ability to precisely yield the magnetic field in the radiating volume. When combined with observationally-deduced plasma density and temperature, these magnetic field measurements offer a straightforward way of tracking evolution of the magnetic and thermal energies in the flare. For the event described here, the magnetic energy density in the radio-emitting volume declines over the flare rise phase, then stays roughly constant during the extended peak phase, but recovers to the original level over the decay phase. At the stage where the magnetic energy density decreases, the thermal energy density increases; however, this increase is insufficient, by roughly an order of magnitude, to compensate for the magnetic energy decrease. When the magnetic energy release is over, the source parameters come back to nearly their original values. We discuss possible scenarios to explain this behavior.

We have undertaken a major enhancement of our IDL-based simulation tools developed earlier for modeling microwave and X-ray emission. The object-based architecture provides an interactive graphical user interface that allows the user to import photospheric magnetic field maps and perform magnetic field extrapolations to almost instantly generate 3D magnetic field models, to investigate the magnetic topology of these models by interactively creating magnetic field lines and associated magnetic flux tubes, to populate the flux tubes with user-defined nonuniform thermal plasma and anisotropic, nonuniform, nonthermal electron distributions; to investigate the spatial and spectral properties of radio and X-ray emission calculated from the model, and to compare the model-derived images and spectra with observational data. The application integrates shared-object libraries containing fast gyrosynchrotron emission codes developed in FORTRAN and C++, soft and hard X-ray codes developed in IDL, a FORTRAN-based potential-field extrapolation routine and an IDL-based linear force free field extrapolation routine. The interactive interface allows users to add any user-defined radiation code that adheres to our interface standards, as well as user-defined magnetic field extrapolation routines. Here we use this tool to analyze a simple single-loop flare and use the model to constrain the 3D structure of the magnetic flaring loop and 3D spatial distribution of the fast electrons inside this loop. We iteratively compute multi-frequency microwave and multi-energy X-ray images from realistic magnetic fluxtubes obtained from an extrapolation of a magnetogram taken prior to the flare, and compare them with imaging data obtained by SDO, NoRH, and RHESSI instruments. We use this event to illustrate use of the tool for general interpretation of solar flares to address disparate problems in solar physics.

Spectra derived from fast Fourier transform (FFT) analysis of time-domain data intrinsically contain statistical fluctuations whose distribution depends on the number of accumulated spectra contributing to a measurement. The tail of this distribution, which is essential for separation of the true signal from the statistical fluctuations, deviates noticeably from the normal distribution for a finite number of the accumulations. In this paper we develop a theory to properly account for the statistical fluctuations when fitting a model to a given accumulated spectrum. The method is implemented in software for the purpose of automatically fitting a large body of such FFT-derived spectra. We apply this tool to analyze a portion of a dense cluster of spikes recorded by our FST instrument during a record-breaking event that occurred on 06 Dec 2006. The outcome of this analysis is briefly discussed.

Based on detailed analysis of radio and X-ray observations of a flare on 2002 April 11 augmented by realistic 3D modeling, we have identified a radio emission component produced directly at the flare acceleration region. This acceleration region radio component has distinctly different (i) spectrum, (ii) light curves, (iii) spatial location, and, thus, (iv) physical parameters from those of the separately identified, trapped or precipitating electron components. To derive evolution of physical parameters of the radio sources we apply forward fitting of the radio spectrum time sequence with the gyrosynchrotron source function with 5 to 6 free parameters. At the stage when the contribution from the acceleration region dominates the radio spectrum, the X-ray- and radio- derived electron energy spectral indices agree well with each other. During this time the maximum energy of the accelerated electron spectrum displays a monotonic increase with time from ~ 300 keV to ~ 2 MeV over roughly one minute duration indicative of an acceleration process in the form of growth of the power-law tail; the fast electron residence time in the acceleration region is about 2 - 4 s, which is much longer than the time of flight and so requires a strong diffusion mode there to inhibit free-streaming propagation. The acceleration region has a relatively strong magnetic field, B ~ 120 G, and a low thermal density, less than 2x10^9 cm^-3. These acceleration region properties are consistent with a stochastic acceleration mechanism.

We have developed a general framework for modeling gyrosynchrotron and free-free emission from solar flaring loops and used it to test the premise that 2D maps of source parameters, particularly magnetic field, can be deduced from spatially resolved microwave spectropolarimetry data. In this paper we show quantitative results for a flaring loop with a realistic magnetic geometry, derived from a magnetic field extrapolation, and containing an electron distribution with typical thermal and nonthermal parameters, after folding through the instrumental profile of a realistic interferometric array. We compare the parameters generated from forward fitting a homogeneous source model to each line of sight through the folded image data cube with both the original parameters used in the model and with parameters generated from forward fitting a homogeneous source model to the original (unfolded) image data cube. We find excellent agreement in general, but with systematic effects that can be understood as due to finite resolution in the folded images and the variation of parameters along the line of sight, which are ignored in the homogeneous source model. We discuss the use of such 2D parameter maps within a larger framework of 3D modeling, and the prospects for applying these methods to data from a new generation of multifrequency radio arrays now or soon to be available.

Microwave emission of solar flares is formed primarily by incoherent gyrosynchrotron radiation generated by accelerated electrons in coronal magnetic loops. The resulting emission depends on many factors, including pitch-angle distribution of the emitting electrons and the source geometry. In this work, we perform systematic simulations of solar microwave emission using recently developed tools (GS Simulator and fast gyrosynchrotron codes) capable of simulating maps of radio brightness and polarization as well as spatially resolved emission spectra. A 3D model of a symmetric dipole magnetic loop is used. We compare the emission from isotropic and anisotropic (of loss-cone type) electron distributions. We also investigate effects caused by inhomogeneous distribution of the emitting particles along the loop. It is found that effect of the adopted moderate electron anisotropy is the most pronounced near the footpoints and it also depends strongly on the loop orientation. Concentration of the emitting particles at the loop top results in a corresponding spatial shift of the radio brightness peak, thus reducing effects of the anisotropy. The high-frequency (around 50 GHz) emission spectral index is specified mainly by the energy spectrum of the emitting electrons; however, at intermediate frequencies (around 10-20 GHz), the spectrum shape is strongly dependent on the electron anisotropy, spatial distribution, and magnetic field nonuniformity. The implications of the obtained results for the diagnostics of the energetic electrons in solar flares are discussed.

We use rotation stereoscopy to estimate the height of a steady-state solar feature relative to the photosphere, based on its apparent motion in the image plane recorded over several days of observation. The stereoscopy algorithm is adapted to work with either one- or two-dimensional data (i.e. from images or from observations that record the projected position of the source along an arbitrary axis). The accuracy of the algorithm is tested on simulated data, and then the algorithm is used to estimate the coronal radio source heights associated with the active region NOAA 10956, based on multifrequency imaging data over 7 days from the Siberian Solar Radio Telescope near 5.7 GHz, the Nobeyama Radio Heliograph at 17 GHz, as well as one-dimensional scans at multiple frequencies spanning the 5.98--15.95 GHz frequency range from the RATAN-600 instrument. The gyroresonance emission mechanism, which is sensitive to the coronal magnetic field strength, is applied to convert the estimated radio source heights at various frequencies, h(f), to information about magnetic field vs. height B(h), and the results are compared to a magnetic field extrapolation derived from photospheric magnetic field observations obtained by Hinode and MDI. We found that the gyroresonant emission comes from the heights exceeding location of the third gyrolayer irrespectively on the magnetic extrapolation method; implications of this finding for the coronal magnetography and coronal plasma physics are discussed.

We report the observation of an unusual cold, tenuous solar flare, which reveals itself via numerous and prominent non-thermal manifestations, while lacking any noticeable thermal emission signature. RHESSI hard X-rays and 0.1-18 GHz radio data from OVSA and Phoenix-2 show copious electron acceleration (10^35 electrons per second above 10 keV) typical for GOES M-class flares with electrons energies up to 100 keV, but GOES temperatures not exceeding 6.1 MK. The imaging, temporal, and spectral characteristics of the flare have led us to a firm conclusion that the bulk of the microwave continuum emission from this flare was produced directly in the acceleration region. The implications of this finding for the flaring energy release and particle acceleration are discussed.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

The coming years will see routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness in the wavelength domain. To capitalize on the soon to be available radio facilities such as the expanded OVSA, SSRT and FASR, and the challenges they present in the visualization and synthesis of the multi-frequency datasets, we propose that realistic, sophisticated 3D active region and flare modeling is timely now and will be a forefront of coronal studies over the coming years. Here we summarize our 3D modeling efforts, aimed at forward fitting of imaging spectroscopy data, and describe currently available 3D modeling tools. We also discuss plans for future generalization of our modeling tools.

Due to its conceptual simplicity and its proven effectiveness in real-time detection and removal of radio frequency interference (RFI) from radio astronomy data, the Spectral Kurtosis (SK) estimator is likely to become a standard tool of a new generation of radio telescopes. However, the SK estimator in its original form must be developed from instantaneous power spectral density (PSD) estimates, and hence cannot be employed as an RFI excision tool downstream of the data pipeline in existing instruments where any time averaging is performed. In this letter, we develop a generalized estimator with wider applicability for both instantaneous and averaged spectral data, which extends its practical use to a much larger pool of radio instruments.

We develop a practical forward fitting method based on the SIMPLEX algorithm with shaking, which allows the derivation of the magnetic field and other parameters along a solar flaring loop using microwave imaging spectroscopy of gyrosynchrotron emission. We illustrate the method using a model loop with spatially varying magnetic field, filled with uniform ambient density and an evenly distributed fast electron population with an isotropic, power-law energy distribution.

Modern observations and models of various astrophysical objects suggest that many of their physical parameters fluctuate substantially at different spatial scales. The rich variety of the emission processes, including Transition Radiation but not limited to it, arising in such turbulent media constitutes the scope of Stochastic Theory of Radiation. We review general approaches applied in the stochastic theory of radiation and specific methods used to calculate the transition radiation produced by fast particles in the magnetized randomly inhomogeneous plasma. The importance of the theory of transition radiation for astrophysics is illustrated by one example of its detailed application to a solar radio burst, including specially designed algorithms of the spectral forward fitting.

Microturbulence, i.e. enhanced fluctuations of plasma density, electric and magnetic fields, is of great interest in astrophysical plasmas, but occurs on spatial scales far too small to resolve by remote sensing, e.g., at ~ 1-100 cm in the solar corona. This paper reports spatially resolved observations that offer strong support for the presence in solar flares of a suspected radio emission mechanism, resonant transition radiation, which is tightly coupled to the level of microturbulence and provides direct diagnostics of the existence and level of fluctuations on decimeter spatial scales. Although the level of the microturbulence derived from the radio data is not particularly high, <\Delta n^2 >/n^2 ~ 10^{-5}$, it is large enough to affect the charged particle diffusion and give rise to effective stochastic acceleration. This finding has exceptionally broad astrophysical implications since modern sophisticated numerical models predict generation of much stronger turbulence in relativistic objects, e.g., in gamma-ray burst sources.