A traveling atmospheric disturbance generated by a soil colour change in a high-resolution climate model experiment

Abstract

The climate model CESM-WACCM is used to study the way a soil colour change of the eastern region of the Sahara affects the dynamics of the troposphere. The soil colour is darkened for 5 days. The difference between the perturbed model run and the control model run is used to isolate the soil colour change-induced atmospheric perturbation from random atmospheric waves which are stronger by an order of magnitude or more. The perturbation generates a circular wave radially propagating away from the Sahara on the first day of the simulation. After nine hours, the wave front reaches the convection zone in Brazil where a secondary wave is generated and can be clearly seen until 23:00 UT. The mean wave velocities of the traveling atmospheric disturbances are 〈v〉=200±50 m/s for the primary wave and 〈v〉=220±40 m/s for the secondary wave. The mean horizontal wavelengths are 〈λ〉=3000±500 km for the primary wave and 〈λ〉=2600±600 km for the secondary wave. The mean wave periods are 〈p〉=4±1 h for the primary wave and 〈p〉=3±1 h for the secondary wave. Since the perturbed model run diverges from the control run with the passage of time, the attribution of cause and effect becomes difficult after a few days. Analysis of the simulation data of the first day leads to a deeper understanding of global teleconnections, radiative transfer and wave-coupling processes between the surface and the atmospheric layers.

Background

In our study we investigate whether a high-resolution global climate model can be used to study the generation, propagation and dissipation of atmospheric waves induced by a soil colour change. We are especially interested to know if a cause-effect study can be performed under realistic atmospheric conditions. The attribution of cause and effects is a fundamental problem of complex systems.

At present two types of numerical models are available for studying the climate processes. The idealized models and the realistic models [1]. Idealized models allow us to conduct detailed studies and hence increase our understanding of simplified or isolated processes. Due to their simplicity these models neglect the role of more complex interactions and hence do not provide a holistic view of the climate system. Realistic models on the other hand, consider almost all processes of the climate system and provide high-resolution climate data, which are rather close to the observations. As a result realistic model data are almost as difficult to understand as the observations. Therefore the recognition of the relevant climate processes, the attribution of climate forcing and their effect on the climate requires advanced methods of data analysis and interpretation (e.g [2]–[4]).

The generation, propagation and dissipation of atmospheric waves are often studied by means of idealized model simulations. Nicholls & Pielke [5] conducted such a study using an idealised three-dimensional, fully compressible atmospheric model. Their goal was to investigate the properties of the atmospheric waves induced by a tropical thunderstorm. Their study showed that the thunderstorm generated not only atmospheric gravity waves, but also thermal compression waves. Gardner & Schunk [6] used a high-resolution global thermosphere-ionosphere model to examine the effect a large scale perturbation (in this case a pulsating geomagnetic storm) has on the atmosphere. They concluded that this type of storm generates multiple Traveling Atmospheric Disturbances (TADs) coincidently travelling from the northern and southern auroral zone towards the equator and into the conjugate hemisphere.

Atmospheric disturbances are also induced by albedo changes, which can be due to changes in the type of the vegetation, the soil colour, the formation of clouds and other causes. Surface albedo itself can include snow cover, surface water and vegetation change too. Only in the case of a desert like the Sahara, surface albedo might be mainly governed by the soil color. Seitz [7] gives a brief survey about the influences of albedo changes on the Earth’s radiation energy balance. According to Seitz [7], the effects of regional albedo changes on the climate system are still not well observed, simulated and understood. Past regional albedo change studies have focused on the effects on long-term climate change [3],[4],[8],[9] and how it can be utilized in geoengineering projects to combat global warming [10]. To our knowledge, there is only one study, which attributed a sudden snow cover in Eurasia to a polar vortex change [3].

Our study is the first to extract and discuss the tropospheric waves induced by a regional soil colour change using a realistic high-resolution 3-D climate model. In the simulation, we change the soil colour of a desert region since this change can be safely implemented in the complex climate model which we use. Furthermore, this simple scenario will allow a clear interpretation of the simulation results.

In the present study, we concentrate on the effects the soil colour change has on the lower troposphere. The generated atmospheric perturbation has a small amplitude so that a small-scale perturbation analysis can be applied to the first 1-2 days of the model simulation. This allows us to study the global propagation of the TADs through the lower troposphere in detail.

Model description

