R-band light curve of Comet 9P/Tempel 1 during the Deep Impact event

We present high-speed CCD photometry of Comet 9P/Tempel 1 during the Deep Impact event on 2005 July 4 UT. Approximately 2 h and 50 min of R-band data were acquired at Mount Laguna Observatory with a temporal resolution of 5.5 s. The flux increased by 9% in the first minute after impact. This was followed by a more gradual two-part linear rise, with a change in slope at 9.2 min post-impact, at which time the rate of brightening increased from ∼ to ∼. An analysis of the light curve obtained with the guide camera on the United Kingdom Infrared Telescope and yields very similar results. These findings are mildly in disagreement with the 3-part linear rise found by Fernández et al. (2007) in that we do not find any evidence for a change at 4 min post-impact. We interpret the linear rise phase as due to solar illumination of the edge of an expanding optically thick dust ejecta plume. After approximately 20 min, the light curves begin to flatten out, perhaps coincident with the start of the transition to becoming optically thin. In the large apertures (>10^″) the light curve continues to gradually rise until the end of the observations. In smaller apertures, the light curves reach a peak at approximately 50 min, then decrease back towards the pre-impact flux level. The drop in flux in the smaller apertures may be caused by the ejecta expanding beyond the edge of the photometric aperture, and if so, we can use this timescale to infer an expansion velocity of ∼, consistent with previous published estimates.     

Identifying storm flow pathways in a rainforest catchment using hydrological and geochemical modelling

The hydrological model TOPMODEL is used to assess the water balance and describe flow paths for the 9·73 ha Lutz Creek Catchment in Central Panama. Monte Carlo results are evaluated based on their fit to the observed hydrograph, catchment‐averaged soil moisture and stream chemistry. TOPMODEL, with a direct‐flow mechanism that is intended to route water through rapid shallow‐soil flow, matched observed chemistry and discharge better than the basic version of TOPMODEL and provided a reasonable fit to observed soil moisture and wet‐season discharge at both 15‐min and daily time‐steps. The improvement of simulations with the implementation of a direct‐flow component indicates that a storm flow path not represented in the original version of TOPMODEL plays a primary role in the response of Lutz Creek Catchment. This flow path may be consistent with the active and abundant pipeflow that is observed or delayed saturation overland flow. The ‘best‐accepted’ simulations from 1991 to 1997 indicate that around 41% of precipitation becomes direct flow and around 10% is saturation overland flow. Other field observations are needed to constrain evaporative and groundwater losses in the model and to characterize chemical end‐members posited in this paper. Published in 2004 by John Wiley & Sons, Ltd.     

Limiting heave response of a wave energy device by draft adjustment with upper surface immersion

Many studies have been published concerning the influence of the immersed shape (in still water) of a floating body on its response and power capture from ocean waves. With a few notable exceptions, much of this analysis has assumed small amplitude motion and linear models have been employed to predict response. The form of the upper surface of such a body has received little attention. Here, we show how the shape of the upper (top) surface of a floating body can be designed to ensure that the response amplitude of the body is within a specified value. This is of considerable importance to the survivability of wave energy devices. The approach used is to achieve a large increase of both natural period and hydrodynamic damping for only a small change of float mass. These two factors impose a hydrodynamic limit on the displacement which may be exploited to avoid the ‘end-stop’ problem often encountered in wave device design. To demonstrate the change of response, experimental measurements are presented of the response of an axisymmetric float with rounded base and conical upper surface with rounded perimeter due to a range of regular, irregular and focused wave conditions. Power extraction is not considered since the mechanically undamped response represents the worst case. In contrast to a simple, straight-sided axisymmetric float, a smaller change of mass is required to satisfy a particular response amplitude limit. Although a significant reduction is not expected, hydrodynamic damping may reduce with increasing physical scale, and this remains to be quantified.     

Assessing plot-scale impacts of land use on overland flow generation in Central Panama

