Radio observations of Comet 9P/Tempel 1 before and after Deep Impact

Comet 9P/Tempel 1 was the target of a multi-wavelength worldwide investigation in 2005. The NASA Deep Impact mission reached the comet on 4.24 July 2005, delivering a 370-kg impactor which hit the comet at 10.3 km s⁻¹. Following this impact, a cloud of gas and dust was excavated from the comet nucleus. The comet was observed in 2005 prior to and after the impact, at 18-cm wavelength with the Nançay radio telescope, in the millimeter range with the IRAM and CSO radio telescopes, and at 557 GHz with the Odin satellite. OH observations at Nançay provided a 4-month monitoring of the outgassing of the comet from March to June, followed by the observation of H₂O with Odin from June to August 2005. The peak of outgassing was found to be around between May and July. Observations conducted with the IRAM 30-m radio telescope in May and July 2005 resulted in detections of HCN, CH₃OH and H₂S with classical abundances relative to water (0.12, 2.7 and 0.5%, respectively). In addition, a variation of the HCN production rate with a period of 1.73±0.10 days was observed in May 2005, consistent with the 1.7-day rotation period of the nucleus. The phase of these variations, as well as those of CN seen in July by Jehin et al. [Jehin, E., Manfroid, J., Hutsemékers, D., Cochran, A.L., Arpigny, C., Jackson, W.M., Rauer, H., Schulz, R., Zucconi, J.-M., 2006. Astrophys. J. 641, L145-L148], is consistent with a rotation period of the nucleus of 1.715 days and a strong variation of the outgassing activity by a factor 3 from minimum to maximum. This also implies that the impact took place on the rising phase of the “natural” outgassing which reached its maximum ≈4 h after the impact. Post-impact observations at IRAM and CSO did not reveal a significant change of the outgassing rates and relative abundances, with the exception of CH₃OH which may have been more abundant by up to one order of magnitude in the ejecta. Most other variations are linked to the intrinsic variability of the comet. The Odin satellite monitored nearly continuously the H₂O line at 557 GHz during the 38 h following the impact on the 4th of July, in addition to weekly monitoring. Once the periodic variations related to the nucleus rotation are removed, a small increase of outgassing related to the impact is present, which corresponds to the release of ≈5000±2000 tons of water. Two other bursts of activity, also observed at other wavelengths, were seen on 23 June and 7 July; they correspond to even larger releases of gas.     

Geology and paleomagnetism of El Potrero pluton, Baja California: Understanding criteria for timing of deformation and evidence of pluton tilt during batholith growth

The 102 Ma El Potrero pluton, in the western foothills of Sierra San Pedro Mártir, in north-central Baja California, was emplaced during a long period of contractional deformation bracketed between 132 and 85 Ma that affected this segment of the Peninsular Ranges Batholith. The pluton records regional and emplacement related deformation manifested by: (1) a solid-state fabric developed on its eastern contact, which is produced by eastward lateral pluton expansion; (2) cleavage triple point zones in the host-rock NW and SE of the pluton; (3) subhorizontal ductile shear zones indicative of top-to-the-east transport; (4) magmatic and tectonic foliations parallel to regional structural trends and regional shear zones; (5) variable axial ratios of microgranitoid enclaves close to pluton–wall rock contacts; (6) evidence of brittle-emplacement mechanisms in the western border of the pluton, which contrast with features indicating mainly ductile mechanisms toward the east; and, (7) markedly discordant paleomagnetic directions that suggest emplacement in an active tectonic setting. The overall mean for 9 accepted paleomagnetic sites is Dec = 34.6°, I = 25.7° (k = 88.3, α₉₅ = 5.5°), and is deviated  35° with respect to the reference cratonic direction. This magnetization is interpreted to indicate a combination of tilt due to initial drag during vertical diapiric ascent (or westward lateral-oblique expansion) of the adjacent San Pedro Mártir pluton and later rotation ( 15°) by Rosarito Fault activity in the southwest; this rotation may have occurred as eastward contraction acted to fill the space emptied by the ascending San Pedro Mártir pluton. The Rosarito fault may have tilted several plutons in the area (Sierra San Pedro Mártir, El Potrero, San José, and Encinosa). Magnetic susceptibility fabrics for 13 sites reflect mostly emplacement-related stress and regional stress. Paleomagnetic data and structural observations lead us to interpret the El Potrero pluton as a syntectonic pluton, emplaced within a regional shear zone delimited by the Main Mártir Thrust and the younger Rosarito Fault.     

Turbulent Pressure and Velocity Perturbations Induced by Gentle Hills Covered with Sparse and Dense Canopies

