A synoptic view of solar transient evolution in the inner heliosphere using the Heliospheric Imagers on STEREO

By exploiting data from the STEREO/heliospheric imagers (HI) we extend a well-established technique developed for coronal analysis by producing time-elongation plots that reveal the nature of solar transient activity over a far more extensive region of the heliosphere than previously possible from coronagraph images. Despite the simplicity of these plots, their power in demonstrating how the plethora of ascending coronal features observed near the Sun evolve as they move antisunward is obvious. The time-elongation profile of a transient tracked by HI can, moreover, be used to establish its angle out of the plane-of-the-sky; an illustration of such analysis reveals coronal mass ejection material that can be clearly observed propagating out to distances beyond 1AU. This work confirms the value of the time-elongation format in identifying/characterising transient activity in the inner heliosphere, whilst also validating the ability of HI to continuously monitor solar ejecta out to and beyond 1A

Joint-bounded crescentic scars formed by subglacial clast-bed contact forces: implications for bedrock failure beneath glaciers

Glaciers and ice sheets are important agents of bedrock erosion, yet the precise processes of bedrock failure beneath glacier ice are incompletely known. Subglacially formed erosional crescentic markings (crescentic gouges, lunate fractures) on bedrock surfaces occur locally in glaciated areas and comprise a conchoidal fracture dipping down-ice and a steep fracture that faces up-ice. Here we report morphologically distinct crescentic scars that are closely associated with preexisting joints, termed here joint-bounded crescentic scars. These hitherto unreported features are ca. 50–200 mm deep and involve considerably more rock removal than previously described crescentic markings. The joint-bounded crescentic scars were found on abraded rhyolite surfaces recently exposed (< 20 years) beneath a retreating glacier in Iceland, as well as on glacially sculpted Precambrian gneisses in NW Scotland and various Precambrian rocks in Ontario, glaciated during the Late Pleistocene. We suggest a common formation mechanism for these contemporary and relict features, whereby a boulder embedded in basal ice produces a continuously migrating clast-bed contact force as it is dragged over the hard (bedrock) bed. As the ice-embedded boulder approaches a preexisting joint in the bedrock, stress concentrations build up in the bed that exceed the intact rock strength, resulting in conchoidal fracturing and detachment of a crescentic wedge-shaped rock fragment. Subsequent removal of the rock fragment probably involves further fracturing or crushing (comminution) under high contact forces. Formation of joint-bounded crescentic scars is favoured by large boulders at the base of the ice, high basal melting rates, and the presence of preexisting subvertical joints in the bedrock bed. We infer that the relative scarcity of crescentic markings in general on deglaciated surfaces shows that fracturing of intact bedrock below ice is difficult, but that preexisting weaknesses such as joints greatly facilitate rock failure. This implies that models of glacial erosion need to take fracture patterns of bedrock into account

Clustering of i-dropout galaxies at z=6 in GOODS and the UDF

We measured the angular clustering at z~6 from a large sample of i-dropout galaxies (293 with z<27.5 from GOODS and 95 with z<29.0 from the UDF). Our largest and most complete subsample (having L>0.5L*) shows the presence of clustering at 94% significance. For this sample we derive a (co-moving) correlation length of r_0=4.5^{+2.1}_{-3.2} h_{72}^{-1} Mpc and bias b=4.1^{+1.5}_{-2.6}, using an accurate model for the redshift distribution. No clustering could be detected in the much deeper but significantly smaller UDF, yielding b<4.4 (1 sigma). We compare our findings to Lyman break galaxies at z=3-5 at a fixed luminosity. Our best estimate of the bias parameter implies that i-dropouts are hosted by dark matter halos having masses of ~10^11 M_sun, similar to that of V-dropouts at z~5. We evaluate a recent claim that at z>5 star formation might have occurred more efficiently compared to that at z=3-4. This may provide an explanation for the very mild evolution observed in the UV luminosity density between z=6 and z=3. Although our results are consistent with such a scenario, the errors are too large to find conclusive evidence for this.Comment: minor changes to match published versio

Interpreting high [O III]/H β ratios with maturing starbursts

Star-forming galaxies at high redshift show ubiquitously high-ionization parameters, as measured by the ratio of optical emission lines. We demonstrate that local (z < 0.2) sources selected as Lyman break analogues also manifest high line ratios with a typical [O III]/Hβ=3.36+0.14−0.04 – comparable to all but the highest ratios seen in star-forming galaxies at z ∼ 2–4. We argue that the stellar population synthesis code BPASS can explain the high-ionization parameters required through the ageing of rapidly formed star populations, without invoking any AGN contribution. Binary stellar evolution pathways prolong the age interval over which a starburst is likely to show elevated line ratios, relative to those predicted by single stellar evolution codes. As a result, model galaxies at near-solar metallicities and with ages of up to ∼100 Myr after a starburst typically have a line ratio [O III]/Hβ ∼ 3, consistent with those seen in Lyman break galaxies and local sources with similar star formation densities. This emphasises the importance of including binary evolution pathways when simulating the nebular line emission of young or bursty stellar populations