Land use in Panama has changed dramatically with ongoing deforestation and conversion to cropland and cattle pastures, potentially altering the soil properties that drive the hydrological processes of infiltration and overland flow. We compared plot-scale overland flow generation between hillslopes in forested and actively cattle-grazed watersheds in Central Panama. Soil physical and hydraulic properties, soil moisture and overland flow data were measured along hillslopes of each land-use type. Soil characteristics and rainfall data were input into a simple, 1-D representative model, HYDRUS-1D, to simulate overland flow that we used to make inferences about overland flow response at forest and pasture sites. Runoff ratios (overland flow/rainfall) were generally higher at the pasture site, although no overall trends were observed between rainfall characteristics and runoff ratios across the two land uses at the plot scale. Saturated hydraulic conductivity (K_s) and bulk density were different between the forest and pasture sites (p < 10⁻⁴). Simulating overland flow in HYDRUS-1D produced more outputs similar to the overland flow recorded at the pasture site than the forest site. Results from our study indicate that, at the plot scale, Hortonian overland flow is the main driver for overland flow generation at the pasture site during storms with high-rainfall totals. We infer that the combination of a leaf litter layer and the activation of shallow preferential flow paths resulting in shallow saturation-excess overland flow are likely the main drivers for plot scale overland flow generation at the forest site. Results from this study contribute to the broader understanding of the delivery of freshwater to streams, which will become increasingly important in the tropics considering freshwater resource scarcity and changing storm intensities.     

Shear zone vs folding origin for spiral inclusion trails in the Canton Schist

The three-dimensional geometry of spiral inclusion trails from the Canton Schist were measured to determine whether the spirals were a product of porphyroblast rotation within a shear zone, or porphyroblast growth during a series of overprinting fold events. The spiral inclusion trails are composed of three separate, sub-planar inclusion trail surfaces occupying texturally distinct parts of the porphyroblasts. These surfaces are correlated across a >10 km² area using textural criteria and relative timing. Measured patterns of inclusion trail orientation and asymmetry suggest they did not form by porphyroblast rotation within a non-coaxial shear zone. Rather, the porphyroblasts grew during three successive overprinting fold events (F₂–F₄), and the spiral inclusion trails represent the accumulated curvature associated with folding of successive axial plane foliations. The data show that spiral garnets are not peculiar to shear zones, and can form by overprinting crenulations and folds. This is consistent with the common occurrence of spiral garnets in multiply-deformed, regionally metamorphosed fold belts.     

Measurement and correlation of microstructures: the case of foliation intersection axes

Recent studies have used the relative rotation axis of sigmoidal and spiral‐shaped inclusion trails, known as Foliation Inflexion/Intersection Axis (FIA), to investigate geological processes such as fold mechanisms and porphyroblast growth. The geological usefulness of this method depends upon the accurate measurement of FIA orientations and correct correlation of temporally related FIAs. This paper uses new data from the Canton Schist to assess the variation in FIA orientations within and between samples, and evaluates criteria for correlating FIAs. For the first time, an entire data set of FIA measurements is published, and data are presented in a way that reflects the variation in FIA orientations within individual samples and provides an indication of the reliability of the data. Analysis of 61 FIA trends determined from the Canton Schist indicate a minimum intrasample range in FIA orientations of 30°. Three competing models are presented for correlation of these FIAs, and each of the models employ different correlation criteria. Correlation of FIAs in Model 1 is based on relative timing and textural criteria, while Model 2 uses relative timing, orientation and patterns of changing FIA orientations, and Model 3 uses relative timing and FIA orientation as correlation criteria. Importantly, the three models differ in the spread of FIA orientations within individual sets, and the number of sets distinguished in the data. Relative timing is the most reliable criterion for correlation, followed by textural criteria and patterns of changing FIA orientations from core to rim of porphyroblasts. It is proposed that within a set of temporally related FIAs, the typical spread of orientations involves clustering of data in a 60° range, but outliers occur at other orientations including near‐normal to the peak distribution. Consequently, in populations of FIA data that contain a wide range of orientations, correlation on the basis of orientation is unreliable in the absence of additional criteria. The results of this study suggest that FIAs are best used as semiquantitative indicators of bulk trends rather than an exact measurement for the purpose of quantitative analyses.     

Variability in the H3 emission of Saturn: Consequences for ionisation rates and temperature

