Capping Inversion Research Articles

Effects of turbulence on warm clouds and precipitation with various aerosol concentrations

Abstract This study investigates the effects of turbulence-induced collision enhancement (TICE) on warm clouds and precipitation by changing the cloud condensation nuclei (CCN) number concentration using a two-dimensional dynamic model with bin microphysics. TICE is determined according to the Taylor microscale Reynolds number and the turbulent dissipation rate. The thermodynamic sounding used in this study is characterized by a warm and humid atmosphere with a capping inversion layer, which is suitable for simulating warm clouds. For all CCN concentrations, TICE slightly reduces the liquid water path during the early stage of cloud development and accelerates the onset of surface precipitation. However, changes in the rainwater path and in the amount of surface precipitation that are caused by TICE depend on the CCN concentrations. For high CCN concentrations, the mean cloud drop number concentration (CDNC) decreases and the mean effective radius increases due to TICE. These changes cause an increase in the amount of surface precipitation. However, for low CCN concentrations, changes in the mean CDNC and in the mean effective radius induced by TICE are small and the amount of surface precipitation decreases slightly due to TICE. A decrease in condensation due to the accelerated coalescence between droplets explains the surface precipitation decrease. In addition, an increase in the CCN concentration can lead to an increase in the amount of surface precipitation, and the relationship between the CCN concentration and the amount of surface precipitation is affected by TICE. It is shown that these results depend on the atmospheric relative humidity.

A Dryline in Southeast Wyoming. Part II: Airborne In Situ and Raman Lidar Observations

Abstract Part I of this study describes the mesoscale structure of a dryline over southeastern Wyoming. This dryline formed just east of the western rim of the high plains on 22 June 2010 and became more defined as it progressed eastward during the afternoon. Part I also describes the numerically simulated structure and evolution of this dryline and the observed initiation of deep convection in the vicinity of the dryline. An instrumented aircraft, the University of Wyoming King Air, repeatedly flew across this dryline, mostly low enough to penetrate the moist-air wedge east of the dryline. Flight-level in situ data along these low-level penetrations indicate relatively high values of convective available potential energy (CAPE; >1500 J kg−1), yet low convective inhibition, within a few kilometers of the dryline. Water vapor transects obtained from a compact nadir-pointing Raman lidar aboard the aircraft reveal an extremely sharp humidity gradient below flight level along the dryline, coinciding with the fineline seen in operational weather radar base reflectivity imagery. They also reveal several plumes of higher specific humidity within the dry elevated mixed layer above the moist-air wedge, possibly precursors of cumulus clouds. The vertical structure of the dryline revealed by Raman lidar and the flight-level data correspond well to that in the high-resolution numerical simulation.

Aerosols in the convective boundary layer: Shortwave radiation effects on the coupled land‐atmosphere system

AbstractBy combining observations and numerical simulations, we investigated the responses of the surface energy budget and the convective boundary layer (CBL) dynamics to the presence of aerosols. A detailed data set containing (thermo)dynamic observations at CESAR (Cabauw Experimental Site for Atmospheric Research) and aerosol information from the European Integrated Project on Aerosol, Cloud, Climate, and Air Quality Interactions was employed to design numerical experiments reproducing two typical clear‐sky days, each characterized by contrasting thermodynamic initial profiles: (i) residual layer above a strong surface inversion and (ii) well‐mixed CBL connected to the free troposphere by a capping inversion, without the residual layer in between. A large‐eddy simulation (LES) model and a mixed‐layer (MXL) model, coupled to a broadband radiative transfer code and a land surface model, were used to study the impacts of aerosols on shortwave radiation. Both the LES model and the MXL model results reproduced satisfactorily the observations for both days. A sensitivity analysis on a wide range of aerosol properties was conducted. Our results showed that higher loads of aerosols decreased irradiance imposing an energy restriction at the surface, delaying the morning onset of the CBL and advancing its afternoon collapse. Moderately to strongly absorbing aerosols increased the heating rate contributing positively to increase the afternoon CBL height and potential temperature and to decrease Bowen ratio. In contrast, scattering aerosols were associated with smaller heating rates and cooler and shallower CBLs. Our findings advocate the need for accounting for the aerosol influence in analyzing surface and CBL dynamics.

