The evolution of magnetic helicity has a close relationship with solar eruptions and is of interest as a predictive diagnostic. In this case study, we analyse the evolution of the normalised emergence, shearing and total magnetic helicity components in the case of three flaring and three non-flaring active regions (ARs) using SHARPs (Spaceweather Helioseismic Magnetic Imager Active Region Patches) vector magnetic field data. The evolution of the three magnetic helicity components is analysed with wavelet transforms, revealing significant common periodicities of the normalised emergence, shearing and total helicity fluxes before flares in the flaring ARs. The three non-flaring ARs do not show such common periodic behaviour. This case study suggests that the presence of significant periodicities in the power spectrum of magnetic helicity components could serve as a valuable precursor for flares.

Flare ribbons are bright manifestations of flare energy dissipation in the lower solar atmosphere. For the first time, we report on high-resolution imaging spectroscopy observations of flare ribbons situated off-limb in the H$\beta$ and Ca II 8542 {\AA} lines and make a detailed comparison with radiative hydrodynamic simulations. Observations of the X8.2-class solar flare SOL2017-09-10T16:06 UT obtained with the Swedish Solar Telescope reveal bright horizontal emission layers in H$\beta$ line wing images located near the footpoints of the flare loops. The apparent separation between the ribbon observed in the H$\beta$ wing and the nominal photospheric limb is about 300 - 500 km. The Ca II 8542 {\AA} line wing images show much fainter ribbon emissions located right on the edge of the limb, without clear separation from the limb. RADYN models are used to investigate synthetic spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. The simulations show that, towards the limb, where the line of sight is substantially oblique with respect to the vertical direction, the flaring atmosphere model reproduces the high contrast of the off-limb H$\beta$ ribbons and their significant elevation above the photosphere. The ribbons in the Ca II 8542 {\AA} line wing images are located deeper in the lower solar atmosphere with a lower contrast. A comparison of the height deposition of electron beam energy and the intensity contribution function shows that the H$\beta$ line wing intensities can be an useful tracer of flare energy deposition in the lower solar atmosphere

Characterising the large-scale structure and plasma properties of the inner corona is crucial to understand the source and subsequent expansion of the solar wind and related space weather effects. Here we apply a new coronal rotational tomography method, along with a method to narrow streamers and refine the density estimate, to COR2A/STEREO observations from a period near solar minimum and maximum, gaining density maps for heights between 4 and 8\Rs. The coronal structure is highly radial at these heights, and the streamers are very narrow, in some regions only a few degrees in width. The mean densities of streamers is almost identical between solar minimum and maximum. However, streamers at solar maximum contain around 50\%\ more total mass due to their larger area. By assuming a constant mass flux, and constraints on proton flux measured by Parker Solar Probe (PSP), we estimate an outflow speed within solar minimum streamers of 50-120\kms\ at 4\Rs, increasing to 90-250\kms\ at 8\Rs. Accelerations of around 6\mss\ are found for streamers at a height of 4\Rs, decreasing with height. The solar maximum slow wind shows a higher acceleration to extended distances compared to solar minimum. To satisfy the solar wind speeds measured by PSP, there must be a mean residual acceleration of around 1-2\mss\ between 8 and 40\Rs. Several aspects of this study strongly suggest that the coronal streamer belt density is highly variable on small scales, and that the tomography can only reveal a local spatial and temporal average.

