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

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

Post-Last Glacial Maximum evolution of a “fjord-type” lake based on high-resolution seismic data: the Lago Roca/Acigami (southern Tierra del Fuego,Argentina/Chile)

Lago Roca/Acigami is a “fjord-type” lake located in the southernmost part of South America, in the proximity of the Beagle Channel. A high-resolution seismic survey was carried out to analyse the seismic stratigraphy of the lake and to shed some light on the post-Last Glacial Maximum history of the area. Six seismic units were recognised, and their nature and depositional context were interpreted using seismic stratigraphy and acoustic facies analysis. A buried large ridge was identified within the glacial unit (SU1), interpreted as a frontal moraine that indicates a stabilisation phase. After retreat of the glacier from the basin, the trough was flooded by meltwater and a lake developed (SU2). The seismic facies, from bottom to top, depict a transition from ice-contact (SU2; SU3) to ice-distal proglacial conditions (SU4). A thick draping unit (SU5) marks a marine transgressive event and the instauration of a fjord environment in the basin. The marine transgression was a rapid event preceded by a fall in the lake level that caused an erosional unconformity. During the fjord phase the sedimentation remained controlled by meltwater discharge. Sea level fall, and subsequent disconnection from the Beagle Channel was accompanied by a progradation of the glaciofluvial deltaic sediments and the occurrence of several mass-wasting deposits (SU6). © 2019 John Wiley & Sons, Ltd.     

Pointing of HAGAR telescope mirrors

An array of seven atmospheric Cherenkov telescopes was commissioned at a high altitude site in Hanle in the Ladakh region of the Himalayas. The array called HAGAR has been designed to observe celestial γ-rays of energy >100 GeV. Each telescope is altitude-azimuth mounted and carries seven parabolic mirrors whose optic axes are co-aligned with the telescope axis. The telescopes point and track a celestial source using a PC-based drive control system. Two important issues in positioning of each HAGAR telescope are pointing accuracy of telescope axis and co-alignment of mirrors’ optic axes with the telescope axis. We have adopted a three pronged strategy to address these issues, namely use of pointing models to improve pointing accuracy of the telescopes, RA-DEC scan technique to measure the pointing offsets of the mirrors and mechanical fine-tuning of off-axis mirrors by sighting a distant stationary light source. This paper discusses our efforts in this regard as well as the current status of pointing and monitoring of HAGAR telescopes.     

A possible 5 km wide impact structure with associated 22 km wide exterior collapse terrain in the Alhabia–Tabernas Basin,southeastern Spain

The Tabernas–Alhabia Basin is a structural depression situated in the province of Almería, southeastern Spain. The basin is filled with Neogene, Pliocene, and Pleistocene sediments resting discordantly on a Paleozoic metamorphic basement. During the marine Tortonian sedimentation, a bed of breccia (Gordo megabed) was formed. It consists of rotated sedimentary megablocks commonly capped and/or surrounded by a polymict breccia composed mainly of up to dm-sized clasts of the crystalline (schist) basement. Previous work has suggested the bed to be a seismite corresponding to events induced by earthquakes. Here, we link the formation of the Gordo megabed with an ∼5 km wide, rimmed depression with exposed breccias on the northern flank of the Sierra de Gádor mountain. This semicircular structure, developed in mainly schists and dolostone of the basement, is delimited to the W, S, and E by an up to 350 m high escarpment with overturned stratigraphy. Toward the north, this crater-like structure opens toward the Gordo megabed of the Tabernas Basin. In the southern sector, the overturned strata transform outward for into a blocky allochthonous breccia with decreasing thickness and clast size. In the interior of the structure, there are occurrences of graded breccia and arenite superposed on a blocky, autochthonous breccia. Based on the presence of mineralogical shock metamorphic evidence, potential shatter cones, and a high Ir anomaly (∼500 ppb) as well as the position of the structure near the town of Alhama de Almería, we propose to call it the Alhama de Almería impact structure.     

Parameter Estimation for Groundwater Models under Uncertain Irrigation Data

The success of modeling groundwater is strongly influenced by the accuracy of the model parameters that are used to characterize the subsurface system. However, the presence of uncertainty and possibly bias in groundwater model source/sink terms may lead to biased estimates of model parameters and model predictions when the standard regression‐based inverse modeling techniques are used. This study first quantifies the levels of bias in groundwater model parameters and predictions due to the presence of errors in irrigation data. Then, a new inverse modeling technique called input uncertainty weighted least‐squares (IUWLS) is presented for unbiased estimation of the parameters when pumping and other source/sink data are uncertain. The approach uses the concept of generalized least‐squares method with the weight of the objective function depending on the level of pumping uncertainty and iteratively adjusted during the parameter optimization process. We have conducted both analytical and numerical experiments, using irrigation pumping data from the Republican River Basin in Nebraska, to evaluate the performance of ordinary least‐squares (OLS) and IUWLS calibration methods under different levels of uncertainty of irrigation data and calibration conditions. The result from the OLS method shows the presence of statistically significant (p < 0.05) bias in estimated parameters and model predictions that persist despite calibrating the models to different calibration data and sample sizes. However, by directly accounting for the irrigation pumping uncertainties during the calibration procedures, the proposed IUWLS is able to minimize the bias effectively without adding significant computational burden to the calibration processes.     

Spectrometric gamma radiation of shale cores applied to sweet spot discrimination in Eastern Pomerania,Poland

This paper describes the application and calculation of hydrocarbon anomalies in two different boreholes located in Eastern Pomerania (northern Poland). Spectrometric data from borehole geophysical probe (borehole 1) and portable gamma logger (borehole 2) were used to analyze shale formations. The results from borehole 1 presented a statistically significant, moderate correlation between calculated hydrocarbon anomalies and hydrocarbon saturation data obtained from well log interpretation. Borehole 2 has been analyzed focusing on the gamma radiation of the core samples, and the positive results of borehole 1. Hydrocarbon anomalies calculated from spectral gamma radiation are reliable indicators of sweet spots, based solely on a cursory evaluation of core measurements. These preliminary information acquired from gamma-ray measurements could help increase sampling precision of further geochemical analysis.     

Chemical characterization of porewaters in an intertidal mudflat of the Seine estuary: relationship to erosion-deposition cycles

A seasonal field study was carried out in the Seine estuary to determine the chemistry of sediment porewaters using the 'peeper' technique and changes in the elevation of the mudflats using the 'Altus' technique. This approach allowed us to evaluate the release of nutrients and to link these releases to the sediment hydrodynamics. Our results show that nutrient and organic matter cycling in a Seine estuary mudflat exhibits a seasonal behaviour, which is mainly influenced by variations in hydrodynamics. Sediments, rich organic matter, were input during floods and they were mineralized during summer and autumn, releasing nutrients and dissolved organic carbon into the sediment porewaters. The nutrient release, including ammonium, is mainly linked to the mineralization of organic matter, while the release of phosphate is delayed. The delay could be the result of phosphate association with organic matter and/or its co-precipitation with calcium and iron.     

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.