Winds in the Lower Cloud Level on the Nightside of Venus from VIRTIS-M (Venus Express) 1.74 µ m Images

The horizontal wind velocity vectors at the lower cloud layer were retrieved by tracking the displacement of cloud features using the 1.74 µ m images of the full Visible and InfraRed Thermal Imaging Spectrometer (VIRTIS-M) dataset. This layer was found to be in a superrotation mode with a westward mean speed of 60–63 m s − 1 in the latitude range of 0–60 ◦ S, with a 1–5 m s − 1 westward deceleration across the nightside. Meridional motion is signiﬁcantly weaker, at 0–2 m s − 1 ; it is equatorward at latitudes higher than 20 ◦ S, and changes its direction to poleward in the equatorial region with a simultaneous increase of wind speed. It was assumed that higher levels of the atmosphere are traced in the equatorial region and a fragment of the poleward branch of the direct lower cloud Hadley cell is observed. The fragment of the equatorward branch reveals itself in the middle latitudes. A diurnal variation of the meridional wind speed was found, as east of 21 h local time, the direction changes from equatorward to poleward in latitudes lower than 20 ◦ S. Signiﬁcant correlation with surface topography was not found, except for a slight decrease of zonal wind speed, which was connected to the volcanic area of Imdr Regio.


Introduction
The thick cloud layer of Venus rotates around the planet in the westward direction with a peak velocity 60 times higher than that of the planet itself at approximately 100 m s −1 for the altitudes of 65-70 km (upper cloud level) above the surface [1,2]; this is a phenomenon known as retrograde superrotation. In the altitude range from the surface up to 90 km, the atmospheric motion is predominantly zonal. The velocity decreases for altitude levels farther from the upper cloud level and for latitudes closer to the poles [3].
Thermal emissions from the hot lower atmosphere and surface are observed on the nightside in the near-infrared range in spectral "windows" between CO 2 bands. The lower clouds have nonhomogeneous distribution of the cloud opacity, which leads to a complex morphology of small-scale features of various observed contrasts [4]. By tracking the displacement of these cloud features with time, one can derive horizontal wind speed at the altitude of remote sensing [5].
In the spectral window of 1.74 µm on the nightside, the thermal radiance comes from the lower clouds and from beneath the main cloud level. The effective level of emissions was estimated at the altitudes of 44-48 km [5], which we will use in this paper. Thus, motion of the observed cloud features is associated with the atmosphere dynamics at these altitudes.
Results of two cloud tracking analyses of the partial dataset of VIRTIS-M (Venus Express) images, acquired in 2006-2008 [5,11] during 18 and 45 orbits, respectively, showed a consistent zonal wind speed of~60 m s −1 between 0 and 70 • S. Those works also analyzed meridional wind speed values, local time (LT), latitude and time dependence of the wind velocity in the lower clouds. Improvements brought by our analysis of the full VIRTIS-M dataset will be shown in the following sections.
More recent lower cloud tracking results came from the IR2 camera onboard JAXA's Akatsuki spacecraft in 2016 [12,13]. Although mean meridional speed values are consistent with the VIRTIS-M results, zonal speed shows a different shape of the latitudinal profile, an "equatorial jet", with zonal wind speed values that reached~75 m s −1 near the equator.

