Saturn’s tropospheric composition and clouds from Cassini/VIMS 4.6-5.1 μm nightside spectroscopy

The latitudinal variation of Saturn’s tropospheric composition (NH₃, PH₃ and AsH₃) and aerosol properties (cloud altitudes and opacities) are derived from Cassini/VIMS 4.6-5.1 μm thermal emission spectroscopy on the planet’s nightside (April 22, 2006). The gaseous and aerosol distributions are used to trace atmospheric circulation and chemistry within and below Saturn’s cloud decks (in the 1- to 4-bar region). Extensive testing of VIMS spectral models is used to assess and minimise the effects of degeneracies between retrieved variables and sensitivity to the choice of aerosol properties. Best fits indicate cloud opacity in two regimes: (a) a compact cloud deck centred in the 2.5-2.8 bar region, symmetric between the northern and southern hemispheres, with small-scale opacity variations responsible for numerous narrow light/dark axisymmetric lanes; and (b) a hemispherically asymmetric population of aerosols at pressures less than 1.4 bar (whose exact altitude and vertical structure is not constrained by nightside spectra) which is 1.5-2.0× more opaque in the summer hemisphere than in the north and shows an equatorial maximum between ±10° (planetocentric).Saturn’s NH₃ spatial variability shows significant enhancement by vertical advection within ±5° of the equator and in axisymmetric bands at 23-25°S and 42-47°N. The latter is consistent with extratropical upwelling in a dark band on the poleward side of the prograde jet at 41°N (planetocentric). PH₃ dominates the morphology of the VIMS spectrum, and high-altitude PH₃ at p < 1.3 bar has an equatorial maximum and a mid-latitude asymmetry (elevated in the summer hemisphere), whereas deep PH₃ is latitudinally-uniform with off-equatorial maxima near ±10°. The spatial distribution of AsH₃ shows similar off-equatorial maxima at ±7° with a global abundance of 2-3 ppb. VIMS appears to be sensitive to both (i) an upper tropospheric circulation (sensed by NH₃ and upper-tropospheric PH₃ and hazes) and (ii) a lower tropospheric circulation (sensed by deep PH₃, AsH₃ and the lower cloud deck).     

Exoplanet atmospheres with EChO: spectral retrievals using EChOSim

We demonstrate the effectiveness of the Exoplanet Characterisation Observatory mission concept for constraining the atmospheric properties of hot and warm gas giants and super Earths. Synthetic primary and secondary transit spectra for a range of planets are passed through EChOSim [13] to obtain the expected level of noise for different observational scenarios; these are then used as inputs for the NEMESIS atmospheric retrieval code and the retrieved atmospheric properties (temperature structure, composition and cloud properties) compared with the known input values, following the method of [1]. To correctly retrieve the temperature structure and composition of the atmosphere to within 2 σ, we find that we require: a single transit or eclipse of a hot Jupiter orbiting a sun-like (G2) star at 35 pc to constrain the terminator and dayside atmospheres; 20 transits or eclipses of a warm Jupiter orbiting a similar star; 10 transits/eclipses of a hot Neptune orbiting an M dwarf at 6 pc; and 30 transits or eclipses of a GJ1214b-like planet.     

A modified peak‐bagging technique for fitting low‐𝓁 solar p‐modes

We introduce a modified version of a standard power spectrum ‘peak‐bagging’ technique which is designed to gain some of the advantages that fitting the entire low‐degree p‐mode power spectrum simultaneously would bring, but without the problems involved in fitting a model incorporating many hundreds of parameters. Employing Monte‐Carlo simulations we show that by using this modified fitting code it is possible to determine the true background level in the vicinity of the p‐mode peaks. In addition to this we show how small biases in other mode parameters, which are related to inaccurate estimates of the true background, are also consequently removed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The visibility of low‐frequency solar acoustic modes