Verification of Convection-Allowing WRF Model Forecasts of the Planetary Boundary Layer Using Sounding Observations

Abstract This study evaluates forecasts of thermodynamic variables from five convection-allowing configurations of the Weather Research and Forecasting Model (WRF) with the Advanced Research core (WRF-ARW). The forecasts vary only in their planetary boundary layer (PBL) scheme, including three “local” schemes [Mellor–Yamada–Janjić (MYJ), quasi-normal scale elimination (QNSE), and Mellor–Yamada–Nakanishi–Niino (MYNN)] and two schemes that include “nonlocal” mixing [the asymmetric cloud model version 2 (ACM2) and the Yonei University (YSU) scheme]. The forecasts are compared to springtime radiosonde observations upstream from deep convection to gain a better understanding of the thermodynamic characteristics of these PBL schemes in this regime. The morning PBLs are all too cool and dry despite having little bias in PBL depth (except for YSU). In the evening, the local schemes produce shallower PBLs that are often too shallow and too moist compared to nonlocal schemes. However, MYNN is nearly unbiased in PBL depth, moisture, and potential temperature, which is comparable to the background North American Mesoscale model (NAM) forecasts. This result gives confidence in the use of the MYNN scheme in convection-allowing configurations of WRF-ARW to alleviate the typical cool, moist bias of the MYJ scheme in convective boundary layers upstream from convection. The morning cool and dry biases lead to an underprediction of mixed-layer CAPE (MLCAPE) and an overprediction of mixed-layer convective inhibition (MLCIN) at that time in all schemes. MLCAPE and MLCIN forecasts improve in the evening, with MYJ, QNSE, and MYNN having small mean errors, but ACM2 and YSU having a somewhat low bias. Strong observed capping inversions tend to be associated with an underprediction of MLCIN in the evening, as the model profiles are too smooth. MLCAPE tends to be overpredicted (underpredicted) by MYJ and QNSE (MYNN, ACM2, and YSU) when the observed MLCAPE is relatively small (large).

Tornadoes in the New York Metropolitan Region: Climatology and Multiscale Analysis of Two Events

Abstract This paper describes the climatology of tornadoes around New York City (NYC) and Long Island (LI), New York, and the structural evolution of two tornadic events that affected NYC on 8 August 2007 and 16 September 2010. Nearly half (18 of 34 events from 1950 to 2010) of NYC–LI tornadoes developed between 0500 and 1300 EDT, and August is the peak tornado month as compared to July for most of the northeast United States. A spatial composite highlights the approaching midlevel trough, moderate most unstable convective available potential energy (MUCAPE), and frontogenesis along a low-level baroclinic zone. Shortly before the early morning tornadoes on 8 August 2007, a mesoscale convective system intensified in the lee of the Appalachians in a region of low-level frontogenesis and moderate MUCAPE (~1500 J kg−1). Warm advection at low levels and evaporative cooling within an elevated mixed layer (EML) ahead of the mesoscale convective system (MCS) helped steepen the low-level lapse rates. Meanwhile, a surface mesolow along a quasi-stationary frontal zone enhanced the warm advection and low-level shear. The late afternoon event on 16 September 2010 was characterized by a quasi-linear convective system (QLCS) that also featured an EML aloft, a surface mesolow just west of NYC, low-level frontogenesis, and a southerly low-level jet ahead of an approaching midlevel trough. The QLCS intensified approaching NYC and generated mesovortices as the QLCS bowed outward. These cases illustrate the benefit of high-density surface observations, terminal Doppler radars, and sounding profiles from commercial aircraft for nowcasting these events.

Vertical Variations of Boundary Layer Potential Buoyancy in Tornadic and Nontornadic Near-Storm Environments

Abstract Despite great strides in understanding the tornadic near-storm environment (NSE), at times it remains difficult to determine why some storms produce significant tornadoes, while others produce none, given similar pretornadic radar reflectivity and velocity signatures. Previous studies have shown that this is likely related to the potential buoyancy (θep) of the rear-flank downdraft (RFD) air. Unfortunately, to date there are few ways to operationally anticipate possible RFD thermodynamic character. Based upon previous research indicating that capping inversions may restrict much of the low-level RFD air to come from within the boundary layer, this study considers the relation of Δθep (vertical change in θep within the boundary layer below the cap) to tornadogenesis potential. This is because when a cap exists above a boundary layer and the descent of lower-θep air from aloft to the surface is potentially limited, then minimal Δθep may indicate more RFD air that has greater potential buoyancy. The Rapid Update Cycle (RUC) soundings used in this study and several observed soundings taken in the vicinity of violent tornadoes suggest that boundary layer Δθep shows promise as an additional means of discriminating between tornadic and nontornadic NSEs.

