Response to “Comment on “Slow-mode shock candidate in the Jovian magnetosheath” by Bebesi et al.”

Hubert and Samsonov addressed our paper published in early 2010 (Bebesi et al., 2010) about a plasma event detected in the magnetosheath of Jupiter by the plasma instruments of the Cassini spacecraft. We proposed that the characteristics of the plasma parameter variations were consistent with a slow mode shock (SMS). As our title indicated, we claimed only that the event was a possible “candidate” for an SMS according to our data, which had some restrictions as discussed in the paper. As to the origin, we proposed the following: since there was a crossing of the then highly tilted Heliospheric Current Sheet in less than two days before the event, it might have played a role in initiating the shock front.We highly appreciate the opinion of the authors, but they do not point out any hard fact that would exclude the possibility of the scenario we suggested.     

A primordial origin for the atmospheric methane of Saturn’s moon Titan

The origin of Titan’s atmospheric methane is a key issue for understanding the origin of the saturnian satellite system. It has been proposed that serpentinization reactions in Titan’s interior could lead to the formation of the observed methane. Meanwhile, alternative scenarios suggest that methane was incorporated in Titan’s planetesimals before its formation. Here, we point out that serpentinization reactions in Titan’s interior are not able to reproduce the deuterium over hydrogen (D/H) ratio observed at present in methane in its atmosphere, and would require a maximum D/H ratio in Titan’s water ice 30% lower than the value likely acquired by the satellite during its formation, based on Cassini observations at Enceladus. Alternatively, production of methane in Titan’s interior via radiolytic reactions with water can be envisaged but the associated production rates remain uncertain. On the other hand, a mechanism that easily explains the presence of large amounts of methane trapped in Titan in a way consistent with its measured atmospheric D/H ratio is its direct capture in the satellite’s planetesimals at the time of their formation in the solar nebula. In this case, the mass of methane trapped in Titan’s interior can be up to ∼1300 times the current mass of atmospheric methane.     

Low-order dynamical behavior in the martian atmosphere: Diagnosis of general circulation model results

The hypothesis of a low dimensional martian climate attractor is investigated by the application of the proper orthogonal decomposition (POD) to a simulation of martian atmospheric circulation using the UK Mars general circulation model (UK-MGCM). In this article we focus on a time series of the interval between autumn and winter in the northern hemisphere, when baroclinic activity is intense. The POD is a statistical technique that allows the attribution of total energy (TE) to particular structures embedded in the UK-MGCM time-evolving circulation. These structures are called empirical orthogonal functions (EOFs). Ordering the EOFs according to their associated energy content, we were able to determine the necessary number to account for a chosen amount of atmospheric TE. We show that for Mars a large fraction of TE is explained by just a few EOFs (with 90% TE in 23 EOFs), which apparently support the initial hypothesis. We also show that the resulting EOFs represent classical types of atmospheric motion, such as thermal tides and transient waves. Thus, POD is shown to be an efficient method for the identification of different classes of atmospheric modes. It also provides insight into the non-linear interaction of these modes.     

The effect of spacecraft radiation sources on electron moments from the Cassini CAPS electron spectrometer

Data from the Cassini plasma spectrometer (CAPS) electron spectrometer (ELS) have been found to be contaminated with an energy-independent background count rate which has been associated with radiation sources on Cassini. In this paper we describe this background radiation and quantitatively assess its impact on numerically integrated electron moments. The general properties of such a background and its effects on numerical moments are derived. The properties of the ELS background are described and a model for the background presented. A model to generate synthetic ELS spectra is presented and used to evaluate the density and temperature of pure noise and then extended to include ambient distributions. It is shown that the presence of noise produces a saturation of the electron density and temperature at quasi-constant values when the instrument is at background, but that these noise level moments are dependent on the floating spacecraft potential and the orientation of the ELS instrument with respect to the spacecraft. When the ambient distribution has a poor signal-to-noise ratio (SNR) the noise determines the density and temperature; however, as the SNR increases (increasing primarily with density) the density and temperature tend to those of the ambient distribution. It is also shown that these noise effects produce highly artificial density-temperature inverse correlations. A method to subtract this noise is presented and shown to correct for the presence of the noise. Simulated error estimates for the density and temperature are also presented. The analysis described in this paper not only applies to weak background noise, but also to more significant penetrating backgrounds such as those in radiation belt regions of planetary magnetospheres.     

