Hard x-ray and Microwave Flux Spectra of the 2 November 1991 Solar Flare

We analysed the hard X-ray and microwave flux spectra of the solar flare (BATSE No. 1791) on 2 November 1991, which started at 16:11:03 UT and ended at 16:56:10 UT. This flare is particularly interesting because of its deep cyclic intensity modulation. Data are available simultaneously from the 16-channel BATSE/LAD hard X-ray and 45-frequency OVRO microwave database. We quantitatively compare the time variations in profiles of the hard X-ray spectral photon index, the 50 keV X-ray flux, and microwave spectral indices (at both high and low frequency ends). As expected, the X-ray photon spectral index decreases as the hard X-ray flux increases. This pattern appears in all the sub-peaks. This is consistent with previous observations that hard X-ray emission hardens at the emission peak. However, the behaviour of the high-frequency microwave index is unexpected. We observe an anti-correlation between the high-frequency microwave index and the hard X-ray photon index during the course of the flare. Finally, we study the arrival time of microwave flux peaks as a function of frequency and find that the microwave peak at a higher frequency comes earlier than that at a lower frequency. A maximum delay of 72 s is found among the main peaks at different frequencies. Shorter delays are found for the other five sub-peaks.     

Galaxy activity in merging binary galaxy clusters

We present the results of a study of galaxy activity in two merging binary clusters (A168 and A1750) using the Sloan Digital Sky Survey (SDSS) data supplemented with the data in the literature. We have investigated the merger histories of A168 and A1750 by combining the results from a two-body dynamical model and X-ray data. In A168, two subclusters appear to have passed each other and to be coming together from the recent maximum separation. In A1750, two major subclusters appear to have started interaction and to be coming together for the first time. We find an enhanced concentration of the galaxies showing star formation (SF) or active galactic nuclei (AGN) activity in the region between two subclusters in A168, which were possibly triggered by the cluster merger. In A1750, we do not find any galaxies with SF/AGN activity in the region between two subclusters, indicating that two major subclusters are in the early stage of merging.     

Spaceborne visible and thermal infrared lithologic mapping of impact‐exposed subsurface lithologies at the Haughton impact structure,Devon Island,Canadian High Arctic: Applications to Mars

Abstract— This study serves as a proof‐of‐concept for the technique of using visible‐near infrared (VNIR), short‐wavelength infrared (SWIR), and thermal infrared (TIR) spectroscopic observations to map impact‐exposed subsurface lithologies and stratigraphy on Earth or Mars. The topmost layer, three subsurface layers and undisturbed outcrops of the target sequence exposed just 10 km to the northeast of the 23 km diameter Haughton impact structure (Devon Island, Nunavut, Canada) were mapped as distinct spectral units using Landsat 7 ETM+ (VNIR/SWIR) and ASTER (VNIR/SWIR/TIR) multispectral images. Spectral mapping was accomplished by using standard image contrast‐stretching algorithms. Both spectral matching and deconvolution algorithms were applied to image‐derived ASTER TIR emissivity spectra using spectra from a library of laboratory‐measured spectra of minerals (Arizona State University) and whole‐rocks (Ward's). These identifications were made without the use of a priori knowledge from the field (i.e., a “blind” analysis). The results from this analysis suggest a sequence of dolomitic rock (in the crater rim), limestone (wall), gypsum‐rich carbonate (floor), and limestone again (central uplift). These matched compositions agree with the lithologic units and the pre‐impact stratigraphic sequence as mapped during recent field studies of the Haughton impact structure by Osinski et al. (2005a). Further conformation of the identity of image‐derived spectra was confirmed by matching these spectra with laboratory‐measured spectra of samples collected from Haughton. The results from the “blind” remote sensing methods used here suggest that these techniques can also be used to understand subsurface lithologies on Mars, where ground truth knowledge may not be generally available.     

Comparison of Five Numerical Codes for Automated Tracing of Coronal Loops

The three-dimensional (3D) modeling of coronal loops and filaments requires algorithms that automatically trace curvilinear features in solar EUV or soft X-ray images. We compare five existing algorithms that have been developed and customized to trace curvilinear features in solar images: i) the oriented-connectivity method (OCM), which is an extension of the Strous pixel-labeling algorithm (developed by Lee, Newman, and Gary); ii) the dynamic aperture-based loop-segmentation method (developed by Lee, Newman, and Gary); iii) unbiased detection of curvilinear structures (developed by Steger, Raghupathy, and Smith); iv) the oriented-direction method (developed by Aschwanden); and v) ridge detection by automated scaling (developed by Inhester). We test the five existing numerical codes with a TRACE image that shows a bipolar active region and contains over 100 discernable loops. We evaluate the performance of the five codes by comparing the cumulative distribution of loop lengths, the median and maximum loop length, the completeness or detection efficiency, the accuracy, and flux sensitivity. These algorithms are useful for the reconstruction of the 3D geometry of coronal loops from stereoscopic observations with the STEREO spacecraft, or for quantitative comparisons of observed EUV loop geometries with (nonlinear force-free) magnetic field extrapolation models.     

