skip to main content


Title: Increasing Daytime Stability Enhances Downslope Moisture Transport in the Subcanopy of an Even‐Aged Conifer Forest in Western Oregon, USA
Abstract

Mountain breezes, including katabatic and anabatic flows, and temperature inversions are common features of forested mountain landscapes. However, the effects of mountain breezes on moisture transport in forests and implications for regional climate change are not well understood. A detailed, instrumented study was conducted from July to September 2012 in an even‐aged conifer forest in the Oregon Cascade Range to investigate how temperature profiles within the forest canopy influenced atmospheric surface layer processes that ventilate the forest. Subcanopy inversion strength has a bimodal relationship to subcanopy wind speed and moisture flux from the forest. On days with relatively modest heating of the top of the canopy and weak subcanopy inversions, above canopy winds more efficiently mix subcanopy air, leading to greater than average vertical moisture flux and weaker than average along‐slope, subcanopy water vapor advection. On days with strong heating of the top of the canopy and a strong subcanopy inversion, vertical moisture flux is suppressed, and daytime downslope winds are stronger than average under the canopy. Increased downslope winds lead to increased downslope transport of water vapor, carbon dioxide, and other scalars under the canopy. Increasing summer vapor pressure deficit in the Pacific Northwest will enhance both processes: vertical moisture transport by mountain breezes when subcanopy inversions are weak and downslope water vapor transport when subcanopy inversions are strong. These mountain breeze dynamics have implications for climate refugia in forested mountains, forest plantations, and other forested regions with a similar canopy structure and regional atmospheric forcings.

 
more » « less
Award ID(s):
1848019 2025755
NSF-PAR ID:
10445911
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Atmospheres
Volume:
127
Issue:
9
ISSN:
2169-897X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process-oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations.   Details about the simulation setup can be found in section 3 of the associated open-source manuscript, "Enhancing understanding of the hydrological cycle via pairing of process‐oriented and isotope ratio tracers." In brief, we conducted two simulations of the atmosphere from 1980-2004 using the isotope-enabled version of the Community Atmosphere Model 6 (iCAM6) at 0.9x1.25° horizontal resolution, and with 30 vertical hybrid layers spanning from the surface to ~3 hPa. In the first simulation, wind and surface pressure fields were "nudged" toward the ERA5 reanalysis dataset by adding a nudging tendency, preventing the model from diverging from observed/reanalysis wind fields. In the second simulation, no additional nudging tendency was included, and the model was allowed to evolve 'freely' with only boundary conditions provided at the top (e.g., incoming solar radiation) and bottom (e.g., observed sea surface temperatures) of the model. In addition to the isotopic variables, our simulation included a suite of 'process-oriented tracers,' which we describe in section 2 of the manuscript. These variables are meant to track a property of water associated with evaporation, condensation, or atmospheric transport. Metadata are provided about each of the files below; moreover, since the attached files are NetCDF data - this information is also provided with the data files. NetCDF metadata can be accessed using standard tools (e.g., ncdump). Each file has 4 variables: the tagged quantity, and the associated coordinate variables (time, latitude, longitude). The latter three are identical across all files, only the tagged quantity changes. Twelve files are provided for the nudged simulation, and an additional three are provided for the free simulations: Nudged simulation files iCAM6_nudged_1980-2004_mon_RHevap: Mass-weighted mean evaporation source property: RH (%) with respect to surface temperature. iCAM6_nudged_1980-2004_mon_Tevap: Mass-weighted mean evaporation source property: surface temperature in Kelvin iCAM6_nudged_1980-2004_mon_Tcond: Mass-weighted mean condensation property: temperature (K) iCAM6_nudged_1980-2004_mon_columnQ: Total (vertically integrated) precipitable water (kg/m2).  Not a tagged quantity, but necessary to calculate depletion times in section 4.3 (e.g., Fig. 11 and 12). iCAM6_nudged_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_nudged_1980-2004_mon_d18Oevap_0: Mass-weighted mean evaporation source property - d18O of the evaporative flux (e.g., the 'initial' isotope ratio prior to condensation), (‰ VSMOW) iCAM6_nudged_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_nudged_1980-2004_mon_dexevap_0: Mass-weighted mean evaporation source property - deuterium excess of the evaporative flux iCAM6_nudged_1980-2004_mon_lnf: Integrated property - ln(f) calculated from the constant-fractionation d18O tracer (see section 3.2). iCAM6_nudged_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. iCAM6_nudged_1980-2004_mon_residencetime: Mean atmospheric water residence time (in days). iCAM6_nudged_1980-2004_mon_transportdistance: Mean atmospheric water transport distance (in km). Free simulation files iCAM6_free_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_free_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_free_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. 
    more » « less
  2. Abstract

    Cold‐air pooling and associated air temperature inversions are important features of mountain landscapes, but incomplete understanding of their controlling factors hinders prediction of how they may mediate potential future climate changes at local scales. We evaluated how topographic and forest canopy effects on insolation and local winds altered the expression of synoptic‐scale meteorological forcing on near‐surface air temperature inversions and how these effects varied by time of day, season, and spatial scale. Using ~13 years of hourly temperature measurements in forest canopy openings and under the forest canopy at the H.J. Andrews Experimental Forest in the western Cascade Range of Oregon (USA), we calculated air temperature gradients at the basin scale (high vs. low elevation) and at the cross‐valley scale for two transects that differed in topography and forest canopy cover. ERA5 and NCEP NCAR R1 reanalysis data were used to evaluate regional‐scale conditions. Basin and cross‐valley temperature inversions were frequent, particularly in winter and often persisted for several days. Nighttime inversions were more frequent at the cross‐valley scale but displayed the same intra‐annual pattern at the basin and regional scales, becoming most frequent in summer. Nighttime temperature gradients at basin and cross‐valley scales responded similarly to regional‐scale controls, particularly free‐air temperature gradients, despite differences in topography and forest cover. In contrast, the intra‐annual pattern of daytime inversions differed between the basin and cross‐valley scales and between the two cross‐valley transects, implying that topographic and canopy effects on insolation and local winds were key controls at these scales.

     
    more » « less
  3. Abstract

    Forest canopies play a critical role in affecting momentum and scalar transfer. Although there have been recent advances in numerical simulations of turbulent flows and scalar transfer across plant canopies and the atmosphere interface, few models have incorporated all important physical and physiological processes in subcanopy layers. Here we describe and evaluate an advanced multiple‐layer canopy module (MCANOPY), which is developed based largely on the Community Land Model version 4.5 and then coupled with the Weather Research and Forecasting model with large‐eddy simulations (WRF‐LES). The MCANOPY includes a suite of subcanopy processes, including radiation transfer, photosynthesis, canopy layer energy balance, momentum drag, and heat, water vapor, and CO2exchange between canopy layers and the canopy atmosphere. Numerical schemes for heat and water transport in soil, ground surface energy balance, and soil respiration are also included. Both the stand‐alone MCANOPY and the coupled system (the WRF‐LES‐MCANOPY) are evaluated against data measured in the Canopy Horizontal Array Turbulence Study field experiment. The MCANOPY performs reasonably well in reproducing vertical profiles of mean and turbulent flows as well as second‐order statistical quantities including heat and scalar fluxes within the canopy under unstable stability conditions. The coupled WRF‐LES‐MCANOPY captures major features of canopy edge flows under both neutral and unstable conditions. Limitations of the MCANOPY are discussed for our further work. Our results suggest that our model can be a promising modeling system for a variety of applications to study canopy flows and scalar transport (e.g., CO2).

     
    more » « less
  4. Abstract

    This study aimed to understand the microphysical processes that affect rapid intensity changes of tropical cyclones (TCs) over the Bay of Bengal (BoB). Four representative TCs were simulated using the Weather Research and Forecasting model with storm tracking nested configuration (at 9‐km and 3‐km resolution). Results indicate that the inner‐core heating strongly correlated (r > 0.85) with the precipitated compared to non‐precipitated hydrometeors. Furthermore, the vertical distribution of hydrometeors and heating is dependent on inner‐core updrafts and relative humidity. A novel composite analysis of microphysical processes indicates that the warmer (2 K) inner core is close to saturation (>90%) with excess water vapor (>2–3 × 10−3 kg·kg−1), which enhances the latent heat release (LHR) through condensation below the freezing level during the rapid intensification (RI) onset. In addition, during RI, strong updrafts transport the water vapor (>2 × 10−3 kg·kg−1) and cloud liquid water (2.5 × 10−4 kg·kg−1) to above freezing level, and enhance the LHR because of deposition and freezing respectively. The increased precipitating particles in the saturated inner core also enhance LHR. The symmetric convection structured by the atmospheric moisture causes the formation of prolonged RI episodes, as seen in TCPhailin. During rapid weakening (RW), asymmetric and relatively fewer hydrometeors are evident, along with the presence of weak updrafts and strong shear. The dry‐air intrusion into the inner core also causes the cooling processes (evaporation and sublimation). The enhancement or reduction of moist static energy and potential vorticity is associated with increased or reduced LHR in the TC rapid intensity changes.

     
    more » « less
  5. Abstract

    Downslope winds are mesoscale mountain meteorological phenomena that contribute to localized temperature extremes and contribute to numerous societal and environmental impacts. Whereas previous studies have examined local downslope winds, no known efforts have attempted to identify and characterize meso‐ to synoptic‐scale downslope winds globally using a common approach. We use a conceptual model for downslope winds that employs cross‐barrier wind speed, near‐mountain top static stability, and downward vertical velocity using thresholds guided by a chronology of local downslope winds and meta‐analysis of downslope wind case studies. This approach was applied to ERA‐5 reanalysis during 1979–2018 to develop a global atlas of downslope winds. Downslope winds adhered to distinct geographic and seasonal patterns, with peak occurrence in north–south oriented midlatitude mountains in the winter hemisphere associated with strong cross‐mountain winds and stability. However, we identify numerous locations from the tropics to the high‐latitudes where downslope winds occur at least 60 days a year as a byproduct of the general circulation and local‐scale circulation interacting with topography. The four‐decade‐long data set is also used to examine statistical relationships between the occurrence of downslope winds and El Niño‐Southern Oscillation as well as long‐term trends in downslope wind occurrence.

     
    more » « less