We make predictions of the detectability of low‐frequency p modes. Estimates of the powers and damping times of these low‐frequency modes are found by extrapolating the observed powers and widths of higher‐frequency modes with large observed signal‐to‐noise ratios. The extrapolations predict that the low‐frequency modes will have small signal‐to‐noise ratios and narrow widths in a frequency‐power spectrum. Monte Carlo simulations were then performed where timeseries containing mode signals and normally distributed Gaussian noise were produced. The mode signals were simulated to have the powers and damping times predicted by the extrapolations. Various statistical tests were then performed on the frequency‐amplitude spectra formed from these timeseries to investigate the fraction of spectra in which the modes could be detected. The results of these simulations were then compared to the number of p‐modes candidates observed in real Sun‐as‐a‐star data at low frequencies. The fraction of simulated spectra in which modes were detected decreases rapidly as the frequency of modes decreases and so the fraction of simulations in which the low‐frequency modes were detected was very small. However, increasing the signal‐to‐noise (S/N) ratio of the low‐frequency modes by a factor of 2 above the extrapolated values led to significantly more detections. Therefore efforts should continue to further improve the quality of solar data that is currently available. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of the Energy Spectra and Number Fluxes From a simple Flare Model to Observations

In this paper, we investigate the energy spectra produced by a simple test particle X-point model of a solar flare for different configurations of the initial electromagnetic field. We find that once the reconnection electric field is larger than 1 Vm^-1 the particle distribution transits from a heated one to a partially accelerated one. As we close the separatrices of the X-point and the angle in the inflow direction widens we find that more particles are accelerated out of the thermal distribution and this power–law component extends to lower energies. When we introduce a guiding magnetic field component we find that more particles are energised, but only up to a maximum energy dictated primarily by the reconnection electric field. Despite being able to accelerate particles to observable energies and demonstrate behaviour in the energy spectra that is consistent with observations, this single X-line model can only deliver the number fluxes required for microflares.     

What is Moss?

TRACE observations of active regions show a peculiar extreme ultraviolet (EUV) emission over certain plage areas. Termed `moss' for its spongy, low-lying, appearance, observations and modeling imply that the phenomenon is caused by thermal conduction from 3–5 MKcoronal loops overlying the plage: moss is the upper transition region emission of hot coronal loops. The spongy appearance is due to the presence of chromospheric jets or `spicules' interspersed with the EUV emission elements. High cadence TRACE observations show that the moss EUV elements interact with the chromospheric jets on 10 s time scales. The location of EUV emission in the moss does not correlate well to the locations of underlying magnetic elements in the chromosphere and photosphere, implying a complex magnetic topology for coronal loop footpoint regions. We summarize here the key observations leading to these conclusions and discuss new implications for understanding the structuring of the outer solar atmosphere. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005286503963     

Vertical profiles of HCN, HC3N, and C2H2 in Titan's atmosphere derived from Cassini/CIRS data

Mid-infrared limb spectra in the range 600-1400 cm⁻¹ taken with the Composite InfraRed Spectrometer (CIRS) on-board the Cassini spacecraft were used to determine vertical profiles of HCN, HC₃N, C₂H₂, and temperature in Titan's atmosphere. Both high (0.5 cm⁻¹) and low (13.5 cm⁻¹) spectral resolution data were used. The 0.5 cm⁻¹ data gave profiles at four latitudes and the 13.5 cm⁻¹ data gave almost complete latitudinal coverage of the atmosphere. Both datasets were found to be consistent with each other. High temperatures in the upper stratosphere and mesosphere were observed at Titan's northern winter pole and were attributed to adiabatic heating in the subsiding branch of a meridional circulation cell. On the other hand, the lower stratosphere was much colder in the north than at the equator, which can be explained by the lack of solar radiation and increased IR emission from volatile enriched air. HC₃N had a vertical profile consistent with previous ground based observations at southern and equatorial latitudes, but was massively enriched near the north pole. This can also be explained in terms of subsidence at the winter pole. A boundary observed at 60° N between enriched and un-enriched air is consistent with a confining polar vortex at 60° N and HC₃N's short lifetime. In the far north, layers were observed in the HC₃N profile that were reminiscent of haze layers observed by Cassini's imaging cameras. HCN was also enriched over the north pole, which gives further evidence for subsidence. However, the atmospheric cross section obtained from 13.5 cm⁻¹ data indicated a HCN enriched layer at 200-250 km, extending into the southern hemisphere. This could be interpreted as advection of polar enriched air towards the south by a meridional circulation cell. This is observed for HCN but not for HC₃N due to HCN's longer photochemical lifetime. C₂H₂ appears to have a uniform abundance with altitude and is not significantly enriched in the north. This is consistent with observations from previous CIRS analysis that show increased abundances of nitriles and hydrocarbons but not C₂H₂ towards the north pole.     