Influence of Low-Level Thermodynamic Structure on the Downdraft Properties of Simulated Supercells

Abstract Idealized simulations using the Weather Research and Forecasting Model (WRF) were performed to examine the role of capping inversions on the near-surface thermodynamic structure of outflow from simulated supercells. Two simulations were performed: one with the traditional noncapped Weisman and Klemp (WK) analytic sounding and the second with a WK sounding modified to contain a capping inversion. Both sounding environments favor splitting storms and a right-moving supercell by 90 min into the simulation. These two supercell simulations evolve in a qualitatively similar fashion, with both storms exhibiting large, quasi-steady updrafts, hook-shaped appendages in the precipitation mixing ratio field, and prominent localized downdrafts. Results show that the supercell simulated in the capped environment has a surface cold pool with larger values of pseudoequivalent potential temperature (θep) than the cold pool of the supercell produced in the noncapped simulation. Parcels in the surface cold pool of the supercell produced in the capped sounding simulation have a lower origin height than those in the surface cold pool of the supercell produced in the noncapped simulation for all times. Although θep values in the surface cold pool are primarily associated with the origin height of downdraft parcels and the environmental θep at that level, it is shown that nonconservation of θep primarily associated with hydrometeor melting can decrease θep values of downdraft parcels as they descend by several degrees.

Controls on Boundary-Layer Thermodynamics and Dynamics in Coastal West Africa During the Rainy Season of 2006

We investigate dominant processes modulating the coastal West African atmospheric boundary layer during August and September 2006. We evaluated boundary-layer attributes using upper air soundings, tower-based observations, and information from the European Centre for Medium-Range Weather Forecasts reanalyses. Boundary-layer thermodynamics exhibited continental and maritime attributes in response to influences from regional onshore (sea to land) flows and local land–atmosphere exchanges of energy and moisture. Onshore flows transported maritime air inland and gave rise to deep (>1 km) nighttime mixed layers whose heat and moisture content resulted in maximum virtual potential temperatures of 306 K and specific humidities up to 20 g kg−1. The presence of the Saharan Air Layer corresponded with capping inversions greater than 4 K and lapse rates exceeding 7 K km−1 above the mixed layer. Mixed layers at these times became deeper than expected (≈1 km) because dust layer events were often concurrent with strong onshore flows. Despite diurnally variable land–atmosphere fluxes of sensible and latent heat that reached maximum values of 200 and 400 W m−2, respectively, the mixed-layer depth exhibited little diurnal variation due to the influences of onshore flows. Daytime heating of the land, the upward transport of moisture, and onshore flows produced boundary layers with high convective available potential energy that often exceeded 3,000 J kg−1. These results demonstrate that the atmospheric boundary-layer thermodynamics in western Senegal can be favorable for storm development during both day and night. Mesoscale and regional models applied in this region should include several processes controlling the boundary-layer attributes to realistically estimate the energy available for storm development.

Effect of a Capping Inversion on the Stability of an Ekman Boundary Layer

Effect of a Capping Inversion on the Stability of an Ekman Boundary Layer

Idealized Mesoscale Model Simulations of Open Cellular Convection Over the Sea

The atmospheric conditions during an observed case of open cellular convection over the North Sea were simulated using the Weather Research and Forecasting (WRF) numerical model. Wind, temperature and water vapour mixing ratio profiles from the WRF simulation were used to initialize an idealized version of the model, which excluded the effects of topography, surface inhomogeneities and large-scale weather forcing. Cells with an average diameter of 17.4 km developed. Simulations both with and without a capping inversion were made, and the cell-scale kinetic energy budget was calculated for each case. By considering all sources of explicit diffusion in the model, the budgets were balanced. In comparison with previous work based on observational studies, the use of three-dimensional, gridded model data afforded the possibility of calculating all terms in the budgets, which showed that the important terms in the budgets were buoyancy, pressure balance and inter-scale transfer to subgrid scales. Cells were also composited to calculate the average cell-scale flow and each of the budget terms on two-dimensional cross-sections through the cells, parallel and perpendicular to the mean wind direction.