Flare loops form an integral part of eruptive events, being detected in the range of temperatures from X-rays down to cool chromospheric-like plasmas. While the hot loops are routinely observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA), cool loops seen off-limb are rare. In this paper we employ unique observations of the SOL2017-09-10T16:06 X8.2-class flare which produced an extended arcade of loops. The Swedish 1-m Solar Telescope (SST) made a series of spectral images of the cool off-limb loops in the Ca II 8542 \r{A} and the hydrogen H$\beta$ lines. Our focus is on the loop apices. Non-LTE spectral inversion is achieved through the construction of extended grids of models covering a realistic range of plasma parameters. The Multilevel Accelerated Lambda Iterations (MALI) code solves the non-LTE radiative-transfer problem in a 1D externally-illuminated slab, approximating the studied loop segment. Inversion of the Ca II 8542 \r{A} and H$\beta$ lines yields two similar solutions, both indicating high electron densities around $2 \times 10^{12}$ cm$^{-3}$ and relatively large microturbulence around 25 kms$^{-1}$. These are in reasonable agreement with other independent studies of the same or similar events. In particular, the high electron densities in the range $10^{12} - 10^{13}$ cm$^{-3}$ are consistent with those derived from the SDO's Helioseismic and Magnetic Imager white-light observations. The presence of such high densities in solar eruptive flares supports the loop interpretation of the optical continuum emission of stars which manifest superflares.

The increasing size of solar datasets demands highly efficient and robust analysis methods. This paper presents an approach that can increase the computational efficiency of differential emission measure (DEM) inversions by an order of magnitude or higher, with the efficiency factor increasing with the size of the input dataset. The method, named the Gridded Solar Iterative Temperature Emission Solver (Grid-SITES) is based on grouping pixels according to the similarity of their intensities in multiple channels, and solving for one DEM per group. This is shown to be a valid approach, given a sufficiently high number of grid bins for each channel. The increase in uncertainty arising from the quantisation of the input data is small compared to the general measurement and calibration uncertainties. In this paper, we use the Solar Iterative Temperature Emission Solver (SITES) as the core method for the DEM inversion, although Grid-SITES provides a general framework which may be used with any DEM inversion method, or indeed any large multi-dimensional data inversion problem. The method is particularly efficient for processing larger images, offering a factor of 30 increase in speed for a 10 megapixel image. For a time series of observations, the gridded results can be passed sequentially to each new image, with new populated bins added as required. This process leads to increasing efficiency with each new image, with potential for a $\approx$100 increase in efficiency dependent on the size of the images.

This is the second of a series of three papers that present a methodology with the aim of creating a set of maps of the coronal density over a period of many years. This paper describes a method for reconstructing the coronal electron density based on spherical harmonics. By assuming a radial structure to the corona at the height of interest, line-of-sight integrations can be made individually on each harmonic basis prior to determining coefficients, i.e. the computationally-expensive integrations are calculated only once during initialization. This approach reduces the problem to finding the set of coefficients which best match the observed brightness using a regularized least-squares approach, and is very efficient. The method is demonstrated on synthetic data created from both a simple and an intricate coronal density model. The quality of reconstruction is found to be reasonable in the presence of noise and large gaps in the data. The method is applied to both LASCO C2 and STEREO COR2 coronagraph observations from 2009/03/20, and the results from both spacecraft compared. Future work will apply the method to large datasets.

Extreme UltraViolet (EUV) images of the optically-thin solar corona in multiple spectral channels give information on the emission as a function of temperature through differential emission measure (DEM) inversions. The aim of this paper is to describe, test, and apply a new DEM method named the Solar Iterative Temperature Emission Solver (SITES). The method creates an initial DEM estimate through a direct redistribution of observed intensities across temperatures according to the temperature response function of the measurement, and iteratively improves on this estimate through calculation of intensity residuals. It is simple in concept and implementation, is non-subjective in the sense that no prior constraints are placed on the solutions other than positivity and smoothness, and can process a thousand DEMs per second on a standard desktop computer. The resulting DEMs replicate model DEMs well in tests on Atmospheric Imaging Assembly (AIA) synthetic data. The same tests show that SITES performs less well on very narrow DEM peaks, and should not be used for temperature diagnostics below ~0.5MK in the case of AIA observations. The SITES accuracy of inversion compares well with two other established methods. A simple yet powerful new method to visualise DEM maps is introduced, based on a fractional emission measure (FEM). Applied to a set of AIA full-disk images, the SITES method and FEM visualisation show very effectively the dominance of certain temperature regimes in different large-scale coronal structures. The method can easily be adapted for any multi-channel observations of optically-thin plasma and, given its simplicity and efficiency, will facilitate the processing of large existing and future datasets.

