Search

The interaction between a Coronal Mass Ejection (CME) and a comet has been observed several times by in-situ observations from the Rosetta Plasma Consortium (RPC), which is designed to investigate the cometary magnetosphere of comet 67P/Churyumov-Gerasimenko (CG). Goetz et al. (2019) reported a magnetic field of up to 300 nT measured in the inner coma, which is among the largest interplanetary magnetic fields observed in the solar system. They suggested the large magnetic field observations in the inner coma come from magnetic field pile-up regions, which are generated by the interaction between a CME and/or corotating interaction region and the cometary magnetosphere. However, the detailed interaction between a CME and the cometary magnetosphere of comet CG in the inner coma has not been investigated by numerical simulations yet. In this manuscript, we will use a numerical model to simulate the interaction between comet CG and a Halloween class CME and investigate its magnetospheric response to the CME. We find that the plasma structures change significantly during the CME event, and the maximum value of the magnetic field strength is more than 500nT close to the nucleus. Virtual satellites at similar distances as Rosetta show that the magnetic field strength can be as large as 250nT, which is slightly less than what Goetz et al. (2019) reported., Comment: 19 pages, 5 figures, submitted to ApJ

ESA's Rosetta spacecraft at comet 67P/Churyumov-Gerasimenko (67P) was the first mission that accompanied a comet over a substantial fraction of its orbit. On board was the ROSINA mass spectrometer suite to measure the local densities of the volatile species sublimating from the ices inside the comet's nucleus. Understanding the nature of these ices was a key goal of Rosetta. We analyzed the primary cometary molecules at 67P, namely H$_2$O and CO$_2$, together with a suite of minor species for almost the entire mission. Our investigation reveals that the local abundances of highly volatile species, such as CH$_4$ and CO, are reproduced by a linear combination of both H$_2$O and CO$_2$ densities. These findings bear similarities to laboratory-based temperature programmed desorption experiments of amorphous ices and imply that highly volatile species are trapped in H$_2$O and CO$_2$ ices. Our results do not show the presence of ices dominated by these highly volatile molecules. Most likely, they were lost due to thermal processing of 67P's interior prior to its deflection to the inner solar system. Deviations in the proportions co-released with H$_2$O and CO$_2$ can only be observed before the inbound equinox, when the comet was still far from the sun and the abundance of highly volatile molecules associated with CO$_2$ outgassing were lower. The corresponding CO$_2$ is likely seasonal frost, which sublimated and lost its trapped highly volatile species before re-freezing during the previous apparition. CO, on the other hand, was elevated during the same time and requires further investigation., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record is available online at: https://doi.org/10.1093/mnras/stad3005

We provide exact analytical solutions for the magnetic field produced by prescribed current distributions located inside a toroidal filament of finite thickness. The solutions are expressed in terms of toroidal functions which are modifications of the Legendre functions. In application to the MHD equilibrium of a twisted toroidal current loop in the solar corona, the Grad-Shafranov equation is decomposed into an analytic solution describing an equilibrium configuration against the pinch-effect from its own current and an approximate solution for an external strapping field to balance the hoop force. Our solutions can be employed in numerical simulations of coronal mass ejections. When superimposed on the background solar coronal magnetic field, the excess magnetic energy of the twisted current loop configuration can be made unstable by applying flux cancellation to reduce the strapping field. Such loss of stability accompanied by the formation of an expanding flux rope is typical for the Titov & D\'emoulin (1999) eruptive event generator. The main new features of the proposed model are: (i) The filament is filled with finite $\beta$ plasma with finite mass and energy, (ii) The model describes an equilibrium solution that will spontaneously erupt due to magnetic reconnection of the strapping magnetic field arcade, and (iii) There are analytic expressions connecting the model parameters to the asymptotic velocity and total mass of the resulting CME, providing a way to connect the simulated CME properties to multipoint coronograph observations., Comment: 20 pages, 7 figures