A bulk cloud parameterization in a Venus General Circulation Model

A condensing cloud parameterization is included in a super-rotating Venus General Circulation Model. A parameterization including condensation, evaporation and sedimentation of mono-modal sulfuric acid cloud particles is described. Saturation vapor pressure of sulfuric acid vapor is used to determine cloud formation through instantaneous condensation and destruction through evaporation, while pressure dependent viscosity of a carbon dioxide atmosphere is used to determine sedimentation rates assuming particles fall at their terminal Stokes velocity. Modifications are described to account for the large range of the Reynolds number seen in the Venus atmosphere.Two GCM experiments initialized with 10 ppm-equivalent of sulfuric acid are integrated for 30 Earth years and the results are discussed with reference to “Y” shaped cloud structures observed on Venus. The GCM is able to produce an analog of the “Y” shaped cloud structure through dynamical processes alone, with contributions from the mean westward wind, the equatorial Kelvin wave, and the mid-latitude/polar Mixed Rossby/Gravity waves. The cloud top height in the GCM decreases from equator to pole and latitudinal gradients of cloud top height are comparable to those observed by Pioneer Venus and Venus Express, and those produced in more complex microphysical models of the sulfur cycle on Venus. Differences between the modeled cloud structures and observations are described and dynamical explanations are suggested for the most prominent differences.     

The roughness of the dark side of Iapetus from the 2004 to 2005 flyby

The roughness of a planetary surface offers clues to its past geologic history. We apply a surface roughness model developed by Buratti and Veverka (Buratti, B.J., Veverka, J. [1985]. Icarus 64, 320-328) to Cassini ISS data from the January 1st, 2005 flyby of Iapetus. This model uses the observed scattering behavior to provide a depth to radius factor q quantifying the size of idealized craters on the surface. Our findings indicate that the surface on the dark side is significantly smoother than the surfaces of other icy low-albedo saturnian satellites. We have found that the average depth to radius on the leading (dark) side is 0.084, corresponding to a Hapke mean slope angle of 6°. As compared to the 13-33° Hapke mean slope angle of other icy satellites (Buratti, B.J., and 10 colleagues [2008]. Icarus 193, 309-322), our results present a clearly different picture for the leading surface of Iapetus, suggesting that the dark deposit contributes to the decrease in macroscopic surface roughness of the leading side. Attempts were made to obtain an average depth to radius value for the trailing (bright) side; however the scans of the bright side from this flyby exhibited large variations in albedo, resulting in results that were physically unrealistic.     

Mass density of the upper atmosphere derived from Starlette’s Precise Orbit Determination with Satellite Laser Ranging

The atmospheric mass density of the upper atmosphere from the spherical Starlette satellite’s Precise Orbit Determination is first derived with Satellite Laser Ranging measurements at 815 to 1115 km during strong solar and geomagnetic activities. Starlette’s orbit is determined using the improved orbit determination techniques combining optimum parameters with a precise empirical drag application to a gravity field. MSIS-86 and NRLMSISE-00 atmospheric density models are compared with the Starlette drag-derived atmospheric density of the upper atmosphere. It is found that the variation in the Starlette’s drag coefficient above 800 km corresponds well with the level of geomagnetic activity. This represents that the satellite orbit is mainly perturbed by the Joule heating from geomagnetic activity at the upper atmosphere. This result concludes that MSIS empirical models strongly underestimate the mass density of the upper atmosphere as compared to the Starlette drag-derived atmospheric density during the geomagnetic storms. We suggest that the atmospheric density models should be analyzed with higher altitude acceleration data for a better understanding of long-term solar and geomagnetic effects.     

The Ashima/MIT Mars GCM and argon in the martian atmosphere

We investigate the ability of modern general circulation models (GCMs) to simulate transport in the martian atmosphere using measurements of argon as a proxy for the transport processes. Argon provides the simplest measure of transport as it is a noble gas with no sinks or sources on seasonal timescales. Variations in argon result solely from ‘freeze distillation’, as the atmosphere condenses at the winter poles, and from atmospheric transport. Comparison of all previously published models when rescaled to a common definition of the argon enhancement factor (EF) suggest that models generally do a poor job in predicting the peak enhancement in southern winter over the winter pole – the time when the capability of the model transport approaches are most severely tested. Despite observed peak EF values of ～6, previously published model predictions peaked at EF values of only 2–3. We introduce a new GCM that provides a better treatment of mass conservation within the dynamical core, includes more sophisticated tracer transport approaches, and utilizes a cube–sphere grid structure thus avoiding the grid-point convergence problem at the pole that exists for most current Mars GCMs. We describe this model – the Ashima Research/Massachusetts Institute of Technology Mars General Circulation Model (Ashima/MIT Mars GCM) and use it to demonstrate the significant sensitivity of peak EF to the choices of transport approach for both tracers and heat. We obtain a peak EF of 4.75 which, while over 50% higher than any prior model, remains well short of the observed value. We show that the polar EF value in winter is primarily determined by the competition between two processes: (1) mean meridional import of lower-latitude air not enriched in argon and (2) the leakage of enriched argon out of the polar column by eddies in the lowest atmospheric levels. We suggest possibilities for improving GCM representation of the CO₂ cycle and the general circulation that may further improve the simulation of the argon cycle. We conclude that current GCMs may be insufficient for detailed simulation of transport-sensitive problems like the water cycle and potentially also the dust cycle.     