The Community Earth System Model (CESM) version 1.04 was used to perform our simulation. It is composed of a coupler (CPL) and five fully coupled geophysical models: atmosphere (ATM), land (LND), ocean (OCN), sea-ice (ICE), land-ice (GLC). The models can be set as fully prognostic, data, or stub and are “state-of-the-art climate prediction and analysis tools” [11] when set in prognostic mode.

The ATM in our simulation is the Whole Atmosphere Community Climate Model (WACCM) version 5 [12]. WACCM is often used for the simulation of circulation, thermal tides, gravity waves, wave-mean flow interaction, and atmospheric composition changes in the middle atmosphere [13]–[20].

It has a fully compressible horizontal discretization, and a quasi-Lagrangian vertical discretization approximation, which ignores the acceleration term in the vertical component of the momentum equation. This approximation is good for scales greater than 10 km [12]. It has 66 vertical levels from the ground up to 5·10⁻⁶ hPa (2.5 – 149 km). The vertical coordinate is purely isobaric above 100 hPa, but is terrain following below that level. The model top is ∼ 150 km. The vertical resolution is 1.1 km in the troposphere, 1.1-1.4 km in the lower stratosphere, 1.75 km at the stratopause and 3.5 km above 65 km. The horizontal resolution of our simulation is 4° × 5° (latitude × longitude), with 72 longitude and 46 latitude grid points. The coupler timestep is Δ t=30 minutes while the time step for the dynamics equations is Δ τ=Δ t/8 [12].

The smallest wavelength that our model experiment can resolve by the finite differencing (that is used by default in WACCM) is twice that of the grid size [21]. Therefore the resolvable waves in the equatorial troposphere have a horizontal wavelength of λ _horizontal >1000 km, a vertical wavelength of λ _vertical >2.2 km, a period of p _horizontal >2 h and a wave velocity of v<Δ x/Δ τ=Δ x/(Δ t/8)=500/(1800/8)≈2200 m/s. Thus the simulation can adequately resolve large scale waves [21].

The land model provides the surface albedo, area-averaged for each atmospheric column, and the upward longwave surface flux, which incorporates the surface emissivity, for input into the radiation scheme. The surface fluxes of momentum, sensible heat, and latent heat serve as the lower flux boundary conditions for the planetary boundary layer parameterization, the vertical diffusion and the gravity wave drag. The atmospheric radiation is calculated using the momentum, sensible heat flux, latent heat flux, land surface albedos and upward longwave radiation. The upward longwave radiation is calculated by taking the difference between the incident and absorbed fluxes. The incident flux values are determined by means of the daily values of the solar radio flux (F10.7) which are provided by the National Oceanic and Atmospheric Administration’s (NOAA) Space Environment Center [22].

The land model used was the Community Land Model (CLM), which has atmosphere-surface coupling, surface colour variability, surface albedo calculation, absorption, reflection, and transmittance of solar radiation, absorption and emission of longwave radiation, sensible heat (ground and canopy) latent heat fluxes and heat transfer in soil and snow [23]. The surface albedo calculation depends on whether the top surface is a vegetated canopy or bare ground.

In the case of bare ground the surface albedo calculation requires only the soil colour, which varies with the colour class (soil colour index). The term bare ground refers to any surface that is not covered by vegetation and can therefore be a glacier, a lake, a wetland, snow covered soil or bare soil. In the case of our simulation the only parameter that affects the surface albedo of the perturbed east Saharan area is the soil colour, which is determined by the value of the soil colour index of the area.

The CLM soil colour indices are prescribed so that they reproduce the observed Moderate-resolution Imaging Spectroradiometer (MODIS) local solar noon surface albedo values at the CLM grid cell [23]. MODIS is a set of spectroradiometers in orbit on board the Terra and Aqua satellites. They provide measurements of large scale global dynamics (e.g changes in the cloud cover, radiation budget, processes on the oceans, the land and lower atmosphere). They capture data in 36 spectral bands (0.4 μ m - 14.4 μ m) at varying spatial resolutions and image the entire Earth every 1 - 2 days [24].

Methods

Simulation setup