Materials and Methods
The infrared channel of the Visible and InfraRed Thermal Imaging Spectrometer (VIRTIS-M) spectrometer onboard the Venus Express covered a wavelength range between 1 and 5 µm with a spectral resolution of 14 nm and obtained images of Venus from a distance range of 40,000-70,000 km with a pixel size of 10-40 km [14]. Images were taken on the approach trajectory, when the spacecraft was moving along the high elliptical orbit towards the pericenter (located above the Northern Pole) [15]. Thus, the field of view of VIRTIS primarily captured Southern Hemisphere.
In this work, we analyzed the entire VIRTIS-M dataset, from which 195 orbits were found to be suitable (see details on analysis methods further) and from which 44,480 horizontal velocity vectors were extracted in turn ( Table 1). The number of analyzed orbits amounts to about 4 times the number used in previous research [11]. It allowed us to decrease the errors of the mean speed values, to provide more coverage in the equatorial latitudes and more spacious longitudinal coverage. The database of retrieved vectors can be found in the Supplementary Materials. VIRTIS-M 1.74 µm data cover mostly (99.8%) Southern Hemisphere. Figure 1 shows how the retrieved velocity vectors are distributed across the nightside. Cloud features are most pronounced between 10 and 40 • S, and as a result, this latitude range contains more than a half-51.2% of all vectors. Vectors are distributed almost symmetrically related to the midnight meridian: 45.9% are located before, 54.1%-after 0 h local time (LT). The number of retrieved vectors decreases towards both terminators. Figure 2 shows distribution of retrieved velocity vectors in the longitude-latitude coordinates. Longitudes between 50-150 • only contain 5.8% of all data. In the Southern Hemisphere, this longitude range includes Aphrodite Terra, a massive geological formation that was previously suggested to influence the atmosphere at the cloud top level on the dayside [16][17][18]. Our research took a step to investigate a possible influence of Venus topography on dynamics in the lower clouds. Unfortunately, data in the aforementioned longitude range are scarce, however some works indicated that other topographical features may also influence the atmosphere-in the upper mesosphere [19] and at the cloud top [20].  Figure 2 shows distribution of retrieved velocity vectors in the longitude-latitude coordinates. Longitudes between 50-150° only contain 5.8% of all data. In the Southern Hemisphere, this longitude range includes Aphrodite Terra, a massive geological formation that was previously suggested to influence the atmosphere at the cloud top level on the dayside [16][17][18]. Our research took a step to investigate a possible influence of Venus topography on dynamics in the lower clouds. Unfortunately, data in the aforementioned longitude range are scarce, however some works indicated that other topographical features may also influence the atmosphere-in the upper mesosphere [19] and at the cloud top [20].    Figure 2 shows distribution of retrieved velocity vectors in the longitude-latitude coordinates. Longitudes between 50-150° only contain 5.8% of all data. In the Southern Hemisphere, this longitude range includes Aphrodite Terra, a massive geological formation that was previously suggested to influence the atmosphere at the cloud top level on the dayside [16][17][18]. Our research took a step to investigate a possible influence of Venus topography on dynamics in the lower clouds. Unfortunately, data in the aforementioned longitude range are scarce, however some works indicated that other topographical features may also influence the atmosphere-in the upper mesosphere [19] and at the cloud top [20].  VIRTIS-M images in the 1.74 µm channel are characterized by complex cloud morphology ( Figure 3), caused by varying opacity, and the abundance of cloud features that are stretched along the zonal direction due to the superrotation. Subsequent images of the same cloud area contain stable cloud features, and tracking of their displacement allows us to retrieve velocity vectors. In this work, we used images separated by 3600-7200 s intervals.
To retrieve the velocity of the displacement of identified cloud features, we used the visual method [21,22]. For a selected pair of images, acquired with a ∆t interval, coordinates of any identified detail were precise enough to calculate zonal (u) and meridional (v) components of the horizontal wind by using the following formulas: where λ 1,2 and ϕ 1,2 are, respectively, longitudes and latitudes of the tracked cloud features in the first and the second images; R = 6051.8 km is the radius of Venus; h = 46 km is the assumed mean altitude of the features; ∆t is the time interval between two images.
To retrieve the velocity of the displacement of identified cloud features, we used the visual method [21,22]. For a selected pair of images, acquired with a Δt interval, coordinates of any identified detail were precise enough to calculate zonal (u) and meridional (v) components of the horizontal wind by using the following formulas: where λ1,2 and φ1,2 are, respectively, longitudes and latitudes of the tracked cloud features in the first and the second images; R = 6051.8 km is the radius of Venus; h = 46 km is the assumed mean altitude of the features; Δt is the time interval between two images.  VIRTIS-M nightside images in 1.74 µm have a low contrast value. To increase contrast of the cloud features, we used 2D-wavelet filtering with the Registax software (http: //www.astronomie.be/registax/) [18,23,24]. An example of this filter application is shown in Figure 4, increasing the Michelson contrast value of the image from 0.75 to 0.99. This technique is especially efficient for cloud details in middle latitudes. It allows us to increase the number of identified cloud features per image pair and provide more uniform spatial coverage through retrieved velocity vectors. In the equatorial region at 0-40 • S latitude, cloud details are small and have high contrast, and therefore can be easily identified (markers 1-4 in Figure 4a). In middle latitudes, cloud details become low-contrast and stretched along the latitudes. Their identification improves after the 2D-wavelet filtering (markers 5, 6 in Figure 4b).
form spatial coverage through retrieved velocity vectors. In the equatorial region at 0-40° S latitude, cloud details are small and have high contrast, and therefore can be easily identified (markers 1-4 in Figure 4a). In middle latitudes, cloud details become low-contrast and stretched along the latitudes. Their identification improves after the 2D-wavelet filtering (markers 5, 6 in Figure 4b). A single measurement error can be attributed to each retrieved velocity vector. It is a positioning error of the coordinates of the identified cloud details, and it is dependent on the pixel size. The errors vary from 5-10 m s −1 closer to the center of the planetary disk (usually in near polar latitudes) to 30-40 m s −1 near the limb (usually in near equatorial latitudes), as was estimated for the VIRTIS-M images previously [19]. Due to high errors and potential cloud feature distortions, the 5-10° area near the limb, and therefore was avoided during the analysis.

Results
Most of the data obtained by VIRTIS-M covers middle and high latitudes of the Southern hemisphere on the night side with the majority of the data cubes around midnight. However, we obtained a disproportionately high number of vectors at the low and middle latitudes because of finer cloud morphology, which resulted in more cloud features to track (Figures 1 and 2). The 1.74 µm VIRTIS-M data, separated in four periods of observations, are described in Table 1.
When obtaining mean values for our research, a two-step approach was made every time. First, the velocity vectors were averaged on a longitude-latitude or local time-latitude grid within each orbit. Second, orbits were averaged on the same grid. Characteristics of the analyzed mean atmospheric motion will be considered in the subsections below. A single measurement error can be attributed to each retrieved velocity vector. It is a positioning error of the coordinates of the identified cloud details, and it is dependent on the pixel size. The errors vary from 5-10 m s −1 closer to the center of the planetary disk (usually in near polar latitudes) to 30-40 m s −1 near the limb (usually in near equatorial latitudes), as was estimated for the VIRTIS-M images previously [19]. Due to high errors and potential cloud feature distortions, the 5-10 • area near the limb, and therefore was avoided during the analysis.

Results
Most of the data obtained by VIRTIS-M covers middle and high latitudes of the Southern hemisphere on the night side with the majority of the data cubes around midnight. However, we obtained a disproportionately high number of vectors at the low and middle latitudes because of finer cloud morphology, which resulted in more cloud features to track (Figures 1 and 2). The 1.74 µm VIRTIS-M data, separated in four periods of observations, are described in Table 1.
When obtaining mean values for our research, a two-step approach was made every time. First, the velocity vectors were averaged on a longitude-latitude or local time-latitude grid within each orbit. Second, orbits were averaged on the same grid. Characteristics of the analyzed mean atmospheric motion will be considered in the subsections below.

Local Time Dependence
The distributions of mean zonal and meridional wind speed in local time-latitude are given in Figure 5. The mean zonal speed in the middle latitudes reaches −60 m s −1 , and this result is in agreement with the superrotation mode and previous measurements [5,11]. There is a slight westward deceleration (1-5 m s −1 , depending on the latitude) of the zonal speed, i.e., towards the morning terminator. The nightside lowest speed (in absolute values), −58 m s −1 , is found at 4-5 h LT, whereas the highest u maxima of −64 m s −1 are located at 19.5-21 and 23-24 h LT. Figure 6 shows local time profiles of zonal speed for cloud features found in three latitude ranges. For a confidence interval of 99% standard error, the evening-to-morning deceleration in middle latitudes is statistically significant. In the equatorial region, two zonal speed minima at 21 and 4 h LT are observed. the zonal speed, i.e., towards the morning terminator. The nightside lowest speed (in absolute values), -58 m s −1 , is found at 4-5 h LT, whereas the highest u maxima of -64 m s −1 are located at 19.5-21 and 23-24 h LT. Figure 6 shows local time profiles of zonal speed for cloud features found in three latitude ranges. For a confidence interval of 99% standard error, the evening-to-morning deceleration in middle latitudes is statistically significant. In the equatorial region, two zonal speed minima at 21 and 4 h LT are observed.  The meridional speed is predominantly northward (equatorward) with the values in the range 0-3 m s −1 . There is a statistically significant westward (towards the morning terminator) acceleration across all night side, more pronounced in the middle and high latitudes than in the low (Figure 7). After the evening terminator, at 19-20 h LT in high latitudes, the mean meridional speed changes direction and becomes southward (poleward). Additionally, the meridional speed has negative values in low latitudes between 19 and 1 h LT. This area coincides with one of the zonal speed maxima. The meridional speed is predominantly northward (equatorward) with the values in the range 0-3 m s −1 . There is a statistically significant westward (towards the morning terminator) acceleration across all night side, more pronounced in the middle and high latitudes than in the low (Figure 7). After the evening terminator, at 19-20 h LT in high

Latitudinal Dependence
Mean latitudinal profiles of the zonal and meridional wind speed are shown in Figure 8. Zonal speed stays unchanging at -60 m s −1 between 20° S and 40° S. It gradually increases both northward and southward of this latitude range. To the south, between 40° S and 60° S, a jet-like increase by 2-4 m s −1 was revealed. Southward of 60° S mean zonal speed decreases to zero at the South Pole.

Latitudinal Dependence
Mean latitudinal profiles of the zonal and meridional wind speed are shown in Figure 8. Zonal speed stays unchanging at −60 m s −1 between 20 • S and 40 • S. It gradually increases both northward and southward of this latitude range. To the south, between 40 • S and 60 • S, a jet-like increase by 2-4 m s −1 was revealed. Southward of 60 • S mean zonal speed decreases to zero at the South Pole.
In most of the Southern Hemisphere, mean meridional speed has a positive sign, indicating that the flow is directed north, towards the equator. This is not the case, however, for latitudes north of 15-20 • S, as the meridional wind likely changes its sign to negative and horizontal flow changes to south-poleward (until about 5 • S). In the equatorial region, the sign of the meridional wind changes abruptly to be positive, and up to 10 • N, the horizontal flow reverts to North-poleward.
We found that the sign change of the meridional wind (from positive to negative) at 20 • S coincides with the increase of the zonal speed (2-5 m s −1 ). From the fact that we detect two flows moving towards each other, we concluded that we could observe levels at different altitudes north and south of 20 • S. Despite sparse data, these variations are significant enough to be contained within a 99% confidence interval. In most of the Southern Hemisphere, mean meridional speed has a positive sign, indicating that the flow is directed north, towards the equator. This is not the case, however, for latitudes north of 15-20° S, as the meridional wind likely changes its sign to negative and horizontal flow changes to south-poleward (until about 5° S). In the equatorial region, the sign of the meridional wind changes abruptly to be positive, and up to 10° N, the horizontal flow reverts to North-poleward.
We found that the sign change of the meridional wind (from positive to negative) at 20° S coincides with the increase of the zonal speed (2-5 m s −1 ). From the fact that we detect two flows moving towards each other, we concluded that we could observe levels at different altitudes north and south of 20° S. Despite sparse data, these variations are significant enough to be contained within a 99% confidence interval.

Longitudinal Dependence
Recent works have shown how Venus topography can influence atmospheric motion on different levels, at the cloud top [16][17][18]20,25] and at 90-110 km [19]. To explain observed variations, aforementioned research suggested that stationary gravity waves emerge near the surface caused by the large highland structures, terrae and regii. Topography-induced influence on the lower cloud motion was not investigated until results of the IR2 camera onboard Akatsuki showed no influence of topography, coupled with the fact that they longitude variations cannot be separated from the local time effects [13].
For VIRTIS-M data uncoupling, the longitudinal variations from the local time variations poses the same challenge. Mean wind speed distribution (Figure 9) shows variations across the Southern Hemisphere.