We investigate the cross-calibration of the Hinode/SOT-SP and SDO/HMI instrument meta-data, specifically the correspondence of the scaling and pointing information. Accurate calibration of these datasets gives the correspondence needed by inter-instrument studies and learning-based magnetogram systems, and is required for physically-meaningful photospheric magnetic field vectors. We approach the problem by robustly fitting geometric models on correspondences between images from each instrument's pipeline. This technique is common in computer vision, but several critical details are required when using scanning slit spectrograph data like Hinode/SOT-SP. We apply this technique to data spanning a decade of the Hinode mission. Our results suggest corrections to the published Level 2 Hinode/SOT-SP data. First, an analysis on approximately 2,700 scans suggests that the reported pixel size in Hinode/SOT-SP Level 2 data is incorrect by around 1%. Second, analysis of over 12,000 scans show that the pointing information is often incorrect by dozens of arcseconds with a strong bias. Regression of these corrections indicates that thermal effects have caused secular and cyclic drift in Hinode/SOT-SP pointing data over its mission. We offer two solutions. First, direct co-alignment with SDO/HMI data via our procedure can improve alignments for many Hinode/SOT-SP scans. Second, since the pointing errors are predictable, simple post-hoc corrections can substantially improve the pointing. We conclude by illustrating the impact of this updated calibration on derived physical data products needed for research and interpretation. Among other things, our results suggest that the pointing errors induce a hemispheric bias in estimates of radial current density., Comment: Under revisions at ApJS

Results. We found that gas-phase silicon was present throughout the Rosetta mission. Furthermore, the presence of sodium and iron atoms near the comet's perihelion confirms that sputtering cannot be the sole release process for refractory elements into the gas phase. Nickel was found to be below the detection limit. The search for parent species of any of the identified gas phase refractories has not been successful. Upper limits for a suite of possible fragment species (SiH, SiC, NaH, etc.) of larger parent and daughter species have been obtained. Furthermore, Si did not exhibit the same drop in signal as do common cometary gases when the spacecraft is pointed away from the nucleus. The combined results suggest that a direct release of elemental species from small grains on the surface of the nucleus or from small grains in the surrounding coma is a more likely explanation than the previous assumption of release via the dissociation of gaseous parent molecules., Comment: Accepted for publication in A&A

We present a reduced magnetohydrodynamic (MHD) mathematical model describing the dynamical behavior of highly conducting plasmas with frozen-in magnetic fields, constrained by the assumption that, there exists a frame of reference, where the magnetic field vector, $\mathbf{B}$, is aligned with the plasma velocity vector, $\mathbf{u}$, at each point. We call this solution "stream-aligned MHD" (SA-MHD). Within the framework of this model, the electric field, $\mathbf{E} = -\mathbf{u} \times \mathbf{B} \equiv 0$, in the induction equation vanishes identically and so does the electromagnetic energy flux (Poynting flux), $\mathbf{E}\times\mathbf{B}\equiv0$, in the energy equation. At the same time, the force effect from the magnetic field on the plasma motion (the Ampere force) is fully taken into account in the momentum equation. Any steady-state solution of the proposed model is a legitimate solution of the full MHD system of equations. However, the converse statement is not true: in an arbitrary steady-state magnetic field the electric field does not have to vanish identically (its curl has to, though). Specifically, realistic tree-dimensional solutions for the steady-state ("ambient") solar atmosphere in the form of so-called Parker spirals, can be efficiently generated within the stream-aligned MHD (SA-MHD) with no loss in generality, Comment: 22 pages, 6 figures. Submitted to The Astrophysical Journal

Both NASA's Solar Dynamics Observatory (SDO) and the JAXA/NASA Hinode mission include spectropolarimetric instruments designed to measure the photospheric magnetic field. SDO's Helioseismic and Magnetic Imager (HMI) emphasizes full-disk high-cadence and good spatial resolution data acquisition while Hinode's Solar Optical Telescope Spectro-Polarimeter (SOT-SP) focuses on high spatial resolution and spectral sampling at the cost of a limited field of view and slower temporal cadence. This work introduces a deep-learning system named SynthIA (Synthetic Inversion Approximation), that can enhance both missions by capturing the best of each instrument's characteristics. We use SynthIA to produce a new magnetogram data product, SynodeP (Synthetic Hinode Pipeline), that mimics magnetograms from the higher spectral resolution Hinode/SOT-SP pipeline, but is derived from full-disk, high-cadence, and lower spectral-resolution SDO/HMI Stokes observations. Results on held-out data show that SynodeP has good agreement with the Hinode/SOT-SP pipeline inversions, including magnetic fill fraction, which is not provided by the current SDO/HMI pipeline. SynodeP further shows a reduction in the magnitude of the 24-hour oscillations present in the SDO/HMI data. To demonstrate SynthIA's generality, we show the use of SDO/AIA data and subsets of the HMI data as inputs, which enables trade-offs between fidelity to the Hinode/SOT-SP inversions, number of observations used, and temporal artifacts. We discuss possible generalizations of SynthIA and its implications for space weather modeling. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan under grant NASA 80NSSC20K0600E, and will be open-sourced.

MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environment Modeling (CSEM) at the University of Michigan, that developed and maintained the Space Weather Modeling Framework (SWMF) and its core element, the BATS-R-US extended MHD code. It took a quarter of a century to develop this capability and reach its present level of maturity that makes it suitable for research use by the space physics community through the Community Coordinated Modeling Center (CCMC) as well as operational use by the NOAA Space Weather Prediction Center (SWPC)., Comment: 105 pages, 36 figures, in press

The Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO) produces estimates of the photospheric magnetic field which are a critical input to many space weather modelling and forecasting systems. The magnetogram products produced by HMI and its analysis pipeline are the result of a per-pixel optimization that estimates solar atmospheric parameters and minimizes disagreement between a synthesized and observed Stokes vector. In this paper, we introduce a deep learning-based approach that can emulate the existing HMI pipeline results two orders of magnitude faster than the current pipeline algorithms. Our system is a U-Net trained on input Stokes vectors and their accompanying optimization-based VFISV inversions. We demonstrate that our system, once trained, can produce high-fidelity estimates of the magnetic field and kinematic and thermodynamic parameters while also producing meaningful confidence intervals. We additionally show that despite penalizing only per-pixel loss terms, our system is able to faithfully reproduce known systematic oscillations in full-disk statistics produced by the pipeline. This emulation system could serve as an initialization for the full Stokes inversion or as an ultra-fast proxy inversion. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan, under grant NASA 80NSSC20K0600E, and has been open sourced.

In 2004, the Rosetta spacecraft was sent to comet 67P/Churyumov-Gerasimenko for the first ever long-term investigation of a comet. After its arrival in 2014, the spacecraft spent more than two years in immediate proximity to the comet. During these two years, the ROSINA Double Focusing Mass Spectrometer (DFMS) onboard Rosetta discovered a coma with an unexpectedly complex chemical composition that included many oxygenated molecules. Determining the exact cometary composition is an essential first step to understanding of the organic rich chemistry in star forming regions and protoplanetary disks that are ultimately conserved in cometary ices. In this study a joint approach of laboratory calibration and space data analysis was used to perform a detailed identification and quantification of CHO-compounds in the coma of 67P/Churyumov-Gerasimenko. The goal was to derive the CHO-compound abundances relative to water for masses up to 100 u. For this study, the May 2015 post-equinox period represent the best bulk abundances of comet 67P/Churyumov-Gerasimenko. A wide variety of CHO-compounds were discovered and their bulk abundances were derived. Finally, these results are compared to abundances of CHO-bearing molecules in other comets, obtained mostly from ground-based observations and modelling., Comment: This document is the unedited Author's version of a Submitted Work that was subsequently accepted for publication in ACS Earth and Space Chemistry, copyright $\copyright$ American Chemical Society after peer review. To access the final edited and published work see https://pubs.acs.org/doi/10.1021/acsearthspacechem.9b00094

A deep learning network, Long-Short Term Memory (LSTM) network, is used in this work to predict whether the maximum flare class an active region (AR) will produce in the next 24 hours is class $\Gamma$. We considered $\Gamma$ are $\ge M$, $\ge C$ and any flare class. The essence of using LSTM, which is a recurrent neural network, is its capability to capture temporal information of the data samples. The input features are time sequences of 20 magnetic parameters from SHARPs - Space-weather HMI Active Region Patches. We analyzed active regions from June 2010 to Dec 2018, using the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalogs and label the data samples with identified ARs in the GOES X-ray flare catalogs. Our results (i) shows consistent skill scores with recently published results using LSTMs and better than the previous work using single time input (eg. DeFN) (ii) The skill scores from the model show essential differences when different years of data was chosen for training and testing.

Comets are considered to be some of the most pristine and unprocessed solar system objects accessible to in-situ exploration. Investigating their molecular and elemental composition takes us on a journey back to the early period of our solar system and possibly even further. In this work, we deduce the bulk abundances of the major volatile species in comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. The basis are measurements obtained with the ROSINA instrument suite on board the Rosetta orbiter during a suitable period of high outgassing near perihelion. The results are combined with both gas and dust composition measurements published in the literature. This provides an integrated inventory of the major elements present in the nucleus of 67P/Churyumov-Gerasimenko. Similar to comet 1P/Halley, which was visited by ESA's Giotto spacecraft in 1986, comet 67P/Churyumov-Gerasimenko also shows near-solar abundances of oxygen and carbon, whereas hydrogen and nitrogen are depleted compared to solar. Still, the degree of devolatilization is lower than that of inner solar system objects, including meteorites and the Earth. This supports the idea that comets are among the most pristine objects in our solar system., Comment: Accepted for publication in MNRAS