Here we report on the unique observation of flaring coronal loops at the solar limb using high resolution imaging spectropolarimetry from the Swedish 1-meter Solar Telescope. The vantage position, orientation and nature of the chromospheric material that filled the flare loops allowed us to determine their magnetic field with unprecedented accuracy using the weak-field approximation method. Our analysis reveals coronal magnetic field strengths as high as 350 Gauss at heights up to 25 Mm above the solar limb. These measurements are substantially higher than a number of previous estimates and may have considerable implications for our current understanding of the extended solar atmosphere.

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

A new, automated method of detecting coronal mass ejections (CMEs) in three dimensions for the LASCO C2 and STEREO COR2 coronagraphs is presented. By triangulating isolated CME signal from the three coronagraphs over a sliding window of five hours, the most likely region through which CMEs pass at 5 solar radii is identified. The centre and size of the region gives the most likely direction of propagation and approximate angular extent. The Automated CME Triangulation (ACT) method is tested extensively using a series of synthetic CME images created using a wireframe flux rope density model, and on a sample of real coronagraph data; including halo CMEs. The accuracy of the angular difference between the detection and true input of the synthetic CMEs is 7.14 degrees, and remains acceptable for a broad range of CME positions relative to the observer, the relative separation of the three observers and even through the loss of one coronagraph. For real data, the method gives results that compare well with the distribution of low coronal sources and results from another instrument and technique made further from the Sun. The true three dimension (3D)-corrected kinematics and mass/density are discussed. The results of the new method will be incorporated into the CORIMP database in the near future, enabling improved space weather diagnostics and forecasting.

Tomography of the solar corona can provide cruicial constraints for models of the low corona, unique information on changes in coronal structure and rotation rates, and a valuable boundary condition for models of the heliospheric solar wind. This is the first of a series of three papers which aim to create a set of maps of the coronal density over an extended period (1996-present). The papers will describe the data processing and calibration (this paper), the tomography method (\paperii) and resulting atlas of coronal electron density at a height of 5\Rs\ between years 1996-2014 (\paperiii). This first paper presents a detailed description of data processing and calibration for the Large-Angle and Spectrometric Coronagraph (LASCO) C2 instrument onboard the Solar and Heliospheric Observatory (SOHO) and the COR2 instruments of the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) package aboard the Solar Terrestial Relations Observatory (STEREO) A \& B spacecraft. The methodology includes noise suppression, background subtraction, separation of large dynamic events, conversion of total brightness to K-coronal brightness and simple functions for crosscalibration between C2/LASCO and COR2/SECCHI. Comparison of the brightness of stars between LASCO C2 total and polarized brightness (\pB) observations provide in-flight calibration factors for the \pB\ observations, resulting in considerable improved agreement between C2 and COR2 A, and elimination of curious artifacts in the C2 \pB\ images. The crosscalibration between LASCO C2 and the STEREO coronagraphs allows, for the first time, the potential use of multi-spacecraft coronagraph data for tomography and for CME analysis.