We present an analysis of observations of the auroral/polar regions of Saturn, carried out in 1999, 2004 and 2005, making use of the facility spectrometer, CGS4, on the United Kingdom Infrared Telescope (UKIRT), Mauna Kea, Hawaii. We obtain temperatures of 380(±70) K in 1999 and 420(±70) K in 2004. (The 2005 data has insufficient spectral resolution for a temperature determination to be made.) Our most probable interpretation is that the temperature of Saturn's auroral/polar H⁺₃ layer should be taken as 400(±50) K. This is lower than the value obtained by Miller et al. [Miller, S., and 10 colleagues, 2000. Philos. Trans. R. Soc. 358, 2485-2502], which is shown to be in error. Our analysis reveals clearly that the line emission due to H⁺₃ varies considerably, showing nearly an order of magnitude increase when one compares the data obtained in 1999 with those of 2004. Our conclusion is that this variability is (mainly) due to the changing H⁺₃ column density. By analogy with modelling results obtained for Jupiter, we estimate that the particle (keV electron) precipitation experienced by Saturn must be ∼20 times greater in 2004 than in 1999, to produce this additional ionisation. The H⁺₃ emission increases, but this is insufficient to offset most of the heating due to the extra particle precipitation, indicating that this ion does not act as a “thermostat” on Saturn, in the same way that it does on Jupiter.     

Saturn's auroral/polar H3 infrared emission: II. A comparison with plasma flow models

We present a detailed analysis of the H⁺₃ intensity and velocity profiles crossing Saturn's auroral/polar region, as described by Stallard et al. [Stallard, T., Miller, S., Melin, H., Lystrup, M., Dougherty, M., Achilleos, N., 2007. Icarus 189, 1-13], with a view to understanding the magnetospheric processes with which they are connected. The data are not consistent with the theory that Saturn's main auroral oval is associated with corotation enforcement currents in the middle magnetosphere. This implies that the main auroral oval can be associated with the open-closed field line boundary [Cowley, S.W.H., Bunce, E.J., O'Rourke, J.M., 2004. J. Geophys. Res. 109. A05212]; a third model, by Sittler et al. [Sittler, E.C., Blanc, M.F., Richardson, J.D., 2006. J. Geophys. Res. 111. A06208] associates the main oval with centrifugal instabilities in the outer magnetosphere, but does not make predictions about ionospheric plasma flows with which we can compare our data. We do, however, tentatively identify emission at latitudes lower than the main auroral oval which may be associated with the corotation enforcement currents in the middle magnetosphere. We also find that at latitudes higher than the main auroral oval there is often a region of the ionosphere that is in rigid corotation with the planet. We suggest that this region corresponds to field lines embedded in the centre of the magnetotail which are shielded from the solar wind such that their rotation is controlled only by the neutral atmosphere.     

Hydrothermal Mn-deposits of the Franciscan Assemblage,II. Isotope and trace element geochemistry,and implications for hydrothermal convection at spreading centers

Major element, trace element and Sr, Nd, Pb and O isotopic data for a Franciscan Mn-deposit suggest an origin by seafloor hydrothermal circulation. Based onQ-mode factor analysis the cherts and Mn-lenses of the Blue Jay mine formed from a combination of 4 components representing 1 biogenic, 1 hydrothermal, and 2 detrital sources. RbSr, UThPb and O isotopic systematics in the Mn-lenses were affected by input from the hydrothermal circulation of material leached from the underlying basalts. Nd isotopic compositions in both cherts and Mn-lenses are identical and within the range measured for Pacific Ocean water suggesting the REE were not mobilized by hydrothermal activity. Correlation of δ¹⁸O with SiO₂ and MnO₂ in the Mn-lenses implies the lenses formed by simple mixing of hydrothermally derived Mn-oxides with seawater and biogenic silica. δ¹⁸O of the cherts is both uniform and depleted relative to DSDP Jurassic cherts but similar to microquartz-bearing cherts of the Monterey Formation: this suggests that diagenetic activity exerted more control on oxygen isotope compositions then hydrothermal alteration or metamorphism. Finally, a well defined RbSr isochron of158 ± 5Myr was obtained for these cherts and opens the possibility of determining absolute radiometric ages for similar cherts throughout the geologic record.