Galaxy shear estimation from stacked images

Statistics of the weak lensing of galaxies can be used to constrain cosmology if the galaxy shear can be estimated accurately. In general, this requires accurate modelling of unlensed galaxy shapes and the point spread function (PSF). I discuss suboptimal but potentially robust methods for estimating galaxy shear by stacking images such that the stacked-image distribution is closely Gaussian by the central limit theorem. The shear can then be determined by radial fitting, requiring only an accurate model of the PSF rather than also needing to model each galaxy accurately. When noise is significant, asymmetric errors in the centroid must be corrected, but the method may ultimately be able to give accurate unbiased results when there is a high galaxy density with constant shear. It provides a useful baseline for more optimal methods, and a test case for estimating biases, though the method is not directly applicable to realistic data. I test stacking methods on the simple toy simulations with constant PSF and shear provided by the GREAT08 project, on which most other existing methods perform significantly more poorly, and briefly discuss generalizations to more realistic cases. In the Appendix , I discuss a simple analytic galaxy population model where stacking gives optimal errors in a perfect ideal case.     

Early Tertiary fluid flow and pressure–temperature conditions in the Shimanto accretionary complex of south-west Japan: constraints from fluid inclusions

The pressure and temperature history of the Tertiary Shimanto belt of south-west Japan has been elucidated by analysing fluids trapped in quartz crystals which grew syn-kinematically along late-stage brittle faults. The samples come from three areas that span the Paleogene exposures on the Muroto Peninsula of Shikoku Island. Applying microthermometric and laser Raman microsampling techniques to coeval water-rich and carbonic fluid inclusions, we have constrained the pressure and temperature conditions that accompanied a widespread and kinematically distinct phase of deformation. The results suggest elevated geothermal gradients during late-stage deformation, conditions that are in disaccord with previous plate reconstructions that have depicted old, thermally mature Pacific crust subducting beneath Eurasia during the early to middle Tertiary. These conditions can most easily be accounted for by including an additional plate boundary in the western Pacific during Paleogene time. Plate reconstructions that include the Kula plate in this region are therefore consistent with our findings. In addition, our results provide clues to the conditions that likely accompany seismogenic deformation at active convergent plate boundaries.     

A finite element solution of a fully coupled implicit formulation for reservoir simulation

An iterative method is presented for solving a fully coupled and implicit formulation of fluid flow in a porous medium. The mathematical model describes a set of fully coupled three-phase flow of compressible and immiscible fluids in a saturated oil reservoir. The finite element method is applied to obtain the simultaneous solution (SS) for the resulting highly non-linear partial differential equations where fluid pressures are the primary unknowns. The final discretized equations are solved iteratively by using a fully implicit numerical scheme. Several examples, illustrating the use of the present model, are described. The increased stability achieved with this scheme has permitted the use of larger time steps with smaller material balance errors.     

Radio sources in the 2dF Galaxy Redshift Survey – II. Local radio luminosity functions for AGN and star-forming galaxies at 1.4 GHz

We have cross-matched the 1.4-GHz NRAO VLA Sky Survey (NVSS) with the first 210 fields observed in the 2dF Galaxy Redshift Survey (2dFGRS), covering an effective area of 325 deg² (about 20 per cent of the final 2dFGRS area). This yields a set of optical spectra of 912 candidate NVSS counterparts, of which we identify 757 as genuine radio identifications – the largest and most homogeneous set of radio source spectra ever obtained. The 2dFGRS radio sources span the redshift range     to 0.438, and are a mixture of active galaxies (60 per cent) and star-forming galaxies (40 per cent). About 25 per cent of the 2dFGRS radio sources are spatially resolved by NVSS, and the sample includes three giant radio galaxies with projected linear size greater than 1 Mpc. The high quality of the 2dF spectra means we can usually distinguish unambiguously between AGN and star-forming galaxies. We make a new determination of the local radio luminosity function at 1.4 GHz for both active and star-forming galaxies, and derive a local star formation density of         .