Laboratory study and predictive modeling for thaw subsidence in deep permafrost

Oil wells on the North Slope of Alaska pass through deep deposits of permafrost. The heat transferred during their operation causes localized thawing, resulting in ground subsidence adjacent to the well casings. This subsidence has a damaging effect, causing the casings to compress, deform, and potentially fail. This paper presents the results of a laboratory study of the thaw consolidation strain of deep permafrost and its predictive modeling. Tests were performed to determine strains due to thaw and post-thaw loading, as well as soil index properties. Results, together with data from an earlier testing program, were used to produce empirical models for predicting strains and ground subsidence. Four distinct strain cases were analyzed with three models by multiple regression analyses, and the best-fitting model was selected for each case. Models were further compared in a ground subsidence prediction using a shared subsurface profile. Laboratory results indicate that strains due to thaw and post-thaw testing in deep core permafrost are insensitive to depth and are more strongly influenced by stress redistributions and the presence of ice lenses and inclusions. Modeling results show that the most statistically valid and useful models were those constructed using moisture content, porosity, and degree of saturation. The applicability of these models was validated by comparison with results from Finite Element modeling.     

Urbanization of the United States over two centuries: an approach based on a long-term database (1790–2010)

This paper presents an original methodology for the construction of a harmonized database for urban areas in the United States from 1790 to 2010 (a period with a census every 10 years, amounting to 23 dates). Upstream, this method is based on a theoretical choice, the use of a spatio-temporal reference frame to construct urban areas over the long term, defined by a maximum range of one hour’s travel for each period. We used an automatic aggregation from a Reilly model. Urban measures (growth, concentration) derived from this harmonized database are then presented, providing a unified picture of urbanization in the United States over more than two centuries. They are compared with results derived from more classical databases established using the official figures from the Census Bureau, or from other constructions differing from the present system. The results diverge considerably in some instances, which confirms the impact of the method of delineation of geographical areas in measures of growth.     

An Innovative Post Grouting Technique for Soil Nails

This study presents an innovative technique of executing soil nails called sectorized post grouting (SPG). The most utilized technique of soil nail grouting is gravity grouting, with the literature reporting advances in pressurized grouting. Although obtaining higher pullout resistance of soil nails, pressurized grouting, mostly done in single-stage grouting, does not compensate for exudation and its use in higher nail lengths is difficult. Thus, a technique has been developed that compensates for exudation, with easier application in lengthier nails. The technique was qualitatively assessed to evaluate its surface roughness and later applied in seven real soil nailing works, where it could be quantitatively assessed. The results show that sectorized post-grouted nails obtained greater pullout resistances than gravity grouted and single-stage grouted nails. Although similar improvement was found in tube-à-manchette (TAM) grouted nails, this method presents lower economic efficiency than sectorized post grouting. The pullout resistance results obtained in this study can be utilized in future soil nailing works executed utilizing SPG.

     

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.     

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 ).     

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₂.     

Fluid-driven low-grade metamorphism in polydeformed rocks of Avalonia (Arisaig Group, Nova Scotia, Canada)

The Lower Silurian??Lower Devonian Arisaig Group (Antigonish Highlands) in the Canadian Appalachians is a sequence of shallow marine strata deposited after the accretion of Avalonia to Baltica during the closure of the Iapetus Ocean. Deformation of the strata is widely attributed to the Devonian Acadian orogeny and produced shallowly plunging regional folds and a cleavage of varying penetrativity. Phyllosilicate minerals from the finest-grained rocks exhibit very low-grade (diagenetic-anchizone) metamorphic conditions. X-ray diffraction study reveals that the sampled rocks contain quartz, K-white mica, chlorite, and feldspars; illite?Csmectite and chlorite?Csmectite mixed-layers are common but Na?CK mica and kaolinite occur only in some samples. The identification of illite?Csmectite mixed-layers in diagenetic samples, with Kübler Index >0.50 ??°2?? and the highly heterogeneous b-cell dimension of the K-white micas are in agreement with the variable chemical composition of dioctahedral micas, which present low illitic substitution and variable phengitic content. The spatial variation in the above crystal-chemical parameters was plotted along a NW?CSE composite cross section across the regional folds. No correlation was found between the metamorphic conditions and either the stratigraphic depth or the strain values measured by phyllosilicates orientation analyses, as a function of the penetrativity of the cleavage. However, the metamorphic grade generally increases towards the Hollow Fault, and is highest in samples located within a 1?km corridor from the fault surface. Incipient cleavage is observed in the anchizonal samples located in the vicinity of the Hollow Fault and in some of the diagenetic samples, indicating cleavage development under low temperatures (<200?oC). These relationships, together with regional syntheses, suggest low-grade metamorphism post-dated regional folding and was coeval with Late Carboniferous dextral movement along the Hollow Fault. Fluid circulation associated with movement along this major fault may be the driving mechanism for the increasing metamorphism towards it.     

Zooplankton and Karenia brevis in the Gulf of Mexico

Blooms of the toxic dinoflagellate Karenia brevis are common in the Gulf of Mexico, yet no in situ studies of zooplankton and K. brevis have been conducted there. Zooplankton abundance and taxonomic composition at non-bloom and K. brevis bloom stations within the Ecology of Harmful Algal Blooms (ECOHAB) study area were compared. At non-bloom stations, the most abundant species of zooplankton were Parvocalanus crassirostris, Oithona colcarva, and Paracalanus quasimodo at the 5-m isobath and P. quasimodo, O. colcarva, and Oikopleura dioica at the 25-m isobath. There was considerable overlap in dominance of zooplankton species between the 5 and 25-m isobaths, with nine species contributing to 90% of abundance at both isobaths. At stations within K. brevis blooms however, Acartia tonsa, Centropages velificatus, Temora turbinata, Evadne tergestina, O. colcarva, O. dioica, and P. crassirostris were dominant. Variations in abundance between non-bloom and bloom assemblages were evident, including the reduction in abundance of three key species within K. brevis blooms.