Longitudinal Dependence
Recent works have shown how Venus topography can influence atmospheric motion on different levels, at the cloud top [16][17][18]20,25] and at 90-110 km [19]. To explain observed variations, aforementioned research suggested that stationary gravity waves emerge near the surface caused by the large highland structures, terrae and regii. Topography-induced influence on the lower cloud motion was not investigated until results of the IR2 camera onboard Akatsuki showed no influence of topography, coupled with the fact that they longitude variations cannot be separated from the local time effects [13]. Strong acceleration up to -70 m s −1 of the zonal wind above the plains south of Aphrodite Terra is the most prominent feature of the zonal wind distribution; however, the data in this area consists of only 10-50 vectors per a 10° × 10° bin. The slowest mean zonal wind of up to -58 m s −1 is located in the middle latitudes between 210° and 260° longitude, corresponding to the area surrounding Imdr Regio on the surface. The meridional speed in the middle latitudes is closer to zero at 220°-350° longitude, whereas at 40°-200° it can reach up to 5-6 m s −1 in the equatorward direction. Figure 10 provides a more detailed look into the behavior of zonal wind in the middle latitudes. Here we consider the changes in wind speed in the 40°-50° S longitude bin. From 180° eastward, zonal speed decelerates from 62-63 m s −1 to 57 m s −1 at 220°-

