Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies

In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (F_T) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating F_T to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (S_T) and F_T have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if S_T and F_T are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.     

Flume Experiments on Turbulent Flows Across Gaps of Permeable and Impermeable Boundaries

Laser Doppler anemometery and laser-induced fluorescence techniques were used to explore the spatial structure of the flow within and above finite cavities created within porous and solid media. The cavities within these two configurations were identical in size and were intended to mimic flow disturbances created by finite gaps and forest clearing. Because flows over permeable boundaries differ from their solid counterparts, the study here addresses how these differences in boundary conditions produce differences in, (i) bulk flow properties including the mean vorticity within and adjacent to the gaps, (ii) second-order statistics such as the standard deviations and turbulent stresses, (iii) the relative importance of advective to turbulent stress terms across various regions within and above the gaps, and (iv) the local imbalance between ejections and sweeps and momentum transport efficiencies of updrafts and downdrafts. Both configurations exhibited a primary recirculation zone of comparable dimensions inside the gap. The mean vorticity spawned at the upstream corner of the gap was more intense for the solid configuration when compared to its porous counterpart. The free-shear layer spawned from the upstream corner-edge deeper into the gap for the porous configuration. The momentum flux at the interface within and above the gap was enhanced by a factor of 1.5–2.0 over its upstream value, and this enhancement zone was much broader in size for the porous configuration. For the turbulent transport terms in the longitudinal and vertical mean momentum balances, these transport terms were significant inside the gap for both boundary configurations when compared to their upstream counterpart. The effectiveness of using incomplete cumulant expansion methods to describe the momentum transport efficiencies, and the relative contributions of ejections and sweeps to turbulent stresses, especially in this zone, were also demonstrated. The flatness factor for both velocity components, often used as a measure of intermittency, was highest in the vicinity of the upstream corner in both configurations. However, immediately following the downstream corner, the flatness factor remained large for the porous configuration, in contrast to its solid configuration counterpart.     

Turbulent Intensities and Velocity Spectra for Bare and Forested Gentle Hills: Flume Experiments

To investigate how velocity variances and spectra are modified by the simultaneous action of topography and canopy, two flume experiments were carried out on a train of gentle cosine hills differing in surface cover. The first experiment was conducted above a bare surface while the second experiment was conducted within and above a densely arrayed rod canopy. The velocity variances and spectra from these two experiments were compared in the middle, inner, and near-surface layers. In the middle layer, and for the canopy surface, longitudinal and vertical velocity variances () were in phase with the hill-induced spatial mean velocity perturbation (Δu) around the so-called background state (taken here as the longitudinal mean at a given height) as predicted by rapid distortion theory (RDT). However, for the bare surface case, and remained out of phase with Δu by about L/2, where L is the hill half-length. In the canopy layer, wake production was a significant source of turbulent energy for , and its action was to re-align velocity variances with Δu in those layers, a mechanism completely absent for the bare surface case. Such a lower ‘boundary condition’ resulted in longitudinal variations of to be nearly in phase with Δu above the canopy surface. In the inner and middle layers, the spectral distortions by the hill remained significant for the background state of the bare surface case but not for the canopy surface case. In particular, in the inner and middle layers of the bare surface case, the effective exponents derived from the locally measured power spectra diverged from their expected  − 5/3 value for inertial subrange scales. These departures spatially correlated with the hill surface. However, for the canopy surface case, the spectral exponents were near  − 5/3 above the canopy though the minor differences from  − 5/3 were also correlated with the hill surface. Inside the canopy, wake production and energy short-circuiting resulted in significant departures from  − 5/3. These departures from  − 5/3 also appeared correlated with the hill surface through the wake production contribution and its alignment with Δu. Moreover, scales commensurate with Von Karman street vorticies well described wake production scales inside the canopy, confirming the important role of the mean flow in producing wakes. The spectra inside the canopy on the lee side of the hill, where a negative mean flow delineated a recirculation zone, suggested that the wake production scales there were ‘broader’ when compared to their counterpart outside the recirculation zone. Inside the recirculation zone, there was significantly more energy at higher frequencies when compared to regions outside the recirculation zone.     

Correction to: Cultural heritage and earthquakes: bridging the gap between geophysics,archaeoseismology and engineering

Journal of Seismology -     

Time–Frequency–Wavenumber Analysis of Surface Waves Using the Continuous Wavelet Transform

A modified approach to surface wave dispersion analysis using active sources is proposed. The method is based on continuous recordings, and uses the continuous wavelet transform to analyze the phase velocity dispersion of surface waves. This gives the possibility to accurately localize the phase information in time, and to isolate the most significant contribution of the surface waves. To extract the dispersion information, then, a hybrid technique is applied to the narrowband filtered seismic recordings. The technique combines the flexibility of the slant stack method in identifying waves that propagate in space and time, with the resolution of f–k approaches. This is particularly beneficial for higher mode identification in cases of high noise levels. To process the continuous wavelet transform, a new mother wavelet is presented and compared to the classical and widely used Morlet type. The proposed wavelet is obtained from a raised-cosine envelope function (Hanning type). The proposed approach is particularly suitable when using continuous recordings (e.g., from seismological-like equipment) since it does not require any hardware-based source triggering. This can be subsequently done with the proposed method. Estimation of the surface wave phase delay is performed in the frequency domain by means of a covariance matrix averaging procedure over successive wave field excitations. Thus, no record stacking is necessary in the time domain and a large number of consecutive shots can be used. This leads to a certain simplification of the field procedures. To demonstrate the effectiveness of the method, we tested it on synthetics as well on real field data. For the real case we also combine dispersion curves from ambient vibrations and active measurements.     

