<p>Biogenic volatile organic compounds (BVOCs) significantly contribute to forest therapy benefits, yet their spatiotemporal emission patterns and environmental drivers in forest air remain inadequately understood. This study aimed to characterize the spatiotemporal dynamics of BVOCs and their environmental drivers in <i>Larix gmelinii</i> var. <i>principis-rupprechtii</i> plantations. Using open-path sampling with GC-MS, we investigated diurnal BVOCs emission profiles from 6:00 to 20:00 in mid-July, mid-August, and mid-September. Vertical and horizontal distribution patterns were characterized at three canopy strata (upper, middle, lower) and microsites (interior, edge, gap). Simultaneous measurements of illuminance, temperature, humidity, wind speed, particulate matter &lt; 10&#xa0;μm (PM₁₀), and negative air ion (NAI) were conducted during diurnal sampling to identify environmental correlates. The study found that total BVOCs emission peaked in the upper and middle canopies—the former enriched in decanal and o-cymene, and the latter in nonanal and monoterpenes—while emission was lowest from the chlorotic leaves in the lower. Horizontally, BVOCs concentrations were highest in the forest interior and decreased outward, being lowest in forest gap. Diurnally, BVOCs decrease during midday photosynthetic depression; at dusk <i>cis</i>-3-hexen-1-ol and 2-methylpentanal emerged, and β-pinene were substantially released; nonanal and decanal show U-shaped patterns. Generalized Additive Models (GAM) identified that nonanal variation was primarily driven by illuminance (<i>F</i>=4.46, <i>p </i>= 0.07), followed by PM₁₀ (<i>F </i>= 3.66, <i>p </i>= 0.06), temperature (<i>F </i>= 2.92,<i> p </i>&lt; 0.05), and humidity (<i>F </i>= 2.46,<i> p </i>&lt; 0.05). Partial Least Squares Path Modeling (PLS-PM) revealed that meteorological factors influenced nonanal through both direct (<i>p </i>&lt; 0.05) and indirect (<i>p </i>&lt; 0.05) pathway. The direct effect was substantial (β = 0.37), while the indirect effect—mediated by PM₁₀ and NAI—was weaker (β = 0.14). These findings established spatiotemporal patterns and key environmental drivers of BVOCs in <i>L. gmelinii</i> var. <i>principis-rupprechtii</i> plantations. They suggest elevating walkways to near mid-canopy height and locating core activity areas inside the forest to optimize exposure to beneficial volatiles. Furthermore, scheduling activities around dawn and dusk during the growing season is recommended. From a management perspective, creating diffuse light conditions and elevating air humidity can foster a more conducive microclimate for therapy.</p>