Large Eddy Simulations of Tall Vegetated Canopies: Heterogeneous Canopies (Patches)

The simulations associated with this dataset form part of a broader investigation into the effects of heterogeneity in tall vegetated canopies, designed to emulate structural and aerodynamic characteristics of the Amazon rainforest. All cases are forced by a geostrophic wind of 9 m s⁻¹ under neutrally stratified conditions. The vegetated canopy configurations feature canopies with circular patches of vegetation surrounded by bare soil that are randomly distributed over the domain. The numerical value in each patch case name corresponds to the diameter of the patches, D, expressed in meters after multiplication by a factor of 100. Thus: the configurations labeled "Patch_4_9mps" and "Patch_4_h_9mps" refer to gaps with D = 400 m, the configuration labeled "Patch_8_9mps" refers to gaps with D = 800 m, and the configuration labeled "Patch_12_9mps" refers to gaps with D = 1200 m. The suffix "_h" in the simulation label "Patch_4_h_9mps" stands for "half," indicating that the number of vegetation patches in this configuration is half that of the corresponding "Patch_4_9mps" simulation. Across all configurations, the canopy height is fixed at hc = 39 m, and the vegetation is characterized by a vertically distributed leaf area density, with a Leaf Area Index (LAI) of 7.05. Simulations were conducted using a Large-Eddy Simulation (LES) domain of size (Lx, Ly, Lz) = (2π, 2π, 1) km, where x, y, and z denote the streamwise, spanwise, and vertical directions. The computational grid consists of (Nx, Ny, Nz) = (256,256,256) points, yielding a spatial resolution of (Δx, Δy, Δz) = (24.5, 24.5, 3.9) m. An aerodynamic roughness length of z0 = 0.01 m is imposed at the surface. Time integration employs a fixed time step of Δt = 0.05 s. Simulations are spun up for 25 h, and data are sampled over the final 5 h of integration. To facilitate data access, a MATLAB script (main_load_snapshots.m) and a supporting function (load_snaps.m) are provided. These scripts load the binary files and organize them into 4D arrays for subsequent analysis.

See README file for data retrieval instructions.