How the spatial perturbations of the first and second moments of the velocity and pressure fields differ for flow over a train of gentle hills covered by either sparse or dense vegetation is explored using large-eddy simulation (LES). Two simulations are investigated where the canopy is composed of uniformly arrayed rods each with a height that is comparable to the hill height. In the first simulation, the rod density is chosen so that much of the momentum is absorbed within the canopy volume yet the canopy is not dense enough to induce separation on the lee side of the hill. In the second simulation, the rod density is large enough to induce recirculation inside the canopy on the lee side of the hill. For this separating flow case, zones of intense shear stress originating near the canopy-atmosphere interface persist all the way up to the middle layer, ‘contaminating’ much of the middle and outer layers with shear stress gradients. The implications of these persistent shear-stress gradients on rapid distortion theory and phase relationships between higher order velocity statistics and hill-induced mean velocity perturbations (Δu) are discussed. Within the inner layer, these intense shear zones improve predictions of the spatial perturbation by K-theory, especially for the phase relationships between the shear stress (~ ?Δu/?z) and the velocity variances, where z is the height. For the upper canopy layers, wake production increases with increasing leaf area density resulting in a vertical velocity variance more in phase with Δu than with ?Δu/?z. However, background turbulence and inactive eddies may have dampened this effect for the longitudinal velocity variance. The increase in leaf area density does not significantly affect the phase relationship between mean surface pressure and topography for the two simulations, though the LES results here confirm earlier findings that the minimum mean pressure shifts downstream from the hill crest. The increase in leaf area density and associated flow separation simply stretches this difference further downstream. This shift increases the pressure drag, the dominant term in the overall drag on the hill surface, by some 15%. With regards to the normalized pressure variance, increasing leaf area density increases ${\sigma_p/u_{*}^{2}}$ near the canopy top, where u _* is the longitudinally averaged friction velocity at the canopy top and σ _p is the standard deviation of the pressure fluctuations. This increase is shown to be consistent with a primitive scaling argument on the leading term describing the mean-flow turbulent interaction. This scaling argument also predicts the spatial variations in σ _p above the canopy reasonably well for both simulations, but not inside the canopy.     

Exploring the Effects of Microscale Structural Heterogeneity of Forest Canopies Using Large-Eddy Simulations

The Regional Atmospheric Modeling System (RAMS)-based Forest Large-Eddy Simulation (RAFLES), developed and evaluated here, is used to explore the effects of three-dimensional canopy heterogeneity, at the individual tree scale, on the statistical properties of turbulence most pertinent to mass and momentum transfer. In RAFLES, the canopy interacts with air by exerting a drag force, by restricting the open volume and apertures available for flow (i.e. finite porosity), and by acting as a heterogeneous source of heat and moisture. The first and second statistical moments of the velocity and flux profiles computed by RAFLES are compared with turbulent velocity and scalar flux measurements collected during spring and winter days. The observations were made at a meteorological tower situated within a southern hardwood canopy at the Duke Forest site, near Durham, North Carolina, U.S.A. Each of the days analyzed is characterized by distinct regimes of atmospheric stability and canopy foliage distribution conditions. RAFLES results agreed with the 30-min averaged flow statistics profiles measured at this single tower. Following this intercomparison, two case studies are numerically considered representing end-members of foliage and midday atmospheric stability conditions: one representing the winter season with strong winds above a sparse canopy and a slightly unstable boundary layer; the other representing the spring season with a dense canopy, calm conditions, and a strongly convective boundary layer. In each case, results from the control canopy, simulating the observed heterogeneous canopy structure at the Duke Forest hardwood stand, are compared with a test case that also includes heterogeneity commensurate in scale to tree-fall gaps. The effects of such tree-scale canopy heterogeneity on the flow are explored at three levels pertinent to biosphere-atmosphere exchange. The first level (zero-dimensional) considers the effects of such heterogeneity on the common representation of the canopy via length scales such as the zero-plane displacement, the aerodynamic roughness length, the surface-layer depth, and the eddy-penetration depth. The second level (one-dimensional) considers the normalized horizontally-averaged profiles of the first and second moments of the flow to assess how tree-scale heterogeneities disturb the entire planar-averaged profiles from their canonical (and well-studied planar-homogeneous) values inside the canopy and in the surface layer. The third level (three-dimensional) considers the effects of such tree-scale heterogeneities on the spatial variability of the ejection-sweep cycle and its propagation to momentum and mass fluxes. From these comparisons, it is shown that such microscale heterogeneity leads to increased spatial correlations between attributes of the ejection-sweep cycle and measures of canopy heterogeneity, resulting in correlated spatial heterogeneity in fluxes. This heterogeneity persisted up to four times the mean height of the canopy (h _c ) for some variables. Interestingly, this estimate is in agreement with the working definition of the thickness of the canopy roughness sublayer (2h _c –5h _c ).     

A Sensitivity Analysis of the Nocturnal Boundary-Layer Properties to Atmospheric Emissivity Formulations