The Magnetic Structure and Generation of EUV Flare Ribbons

The `ribbons' of two-ribbon flares show complicated patterns reflecting the linkages of coronal magnetic field lines through the lower solar atmosphere. We describe the morphology of the EUV ribbons of the July 14, 2000 flare, as seen in SOHO, TRACE, and Yohkoh data, from this point of view. A successful co-alignment of the TRACE, SOHO/MDI and Yohkoh/HXT data has allowed us to locate the EUV ribbon positions on the underlying field to within ∼ 2′′, and thus to investigate the relationship between the ribbons and the field, and also the sites of electron precipitation. We have also made a determination of the longitudinal magnetic flux involved in the flare reconnection event, an important parameter in flare energetic considerations. There are several respects in which the observations differ from what would be expected in the commonly-adopted models for flares. Firstly, the flare ribbons differ in fine structure from the (line-of-sight) magnetic field patterns underlying them, apparently propagating through regions of very weak and probably mixed polarity. Secondly, the ribbons split or bifurcate. Thirdly, the amount of line-of-sight flux passed over by the ribbons in the negative and positive fields is not equal. Fourthly, the strongest hard X-ray sources are observed to originate in stronger field regions. Based on a comparison between HXT and EUV time-profiles we suggest that emission in the EUV ribbons is caused by electron bombardment of the lower atmosphere, supporting the hypothesis that flare ribbons map out the chromospheric footpoints of magnetic field lines newly linked by reconnection. We describe the interpretation of our observations within the standard model, and the implications for the distribution of magnetic fields in this active region.     

A Relationship Between Transition Region Brightenings,Abundances, and Magnetic Topology

We present multi-instrument observations of active region (AR) 8048, made between 3 June and 5 June 1997, as part of the SOHO Joint Observing Program 33. This AR has a sigmoid-like global shape and undergoes transient brightenings in both soft X-rays and transition region (TR) lines. We compute a magneto-hydrostatic model of the AR magnetic field, using as boundary condition the photospheric observations of SOHO/MDI. The computed large-scale magnetic field lines show that the large-scale sigmoid is formed by two sets of coronal loops. Shorter loops, associated with the core of the SXT emission, coincide with the loops observed in the hotter CDS lines. These loops reveal a gradient of temperature, from 2 MK at the top to 1 MK at the ends. The field lines most closely matching these hot loops extend along the quasi-separatrix layers (QSLs) of the computed coronal field. The TR brightenings observed with SOHO/CDS can also be associated with the magnetic field topology, both QSL intersections with the photosphere, and places where separatrices issuing from bald patches (sites where field lines coming from the corona are tangent to the photosphere) intersect the photosphere. There are, furthermore, suggestions that the element abundances measured in the TR may depend on the type of topological structure present. Typically, the TR brightenings associated with QSLs have coronal abundances, while those associated with BP separatrices have abundances closer to photospheric values. We suggest that this difference is due to the location and manner in which magnetic reconnection occurs in two different topological structures. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1013302317042     

Precise measurement of seismic travel times from an air gun source

Summary. Seismic travel times for extrema, zero-crossings, or entire body waves need to be determined precisely to one part in 10³ or better in several varieties of seismic studies employing an impulsive artificial source. Examples are crosshole surveys which delineate rock crack distribution separating the holes and monitoring of crustal seismic travel times in earthquake precursor studies. A timing resolution of one part in 10³ has been achieved previously using digitally recorded seismic data. These methods, however, do not use interpolation between digitized data points as a method to increase the timing resolution. We report travel-time determinations based on interpolation between digitized points which achieve a precision of two parts in 10⁴, a five-fold improvement over the existing methods. In addition, the effects of seismic noise on travel-time measurement have been compared for the extremum location, the unnormalized correlation, and the normalized correlation method. The following conclusions are drawn from this comparison: (1) the normalized correlation method provides an 18–55 per cent improvement in the standard deviation of the mean over the extremum location method, and (2) results as accurate as those by the normalized correlation can be obtained by the unnormalized correlation if a complete up-and-down swing of the waveform is used as the master trace and if the master trace is close to being sinusoidal. The advantage of the unnormalized correlation over the normalized correlation is speed; the unnormalized correlation is faster by a factor of 28 in computing time.