Climatological aspects of seasonal variation of aerosol vertical distribution over central Indo-Gangetic belt (IGB) inferred by the space-borne lidar CALIOP

Abstract Climatological information on the vertical distribution of aerosol in the lower atmosphere is needed to assess the effects of aerosols on climate. Altitude profiles of the backscatter coefficient (at 532 nm), color ratio (CR) and particle depolarization ratio (PDR) derived from space-borne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite were used to investigate the climatological aspects of seasonal variation of aerosol vertical distribution over central Indo-Gangetic belt (IGB). Strong convective activity coupled with capping inversion lead to vertical extension (4–5 km) of aerosol distribution during the pre-monsoon (PrM) season. Trapping of pollution at low altitudes (below 1.5 km) by subsidence have been found during the late post-monsoon (PoM) and winter seasons. Combined results of CR and PDR suggest presence of elevated mineral dust particles during the PrM season, whereas fine mode and spherical (biomass burning/industrial pollution) particles are dominant during the PoM and winter seasons. The observations indicate a seasonal dependency of the aerosol vertical distribution over central IGB.

South East Pacific atmospheric composition and variability sampled along 20° S during VOCALS-REx

Abstract. The VAMOS Ocean-Cloud-Atmosphere-Land Regional Experiment (VOCALS-REx) was conducted from 15 October to 15 November 2008 in the South East Pacific (SEP) region to investigate interactions between land, sea and atmosphere in this unique tropical eastern ocean environment and to improve the skill of global and regional models in representing the region. This study synthesises selected aircraft, ship and surface site observations from VOCALS-REx to statistically summarise and characterise the atmospheric composition and variability of the Marine Boundary Layer (MBL) and Free Troposphere (FT) along the 20° S parallel between 70° W and 85° W. Significant zonal gradients in mean MBL sub-micron aerosol particle size and composition, carbon monoxide, sulphur dioxide and ozone were seen over the campaign, with a generally more variable and polluted coastal environment and a less variable, more pristine remote maritime regime. Gradients in aerosol and trace gas concentrations were observed to be associated with strong gradients in cloud droplet number. The FT was often more polluted in terms of trace gases than the MBL in the mean; however increased variability in the FT composition suggests an episodic nature to elevated concentrations. This is consistent with a complex vertical interleaving of airmasses with diverse sources and hence pollutant concentrations as seen by generalised back trajectory analysis, which suggests contributions from both local and long-range sources. Furthermore, back trajectory analysis demonstrates that the observed zonal gradients both in the boundary layer and the free troposphere are characteristic of marked changes in airmass history with distance offshore – coastal boundary layer airmasses having been in recent contact with the local land surface and remote maritime airmasses having resided over ocean for in excess of ten days. Boundary layer composition to the east of 75° W was observed to be dominated by coastal emissions from sources to the west of the Andes, with evidence for diurnal pumping of the Andean boundary layer above the height of the marine capping inversion. Analysis of intra-campaign variability in atmospheric composition was not found to be significantly correlated with observed low-frequency variability in the large scale flow pattern; campaign-average interquartile ranges of CO, SO2 and O3 concentrations at all longitudes were observed to dominate over much smaller differences in median concentrations calculated between periods of different flow regimes. The campaign climatology presented here aims to provide a valuable dataset to inform model simulation and future process studies, particularly in the context of aerosol-cloud interaction and further evaluation of dynamical processes in the SEP region for conditions analogous to those during VOCALS-REx. To this end, our results are discussed in terms of coastal, transitional and remote spatial regimes in the MBL and FT and a gridded dataset are provided as a resource.

High-Resolution Numerical Modeling of Thermally Driven Slope Winds in a Valley with Strong Capping