In the simulation setup we used the perpetual year 2000 component set (F_2000_WACCM) [11]. A component set is the assemble of a particular mix of geophysical models along with geophysical model-specific configurations and namelist settings. In this case it is a set of two fully prognostic, present day, coupled atmosphere and land geophysical models (WACCM,CLM), a prescribed data ocean geophysical model (docn), a prescribed sea-ice (CICE) geophysical model and no land-ice geophysical model ([11]). With it, we conduct two simulations. In the first one, from now on referred to as control run, all the input parameter fields remain unchanged. In the second one, from now on referred to as perturbed run, the colour of a small region of the east Sahara, that spans 3 × 3 pixels (longitude: 10°-20°, latitude: 18°-26°), is darkened from beige (bright sand) to very dark green (the colour of the darkest forests). To avoid a discontinuity in the soil colour map, the soil colour change is gradually performed with a small change of 50% at the edge of the Sahara array (soil colour index = 15) and a maximal change in the center of the array (soil colour index = 20) as shown in Figure 1. The simulation spans 5 days starting at 0:00 UT on 01/01/2000 and ending at 23:00 UT on 05/01/2000.

Figure 1

Simulation setup. (a) soil colour index of the control run, (b) soil colour index of the perturbed run, (c) difference between the soil colour indices of the perturbed and the control run.

Data analysis

The surface albedo, the surface temperature, the atmospheric pressure and the vertical wind are extracted from the output datasets for the control and the perturbed run. In the present letter, we only show the results obtained for the vertical wind at 2 km altitude. A comprehensive study of the disturbances in all parameters at all altitudes is planed as a follow-up-study. In addition we have to design advanced algorithms for the analysis of different wave modes. However the initial data analysis presented here already provides many new results.

As a first step the above parameters are interpolated to the altitude of 2 km from 0:00 UT to 23:00 UT. Then the control run is subtracted from the perturbed run at each timestep. The mean standard deviation σ=2·10⁻⁴ Pa/s for the first day at 23:00 UT is taken as reference for the normalization of the vertical wind fluctuations. To calculate it we derive at first the zonal means of σ as a function of latitude for 23:00 UT. Then we obtain the global mean of sigma σ from the zonal means of sigma σ _z by calculation of the surface area preserving mean. For recognition of significant atmospheric waves in the global plots, we divide the vertical wind fluctuations d Ω by σ. Variations > 2 σ have a significance (confidence level) of 90%.

Results and discussion

The change in surface colour resulted in the appearance of fast propagating primary and secondary waves at 2 km altitude.

Primary perturbation

As can be seen in Figure 2a, the primary perturbation, appears over the Sahara at 9:00 UT (four hours after sunrise). The surface colour change causes a convective perturbation in the lower troposphere, a buoyancy oscillation that rises and falls as the day progresses (Figure 2). The warm colours in the figure indicate downward air motion and the cold colours indicate upward motion since WACCM provides the vertical wind in [Pa/s].

Figure 2

Evolution of the atmospheric perturbation due to the soil colour change in Eastern Sahara during day 1. The difference between the vertical wind at 2 km altitude of the perturbed run minus the control run is shown at (a) 9:00 UT, (b) 12:00 UT, (c) 15:00 UT, (d) 23:00 UT. Shades of blue indicate upward motion of the air, while shades of red indicate downward motion.

As can be seen in Figure 2b,c,d the perturbation amplitude increases with time to a value of A>5σ. We have therefore a 5σ confidence (σ=2·10⁻⁴ Pa/s at this altitude).

Figure 2b shows a nearly circular wave ring outgoing from the surface colour change region in the eastern Sahara. The wave ring is also visible in Figure 2c (15:00 UT). The mean horizontal wavelength 〈λ〉 of the wave is calculated as follows:

• i. The center of the disturbance is known for the primary wave (22° N, 15° E). In the case of the secondary wave, its center is determined by an educated guess (10° S, 55° W).

• ii. The positions of the inner and outer ring of the wave (R _inner and R _outer ) are visible in the global maps of the vertical wind fluctuations at 12:00 UT). Then the ring positions are manually determined, from the center towards the North, East, South and West direction. This gives us four values for the inner and four values for the outer ring.

• iii. Then the mean horizontal wavelength is calculated: 〈λ〉 = <λ>=2·(<R _outer >−<R _inner >).

The calculated values for <R _inner > and <R _outer > together with their standard deviations can be seen in Figure 3. The calculated mean horizontal wavelength of the primary wave is 〈λ〉=3000 km.

Figure 3

Evolution of the inner and outer radius of the primary wave ring. The three point pairs of the inner and outer radius are at 10:00 UT, 11:00 UT and 12:00 UT. The horizontal wavelength is two times the mean distance between the two linear regression lines.