Discussion
The results of our research expand on the previous partial analysis of VIRTIS-M data. Main wind speed values correspond with [5,11] within the standard error limits. A similar east-to-west increase of the meridional speed was reported, as well as the acceleration of the zonal speed in middle latitudes. However, the expansion of analyzed data allowed us to take a more detailed look at the lower cloud motion, discerning local time, latitude and longitude-dependent trends that previously were not detected. Comparison with NIMS and ground-based observations [6][7][8][9][10] also shows agreement of the mean values within their error limits.
The meridional circulation is organized as Hadley cells. A theory suggests the existence of several direct and backward cells [2]; however, experimentally only one direct cell was confirmed, in the upper clouds. General circulation models (GCM) have shown that several meridional circulation cells can exist in the lower clouds [26]. Detailed observations from UV VMC, Pioneer Venus and ground-based observations showed the South pole directed branch at the upper boundary of clouds (68-72 km), and the backward branch directed to the equator from IR VMC near the upper boundary of middle clouds [18,27,28].
As was discussed (in Section 3.2), the fact that we see two flows moving towards each other northward and southward from 20° S, signifies that we observe levels at different altitudes in the clouds. Additional insight can be achieved by analyzing the change in zonal wind speed when the meridional speed changes direction. Based on the latitudinal profile (Figure 8a) we selected two neighboring latitude ranges, where the mean meridional speed is equatorward: 20° S-40° S and poleward: 3° S-20° S. We then calculated the mean zonal speed for both ranges, which turned out to be −60.3 and −61.8 m s −1 , respectively, with the standard error of each value not exceeding 0.9 m s −1 . Taking into account that the estimated vertical shear in the lower clouds is <1 m s −1 km −1 [5], we suggest that the layer in the equatorial region is located about 2 km higher (48-50 km) than the layer with the positive meridional speed (44-46 km). The higher level has a negative meridional speed and horizontal flow there moves towards the South Pole, whereas in the latitudes of 20° S-40° S, at the lower level, the meridional speed is positive and the horizontal flow moves towards the equator. We assume that we observed two fragments of the lower cloud direct Hadley cell: the South pole directed branch is located at 48-50