Abstract The complete day–night cycle of the circulation over a slope under simplified idealized boundary conditions is investigated by means of large-eddy simulations (LES). The thermal forcing is given with a time-varying law for the surface temperature. A surface layer parameterization based on the Monin–Obukhov similarity theory is used as a wall layer model. The domain geometry is symmetric, having an infinitely long straight valley in the y direction. Since the depth of the katabatic flow in midlatitude climates is limited to 5–30 m, the authors introduced a vertically stretched grid to obtain a finer mesh near the ground. The length scale for the calculation of eddy viscosities is modified to take into account the grid anisotropy. A preintegration of 24 h is made to obtain a capping inversion over the valley. Results show that the model is able to reproduce microscale circulation dynamics driven by thermal forcing over sloping terrain. The diurnal growth of the convective boundary layer leading to the development of the anabatic wind as well as the evolution of the cold pool in the valley during the night and its interaction with the katabatic flow are shown. Waves develop at the interface between the anabatic current and the return flow. During the day, as a combined effect of the geometry and the forcing, a horizontal breeze develops directed from the middle of the valley toward the ridges. The impact of the gravity current on the quiescent atmosphere in the valley generates a weak hydraulic jump during the night.

The Association of the Elevated Mixed Layer with Significant Severe Weather Events in the Northeastern United States*

AbstractThe occurrence of rare but significant severe weather events associated with elevated mixed-layer (EML) air in the northeastern United States is investigated herein. A total of 447 convective event days with one or more significant severe weather report [where significant is defined as hail 2 in. (5.1 cm) in diameter or greater, a convective gust of 65 kt (33 m s−1) or greater, and/or a tornado of F2 or greater intensity] were identified from 1970 through 2006 during the warm season (1 May–30 September). Of these, 34 event days (7.6%) were associated with identifiable EML air in regional rawinsondes preceding the event. Taken with two other noteworthy events in 1953 and 1969, a total of 36 significant severe weather events associated with EML air were studied via composite and trajectory analysis. Though a small percentage of the total, these 36 events compose a noteworthy list of historically significant derechos and tornadic events to affect the northeastern United States. It is demonstrated that plumes of EML air emanating from the Intermountain West in subsiding, anticyclonically curved flows can reinforce the capping inversion and maintain the integrity of the EML across the central United States over a few days. The EML plume can ultimately become entrained into a moderately fast westerly to northwesterly midtropospheric flow allowing for the plume’s advection into the northeastern United States. Resultant thermodynamic conditions in the convective storm environment are similar to those more typically observed closer to the EML source region in the Great Plains of the United States. In addition to composite and trajectory analysis, two case studies are employed to demonstrate salient and evolutionary aspects of the EML in such events. A lapse rate tendency equation is explored to put EML advection in context with other processes affecting lapse rate.

Observations of the Relationship between 700-mb Temperatures and Severe Weather Reports across the Contiguous United States

Abstract The relationship between 700-mb temperatures and convective severe storm reports is examined using data from 1993 to 2006 for the contiguous United States. Severe storm reports are used as a rough “proxy” for the occurrence of deep moist convection, and spatial and temporal distributions of 700-mb temperatures associated with these reports are analyzed. Secondarily, the distributions are assessed by individual severe storm report type, and convective inhibition also is evaluated. The motivation for this study is derived from the occasionally used 10°–12°C at 700 mb rule of thumb for estimating the extent and strength of the capping inversion. Whereas there is a semblance of merit for using this rule at times, its utility is shown to be strongly dependent on 1) geographic location, particularly with respect to surface elevation and the frequency of elevated mixed layers, and 2) the time of year. Calculation of convective inhibition, careful examination of the sounding, and assessment of lifting mechanisms likely are more valuable than 700-mb temperatures when forecasting the potential for deep moist convection and severe storms.

The Echo Size Distribution of Precipitating Shallow Cumuli

AbstractS-band radar surveillance scans of precipitating shallow convection are analyzed. The scans are complied from 52 days of near-continuous measurements in the winter trades of the North Atlantic during the Rain in Cumulus over the Ocean (RICO) field campaign. After being analyzed and filtered to exclude spurious returns, the scans are segmented to identify contiguous returns, or echoes. The echo size and reflectivity statistics are then analyzed. A new normalization method is developed to account for biases associated with the nonuniformity in the native radar grid. The echo area distribution robustly exhibits power-law scaling up until sizes of about 10 km2, with a scaling exponent of about −1.1. At larger sizes the scaling behavior breaks down and varies more markedly across subsamples of the data. Conditional sampling suggests that the scaling behavior of the larger echoes does, however, approach that of the smaller echoes as echo coverage increases, which supports the idea of a limiting distribution. Departures from this limiting distribution are argued to reflect finite size effects, modulated by the presence of a capping inversion whose height and strength varies across the samples.

