3.1. Free Decaying Convective Boundary Layer
The simulation without backgound wind shows developing and decaying phases of the free CBL and its transition to a residual layer. To illustrate the PBL depth and its time evolution, we calculate the inversion height
, defined as the level of minimum sensible heat flux, and the level of maximum vertical gradient of potential temperature
(
Figure 1a) [
24]. The periods when the magnitude of minimum heat flux is smaller than 0.005 K m
and when the maximum heat flux is negative are excluded from the calculation of
. The inversion height increases until
and then stays almost constant for ∼540 s. Then, the heat flux becomes negative at all the vertical levels (see later discussions and Figure 5c), implying that the vertical structure of CBL is broken. The time interval between
and the CBL breakage is comparable to
(648 s at
). The level
starts from around the bottom of initial capping inversion (0.7 km) and increases to 1010 m at
t =
, and then it is almost constant even after the CBL breakage. Note that
is higher than
by 30–70 m because of their different definitions [
24]. Actually,
represents the entrainment interface, irrespective of heat flux profile below, and it is close to the top of entrainment zone or the top of residual layer while
is usually located in the middle of the entrainment zone.
Figure 1b shows the time series of the TKE averaged below
. After 45-min spinup, the vertically averaged TKE increases slowly until
t =
, followed by its abrupt decrease after ∼200 s. The time evolution of the minimum heat flux is presented in
Figure 1c. The minimum heat flux represents the strength of top entrainment, proportional to the activity of overshooting thermals and downward heat transport [
24]. After the spinup, its magnitude increases slightly until
t =
and then decreases immediately.
Figure 2 shows the vertical velocity fields at
and
z = 40, 460 (
), and 940 m (
) and the same fields 604, 1204, and 3603 s after
. Convection cells with narrow branches of near-surface updrafts, downdrafts and narrow circumferential updrafts in the mixed layer, so-called cellular convection, and overshooting thermals (or updrafts) across the entrainment zone are clearly visible at
(
Figure 2a,e,i). Turbulent eddies seem to decay quickly from the bottom up. For example, the narrow branches of updrafts near the surface become weak and broad within a period of
(
Figure 2i,j). The number of overshooting updrafts decreases, too, corresponding to the decrease of minimum heat flux after
(
Figure 1c). In contrast, the cellular up- and downdrafts in the middle of CBL maintain their strength and position for more than 604 s. For example, updrafts on the rising branches are still faster than 1 m
at
+ 604 s. This illustrates that boundary layer thermals survive for a time period in the order of
whereas near-surface eddies decay quickly. The cellular up- and downdrafts decay and overshooting updrafts disappear completely afterwards. One hour after
, the up- and downdrafts become much weaker than before and they lost their cellular structure at
z =
, now in the residual layer.
The vertical velocity fields in the
y-
z plane (
Figure 3) and the potential temperature fields in the same
y-
z plane (
Figure 4) illustrate the vertical structure of decaying convection cells. While the near-surface updrafts decay quickly, the convection cells seem to survive longer. At the initial decaying phase, the individual updrafts of the cells become weak but diffused, filling gaps between updrafts and making wider convection cells (
Figure 3b,c). During the decay, small downdrafts penetrate down through cellular updrafts but convective circulation and diffusing motions seem to maintain the boundary layer circulations until
+ 905 s. The large convection cells, however, break into smaller eddies 1202 s after
(
Figure 3e). The black contours of 301.51 K potential temperature in
Figure 3c–g illustrate thermal stratification in the middle of PBL briefly, e.g., the air below and above the contours are cooler and warmer than 301.51 K, respectively, as can be confirmed in
Figure 4. The value of 301.51 K is selected to visualize the thermal stratification as best as possible. Cellular up- and downdrafts are shown to be located in cooler and warmer regions, respectively, at
= 604, 905, and 1202 s (
Figure 3c–e). Thus, updrafts push cool air upward and downdrafts drag warm air downward during the decay. The potential temperature fields in
Figure 4 confirms that cellular up- and downdrafts stack cooler and warmer air, respectively, inducing the undulating distribution of potential temperature at
= 604, 905, and 1202 s. This kind of undulating distribution appears several hundred seconds more but becomes less and less distinct over time (not shown). Along with the boundary layer scale circulations, local downdrafts touch the bottom surface as the cold pools from deep convective clouds do [
25], but they are not spreading out at the bottom surface because they are warmer than near-surface air (not shown). This kind of abnormal and unsustainable circulations weaken gradually with the contours getting flatter over time. Eddies, with scales of several hundred meters, still exist ∼1 h after
but they disappear completely two more hours later with a stably stratified residual layer (
Figure 3f,g).
The vertical structure of the decaying CBL is presented in the time series of vertical profiles of TKE, vertical velocity variance
, and vertical heat flux
(
Figure 5). In this study, overbars and angle brackets denote temporal (60 s) and horizontal averages and primes denote perturbations from the horizontal averages. The TKE in the mixed layer maintains the initial level for several hundred seconds, for example, TKE at
lost 5% in 360 s, and then decays from the bottom up. The
e-folding decay time of TKE also increases with height, confirming the bottom-up decay. The variance of vertical velocity, representing vertical turbulence, is larger in the middle of CBL than near the surface or the CBL top. The
e-folding decay time increases with height, too, but it is shorter than that of TKE especially near the surface. This indicates that the near-surface TKE is maintained by horizontal diverging and converging motions induced by subsiding downdrafts. The profiles of vertical heat flux illustrate that the vertical thermal structure of CBL, positive and negative heat flux in and above the mixed layer, is maintained until
+ 540 s, then negative heat flux propagates down to the bottom. After the negative heat flux dominates throughout the boundary layer, an oscillation of heat flux occurs despite of its small amplitude. This downward propagation and the following oscillation were simulated in previous numerical studies [
10,
12] and observed in the real atmosphere [
5,
7]. These are known to be induced by demixing of air parcels, entrained from above the PBL top, and their returning to equilibrium levels [
12]. The
e-folding time of vertical heat flux increases vertically in the mixed layer but the time scale is shorter than that of TKE. This is attributable to that potential temperature perturbations decay faster than velocity perturbations [
10].
The downward propagating negative heat flux and the following heat flux oscillation have been explained by circulations of demixed (or non-mixed) air parcels after small-scale turbulence decays [
10,
12] but the explanation has not been proven yet. For this, a quadrant analysis of
and
is done (
Figure 6). Perturbations at every grid point and every ∼60 s are classified into the four quadrants—warm air rising, cool air rising, cool air sinking, and warm air sinking events—and the classified events are averaged to show the contribution of the individual quadrants to vertical heat transport [
24,
26]. The first and third quadrants represent the contributions by rising thermals and subsidence, respectively. They are the main upward heat transporters in the mixed layer. The other two quadrants do a minor role in the mixed layer, but they are important in the entrainment zone. Overshooting thermals become cooler than the stably stratified environment and thus their contribution to downward heat flux is represented by the second quadrant, cool air rising. The fourth quadrant, representing warm air sinking events, is weaker than the second quadrant not only in the entrainment zone but also in the lower mixed layer. However, the contribution by the fourth quadrant events, especially in the lower mixed layer, becomes stronger several hundred seconds after
. Thus, warm air sinking events appear to induce the downward propagation of negative heat flux. After the cutoff of surface heat flux, near-surface air is a little cooler than previously heated and already lifted air and thus the lower mixed layer becomes weakly stable. We speculate that demixed downdrafts from the entrainment zone or from the middle of CBL are now warmer than the environment (
Figure 3c–e) and they contribute to the downward propagating negative heat flux. However, the histories of demixed downdrafts are still questionable and a more sophisticated method (e.g., a Lagrangian tracking) is required to reveal the mechanism of this downward propagation completely.
Time series of the pattern correlation between the perturbations at
and later perturbations illustrate the decaying characteristics in a different perspective (
Figure 7). This kind of analysis is meaningful only for the free CBL because convection cells develop and decay nearly at the same place. The pattern correlation of the vertical velocity shows that convective turbulence in the middle of CBL decays slower than that near the surface or at the CBL top, as the vertical velocity variance decays (
Figure 5b). The pattern correlation at
, for instance, is higher than 0.4 at 600 s (
) after
. The pattern correlation of potential temperature decreases more quickly than that of vertical velocity. In contrast to potential temperature, the pattern correlation of passive scalar concentration decreases slower than that of vertical velocity, implying that passive scalar has a longer memory. It is also notable that the pattern correlation of passive scalar has maximum peaks at higher levels than that of vertical velocity. This may be due to the fact that scalar perturbations have less resistance (or more inertia) in the entrainment zone than thermal and momentum perturbations, but this needs a further investigation.
To demonstrate the height- and scale-dependent characteristics of convective turbulence, two-dimensional spectra of vertical velocity and potential temperature near the surface and at the middle of CBL are calculated and plotted in
Figure 8. In this study, Fourier coefficients in a two-dimensional wavenumber space are calculated first, and then the radial averages of the coefficients are computed. While one-dimensional spectra can show only one directional spectral energy, this kind of two-dimensional spectra can represent spectral energy in all horizontal directions, thus illustrating spectral energy of convection cells better than its one-dimensional counterpart. The spectra of vertical velocity have maximum peaks in the wavelength range of 1280–2560 m and monotonically decrease with increasing frequency (decreasing wavelength) at the middle of CBL (
Figure 8b). The spectra near the surface are flat in the wavelength range of 320–1280 m with cascading in the smaller wavelength range (
Figure 8a). The flat spectral range is attributed to the near-surface local eddies and the near-surface spectral energy decays quickly shortly after the surface heat flux is stopped. In contrast to the fast decay near the surface, the spectrum of vertical velocity at the middle of CBL changes very little for the initial 362 s. This demonstrates that nonlocal eddies are active for the initial several hundred seconds. Then, the nonlocal eddies decay over time along with the decaying near-surface local eddies. The spectral energy of potential temperature near the surface decays quickly over time (
Figure 8c). For example, the sum of spectral energies (=variance) of potential temperature decreases by a factor of 67.2 for the initial 362 s. At the middle of CBL, the spectral energy of potential temperature decays more slowly than that near the surface, too (
Figure 8d).
3.2. Advective Decaying Convective Boundary Layer
The simulation of advective decaying CBL reveals the role of mechanically generated eddies near the surface.
Figure 9 shows vertical velocity fields in the
y-
z plane, illustrating the vertical structure of convective rolls. Convective rolls, pairs of linearly aligned and circulating up- and downdrafts, are being advected in the streamwise direction during their developing and decaying phases and thus tracking of individual rolls is difficult at a fixed position. Nonetheless, decaying characteristics in this advective CBL are identifiable in the series of vertical velocity fields in the fixed
y-
z plane (
Figure 9) and the characteristics are similar to those in the free CBL. For example, convective rolls become weak but a little wider and then break into much smaller eddies (
Figure 9a–e). The abnormal circulations observed in
Figure 3, in which cool air rises and warm air sinks, are less distinct in this advective CBL. Moreover, up- and downdrafts in the convective rolls seem to decay more slowly than those in the free CBL. It is also notable that local eddies are generated constantly near the surface even after nonlocal eddies disappear and they survive until the end of the 6-h simulation. However, in the absence of surface cooling and strong stratification, the surface eddies grow up to several hundred meters, thus the eddies do not seem to be local at the end of the simulation.
The potential temperature fields in
Figure 10 confirm that the undulating potential temperature distribution in the free decaying CBL (
Figure 4) is not very distinct in the advective decaying CBL. It is also remarkable that the regions cooler than 301.4 K in the advective decaying CBL are wider than those in the free decaying CBL. Another remarkable feature is that potential temperature at the end of the simulation seems to be less stratified than in the free decaying CBL (
Figure 10g) due to the more active near-surface eddies in the advective decaying CBL (
Figure 9g).
The time series of vertical profiles of TKE, vertical velocity variance, vertical heat flux, vertical momentum flux
, and horizontally and temporally averaged streamwise velocity
are plotted in
Figure 11. While the vertical velocity variance and the heat flux decay from the bottom up as in the free decaying CBL (
Figure 5b,c), near-surface TKE is higher and decays slower than TKE in the middle of CBL (
Figure 11a). This slower decay is attributable to the near-surface eddies such as sweeps and ejections interacting with the convective rolls (not shown). Note that the
e-folding times of TKE, vertical velocity variance, and vertical heat flux in this advective CBL are longer than those in the free convective boundary layer, as shown by the black and gray dots in
Figure 11. For instance, the
e-folding time of TKE at
is 1261 s in the advective CBL and it is 966 s in the free CBL. We speculate that surface wind shear induces more lasting vertical circulations. Momentum flux decays quickly, too, except near the surface. Downward momentum flux is maintained to a certain degree and the vertical range of negative momentum flux extends upward over time (
Figure 11d). This is consistent with the upward growth of the surface eddies seen in
Figure 9. As a result, mean flow in the lower CBL or in the lower residual layer (after the CBL breakage) is decelerated over time until the end of the simulation (
Figure 11e). In the real atmosphere, however, the growth of the surface eddies is more suppressed in the shallow and stable boundary layer over the radiatively cooled surface.
Two-dimensional spectra of vertical velocity and streamwise velocity in the advective CBL are presented in
Figure 12. Near the surface and in the wavelength range of 320–1280 m, the spectra of vertical velocity are quite flat for the initial 181 s but large eddies decay faster than small eddies, and parabolic spectra of vertical velocity are left finally (
Figure 12a). The wavelength of the spectral peak ranges between 160 and 320 m, matching with the scale of the near-surface eddies (
Figure 9). The spectrum of vertical velocity at the middle of CBL stays almost the same for the initial 361 s as in the free CBL. Since then, the spectra decrease and the gap between two successive spectra becomes narrower over time than the gap in the free CBL (
Figure 8b), indicating the slower decay in the advective CBL. The slower decay, distinct in the wavelengths smaller than 640 m, is related to the surface eddies growing up to several hundred meters (
Figure 9g). The spectra of streamwise velocity near the surface and in the wavelengths smaller than 320 m decay very slowly for the three hours after
but the spectra of larger eddies decay faster as the spectra of vertical velocity (
Figure 12a,c). At the middle of CBL, the spectral energies of streamwise velocity stay almost the same for the initial 721 s and then decay but the impact of the surface eddies is hardly detectable (
Figure 12d).