The mean wave speed v is calculated as follows. First the mean horizontal speeds of the inner and outer radii are taken from the inclination of the above mentioned linear regression lines: v _inner =d R _inner /d t=240 m/s, v _outer =d R _outer /d t= 170 m/s. Then the mean of the above calculated horizontal speeds is taken 〈v〉=(v _inner +v _outer )/2≈200 m/s. Its standard deviation is σ _v =50 m/s.

The mean wave period is calculated taking the ratio of the mean horizontal wavelength and the mean horizontal velocity of the wave 〈p〉=〈λ〉/〈v〉=4 h. Its uncertainty is calculated using the Gaussian error propagation law${\sigma }_p=\sqrt{{\left({\sigma }_{\lambda }/〈v〉\right)}^2+{\left(〈\lambda 〉·{\sigma }_v/{〈v〉}^2\right)}^2}=1h$.

The derived uncertainties mainly reflect the azimuthal variations of the wave parameters. These variations are possibly due to azimuthal changes in the background wind flow, topography, convective activity, Coriolis force and other factors which modulate the radial propagation of the wave rings.

In summary, the soil colour change induced a circular large scale wave with a wave speed close to the speed of sound. The wave characteristics are quite similar to those in the numerical simulation of Nicholls & Pielke [5]. They interpreted the horizontally propagating circular wave outgoing from a thunderstorm region with the speed of sound as a Lamb wave mode n1 and n2. The n1 mode moved faster and resulted in deep subsidence warming, while the second mode propagated at half the speed and resulted in weak-low level uplift. The appearance of an uplift as a result of a thermally induced gravity wave is in accordance with our own results (Figure 2a). Parallel to Nicholls & Pielke, the primary wave in our simulation is mostly soliton-like which might be interpreted as a shock wave front (as done by Nicholls & Pielke). On the other hand large scale waves with speeds close to the speed of sound are often found in the thermosphere. They have periods of several hours and horizontal wavelengths of a few thousands kilometres. Gardner & Shunk [6] as well as Vadas [25] classify these large-scale waves as atmospheric gravity waves.

In the case of our simulation, the surface colour change resulted in an increase in the fraction of solar radiation absorbed by the region. This resulted in an increase in the surface temperature and therefore in an increase in the temperature of the air located directly over the perturbed region. The observed buoyancy oscillation at 2 km altitude could be a result of either convective or radiative heating. The convective heating would be consistent with the evolution of the central perturbation, as the air seems to rise over the perturbation and then fall as the day progresses, with the opposite happening to the air around it (Figure 2).

Secondary perturbations

Surprisingly, two secondary perturbations (or secondary waves) appear in Figure 4. The first appears over Indonesia at 14:00 UT. The pathway between this perturbation and the source region, in the Sahara, is unclear because the primary wavefront has not arrived in Indonesia at the time the perturbation appears.

Figure 4

Evolution of the secondary wave above Brazil during day 1. The difference between the vertical wind at 2 km altitude of the perturbed run minus the control run is shown at (a) 16:00 UT, (b) 17:00 UT, (c) 19:00 UT, (d) 23:00 UT.

Later the wave front of the primary wave approaches Brazil at 15:00 UT (Figure 2c) and reaches its tropical convection zone at 16:00 UT (Figure 4a).

As the wave front comes in contact with the tropical convection zone, it scatters and a secondary wave is generated. The secondary wave is clearly visible as concentric wave rings outgoing from the center of the tropical convection zone (10° S, 55° W). The vertical wind in the tropical convection zone also has a periodic oscillation which can be seen in Figure 4a. The outward propagation of the generated secondary wave can be seen in Figure 4a,b.

As can be seen in Figure 4c at 19:00 UT another wave front generated by the perturbation over the Sahara, scatters on the Brazilian tropical convection zone. Besides wave scattering based on the Huygens-Fresnel principle, we think that the latent heat release by periodic, vertical advection of moist air masses is another reason for the generation and amplification of a secondary wave over Brazil (Figure 4a,c).

The parameters of the secondary wave are estimated for the time interval between 19:00 UT and 23:00 UT. The wave parameters are derived in the same manner as for the primary wave. The point pairs of the inner and outer rings and their linear regression lines are shown in Figure 5. The estimated wave parameters are: 〈λ〉=2600±600 km, 〈v〉=220±40 m/s, 〈p〉=3±1 h.