Jets on daytime boundary layer capping inversions

Jets on daytime boundary layer capping inversions

Characterising the background errors for the boundary-layer capping inversion

The one-dimensional variational assimilation of vertical temperature information in the presence of a boundary-layer capping inversion is studied. For an optimal analysis of the vertical temperature profile, an accurate representation of the background error covariances is essential. The background error covariances are highly flow-dependent due to the variability in the presence, structure and height of the boundary-layer capping inversion. Flow-dependent estimates of the background error covariances are shown by studying the spread in an ensemble of forecasts. A forecast of the temperature profile (used as a background state) may have a significant error in the position of the capping inversion with respect to observations. It is shown that the assimilation of observations may weaken the inversion structure in the analysis if only magnitude errors are accounted for as is the case for traditional data assimilation methods used for operational weather prediction. The positional error is treated explicitly here in a new data assimilation scheme to reduce positional error, in addition to the traditional framework to reduce magnitude error. The distribution of the positional error of the background inversion is estimated for use with the new scheme.

A comparison between latent heat fluxes over grass using a weighing lysimeter and surface renewal analysis

An experiment to study the performance of the hourly sensible ( H ) and latent heat ( LE ) estimates over grass using surface renewal, SR , analysis was carried out at the University of California West Side Field Station near Five Points (San Joaquin Valley, CA). Hourly eddy covariance, EC , and SR sensible heat fluxes were averaged from their respective half-hourly fluxes. They were used to estimate hourly latent heat flux, LE , by calculating the residual of the energy balance equation, i.e., LE EC = R n − G − H EC and LE SR = R n − G − H SR , where R n is the net radiation; G is the soil heat flux; H is the buoyant sensible heat flux and the H subscripts denote measurements values from the EC and SR methods. Regardless of the weather conditions and time of day, H SR was comparable to H EC ; especially from sunrise until the afternoon formation of a capping inversion due to regional advection of sensible heat flux. The main source of energy contributing to LE was ( R n − G ). When ( R n − G ) was positive, on average, H contributed to 9% of ( R n − G ) when data were collected under unstable conditions, but it accounted for 21% under stable conditions . In general, correlation between measured LE using a weighing lysimeter, LE Lys , and LE SR was R 2 = 0.97 which was substantially higher than between LE Lys and ( R n − G ), R 2 = 0.93. In principle, the observed performance can be extrapolated to all irrigated, short, and dense-canopy vegetation that is mostly decoupled from the environment. Because for these surfaces the zero-plane displacement and the aerodynamic surface roughness length can be estimated as a portion (about 0.7 and 0.12, respectively) of the canopy height, using a two dimensional sonic anemometer capable to record high-frequency sonic temperature provides all the input to determine H SR which in conjunction with ( R n − G ) is required to estimate LE .

Observations of Marine Atmospheric Boundary Layer Transitions across the Summer Kuroshio Extension*

Abstract The baiu and Kuroshio Extension (KE) fronts, both zonally oriented and nearly collocated east of Japan, are the dominant summertime features of the atmosphere and ocean, respectively, over the midlatitude northwest Pacific. An atmospheric sounding campaign was conducted on board the R/V Roger Revelle during the 2005 summer. Transects of soundings across the KE front are analyzed to study its effects on the atmosphere, along with continuous surface meteorological and ceilometer cloud-base observations. While the KE front remained nearly stationary during the cruise, the baiu front displayed large meridional displacements that changed wind direction across the KE front. The presence of sharp sea surface temperature (SST) gradients anchored by the KE enhanced the thermal and moisture advection, causing substantial changes in the marine atmospheric boundary layer (MABL) structure. When the baiu front was displaced north of the KE front, southwesterly winds advected warm, humid air from the subtropics over the cold water, producing a surface inversion favorable to fog formation. When the baiu front was to the south, on the other hand, northerly winds across the KE front destabilized the MABL, leading to the formation of a solid low-cloud deck beneath a strong capping inversion. The wind changes with the meridional displacement of the baiu front thus caused large variations in near-surface atmospheric stability and surface turbulent heat fluxes, with potential feedback on deep convection and fog/low-cloud formation around the front.