The solar atmosphere is a hot (about 1MK), magnetised plasma of great interest to physicists. There have been many previous studies of the temperature of the Sun's atmosphere (Plowman2012, Wit2012, Hannah2012, Aschwanden2013, etc.). Almost all of these studies use the SolarSoft software package written in the commercial Interactive Data Language (IDL), which has been the standard language for solar physics. The SunPy project aims to provide an open-source library for solar physics. This work presents (to the authors' knowledge) the first study of its type to use SunPy rather than SolarSoft. This work uses SunPy to process multi-wavelength solar observations made by the Atmospheric Imaging Assembly (AIA) instrument aboard the Solar Dynamics Observatory (SDO) and produce temperature maps of the Sun's atmosphere. The method uses SunPy's utilities for querying databases of solar events, downloading solar image data, storing and processing images as spatially aware Map objects, and tracking solar features as the Sun rotates. An essential consideration in developing this software is computational efficiency due to the large amount of data collected by AIA/SDO, and in anticipating new solar missions which will result in even larger sets of data. An overview of the method and implementation is given, along with tests involving synthetic data and examples of results using real data for various regions in the Sun's atmosphere.

The initiation phase of CMEs is a very important aspect of solar physics, as these phenomena ultimately drive space weather in the heliosphere. This phase is known to occur between the photosphere and low corona, where many models introduce an instability and/or magnetic reconnection that triggers a CME, often with associated flaring activity. To this end, it is important to obtain a variety of observations of the low corona in order to build as clear a picture as possible of the dynamics that occur therein. Here, we combine the EUV imagery of the SWAP instrument on board PROBA2 with the white-light imagery of the ground-based Mk4 coronameter at MLSO in order to bridge the observational gap that exists between the disk imagery of AIA on board SDO and the coronal imagery of LASCO on board SOHO. Methods of multiscale image analysis were applied to the observations to better reveal the coronal signal while suppressing noise and other features. This allowed an investigation into the initiation phase of a CME that was driven by a rising flux rope structure from a "two-stage" flaring active region underlying an extended helmet streamer. It was found that the initial outward motion of the erupting loop system in the EUV observations coincided with the first X-ray flare peak, and led to a plasma pile-up of the white-light CME core material. The characterized CME core then underwent a strong jerk in its motion, as the early acceleration increased abruptly, simultaneous with the second X-ray flare peak. The overall system expanded into the helmet streamer to become the larger CME structure observed in the LASCO coronagraph images, which later became concave-outward in shape. Theoretical models for the event are discussed in light of these unique observations, and it is concluded that the formation of either a kink-unstable or torus-unstable flux rope may be the likeliest scenario.

Extreme UltraViolet images of the corona contain information over a large range of spatial scales, and different structures such as active regions, quiet Sun and filament channels contain information at very different brightness regimes. Processing of these images is important to reveal information, often hidden within the data, without introducing artifacts or bias. It is also important that any process be computationally efficient, particularly given the fine spatial and temporal resolution of Atmospheric Imaging Assembly on the Solar Dynamics Observatory (AIA/SDO) , and consideration of future higher-resolution observations. A very efficient process is described here which is based on localized normalizing of the data at many different spatial scales. The method reveals information at the finest scales, whilst maintaining enough of the larger-scale information to provide context. It also intrinsically flattens noisy regions and can reveal structure in off-limb regions out to the edge of the field of view. The method is also successfully applied to a white light coronagraph observation.

The study of solar eruptive events and associated phenomena is of great importance in the context of solar and heliophysics. Coronal mass ejections (CMEs) and coronal waves are energetic manifestations of the restructuring of the solar magnetic field and mass motion of the plasma. Characterising this motion is vital for deriving the dynamics of these events and thus understanding the physics driving their initiation and propagation. The development and use of appropriate methods for measuring event kinematics is therefore imperative. Traditional approaches to the study of CME and coronal wave kinematics do not return wholly accurate nor robust estimates of the true event kinematics and associated uncertainties. We highlight the drawbacks of these approaches, and demonstrate improved methods for accurate and reliable determination of the kinematics. The Savitzky-Golay filter is demonstrated as a more appropriate fitting technique for CME and coronal wave studies, and a residual resampling bootstrap technique is demonstrated as a statistically rigorous method for the determination of kinematic error estimates and goodness-of-fit tests. It is shown that the scatter on distance-time measurements of small sample size can significantly limit the ability to derive accurate and reliable kinematics. This may be overcome by (i) increasing measurement precision and sampling cadence, and (ii) applying robust methods for deriving the kinematics and reliably determining their associated uncertainties. If a priori knowledge exists and a pre-determined model form for the kinematics is available (or indeed any justified fitting-form to be tested against the data), then its precision can be examined using a bootstrapping technique to determine the confidence interval associated with the model/fitting parameters.