Interpretation of combined infrared,submillimeter, and millimeter thermal flux data obtained during the Rosetta fly-by of Asteroid (21) Lutetia

The European Space Agency’s Rosetta spacecraft is the first Solar System mission to include instrumentation capable of measuring planetary thermal fluxes at both near-IR (VIRTIS) and submillimeter–millimeter (smm–mm, MIRO) wavelengths. Its primary mission is a 1 year reconnaissance of Comet 67P/Churyumov–Gerasimenko beginning in 2014. During a 2010 close fly-by of Asteroid 21 Lutetia, the VIRTIS and MIRO instruments provided complementary data that have been analyzed to produce a consistent model of Lutetia’s surface layer thermal and electrical properties, including a physical model of self-heating. VIRTIS dayside measurements provided highly resolved 1 K accuracy surface temperatures that required a low thermal inertia, I < 30 J/(K m² s^0.5). MIRO smm and mm measurements of polar night thermal fluxes produced constraints on Lutetia’s subsurface thermal properties to depths comparable to the seasonal thermal wave, yielding a model of I < 20 J/(K m² s^0.5) in the upper few centimeters, increasing with depth in a manner very similar to that of Earth’s Moon. Subsequent MIRO-based model predictions of the dayside surface temperatures reveal negative offsets of ～5–30 K from the higher VIRTIS-measurements. By adding surface roughness in the form of 50% fractional coverage of hemispherical mini-craters to the MIRO-based thermal model, sufficient self-heating is produced to largely remove the offsets relative to the VIRTIS measurements and also reproduce the thermal limb brightening features (relative to a smooth surface model) seen by VIRTIS. The Lutetia physical property constraints provided by the VIRTIS and MIRO data sets demonstrate the unique diagnostic capabilities of combined infrared and submillimeter/millimeter thermal flux measurements.     

Volcanic hazards from Bezymianny- and Bandai-type eruptions

Major slope failures are a significant degradational process at volcanoes. Slope failures and associated explosive eruptions have resulted in more than 20 000 fatalities in the past 400 years; the historic record provides evidence for at least six of these events in the past century. Several historic debris avalanches exceed 1 km³ in volume. Holocene avalanches an order of magnitude larger have traveled 50–100 km from the source volcano and affected areas of 500–1500 km². Historic eruptions associated with major slope failures include those with a magmatic component (Bezymianny type) and those solely phreatic (Bandai type). The associated gravitational failures remove major segments of the volcanoes, creating massive horseshoe-shaped depressions commonly of caldera size. The paroxysmal phase of a Bezymianny-type eruption may include powerful lateral explosions and pumiceous pyroclastic flows; it is often followed by construction of lava dome or pyroclastic cone in the new crater. Bandai-type eruptions begin and end with the paroxysmal phase, during which slope failure removes a portion of the edifice. Massive volcanic landslides can also occur without related explosive eruptions, as at the Unzen volcano in 1792.The main potential hazards from these events derive from lateral blasts, the debris avalanche itself, and avalanche-induced tsunamis. Lateral blasts produced by sudden decompression of hydrothermal and/or magmatic systems can devastate areas in excess of 500km² at velocities exceeding 100 m s^–1. The ratio of area covered to distance traveled for the Mount St. Helens and Bezymianny lateral blasts exceeds that of many pyroclastic flows or surges of comparable volume. The potential for large-scale lateral blasts is likely related to the location of magma at the time of slope failure and appears highest when magma has intruded into the upper edifice, as at Mount St. Helens and Bezymianny.Debris avalanches can move faster than 100 ms^–1 and travel tens of kilometers. When not confined by valley walls, avalanches can affect wide areas beyond the volcano's flanks. Tsunamis from debris avalanches at coastal volcanoes have caused more fatalities than have the landslides themselves or associated eruptions. The probable travel distance (L) of avalanches can be estimated by considering the potential vertical drop (H). Data from a catalog of around 200 debris avalanches indicates that the H/L rations for avalanches with volumes of 0.1–1 km³ average 0.13 and range 0.09–0.18; for avalanches exceeding 1 km³, H/L ratios average 0.09 and range 0.5–0.13.Large-scale deformation of the volcanic edefice and intense local seismicity precede many slope failures and can indicate the likely failure direction and orientation of potential lateral blasts. The nature and duration of precursory activity vary widely, and the timing of slope faliure greatly affects the type of associated eruption. Bandai-type eruptions are particularly difficult to anticipate because they typically climax suddenly without precursory eruptions and may be preceded by only short periods of seismicity.