A one-dimensional model for the mean potential temperature within the nocturnal boundary layer (NBL) was used to assess the sensitivity of three NBL properties (height, thermal stratification strength, and near-surface cooling) to three widely used atmospheric emissivity formulations. The calculations revealed that the NBL height is robust to the choice of the emissivity function, though this is not the case for NBL Richardson number and near-surface cooling rate. Rather than endorse one formulation, our analysis highlights the importance of atmospheric emissivity in modelling the radiative properties of the NBL especially for clear-sky conditions.     

The quest for the solar g modes

Solar gravity modes (or g modes)—oscillations of the solar interior on which buoyancy acts as the restoring force—have the potential to provide unprecedented inference on the structure and dynamics of the solar core, inference that is not possible with the well-observed acoustic modes (or p modes). The relative high amplitude of the g-mode eigenfunctions in the core and the evanesence of the modes in the convection zone make the modes particularly sensitive to the physical and dynamical conditions in the core. Owing to the existence of the convection zone, the g modes have very low amplitudes at photospheric levels, which makes the modes extremely hard to detect. In this article, we review the current state of play regarding attempts to detect g modes. We review the theory of g modes, including theoretical estimation of the g-mode frequencies, amplitudes and damping rates. Then we go on to discuss the techniques that have been used to try to detect g modes. We review results in the literature, and finish by looking to the future, and the potential advances that can be made—from both data and data-analysis perspectives—to give unambiguous detections of individual g modes. The review ends by concluding that, at the time of writing, there is indeed a consensus amongst the authors that there is currently no undisputed detection of solar g modes.     

Potential for deep geological sequestration of CO2 in Switzerland: a first appraisal

Possibilities to sequester anthropogenic CO₂ in deep geological formations are being investigated worldwide, but the potential within Switzerland has not yet been evaluated. This study presents a first-order appraisal based solely on geological criteria collated from the literature. The Swiss Molasse Basin (SMB) and the adjacent Folded Jura are the only realms of the country where CO₂ could conceivably be stored in saline aquifers. Evaluation of geological criteria at the basin-wide scale shows that the SMB–Jura has moderate potential (score of 0.6 on a scale from 0 to 1) when compared to basins elsewhere. At the intrabasinal scale, inspection of the stratigraphy reveals four regional candidate aquifers that are sealed by suitable caprocks: top Basement plus basal Mesozoic sandstones, all sealed by the Anhydrite Group; Upper Muschelkalk sealed by the Gipskeuper; Hauptrogenstein sealed by the Effinger Member, and Upper Malm plus Lower Cretaceous sealed by the Lower Freshwater Molasse. Nine geological criteria are defined to evaluate the storage potential of these and other smaller scale candidates. A numerical scoring and weighting scheme allows the criteria to be assessed simultaneously, permitting the storage potential to be depicted using the 0–1 scale in contoured maps. Approximately 5,000 km² of the central SMB exhibits potentials between 0.6 and 0.96. The Fribourg–Olten–Luzern area is the most favoured owing to the presence of several sealed aquifers within the preferred 800–2,500 m depth interval, and to its low seismicity, low geothermal gradient, low fault density, and long groundwater residence times. Smaller areas with good potential lie between Zürich and St. Gallen. In contrast, western Switzerland, the Jura and the southern SMB have markedly poorer potential. Considering only the portions of the aquifers with potential above 0.6, the theoretical, effective storage capacity of the basin is estimated to be 2,680 million tonnes of CO₂.     

Multiscale analysis of vegetation surface fluxes: from seconds to years

The variability in land surface heat (H), water vapor (LE), and CO₂ (or net ecosystem exchange, NEE) fluxes was investigated at scales ranging from fractions of seconds to years using eddy-covariance flux measurements above a pine forest. Because these fluxes significantly vary at all these time scales and because large gaps in the record are unavoidable in such experiments, standard Fourier expansion methods for computing the spectral and cospectral statistical properties were not possible. Instead, orthonormal wavelet transformations are proposed and used. The are ideal at resolving process variability with respect to both scale and time and are able to isolate and remove the effects of missing data (or gaps) from spectral and cospectral calculations. Using the spectra, we demonstrated unique aspects in three appropriate ranges of time scales: turbulent time scales (fractions of seconds to minutes), meteorological time scales (hour to weeks), and seasonal to interannual time scales corresponding to climate and vegetation dynamics. We have shown that: (1) existing turbulence theories describe the short time scales well, (2) coupled physiological and transport models (e.g. CANVEG) reproduce the wavelet spectral characteristics of all three land surface fluxes for meteorological time scales, and (3) seasonal dynamics in vegetation physiology and structure inject strong correlations between land surface fluxes and forcing variables at monthly to seasonal time scales. The broad implications of this study center on the possibility of developing low-dimensional models of land surface water, energy, and carbon exchange. If the bulk of the flux variability is dominated by a narrow band or bands of modes, and these modes “resonate” with key state and forcing variables, then low-dimensional models may relate these forcing and state variables to NEE and LE.