Optical to mid-IR observations of Lyman-alpha galaxies at redshift 5 in the HUDF: a young and low mass population

High redshift galaxies selected on the basis of their strong Lyman-alpha emission tend to be young ages and small physical sizes. We show this by analyzing the spectral energy distribution (SED) of 9 Lyman-alpha emitting (LAE) galaxies at 4.0 < z < 5.7 n the Hubble Ultra Deep Field (HUDF). Rest-frame UV to optical 700A < lambda < 7500A luminosities, or upper limits, are used to constrain old stellar populations. We derive best fit, as well as maximally massive and maximally old, properties of all 9 objects. We show that these faint and distant objects are all very young, being most likely only a few millions years old, and not massive, the mass in stars being ~10^6-10^8 M_sun. Deep Spitzer Infrared Array Camera (IRAC) observations of these objects, even in cases where objects were not detected, were crucial in constraining the masses of these objects. The space density of these objects, ~1.25x10^-4 Mpc^-3 is comparable to previously reported space density of LAEs at moderate to high redshifts. These Lyman-alpha galaxies show modest star formation rates of ~8 M_sun yr^-1, which is nevertheless strong enough to have allowed these galaxies to assemble their stellar mass in less than a few x10^6 years. These sources appear to have small physical sizes, usually smaller than 1 Kpc, and are also rather concentrated. They are likely to be some of the least massive and youngest high redshift galaxies observed to date.Comment: 27 pages. Accepted. To appear in The Astrophysical Journa

Speeds and arrival times of solar transients approximated by self-similar expanding circular fronts

The NASA STEREO mission opened up the possibility to forecast the arrival times, speeds and directions of solar transients from outside the Sun-Earth line. In particular, we are interested in predicting potentially geo-effective Interplanetary Coronal Mass Ejections (ICMEs) from observations of density structures at large observation angles from the Sun (with the STEREO Heliospheric Imager instrument). We contribute to this endeavor by deriving analytical formulas concerning a geometric correction for the ICME speed and arrival time for the technique introduced by Davies et al. (2012, ApJ, in press) called Self-Similar Expansion Fitting (SSEF). This model assumes that a circle propagates outward, along a plane specified by a position angle (e.g. the ecliptic), with constant angular half width (lambda). This is an extension to earlier, more simple models: Fixed-Phi-Fitting (lambda = 0 degree) and Harmonic Mean Fitting (lambda = 90 degree). This approach has the advantage that it is possible to assess clearly, in contrast to previous models, if a particular location in the heliosphere, such as a planet or spacecraft, might be expected to be hit by the ICME front. Our correction formulas are especially significant for glancing hits, where small differences in the direction greatly influence the expected speeds (up to 100-200 km/s) and arrival times (up to two days later than the apex). For very wide ICMEs (2 lambda > 120 degree), the geometric correction becomes very similar to the one derived by M\"ostl et al. (2011, ApJ, 741, id. 34) for the Harmonic Mean model. These analytic expressions can also be used for empirical or analytical models to predict the 1 AU arrival time of an ICME by correcting for effects of hits by the flank rather than the apex, if the width and direction of the ICME in a plane are known and a circular geometry of the ICME front is assumed.Comment: 15 pages, 5 figures, accepted for publication in "Solar Physics

Debris cover and surface melt at a temperate maritime alpine glacier: Franz Josef Glacier, New Zealand

Melt rates on glaciers are strongly influenced by the presence of supraglacial debris, which can either enhance or reduce ablation relative to bare ice. Most recently, Franz Josef Glacier has entered into a phase of strong retreat and downwasting, with the increasing emergence of debris on the surface in the ablation zone. Previously at Franz Josef Glacier, melt has only been measured on bare ice. During February 2012, a network of 11 ablation stakes was drilled into locations of varying supraglacial debris thickness on the lower glacier. Mean ablation rates over 9 days varied over the range 1.2–10.1 cm d−1, and were closely related to debris thickness. Concomitant observations of air temperature allowed the application of a degree-day approach to the calculation of melt rates, with air temperature providing a strong indicator of melt. Degree-day factors (d f) varied over the range 1.1–8.1 mm d−1 °C−1 (mean of 4.4 mm d−1 °C−1), comparable with rates reported in other studies. Mapping of the current debris cover revealed 0.7 km2 of the 4.9 km2 ablation zone surface was debris-covered, with thicknesses ranging 1–50 cm. Based on measured debris thicknesses and d f, ablation on debris-covered areas of the glacier is reduced by a total of 41% which equates to a 6% reduction in melt overall across the entire ablation zone. This study highlights the usefulness of a short-term survey to gather representative ablation data, consistent with numerous overseas ablation studies on debris-covered glaciers

