Homologous Flares and Magnetic Field Topology in Active Region NOAA 10501 on 20 November 2003

We present and interpret observations of two morphologically homologous flares that occurred in active region (AR) NOAA 10501 on 20 November 2003. Both flares displayed four homologous Hα ribbons and were both accompanied by coronal mass ejections (CMEs). The central flare ribbons were located at the site of an emerging bipole in the centre of the active region. The negative polarity of this bipole fragmented in two main pieces, one rotating around the positive polarity by ≈ 110° within 32 hours. We model the coronal magnetic field and compute its topology, using as boundary condition the magnetogram closest in time to each flare. In particular, we calculate the location of quasi-separatrix layers (QSLs) in order to understand the connectivity between the flare ribbons. Though several polarities were present in AR 10501, the global magnetic field topology corresponds to a quadrupolar magnetic field distribution without magnetic null points. For both flares, the photospheric traces of QSLs are similar and match well the locations of the four Hα ribbons. This globally unchanged topology and the continuous shearing by the rotating bipole are two key factors responsible for the flare homology. However, our analyses also indicate that different magnetic connectivity domains of the quadrupolar configuration become unstable during each flare, so that magnetic reconnection proceeds differently in both events.     

Erratum to: Width of Sunspot Generating Zone and Reconstruction of Butterfly Diagram

Based on the extended Greenwich – NOAA/USAF catalogue of sunspot groups, it is demonstrated that the parameters describing the latitudinal width of the sunspot generating zone (SGZ) are closely related to the current level of solar activity, and the growth of the activity leads to the expansion of the SGZ. The ratio of the sunspot number to the width of the SGZ shows saturation at a certain level of the sunspot number, and above this level the increase of the activity takes place mostly due to the expansion of the SGZ. It is shown that the mean latitudes of sunspots can be reconstructed from the amplitudes of solar activity. Using the obtained relations and the group sunspot numbers by Hoyt and Schatten (Solar Phys. 179, 189, 1998), the latitude distribution of sunspot groups (“the Maunder butterfly diagram”) for the eighteenth and the first half of the nineteenth centuries is reconstructed and compared with historical sunspot observations.     

Constraints on Mercury’s Na exosphere: Combined MESSENGER and ground-based data

We have used observations of sodium emission obtained with the McMath-Pierce solar telescope and MESSENGER’s Mercury Atmospheric and Surface Composition Spectrometer (MASCS) to constrain models of Mercury’s sodium exosphere. The distribution of sodium in Mercury’s exosphere during the period January 12-15, 2008, was mapped using the McMath-Pierce solar telescope with the 5″ × 5″ image slicer to observe the D-line emission. On January 14, 2008, the Ultraviolet and Visible Spectrometer (UVVS) channel on MASCS sampled the sodium in Mercury’s anti-sunward tail region. We find that the bound exosphere has an equivalent temperature of 900-1200 K, and that this temperature can be achieved if the sodium is ejected either by photon-stimulated desorption (PSD) with a 1200 K Maxwellian velocity distribution, or by thermal accommodation of a hotter source. We were not able to discriminate between the two assumed velocity distributions of the ejected particles for the PSD, but the velocity distributions require different values of the thermal accommodation coefficient and result in different upper limits on impact vaporization. We were able to place a strong constraint on the impact vaporization rate that results in the release of neutral Na atoms with an upper limit of 2.1 × 10⁶ cm⁻² s⁻¹. The variability of the week-long ground-based observations can be explained by variations in the sources, including both PSD and ion-enhanced PSD, as well as possible temporal enhancements in meteoroid vaporization. Knowledge of both dayside and anti-sunward tail morphologies and radiances are necessary to correctly deduce the exospheric source rates, processes, velocity distribution, and surface interaction.     

Production and detection of carbon dioxide on Iapetus