Discussion
The results of our research expand on the previous partial analysis of VIRTIS-M data. Main wind speed values correspond with [5,11] within the standard error limits. A similar east-to-west increase of the meridional speed was reported, as well as the acceleration of the zonal speed in middle latitudes. However, the expansion of analyzed data allowed us to take a more detailed look at the lower cloud motion, discerning local time, latitude and longitude-dependent trends that previously were not detected. Comparison with NIMS and ground-based observations [6][7][8][9][10] also shows agreement of the mean values within their error limits.
The meridional circulation is organized as Hadley cells. A theory suggests the existence of several direct and backward cells [2]; however, experimentally only one direct cell was confirmed, in the upper clouds. General circulation models (GCM) have shown that several meridional circulation cells can exist in the lower clouds [26]. Detailed observations from UV VMC, Pioneer Venus and ground-based observations showed the South pole directed branch at the upper boundary of clouds (68-72 km), and the backward branch directed to the equator from IR VMC near the upper boundary of middle clouds [18,27,28].
As was discussed (in Section 3.2), the fact that we see two flows moving towards each other northward and southward from 20 • S, signifies that we observe levels at different altitudes in the clouds. Additional insight can be achieved by analyzing the change in zonal wind speed when the meridional speed changes direction. Based on the latitudinal profile (Figure 8a) we selected two neighboring latitude ranges, where the mean meridional speed is equatorward: 20 • S-40 • S and poleward: 3 • S-20 • S. We then calculated the mean zonal speed for both ranges, which turned out to be −60.3 and −61.8 m s −1 , respectively, with the standard error of each value not exceeding 0.9 m s −1 . Taking into account that the estimated vertical shear in the lower clouds is <1 m s −1 km −1 [5], we suggest that the layer in the equatorial region is located about 2 km higher (48-50 km) than the layer with the positive meridional speed (44-46 km). The higher level has a negative meridional speed and horizontal flow there moves towards the South Pole, whereas in the latitudes of 20 • S-40 • S, at the lower level, the meridional speed is positive and the horizontal flow moves towards the equator. We assume that we observed two fragments of the lower cloud direct Hadley cell: the South pole directed branch is located at 48-50 km, and the equatorward return branch is at 44-46 km. Thus, these two fragments can be considered as the fragments of the second direct Hadley cell in the lower clouds observed at 1.74 µm. This proposition is also supported by the existence of the mid-latitude jet at 40-60 • S (Figure 8a), which is most likely connected to the lower cloud direct Hadley cell.
We compared our results with previous cloud tracking in other wavelengths on the dayside from UV (VMC/VEX [22] and two Pioneer Venus profiles for two different years of observation [29]) and IR (VMC/VEX [18]) ( Figure 11). Blue arrows in Figure 11a  km, and the equatorward return branch is at 44-46 km. Thus, these two fragments can be considered as the fragments of the second direct Hadley cell in the lower clouds observed at 1.74 µm. This proposition is also supported by the existence of the mid-latitude jet at 40-60° S (Figure 8a), which is most likely connected to the lower cloud direct Hadley cell. We compared our results with previous cloud tracking in other wavelengths on the dayside from UV (VMC/VEX [22] and two Pioneer Venus profiles for two different years of observation [29]) and IR (VMC/VEX [18]) ( Figure 11). Blue arrows in Figure 11a   (a) Mean direction of the meridional wind at different altitudes (arrows) overlaid on the thermal structure of the atmosphere in latitude-altitude coordinates from combined Venus Express and Akatsuki radio occultation experiment results, averaged over both hemispheres [30]. We changed the latitude sign to negative to adopt the figure to the Southern hemisphere. Numbers at the contours indicate the temperature in Kelvin. Latitudes represent both hemispheres assuming north-south symmetry. The thermal structure is used for illustration to designate different levels in the Southern Hemisphere. Blue arrows represent UV cloud tracking results (PV and VMC/VEX), red arrow: IR cloud tracking from VMC/VEX, green arrows:  (a) Mean direction of the meridional wind at different altitudes (arrows) overlaid on the thermal structure of the atmosphere in latitude-altitude coordinates from combined Venus Express and Akatsuki radio occultation experiment results, averaged over both hemispheres [30]. We changed the latitude sign to negative to adopt the figure to the Southern hemisphere. Numbers at the contours indicate the temperature in Kelvin. Latitudes represent both hemispheres assuming north-south symmetry. The thermal structure is used for illustration to designate different levels in the Southern Hemisphere. Observed deceleration of the zonal speed from evening to morning terminator is the most notable feature of the local time-dependent variations. A similar deceleration has been suggested for the middle clouds based on near-IR cloud tracking [18], although it was less pronounced than in our case at only up to 2 m s −1 , close to the error margins.
Emergence of the mid-latitude jet observed at the cloud top can be connected to a Hadley circulation cell [31]. In the lower clouds the jet found between 40 • S and 60 • S can as well be driven by Hadley circulation. This can also indirectly indicate that we observe fragments of a direct Hadley cell. Therefore, features described above repeat those of the upper cloud cell.
Discovered longitudinal variations are in many instances comparable to the diurnal variations (up to ±5 m s −1 ). Moreover, searching for the influence of topography is challenging because VIRTIS nadir data do not cover any prominent highlands, such as Aphrodite Terra, instead covering several smaller ones, the tallest of which is Lada Terra (tallest areas are up to~3 km peak altitude, located at 65-75 • S/320-30 • E). The notable deceleration of zonal speed, roughly in the coordinates of Imdr Regio, is a volcanic area with young lava flows and excess heat radiation [32], which as shown in Section 3.3, is stronger than any possible local time bias; however, other highlands such as Lada Terra (center at 60 • S/20 • E), Themis Regio (center at 40 • S/280 • E) and Alpha Regio (center at 22 • S/5 • E) do not show any coherent signs of influence on the circulation.
We can compare our results to the general circulation models. The AFES-Venus model [33] produced the same deceleration of the zonal wind from evening to morning terminator at the altitude of 50 km, as in our results ( Figure 5). The local time dependency for the meridional speed also shows similarities, as the poleward motion is found before midnight, and equatorward-after midnight. Jet-like increases of zonal speed at the lower cloud level near equator and in middle latitudes has been shown by the T63 GCM [34]. Longitudinal variations have also been simulated by the GCM at~54 km [35]. Although the general behavior of wind velocity does not coincide with our results (Figure 9), model meridional speed was found to show more variability compared to the zonal speed.