Studying CMEs in coronagraph data can be challenging due to their diffuse structure and transient nature, and user-specific biases may be introduced through visual inspection of the images. The large amount of data available from the SOHO, STEREO, and future coronagraph missions, also makes manual cataloguing of CMEs tedious, and so a robust method of detection and analysis is required. This has led to the development of automated CME detection and cata- loguing packages such as CACTus, SEEDS and ARTEMIS. Here we present the development of a new CORIMP (coronal image processing) CME detection and tracking technique that overcomes many of the drawbacks of current catalogues. It works by first employing the dynamic CME separation technique outlined in a companion paper, and then characterising CME structure via a multiscale edge-detection algorithm. The detections are chained through time to determine the CME kinematics and morphological changes as it propagates across the plane-of-sky. The effectiveness of the method is demonstrated by its application to a selection of SOHO/LASCO and STEREO/SECCHI images, as well as to synthetic coronagraph images created from a model corona with a variety of CMEs. The algorithms described in this article are being applied to the whole LASCO and SECCHI datasets, and a catalogue of results will soon be available to the public.

During 2011/09/24, as observed by the Atmospheric Imaging Assembly (AIA) instrument of the Solar Dynamic Observatory (SDO) and ground-based \Ha\ telescopes, a prominence and associated cavity appeared above the southwest limb. On 2011/09/25 8:00UT material flows upwards from the prominence core along a narrow loop-like structure, accompanied by a rise ($\geq$50,000km) of the prominence core and the loop. As the loop fades by 10:00, small blobs and streaks of varying brightness rotate around the top part of the prominence and cavity, mimicking a cyclone. The most intense and coherent rotation lasts for over three hours, with emission in both hot ($\sim$1MK) and cold (hydrogen and helium) lines. We suggest that the cyclonic appearance and overall evolution of the structure can be interpreted in terms of the expansion of helical structures into the cavity, and the movement of plasma along helical structures which appears as a rotation when viewed along the helix axis. The coordinated movement of material between prominence and cavity suggest that they are structurally linked. Complexity is great due to the combined effect of these actions and the line-of-sight integration through the structure which contains tangled fields.

This white paper is a call for a concerted effort to support total solar eclipse observations over the next decade, in particular for the 21 August 2017 eclipse which will traverse the US continent. With the recent advances in image processing techniques and detector technology, the time is ripe to capitalize on the unique diagnostic tools available in the visible and near infrared wavelength range to explore the physics of the corona. The advantage of coronal emission lines in this wavelength range, over their extreme ultraviolet counterparts, is (1) the significant radiative component in their excitation process (in addition to the collisional excitation), which allows for observations out to a few solar radii, (2) the higher spectral selectivity available for imaging, giving well-defined temperature responses for each bandpass (one line as opposed to many), and (3) the capability of polarization measurements in a number of spectral lines. Consequently, the evolution of the thermodynamic and magnetic properties of the coronal plasma can be explored starting from the solar surface out to a few solar radii, namely the most important region of the corona where the expansion of the solar magnetic field and the acceleration of the solar wind occur. Since the planning of eclipse observations will not be possible without the invaluable NASA-published total solar eclipse bulletins by Espenak and Andersen, a call is also made to ensure continued support for these efforts.