Experimental assessment of sliding materials for seismic isolation systems

Novel self-lubricating materials are currently being proposed to be used in sliding isolation systems such as the curved surface sliding isolator system, or pendulum bearing system, for the protection of buildings and structures. The current codes for anti-seismic devices are focused on the evaluation of the performance of the whole isolation system; as a consequence, a reliable procedure for the pre-assessment of the material combinations of sliding interfaces is currently missing. Therefore in this paper, an experimental methodology is proposed for the characterization of self-lubricating materials through tests on small-scale specimens performed using customized equipment able to reproduce the operational conditions of real isolation systems as per contact pressure, sliding velocity, temperature and slide path. The testing sequence has been designed in order to evaluate the sliding properties of the material in terms of static and dynamic coefficient of friction and wear resistance. Examples are reported for the assessment of two self-lubricating materials with different sliding characteristics. In order to validate the method and to confirm the reliability of extrapolating the results to real working conditions, prototypes of pendulum isolation systems incorporating the assessed materials have been tested according to the AASHTO (Guide specifications for seismic isolation design, 2nd edn. American Association of State Highway and Transportation Officials, Washington, DC, 1999) specifications and the relevant dynamic properties assessed from the Horizontal Load—Displacement loops. The experimental outcomes confirmed that the frictional characteristics provided by the proposed procedure can be reliably used in the design of seismically isolated structures.     

Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project

Surface wave methods gained in the past decades a primary role in many seismic projects. Specifically, they are often used to retrieve a 1D shear wave velocity model or to estimate the V_S,30 at a site. The complexity of the interpretation process and the variety of possible approaches to surface wave analysis make it very hard to set a fixed standard to assure quality and reliability of the results. The present guidelines provide practical information on the acquisition and analysis of surface wave data by giving some basic principles and specific suggestions related to the most common situations. They are primarily targeted to non-expert users approaching surface wave testing, but can be useful to specialists in the field as a general reference. The guidelines are based on the experience gained within the InterPACIFIC project and on the expertise of the participants in acquisition and analysis of surface wave data.     

Technical guidelines for the assessment of earthquake induced liquefaction hazard at urban scale

Bulletin of Earthquake Engineering - Microzonation for earthquake-induced liquefaction hazard is the subdivision of a territory at a municipal or submunicipal scale in areas characterized by the...     

Two-Dimensional Scalar Spectra in the Deeper Layers of a Dense and Uniform Model Canopy

The turbulent flow inside dense canopies is characterized by wake production and short-circuiting of the energy cascade. How these processes affect passive scalar concentration variability in general and their spectral properties in particular remains a vexing problem. Progress on this problem is frustrated by the shortage of high resolution spatial concentration measurements, and by the lack of simplified analytical models that connect spectral modulations in the turbulent kinetic energy (TKE) cascade to scalar spectra. Here, we report the first planar two-dimensional scalar concentration spectra (ϕ _cc ) inside tall canopies derived from flow visualization experiments. These experiments were conducted within the deeper layers of a model canopy composed of densely arrayed cylinders welded to the bottom of a large recirculating water channel. We found that in the spectral region experiencing wake production, the ϕ _cc exhibits directional scaling power laws. In the longitudinal direction (x), or the direction experiencing the largest drag force, the ϕ _cc (k _x ) was steeper than and followed an approximate at wavenumbers larger than the injection scale of wake energy, where k _x is the longitudinal wavenumber. In the lateral direction (y), the spectra scaled as up to the injection scale, and then decayed at an approximate power law. This departure from the classical inertial subrange scaling (i.e., k ^−5/3) was reproduced using a newly proposed analytical solution to a simplified scalar spectral budget equation. Near the velocity viscous dissipation range, the scalar spectra appear to approach an approximate k ⁻³, a tantalizing result consistent with dimensional analysis used in the inertial-diffusive range. Implications to subgrid modelling for large-eddy simulations (LES) inside canopies are briefly discussed.     

Momentum Transfer and Turbulent Kinetic Energy Budgets within a Dense Model Canopy

Second-order closure models for the canopy sublayer (CSL) employ aset of closure schemes developed for `free-air' flow equations andthen add extra terms to account for canopy related processes. Muchof the current research thrust in CSL closure has focused on thesecanopy modifications. Instead of offering new closure formulationshere, we propose a new mixing length model that accounts for basicenergetic modes within the CSL. Detailed flume experiments withcylindrical rods in dense arrays to represent a rigid canopy areconducted to test the closure model. We show that when this lengthscale model is combined with standard second-order closureschemes, first and second moments, triple velocity correlations,the mean turbulent kinetic energy dissipation rate, and the wakeproduction are all well reproduced within the CSL provided thedrag coefficient (C_D) is well parameterized. The maintheoretical novelty here is the analytical linkage betweengradient-diffusion closure schemes for the triple velocitycorrelation and non-local momentum transfer via cumulant expansionmethods. We showed that second-order closure models reproducereasonably well the relative importance of ejections and sweeps onmomentum transfer despite their local closure approximations.Hence, it is demonstrated that for simple canopy morphology (e.g.,cylindrical rods) with well-defined length scales, standard closureschemes can reproduce key flow statistics without much revision.When all these results are taken together, it appears that thepredictive skills of second-order closure models are not limitedby closure formulations; rather, they are limited by our abilityto independently connect the drag coefficient and the effectivemixing length to the canopy roughness density. With rapidadvancements in laser altimetry, the canopy roughness densitydistribution will become available for many terrestrialecosystems. Quantifying the sheltering effect, the homogeneity andisotropy of the drag coefficient, and more importantly, thecanonical mixing length, for such variable roughness density isstill lacking.