Cassini VIMS detected carbon dioxide on the surface of Iapetus during its insertion orbit. We evaluated the CO₂ distribution on Iapetus and determined that it is concentrated almost exclusively on Iapetus’ dark material. VIMS spectra show a 4.27-μm feature with an absorption depth of 24%, which, if it were in the form of free ice, requires a layer 31 nm thick. Extrapolating for all dark material on Iapetus, the total observable CO₂ would be 2.3 × 10⁸ kg.Previous studies note that free CO₂ is unstable at 10 AU over geologic timescales. Carbon dioxide could, however, be stable if trapped or complexed, such as in inclusions or clathrates. While complexed CO₂ has a lower thermal volatility, loss due to photodissociation by UV radiation and gravitational escape would occur at a rate of 2.6 × 10⁷ kg year⁻¹. Thus, Iapetus’ entire inventory of surface CO₂ could be lost within a few decades.The high loss/destruction rate of CO₂ requires an active source. We conducted experiments that generated CO₂ by UV radiation of simulated icy regolith under Iapetus-like conditions. The simulated regolith was created by flash-freezing degassed water, crushing it into sub-millimeter sized particles, and then mixing it with isotopically labeled amorphous carbon (¹³C) dust. These samples were placed in a vacuum chamber and cooled to temperatures between 50 K and 160 K. The samples were irradiated with UV light, and the products were measured using a mass spectrometer, from which we measured ¹³CO₂ production at a rate of 2.0 × 10¹² mol s⁻¹. Extrapolating to Iapetus and adjusting for the solar UV intensity and Iapetus’ surface area, we calculated that CO₂ production for the entire surface would be 1.1 × 10⁷ kg year⁻¹, which is only a factor of two less than the loss rate. As such, UV photochemical generation of CO₂ is a plausible source of the detected CO₂.     

Modeling the morphological diversity of impact craters on icy satellites

Impact craters on icy satellites display a wide range of morphologies, some of which have no counterpart on rocky bodies. Numerical simulation studies have struggled to reproduce the diversity of features, such as central pits and transitions in crater depth with increasing diameter, observed on the icy Galilean satellites. The transitions in crater depth (at diameters of about 26 and 150 km on Ganymede and Callisto) have been interpreted as reflecting subsurface structure. Using the CTH shock physics code, we model the formation of craters with diameters between 400 m and about 200 km on Ganymede using different subsurface temperature profiles. Our calculations include recent improvements in the model equation of state for H₂O and quasi-static strength parameters for ice. We find that the shock-induced formation of dense high-pressure polymorphs (ices VI and VII) creates a gap in the crater excavation flow, which we call discontinuous excavation. For craters larger than about 20 km, discontinuous excavation concentrates a hot plug of material (>270 K and mostly on the melting curve) in the center of the crater floor. The size and occurrence of the hot plug are in good agreement with the observed characteristics of central pit craters, and we propose that a genetic link exists between them. We also derive depth versus diameter curves for different internal temperature profiles. In a 120 K isothermal crust, calculated craters larger than about 30 km diameter are deeper than observed and do not reproduce the transition at about 26 km diameter. Calculated crater depths are shallower and in good agreement with observations between about 30 and 150 km diameter using a warm thermal gradient representing a convective interior. Hence, the depth-to-diameter transition at about 26 km reflects thermal weakening of ice. Finally, simulation results generally support the hypothesis that the anomalous interior morphologies for craters larger than 100 km are related to the presence of a subsurface ocean.     

Mid-infrared lightcurve of Vesta

We present a thermal mid-infrared lightcurve of Asteroid 4 Vesta and use this to infer variations in thermophysical properties over the surface. Vesta was observed over three nights during the May 2007 opposition with the Infrared Telescope Facility on Mauna Kea. Mid-infrared observations are compared to a model based on the Standard Thermal Model which is draped over a Vesta shape model derived from Hubble Space Telescope observations.A visible lightcurve with similar aspect was used to estimate the albedo as Vesta rotates. Shape and albedo can explain some of the features observed in the mid-infrared lightcurve. However, variations in the thermophysical properties, such as the “beaming parameter,” over Vesta’s surface are required to completely explain the observations.In order to match the mid-infrared magnitudes observed of Vesta, a beaming parameter of ∼0.862 is required which is higher than other Main Belt Asteroids such as Ceres and Pallas (0.756), indicating a smoother and/or rockier surface on Vesta. Variations in the beaming parameter with longitude are invoked to reproduce the observed thermal variations. Surface materials with relatively high beaming values, indicating a smoother and/or rockier surface, in the eastern hemisphere of Vesta coincide with locations where impact excavations may have produced surfaces that are younger and brighter relative to the western hemisphere.     

High kinetic energy density jets in the Earth’s magnetosheath: A case study

We study several high kinetic energy density jets observed during a traversal of the dayside magnetosheath by the Cluster spacecraft on March 17, 2001, at various distances from the magnetopause, generally characterised by anomalously high values of the local magnetosonic Mach number. We concentrate on two jets observed just outside the magnetopause, the first almost parallel to the GSM x axis and the second directed northward-tailward along the nominal magnetopause surface. We present evidence that none of them can be ascribed to magnetic reconnection at the magnetopause and show that the magnetopause is severely deformed by the jets, so that its local normal forms an angle of 97° with the quiet time magnetopause normal. On these grounds, we suggest that the indentation of the magnetopause is caused by an anti-sunward jet ramming into the magnetopause slightly equatorward of the northern cusp and that the northward-tailward jet is the result of its reflection at the deformed magnetopause. Finally, we briefly discuss our results by comparing them with past studies of events which in some way recall the one analysed herein.     