Figure 5

Evolution of the inner and outer radius of the secondary wave ring outgoing from Brazil. The five point pairs of the inner and outer radius are at 19:00 UT, 20:00 UT, 21:00 UT, 22:00 UT and 23:00 UT. The horizontal wavelength is two times the mean distance between the two linear regression lines.

In summary, the soil colour change induced a secondary circular large scale wave with a wave speed ∼ 70% of the speed of sound.

The generation of secondary waves from convectively generated gravity waves was simulated by Vadas [26], with the wave characteristics being quite similar to those in our simulation. The generation of secondary oceanic waves through interaction with the topography in a manner similar to the interaction of our wave with the tropical convection zone, was reported by Vlasenko [27].

The scattering of the primary wave front results in a change in the vertical motion. Updrafts might produce condensation of moist air resulting in latent heat release which may further amplify the secondary wave. Indeed the strong amplitude of the secondary wave points to such an amplification process which would affect the precipitation rate and the regional climate.

An ensemble simulation was also performed, between 01.01.2000 and 01.06.2000, with a 5 day gap between each run. The main features (primary wave and excitation of secondary waves) are present in both the ensemble average and in each individual run. However the shape of the primary wave ring can differ indicating a seasonal variation in the azimuthal wave speeds, which is likely due to a seasonal variation of the mean flow and the thermal structure of the lower troposphere. For example, in June the southward propagation of the primary wave is enhanced compared to January. In March, a secondary wave is generated in the vicinity of Lake Victoria (Tanzania), where a tropical convention zone is present during that time of the year (rainy season) [28]. These examples show the potential of wave propagation studies with a realistic high-resolution climate model.

Evolution five days later

Five days after the initiation of the simulation, no circular wave patterns are visible. And while enhanced fluctuations are visible in the vicinity of the Sahara and the convective zones over Brazil and Indonesia, disturbances are also visible in a seemingly random distribution all over the globe. Furthermore the global mean standard deviation σ linearly increases over the 5 days of our simulation from σ=2·10⁻⁴ on the first day to σ=27·10⁻⁴ on the fifth day. This indicates that the perturbed run diverges more and more from the control run with increase of time, making the attribution of cause and effect difficult. Our small perturbation analysis is only applicable for the first two days when the perturbation waves can be clearly separated from the random atmospheric waves. The coupling of the primary and secondary waves in our study could be exemplary to a similar coupling between migrating and non-migrating tides.

Conclusions

The soil colour change of a surface area in the eastern Sahara generates a traveling atmospheric disturbance. The advantage of our high-resolution climate model experiment is that we can study, for the first time, the generation and propagation of the TAD under realistic atmospheric conditions. This means that all the interactions between the TAD and the atmospheric jets, circulation cells, solar tides, planetary waves, random gravity waves, orography and tropical convection zones are included in the simulation. We find that the interaction of the TAD with the tropical convection zones over Brazil is most obvious.

The soil colour perturbation generates a primary atmospheric wave outgoing from the source region in the Sahara. It has an amplitude of A>5σ, a mean horizontal wavelength of 〈λ〉=3000±500 km, a mean wave velocity of 〈v〉=200±50 m/s and a mean period of 〈p〉=4±1 h.

As the primary wave reaches Brazil, the wave front scatters at the tropical convection zone and generates secondary atmospheric outgoing waves. The latent heat release by moist air amplify the secondary wave to an amplitude of A>5σ. The wave has a mean horizontal wavelength of 〈λ〉=2600±600 km, a mean wave velocity of 〈v〉=220±40 m/s and a mean period of 〈p〉=3±1 h.

In the present letter we only discussed the dynamical response of the global troposphere at 2 km height to the soil colour change in Sahara. This brief analysis of the numerical experiment already gave new results on the land-atmosphere interaction, global teleconnections and wave dynamics. We expect that a future analysis of the simulation data at different altitudes will provide us with a better understanding of how regional soil colour changes affect the atmospheric layers above. For example how much energy and momentum are transported by wave coupling and radiative transfer, which wave modes are most important and which role do the soil colour changes play on the generation of non-migrating atmospheric tides.

Another topic would be the investigation of the dependences of the primary and secondary TAD parameters on the location, size and strength of the sources. Such studies are relevant for geoengineering as well as for understanding of the climate processes and wave-wave interactions.