The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios., Comment: 10 pages, 6 figures. Submitted to MNRAS

We present several methods towards construction of precursors, which show great promise towards early predictions, of solar flare events in this paper. A data pre-processing pipeline is built to extract useful data from multiple sources, Geostationary Operational Environmental Satellites (GOES) and Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI), to prepare inputs for machine learning algorithms. Two classification models are presented: classification of flares from quiet times for active regions and classification of strong versus weak flare events. We adopt deep learning algorithms to capture both the spatial and temporal information from HMI magnetogram data. Effective feature extraction and feature selection with raw magnetogram data using deep learning and statistical algorithms enable us to train classification models to achieve almost as good performance as using active region parameters provided in HMI/Space-Weather HMI-Active Region Patch (SHARP) data files. Case studies show a significant increase in the prediction score around 20 hours before strong solar flare events.

The European Space Agency (ESA) spacecraft Rosetta accompanied the Jupiter-family comet (JFC) 67P / Churyumov-Gerasimenko for over two years along its trajectory through the inner solar system. Between 2014 and 2016, it performed almost continuous in-situ measurements of the comet's gaseous atmosphere in close proximity to its nucleus. In this study, the $^{16}$O/$^{18}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko, as measured by the ROSINA DFMS mass spectrometer on board Rosetta, was determined from the ratio of H${_2}^{16}$O / H${_2}^{18}$O and $^{16}$OH / $^{18}$OH. The value of 445 $\pm$ 35 represents an $\sim$ 11% enrichment of $^{18}$O compared with the terrestrial ratio of 498.7 $\pm$ 0.1. This cometary value is consistent with the comet containing primordial water, in accordance with leading self-shielding models. These models predict primordial water to be between 5% to 20% enriched in heavier oxygen isotopes compared to terrestrial water. Addendum: The $^{16}$O/$^{17}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko., Comment: 14 pages, 10 figures, 4 tables + Addendum. Accepted for publication in Astronomy and Astrophysics (A&A) Journal, Rosetta Special Issue 2019

The ROSINA mass spectrometer DFMS on board ESA's Rosetta spacecraft detected the major isotopes of the noble gases argon, krypton, and xenon in the coma of comet 67P/Churyumov-Gerasimenko. Earlier, it has been shown that xenon exhibits an isotopic composition distinct from anywhere else in the solar system. However, argon isotopes, within error, were shown to be consistent with solar isotope abundances. This discrepancy suggested an additional exotic component of xenon in comet 67P/Churyumov-Gerasimenko. Here we show that also krypton exhibits an isotopic composition close to solar. Furthermore, we found a depletion compared to solar of argon with respect to krypton and of krypton with respect to xenon, which is a necessity to postulate an addition of exotic xenon in the comet., Comment: http://advances.sciencemag.org/content/4/7/eaar6297

The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind models are described in detail. This historical perspective is also applied to the present suite of global corona and solar wind models. We discuss the evolution of new ideas and their implementation into numerical simulation codes. We point out the scientific and computational challenges facing these models and discuss the ways various groups tried to overcome these challenges. Next, we discuss the latest, state-of-the art models and point to the expected next steps in modeling the corona and the interplanetary medium., Comment: Accepted for publication, Living Reviews in Solar Physics, 2018

Coronal Mass Ejections (CMEs), a kind of energetic solar eruptions, are an integral subject of space weather research. Numerical magnetohydrodynamic (MHD) modeling, which requires powerful computational resources, is one of the primary means of studying the phenomenon. With increasing accessibility of such resources, grows the demand for user-friendly tools that would facilitate the process of simulating CMEs for scientific and operational purposes. The Eruptive Event Generator based on Gibson-Low flux rope (EEGGL), a new publicly available computational model presented in this paper, is an effort to meet this demand. EEGGL allows one to compute the parameters of a model flux rope driving a CME via an intuitive graphical user interface (GUI). We provide a brief overview of the physical principles behind EEGGL and its functionality. Ways towards future improvements of the tool are outlined.