Conclusions
In this work, the full set of VIRTIS data was processed and analyzed. It allowed us to study four times more data than in the previous works and to obtain some new findings.
An experimental indication for the existence of the lower cloud direct Hadley cell was found for the first time. Although GCMs allow the existence of one or more cells in the lower atmosphere, they were not experimentally found before.
Dynamics of the Venusian atmosphere is very variable (if we compare VIRTIS and Akatsuki measurements), and we were lucky that during VIRTIS observations, a lower cloud layer was arranged in such a way that in the same wavelength (1.74 µm), we could observe the levels at different altitudes in the atmosphere depending on latitude. Parallel with the change in the value of the zonal speed, we also observed the change of the sign of the meridional speed, which leads to a conclusion that there most likely are different circulation cells in the considered altitude range of the atmosphere.
It is known that in the upper cloud meridional circulation is organized as a direct Hadley cell, with the upper branch directed to the pole at the altitude of 68-70 km and the lower one directed to the equator more than 10 km below. The mid-latitude jet is considered as being connected to the Hadley cell.
Meridional flow directions of the upper cloud Hadley cell branches and the observed lower cloud cell branches show a correlation. The mid-latitude jet observed at 40-60 • S in the lower clouds, most likely connected to the lower cloud cell, also resembles the jet in the upper cloud layer. Thus, by analogy with the direct Hadley cell in the upper cloud layer, we came to a conclusion that in 1.74 µm, we observe fragments of the direct Hadley cell in the lower cloud layer.