Reliable non-destructive testing (NDT) ultrasound systems for timber composite structures require quantitative understanding of the propagation of ultrasound beams in wood. A finite-difference time-domain (FDTD) model is described, which incorporates local anisotropy variations of stiffness, damping and density in timber elements. The propagation of pulsed air-coupled ultrasound (ACU) beams in normal and slanted incidence configurations is reproduced by direct definition of material properties (gas, solid) at each model pixel. First, the model was quantitatively validated against analytical derivations. Time-varying wavefronts in unbounded timber with curved growth rings were accurately reproduced, as well as the acoustic properties (velocity, attenuation, beam skewing) of ACU beams transmitted through timber lamellas. An experimental sound field imaging (SFI) setup was implemented at NDT frequencies (120kHz), which for specific beam incidence positions allows spatially resolved ACU field characterization at the receiver side. The good agreement of experimental and modeled beam shifts across timber laminates allowed extrapolation of the inner propagation paths. The modeling base is an orthotropic stiffness dataset for the desired wood species. In cross-grain planes, beam skewing leads to position-dependent wave paths. They are well-described in terms of the growth ring curvature, which is obtained by visual observation of the laminate. Extraordinary refraction phenomena were observed, which lead to well-collimated quasi-shear wave coupling at grazing beam incidence angles. The anisotropic damping in cross-grain planes is satisfactorily explained in terms of the known anisotropic stiffness dataset and a constant loss tangent. The incorporation of high-resolution density maps (X-ray computed tomography) provided insight into ultrasound scattering effects in the layered growth ring structure. Finally, the combined potential of the FDTD model and the SFI setup for material property and defect inversion in anisotropic materials was demonstrated. A portable SFI demonstrator was implemented with a multi-sensor MEMs receiver array that captures and compensates for variable wave propagation paths in glued laminated timber, and improves the imaging of lamination defects.