Radio remote sensing of the heliosphere using spacecraft radio signals has been used to study the near-sun plasma in and out of the ecliptic, close to the sun, and on spatial and temporal scales not accessible with other techniques. Studies of space-time variations in the inner solar wind are particularly timely because of the desire to understand and predict space weather, which can disturb satellites and systems at 1AU and affect human space exploration. Here we demonstrate proof-of-concept of a new radio science application for spacecraft radio science links. The differing transfer functions of plasma irregularities to spacecraft radio up- and downlinks can be exploited to localize plasma scattering along the line of sight. We demonstrate the utility of this idea using Cassini radio data taken in 2001-2002. Under favorable circumstances we demonstrate how this technique, unlike other remote sensing methods, can determine center-of-scattering position to within a few thousandths of an AU and thickness of scattering region to less than about 0.02 AU. This method, applied to large data sets and used in conjunction with other solar remote sensing data such as white light data, has space weather application in studies of inhomogeneity and nonstationarity in the near-sun solar wind.

Sunspots are locations on the Sun where unique atmospheric conditions prevail. In particular, the very low temperatures found above sunspots allow the emission of H_2 lines. In this study we are interested in the radiation emitted by sunspots in the O VI lines at 1031.96 A and 1037.60 A. We use SOHO/SUMER observations of a sunspot performed in March 1999 and investigate the interaction between the O VI lines and a H_2 line at 1031.87 A found in the Werner band. The unique features of sunspots atmospheres may very well have important implications regarding the illumination of coronal O+5 ions in the low corona, affecting our interpretation of Doppler dimming diagnostics.

Despite the wealth of solar data currently available, the explicit connection between coronal streamers and features on the solar disk remains unresolved. An empirical three-dimensional model, which reproduces the evolution of the large-scale coronal structure starting from the solar surface, is presented. The model is based on the view that the source of the large-scale coronal structure, namely streamers, is a consequence of the evolution of twisted sheet-like structures originating from prominences (or, equivalently, filaments) at the base of the corona. The high-density sheets evolve and merge with height into a final radial configuration constrained by the oberved position of streamers stalks higher up in the corona. The observational constraints are provided by white light observations from the LASCO/C2 data during the declining phase of solar activity, spanning the end of Carrington Rotation (CR) 2005 and the start of CR 2006, i.e. July-August 2003, and the position of filaments from the corresponding H$\alpha$ synoptic maps of the Paris-Meudon Observatory. The 3D model thus derived yields a reasonable agreement with the observed large-scale coronal structure, in particular the shape of large helmet streamers. These results give confidence in the underlying assumption that large helmet streamers can be the result of the convergence of two or more sheet-like structures originating from a distribution of prominences on the solar disk. The model supports the view that streamers, during that time of the solar cycle, are often associated with multiple current sheets.

The very sharp decrease of density with heliocentric distance makes imaging of coronal density structures out to a few solar radii challenging. The radial gradient in brightness can be reduced using numerous image processing techniques, thus quantitative data are manipulated to provide qualitative images. Introduced in this study is a new normalizing radial graded filter (NRGF), a simple filter for removing the radial gradient to reveal coronal structure. Applied to polarized brightness observations of the corona, the NRGF produces images which are striking in their detail. Total brightness white light images include contributions from the F corona, stray light and other instrumental contributions which need to be removed as effectively as possible to properly reveal the electron corona structure. A new procedure for subtracting this background from LASCO C2 white light total brightness images is introduced. The background is created from the unpolarized component of total brightness images and is found to be remarkably time-invariant, remaining virtually unchanged over the solar cycle. By direct comparison with polarized brightness data, we show that the new background subtracting procedure is superior in depicting coronal structure accurately, particularly when used in conjunction with the NRGF. The effectiveness of the procedures is demonstrated on a series of LASCO C2 observations of a coronal mass ejection (CME).