Accuracy and Limitations of Fitting and Stereoscopic Methods to Determine the Direction of Coronal Mass Ejections from Heliospheric Imagers Observations

Using data from the Heliospheric Imagers (HIs) onboard STEREO, it is possible to derive the direction of propagation of coronal mass ejections (CMEs) in addition to their speed with a variety of methods. For CMEs observed by both STEREO spacecraft, it is possible to derive their direction using simultaneous observations from the twin spacecraft and also, using observations from only one spacecraft with fitting methods. This makes it possible to test and compare different analyses techniques. In this article, we propose a new fitting method based on observations from one spacecraft, which we compare to the commonly used fitting method of Sheeley et al. (1999). We also compare the results from these two fitting methods with those from two stereoscopic methods, focusing on 12 CMEs observed simultaneously by the two STEREO spacecraft in 2008 and 2009. We find evidence that the fitting method of Sheeley et al. (1999) can result in significant errors in the determination of the CME direction when the CME propagates outside of 60deg \pm 20 deg from the Sun-spacecraft line. We expect our new fitting method to be better adapted to the analysis of halo or limb CMEs with respect to the observing spacecraft. We also find some evidence that direct triangulation in the HI fields-of-view should only be applied to CMEs propagating approximatively towards Earth (\pm 20deg from the Sun-Earth line). Last, we address one of the possible sources of errors of fitting methods: the assumption of radial propagation. Using stereoscopic methods, we find that at least seven of the 12 studied CMEs had an heliospheric deflection of less than 20deg as they propagated in the HI fields-of-view, which, we believe, validates this approximation.Comment: 17 pages, 6 figures, 2 tables, accepted to Solar Physic

Effect of Solar Wind Drag on the Determination of the Properties of Coronal Mass Ejections from Heliospheric Images

The Fixed-\Phi (F\Phi) and Harmonic Mean (HM) fitting methods are two methods to determine the average direction and velocity of coronal mass ejections (CMEs) from time-elongation tracks produced by Heliospheric Imagers (HIs), such as the HIs onboard the STEREO spacecraft. Both methods assume a constant velocity in their descriptions of the time-elongation profiles of CMEs, which are used to fit the observed time-elongation data. Here, we analyze the effect of aerodynamic drag on CMEs propagating through interplanetary space, and how this drag affects the result of the F\Phi and HM fitting methods. A simple drag model is used to analytically construct time-elongation profiles which are then fitted with the two methods. It is found that higher angles and velocities give rise to greater error in both methods, reaching errors in the direction of propagation of up to 15 deg and 30 deg for the F\Phi and HM fitting methods, respectively. This is due to the physical accelerations of the CMEs being interpreted as geometrical accelerations by the fitting methods. Because of the geometrical definition of the HM fitting method, it is affected by the acceleration more greatly than the F\Phi fitting method. Overall, we find that both techniques overestimate the initial (and final) velocity and direction for fast CMEs propagating beyond 90 deg from the Sun-spacecraft line, meaning that arrival times at 1 AU would be predicted early (by up to 12 hours). We also find that the direction and arrival time of a wide and decelerating CME can be better reproduced by the F\Phi due to the cancellation of two errors: neglecting the CME width and neglecting the CME deceleration. Overall, the inaccuracies of the two fitting methods are expected to play an important role in the prediction of CME hit and arrival times as we head towards solar maximum and the STEREO spacecraft further move behind the Sun.Comment: Solar Physics, Online First, 17 page

Collisional depolarization of NO(A) by He and Ar studied by quantum beat spectroscopy

Zeeman and hyperfine quantum beat spectroscopies have been used to measure the total elastic plus inelastic angular momentum depolarization rate constants at 300 K for NO (A 2 σ+) in the presence of He and Ar. In the case of Zeeman quantum beats it is shown how the applied magnetic field can be used to allow measurement of depolarization rates for both angular momentum orientation and alignment. For the systems studied here, collisional loss of alignment is more efficient than loss of orientation. In the case of NO (A) with He, and to a lesser extent NO (A) with Ar, collisional depolarization is found to be a relatively minor process compared to rotational energy transfer, reflecting the very weak long-range forces in these systems. Detailed comparisons are made with quantum mechanical and quasiclassical trajectory calculations performed on recently developed potential energy surfaces. For both systems, the agreement between the calculated depolarization cross sections and the present measurements is found to be very good, suggesting that it is reasonable to consider the NO (A) bond as frozen during these angular momentum transferring collisions. A combination of kinematic effects and differences in the potential energy surfaces are shown to be responsible for the differences observed in depolarization cross section with He and Ar as a collider. © 2009 American Institute of Physics