Starting with the description of meteorological variables in forest canopies and its parameter variations, a numerical three-dimentional model is developed. Its applicability is demonstrated, first, by wind sheltering effects of hedges and, second, by the effects of deforestation on local climate in complex terrain. Scientists in ecology, agricultural botany and meteorology, but also urban and regional lanners will profit from this study finding the most effective solution for their specific problems.

Provides an essential introduction to modeling terrestrial ecosystems in Earth system models for graduate students and researchers.

The turbulent flow within and above plant canopies is responsible for the exchange of momentum, heat, gases and particles between vegetation and the atmosphere. Turbulence is also responsible for the mixing of air inside the canopy, playing an important role in chemical and biophysical processes occurring in the plants environment. In the last fifty years, research has significantly advanced the un- derstanding of and ability to model the flow field within and above the canopy, but important issues remain unsolved. In this work, we focus on (i) the estimation of turbulent mixing timescales within the canopy from field data; and (ii) the development of new computationally efficient modeling approaches for the coupled canopy-atmosphere flow field.The turbulent mixing timescale represents how quickly turbulence creates a well- mixed environment within the canopy. When the mixing timescale is much smaller than the timescale of other relevant processes (e.g. chemical reactions, deposition), the system can be assumed to be well-mixed and detailed modeling of turbulence is not critical to predict the system evolution. Conversely, if the mixing timescale is comparable or larger than the other timescales, turbulence becomes a controlling factor for the concentration of the variables involved; hence, turbulence needs to be taken into account when studying and modeling such processes. In this work, we used a combination of ozone concentration and high-frequency velocity data measured within and above the canopy in the Amazon rainforest to characterize turbulent mixing. The eddy diffusivity parameter (used as a proxy for mixing efficiency) was applied in a simple theoretical model of one-dimensional diffusion, providing an estimate of turbulent mixing timescales as a function of height within the canopy and time-of-day. Results showed that, during the day, the Amazon rainforest is characterized by well-mixed conditions with mixing timescales smaller than thirty minutes in the upper-half of the canopy, and partially mixed conditions in the lower half of the canopy. During the night, most of the canopy (except for the upper 20%) is either partially or poorly mixed, resulting in mixing timescales of up to several hours. For the specific case of ozone, the mixing timescales observed during the day are much lower than the chemical and deposition timescales, whereas chemical processes and turbulence have comparable timescales during the night. In addition, the high day-to-day variability in mixing conditions and the fast increase in mixing during the morning transition period indicate that turbulence within the canopy needs to be properly investigated and modeled in many studies involving plant-atmosphere interactions.Motivated by the findings described above, this work proposes and tests a new approach for modeling canopy flows. Typically, vertical profiles of flow statistics are needed to represent canopy-atmosphere exchanges in chemical and biophysical processes happening within the canopy. Current single-column models provide only steady-state (equilibrium) profiles, and rely on closure assumptions that do not represent the dominant non-local turbulent fluxes present in canopy flows. We overcome these issues by adapting the one-dimensional turbulent (ODT) model to represent atmospheric flows from the ground up to the top of the atmospheric boundary layer (ABL). The ODT model numerically resolves the one-dimensional diffusion equation along a vertical line (representing a horizontally homogeneous ABL column), and the presence of three-dimensional turbulence is added through the effect of stochastic eddies. Simulations of ABL without canopy were performed for different atmospheric stabilities and a diurnal cycle, to test the capabilities of this modeling approach in representing unsteady flows with strong non-local transport. In addition, four different types of canopies were simulated, one of them including the transport of scalar with a point source located inside the canopy. The comparison of all simulations with theory and field data provided satisfactory results. The main advantages of using ODT compared to typical 1D canopy-flow models are the ability to represent the coupled canopy-ABL flow with one single modeling approach, the presence of non-local turbulent fluxes, the ability to simulate transient conditions, the straightforward representation of multiple scalar fields, and the presence of only one adjustable parameter (as opposed to the several adjustable constants and boundary conditions needed for other modeling approaches).The results obtained with ODT as a stand-alone model motivated its use as a surface parameterization for Large-Eddy Simulation (LES). In this two-way coupling between LES and ODT, the former is used to simulate the ABL in a case where a canopy is present but cannot be resolved by the LES (i.e., the LES first vertical grid point is above the canopy). ODT is used to represent the flow field between the ground and the first LES grid point, including the region within and just above the canopy. In this work, we tested the ODT-LES model for three different types of canopies and obtained promising results. Although more work is needed in order to improve first and second-order statistics within the canopy (i.e. in the ODT domain), the results obtained for the flow statistics in the LES domain and for the third order statistics in the ODT domain demonstrate that the ODT-LES model is capable of capturing some important features of the canopy-atmosphere interaction. This new surface superparameterization approach using ODT provides a new alternative for simulations that require complex interactions between the flow field and near-surface processes (e.g. sand and snow drift, waves over water surfaces) and can potentially be extended to other large-scale models, such as mesoscale and global circulation models.

Forest thinning has been going on for decades for reasons such as removing dead trees, controlling potential wildfires and controlling invasive bark beetle infestations. However more needs to be done in the way of examining the effects of the thinning on canopy flow. This research presents the results of a large-eddy simulation (LES) analysis of flow within and above a loblolly pine forest. The simulations represent an actual forest where the US Forest Service engaged in a project to measure canopy flow properties while progressively thinning the forest vegetation. The LES model is compared to an eddy-viscosity model as well as data from the forest thinning. It was found that the number of large coherent structures decreased as the canopy density decreased. The spacing between the structures increased with decreasing canopy density. As the forest was thinned the velocity in the lower canopy increased but the average velocity profile continued to show an inflected shape. The height of the inflection remaining unchanged with the change in canopy densities examined. Zero displacement height tracked linearly with the percentage of density reduction and it also proved to be a good reference point for validating reestablishment of transitioning flow. Studies of canopy edge flow where the flow transitions from dense forest to open ground indicate a distance of approximately 15 canopy heights before the flow completes its transition and behaves as if there had been no canopy present. When the flow transitions to a thinned canopy instead of bare ground the presence of the canopy has a measureable effect on the distance required to reestablish steady state flow. The denser the canopy the more quickly the flow is reestablished. The thinner the canopy, the higher the downward sweeps of velocity occur to reestablish the flow. Research into flow transitioning from canopy to bare ground found a recirculating zone at the transition edge for some canopy densities. No recirculation zone was found at the transition edge from the fully dense canopy to the three different thinnings examined here.

This book is an introduction to the theory, practice, and implementation of the Lattice Boltzmann (LB) method, a powerful computational fluid dynamics method that is steadily gaining attention due to its simplicity, scalability, extensibility, and simple handling of complex geometries. The book contains chapters on the method's background, fundamental theory, advanced extensions, and implementation. To aid beginners, the most essential paragraphs in each chapter are highlighted, and the introductory chapters on various LB topics are front-loaded with special "in a nutshell" sections that condense the chapter's most important practical results. Together, these sections can be used to quickly get up and running with the method. Exercises are integrated throughout the text, and frequently asked questions about the method are dealt with in a special section at the beginning. In the book itself and through its web page, readers can find example codes showing how the LB method can be implemented efficiently on a variety of hardware platforms, including multi-core processors, clusters, and graphics processing units. Students and scientists learning and using the LB method will appreciate the wealth of clearly presented and structured information in this volume.