Abstract. Airborne LiDAR is increasingly used in forest carbon, ecosystem, and resource monitoring. For practical design and manufacture reasons, the 1064 nm near-infrared (NIR) wavelength has been the most commonly adopted, and most literature in this field represents sampling characteristics in this wavelength. However, due to eye-safety and application-specific needs, other common wavelengths are 1550 nm and 532 nm. All provide canopy structure reconstructions that can be integrated or compared through space and time but the consistency or complementarity of 3D airborne LiDAR data sampled at multiple wavelengths is poorly understood. Here, we report on multispectral LiDAR missions carried out in 2013 and 2015 over a managed forest research site. The 1st used 3 independent sensors, and the 2nd used a single sensor carrying 3 lasers. The experiment revealed differences in proportions of returns at ground level, vertical foliage distributions, and gap probability across wavelengths. Canopy attenuation was greatest at 532 nm, presumably due to leaf tissue absorption. Relative to 1064 nm, foliage was undersampled at midheight percentiles at 1550 nm and 532 nm. Multisensor data demonstrated differences in foliage characterization due to combined influences of wavelength and acquisition configuration. Single-sensor multispectral data were more stable but demonstrated clear wavelength-dependent variation that could be exploited in intensity-based land cover classification without the aid of 3D derivatives. This work sets the stage for improvements in land surface classification and vertical foliage partitioning through the integration of active spectral and structural laser return information.