Due to their high intensity of emission in the O VI 1031.9 and 1037.6A lines, even small sunspots on the solar disk can strongly influence the intensity of the radiative scattering component of O VI lines in the corona. Observations of O VI disk spectra show a 1032/1038 line intensity ratio of >2.6 in a sunspot compared to quiet disk values of 2. The enhancement of the 1032 line in comparison to the 1038 is likely due to interaction between molecular hydrogen emission from the sunspot and the chromospheric O5+. Modeling shows that a contribution from sunspots increases the coronal O VI 1032/1038 intensity ratio to values considerably higher than those achieved with a quiet disk or coronal hole spectrum. Therefore a re-examination of flow velocities derived from UVCS/SOHO streamer observations must be made. This modeling demonstrates that the inclusion of sunspots, when present, may lead to non-zero outflow velocities at lower heights in streamer cores in contrast to some existing model results.

The low ionization state in parts of a sunspot may play an important role in its evolution and dynamical state. The cool magnetic interior region of the sunspot develops a substantial neutral atomic and molecular hydrogen osmotic pressure which can drive a wind outward from the umbra. Ambipolar diffusion against the magnetically pinned ionized plasma component can also distort the umbral magnetic field into a filamentary penumbral structure. This may be important for explaining the development of the sunspot penumbra and the Evershed flow. This fractionation process may also be important for the ``First Ionization Potential'' (FIP) effect seen in the solar wind. In support of this mechanism we find evidence for such ionization fractionization in UV observations of molecular hydrogen in a sunspot umbra and penumbra.

The timescales to replenish dust from the cool, dense winds of Asymptotic Giant Branch stars are believed to be greater than the timescales for dust destruction. In high redshift galaxies, this problem is further compounded as the stars take longer than the age of the Universe to evolve into the dust production stages. To explain these discrepancies, dust formation in supernovae (SNe) is required to be an important process but until very recently dust in supernova remnants has only been detected in very small quantities. We present the first submillimeter observations of cold dust in Kepler's supernova remnant (SNR) using SCUBA. A two component dust temperature model is required to fit the Spectral Energy Distribution (SED) with $T_{warm} \sim 102$K and $T_{cold} \sim 17$K. The total mass of dust implied for Kepler is $\sim 1M_{\odot}$ - 1000 times greater than previous estimates. Thus SNe, or their progenitors may be important dust formation sites.

Large amounts of dust have recently been discovered in high-z galaxies and QSOs. The stellar winds produced by AGB stars are thought to be the main source of dust in galaxies, but they cannot produce that dust on a short enough timescale (<1 Gyr) to explain the results in the high-z galaxies. Supernova explosions of massive stars (TypeII) are also a potential source, with models predicting 0.2-4 Msun of dust from each. As massive stars evolve rapidly, on timescales of a few Myr, these supernovae could be responsible for the high-z dust. Observations of supernova remnants in the Milky Way have previously only revealed 10^-7 - 10^-3 Msun of dust each, which is insufficient to explain the high-z data. Here we report the detection of ~2-4 Msun of cold dust in Cas A. This implies that supernovae are at least as important as stellar winds in producing dust in our Galaxy and would have been the dominant source of dust at high-z.

We investigate the sources and amount of dust in early galaxies. We discuss dust nucleation in stellar atmospheres using published extended atmosphere models, stellar evolution tracks and nucleation conditions and conclude that the (TPAGB) phase of intermediate mass stars is likely to be the most promising site for dust formation in stellar winds. The implications of chemical evolution models for high redshift galaxies are investigated and we show there is no difficulty in producing dusty galaxies at redshifts above 5 if supernovae are a dominant source of interstellar dust. If dust does not condense efficiently in SNe then significant dust masses can only be generated at by galaxies with a high star formation efficiency. We find the visual optical depth for individual star forming clouds can reach values greater than 1 at very low metallicity (1/100 solar) provided that the mass-radius exponent of molecular clouds is less than two. Most of the radiation from star formation will emerge at IR wavelengths in the early universe provided that dust is present. The (patchy) visual optical depth through a typical early galaxy will however, remain less than 1 on average until a metallicity of 1/10 solar is reached.