Aerial electromagnetic sounding of the lithosphere of Venus

Electromagnetic (EM) investigation depths are larger on Venus than Earth due to the dearth of water in rocks, in spite of higher temperatures. Whistlers detected by Venus Express proved that lightning is present, so the Schumann resonances ～10–40 Hz may provide a global source of electromagnetic energy that penetrates ～10–100 km. Electrical conductivity will be sensitive at these depths to temperature structure and hence thermal lithospheric thickness. Using 1D analytic and 2D numerical models, we demonstrate that the Schumann resonances—transverse EM waves in the ground-ionosphere waveguide—remain sensitive at all altitudes to the properties of the boundaries. This is in marked contrast to other EM methods in which sensitivity to the ground falls off sharply with altitude. We develop a 1D analytical model for aerial EM sounding that treats the electrical properties of the subsurface (thermal gradient, water content, and presence of conductive crust) and ionosphere, and the effects of both random errors and biases that can influence the measurements. We initially consider specified 1D lithospheric thicknesses 100–500 km, but we turn to 2D convection models with Newtonian temperature-dependent viscosity to provide representative vertical and lateral temperature variations. We invert for the conductivity-depth structure and then temperature gradient. For a dry Venus, we find that the error on temperature gradient obtained from any single local measurement is ～100%—perhaps enough to distinguish “thick” vs. “thin” lithospheres. When averaging over thousands of kilometers, however, the standard deviation of the recovered thermal gradient is within the natural variability of the convection models, <25%. A “wet” interior (hundreds of ppm H₂O) limits EM sounding depths using the Schumann resonances to <20 km, and errors are too large to estimate lithospheric properties. A 30-km conductive crust has little influence on the dry-interior models because the Schumann penetration depths are significantly larger. We conclude that EM sounding of the interior of Venus is feasible from a 55-km high balloon. Lithospheric thickness can be measured if the upper-mantle water content is low. If H₂O at hundreds of ppm is present, the deeper, temperature-sensitive structure is screened, but the “wet” nature of the upper mantle, as well as structure of the upper crust, is revealed.     

Dynamical and collisional evolution of Halley-type comets

The number of observed Halley-type comets is hundreds of times less than predicted by models (Levison, H.F., Dones, L., Duncan, M.J. [2001]. Astron. J. 121, 2253–2267). In this paper we investigate the impact of collisions with planetesimals on the evolution of Halley-type comets. First we compute the dynamical evolution of a sub-set of 21 comets using the Mercury integrator package over 100 Myr. The dynamical lifetime is determined to be of the order of 10⁵–10⁶ years in agreement with previous work. The collisional probability of Halley-type comets colliding with known asteroids, a simulated population of Kuiper-belt objects, and planets, is calculated using a modified, Öpik-based collision code. Our results show that the catastrophic disruption of the cometary nucleus has a very low probability of occurring, and disruption through cumulative minor impacts is concluded to be negligible. The dust mantle formed from ejected material falling back to the comet’s surface is calculated to be less than a few centimeters thick, which is insignificant compared to the mantle formed by volatile depletion, while planetary encounters were found to be a negligible disruption mechanism.     

Extreme Ultraviolet Variability Experiment (EVE) Multiple?EUV Grating Spectrographs (MEGS): Radiometric Calibrations and Results

The NASA Solar Dynamics Observatory (SDO), scheduled for launch in early 2010, incorporates a suite of instruments including the Extreme Ultraviolet Variability Experiment (EVE). EVE has multiple instruments including the Multiple Extreme ultraviolet Grating Spectrographs (MEGS) A, B, and P instruments, the Solar Aspect Monitor (SAM), and the Extreme ultraviolet SpectroPhotometer (ESP). The radiometric calibration of EVE, necessary to convert the instrument counts to physical units, was performed at the National Institute of Standards and Technology (NIST) Synchrotron Ultraviolet Radiation Facility (SURF III) located in Gaithersburg, Maryland. This paper presents the results and derived accuracy of this radiometric calibration for the MEGS A, B, P, and SAM instruments, while the calibration of the ESP instrument is addressed by Didkovsky et?al. (Solar Phys., 2010, doi: 10.1007/s11207-009-9485-8 ). In addition, solar measurements that were taken on 14 April 2008, during the NASA 36.240 sounding-rocket flight, are shown for the prototype EVE instruments.