We perform a three-dimensional (3D) global simulation of Earth's magnetosphere with kinetic reconnection physics to study the flux transfer events (FTEs) and dayside magnetic reconnection with the recently developed magnetohydrodynamics with embedded particle-in-cell model (MHD-EPIC). During the one-hour long simulation, the FTEs are generated quasi-periodically near the subsolar point and move toward the poles. We find the magnetic field signature of FTEs at their early formation stage is similar to a `crater FTE', which is characterized by a magnetic field strength dip at the FTE center. After the FTE core field grows to a significant value, it becomes an FTE with typical flux rope structure. When an FTE moves across the cusp, reconnection between the FTE field lines and the cusp field lines can dissipate the FTE. The kinetic features are also captured by our model. A crescent electron phase space distribution is found near the reconnection site. A similar distribution is found for ions at the location where the Larmor electric field appears. The lower hybrid drift instability (LHDI) along the current sheet direction also arises at the interface of magnetosheath and magnetosphere plasma. The LHDI electric field is about 8 mV/m and its dominant wavelength relative to the electron gyroradius agrees reasonably with MMS observations.

We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfv\'en wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfv\'en wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}

The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSoM, a global magnetohydrodynamic model driven by Alfv{\'e}n waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge states calculation covering all latitudes in a realistic magnetic field. The ratios $O^{+7}/O^{+6}$ and $C^{+6}/C^{+5}$ are compared to in-situ Ulysses observations, and are found to be higher in the slow wind, as observed; however, they are under-predicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to EIS observations above a coronal hole. The agreement is partial, and suggests that all ionization rates are under-predicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows it originates from coronal hole boundaries (CHB), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global, and observationally supported by EUV tomography., Comment: Submitted to the Astrophysical Journal

We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently., Comment: 26 pages, 5 figures; submitted to the Astrophysical Journal

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107., Comment: Submitted to the Astrophysical Journal

Magnetic loops are building blocks of the closed-field corona. While active region loops are readily seen in images taken at EUV and X-ray wavelengths, quiet Sun loops are seldom identifiable and therefore difficult to study on an individual basis. The first analysis of solar minimum (Carrington Rotation 2077) quiet Sun (QS) coronal loops utilizing a novel technique called the Michigan Loop Diagnostic Technique (MLDT) is presented. This technique combines Differential Emission Measure Tomography (DEMT) and a potential field source surface (PFSS) model, and consists of tracing PFSS field lines through the tomographic grid on which the Local Differential Emission Measure (LDEM) is determined. As a result, the electron temperature Te and density Ne at each point along each individual field line can be obtained. Using data from STEREO/EUVI and SOHO/MDI, the MLDT identifies two types of QS loops in the corona: so-called "up" loops in which the temperature increases with height, and so-called "down" loops in which the temperature decreases with height. Up loops are expected, however, down loops are a surprise, and furthermore, they are ubiquitous in the low-latitude corona. Up loops dominate the QS at higher latitudes. The MLDT allows independent determination of the empirical pressure and density scale heights, and the differences between the two remain to be explained. The down loops appear to be a newly discovered property of the solar minimum corona that may shed light on the physics of coronal heating. The results are shown to be robust to the calibration uncertainties of the EUVI instrument., Comment: Accepted for publication in The Astrophysical Journal, waiting for the full biblio info

We present a focused parameter study of solar wind - magnetosphere interaction for the young Sun and Earth, $~3.5$ Ga ago, that relies on magnetohydrodynamic (MHD) simulations for both the solar wind and the magnetosphere. By simulating the quiescent young Sun and its wind we are able to propagate the MHD simulations up to Earth's magnetosphere and obtain a physically realistic solar forcing of it. We assess how sensitive the young solar wind is to changes in the coronal base density, sunspot placement and magnetic field strength, dipole magnetic field strength and the Sun's rotation period. From this analysis we obtain a range of plausible solar wind conditions the paleomagnetosphere may have been subject to. Scaling relationships from the literature suggest that a young Sun would have had a mass flux different from the present Sun. We evaluate how the mass flux changes with the aforementioned factors and determine the importance of this and several other key solar and magnetospheric variables with respect to their impact on the paleomagnetosphere. We vary the solar wind speed, density, interplanetary magnetic field strength and orientation as well as Earth's dipole magnetic field strength and tilt in a number of steady-state scenarios that are representative of young Sun-Earth interaction. This study is done as a first step of a more comprehensive effort towards understanding the implications of Sun-Earth interaction for planetary atmospheric evolution., Comment: 16 pages, 7 figures

In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona., Comment: 17 pages, 11 figures, Submitted to ApJ on 12/08/2009