Venus’ Cloud-Tracked Winds Using Ground- and Space-Based Observations with TNG/NICS and VEx/VIRTIS †

: Characterizing the wind speeds of Venus and their variability at multiple vertical levels is essential for a better understanding of the atmospheric superrotation, constraining the role of large-scale planetary waves in the maintenance of this superrotation, and in studying how the wind ﬁeld affects clouds’ distribution. Here, we present cloud-tracked wind results of the Venus nightside, obtained with unprecedented quality using ground-based observations during July 2012 with the near-infrared camera and spectrograph (NICS) of the Telescopio Nazionale Galileo (TNG) in La Palma. These observations were performed during 3 consecutive days for periods of 2.5 h starting just before dawn, sensing the nightside lower clouds of Venus close to 48 km of altitude with images taken at continuum K ﬁlter at 2.28 µ m. Our observations cover a period of time when ESA’s Venus Express was not able to observe these deeper clouds of Venus due to a failure in the infrared channel of its imaging spectrometer, VIRTIS-M, and the dates were chosen to coordinate these ground-based observations with Venus Express’ observations of the dayside cloud tops (at about 70 km) with images at 380 nm acquired with the imaging spectrometer VIRTIS-M. Thanks to the quality and spatial resolution of TNG/NICS images and the use of an accurate technique of template matching to perform cloud tracking, we present the most detailed and complete proﬁle of wind speeds ever performed using ground-based observations of Venus. The vertical shear of the wind was also obtained for the ﬁrst time, obtained by the combination of ground-based and space-based observations, during the Venus Express mission since the year 2008, when the infrared channel of VIRTIS-M stopped working. Our observations exhibit day-to-day changes in the nightside lower clouds, the probable manifestation of the cloud discontinuity, no relevant variations in the zonal winds, and an accurate characterization of their decay towards the poles, along with the meridional circulation. Finally, we also present the latitudinal proﬁles of zonal winds, meridional winds, and vertical shear of the zonal wind between the upper clouds’ top and lower clouds, conﬁrming previous ﬁndings by Venus Express.


Introduction
The atmosphere of Venus is in superrotation, a state in which its averaged angular momentum is much greater than that corresponding to corotation with the surface. The circulation up to the cloud tops is characterised by an increasing zonal retrograde wind (in the east-west direction). The retrograde zonal superrotation (RZS) wind starts to build up From 2006, the southern hemisphere lower cloud deck was thoroughly studied with Venus Express (VEx) observations and the related wind field was explored with unprecedented precision [5,20,21].
Recently, with the advent of the Akatsuki space probe and its equatorial orbit, it was possible to analyse with great detail the clouds and dynamics of the lower-to-middle Venus' atmosphere layers, in both hemispheres [22,23]. The accuracy of Akatsuki IR2 camera images were aligned with the long-term observations of the deep Venusian atmosphere, allowing researchers to retrieve striking new features, as this is the case in the planetaryscale, bow-shaped equatorial discontinuity [9] and the profusion of mesoscale vortices [23]. Since the IR1 and IR2 cameras stopped working in December 2016 [24], ground-based NIR observations constitute the only way of resuming the outstanding Akatsuki's infrared atmospheric research.
Recent developments in ground-based techniques in enhancing the precision of the images obtained and improving accuracy in the images' navigation (i.e., extreme precision in the attribution of latitude/longitude coordinates in the observed images from ground telescopes) allowed researchers to improve the contribution of ground-based NIR observations on Venus' nightside lower cloud deck investigation. In fact, the use of an improved cloud tracking method, based on phase correlation between images [22], and the use of NASA SPICE kernels [25], combined with a more accurate adjustment of the planetary limb to the coordinates grid, increased the precision in the identification of telescope observation images' pixel-related planetary coordinates.
The study of the lower cloud wind field (∼48-55 km), on both Venus' hemispheres, could provide crucial hints in understanding the source of the atmospheric angular momentum excess that leads to the superrotation phenomena observed on Venus' atmosphere. We note that, in the period between the failure of the visible and infrared thermal imaging spectrometer (VIRTIS) on board the Venus Express [26] infrared channel and the Akatsuki Venus' orbital insertion, and after the Akatsuki/IR2 camera stopped working in 2016 [24], it was not possible to sense and study Venus' dynamics at the bottom of the cloud level using space-based observations. We also note that the ground-based NIR dataset presented in this work is a unique way to study the Venus' atmosphere dynamics at the cloud bottom layer for 2012.
Taking in account the relevance of presenting results from coordinated observations, coming from different instruments and wavelength ranges, we took advantage of a coordinated observational campaign between space-based VEx and ground-based observations with TNG/NICS. This multi-instrument and multi-technique observation method has largely proved its efficiency in previous studies [10,11,[27][28][29].
In this work, we present nightside wind velocity results, these results were retrieved from ground-based TNG/NICS observations (see Figure 1) in the 2.28 µm NIR windows sounding the bottom of the cloud level (∼48 km). Our observational efforts were coordinated with VEx/VIRTIS-M cloud tops dayside observations; therefore, we will also present simultaneous dayside winds measured at the cloud top level (∼70 km) and, finally, we will present a first approach to the vertical wind shear between the altitudes, sensed in the context of the present work. We will also compare the obtained results with GCM modelling predictions (LMD). The fact that the zonal wind at the bottom of the cloud layer is highly stable [30] allowed the researchers to present, here, an estimate of the vertical shear of the wind between the two altitudes sensed in the context of this study (the lower and top cloud layers). Venus' lower cloud deck as seen from the Earth's point of view with the Telescopio Nazionale Galileo (TNG). These images were obtained in the near-infrared Earth transparency window at 2.28 µm with the camera NICS after daylight correction. Given the diversity of cloud's morphological features and their high variability, we note that these images were obtained over three consecutive days.

Observations and Image Processing
In this work, we present results based on observations carried out with the nearinfrared camera and spectrograph (NICS) of the Telescopio Nazionale Galileo (TNG), in La Palma, between 10 and 12 July 2012 ( Figure 1). We observed for periods of 2.5 h starting just before dawn, for 3 consecutive nights ( Figure 2 shows the observation's geometry for each one of the observing days). We acquired a set of images of the nightside of Venus with the continuum K filter at 2.28 µm (k cont ), which allowed us to monitor motions at the lower cloud level of the atmosphere of Venus, close to 48 km altitude. We acquired a series of short exposures of the Venus disk (exposure time of 3 s). These observations were part of the network of ground-based observations of Venus coordinated with ESA's Venus Express orbiter for the 2012 Venus transit campaign [31].

Ground-Based Observations with TNG/NICS
Venus' apparent diameter at observational dates was greater than 36 arcsec, allowing a high spatial resolution. Even though the pixel scale of the NICS narrow field camera is 0.13 arcsec, the atmospheric seeing during our observations (ranging 0.8-1.0 arcsec) limited the effective spatial resolution of our images. Super-resolution imaging for Venus can be attained using the lucky imaging technique [32][33][34], and, in previous works, the lucky imaging technique allowed researchers to obtain stacked images of Venus that improved, by a factor of 3, the spatial resolution imposed by the atmospheric seeing [29]. In this work, we selected the best-aligned frames where turbulence was not present, stacked those frames, and obtained an ultra-high-resolution image below the typical seeing, though it was limited by the pixel size of the NICS camera. Lucky imaging has been proven to improve the spatial resolution under a varied set of seeing conditions [35], and although resolutions as small as 0.1 have been reported in some cases [36], we assessed that a more conservative resolution of 0.2 arcsec per pixel matched the resolution of our stacked images better. Note that this value is typical of the resolutions achieved at very high altitude observatories or with adaptive optics [37].

Figure 2. Left:
we can see the geometry of the ground-based observations. It also shows the apparent size of Venus on 11th of June 2012, as it was seen at the Telescopio Nazionale Galileo (TNG) site at the La Palma's observatory. Right: example of observations of the nightside of Venus where our use of SPICE kernels provided accurate positional referencing of the images. In the observations that we performed, we note that by using SPICE kernels for obtaining an accurate positional reference, we were able to very accurately assign to each pixel in the image their latitude and longitude data (besides the presented grid).
The absolute spatial resolution, averaged across the 3 days of observations, on the disk was ∼106 km/px at disk centre. Since the observing nights, we obtained two stacked images, obtained from sets of images separated by a time interval of about 60 min, assuming wind speeds of the order of 60 m/s [22] we expect that clouds would move approximately 216 km (above 2 pixels) between the pair of images.
In order to correct for scattered light from the (saturated) dayside crescent into the nightside, a set of observations with a Bracket gamma (Brγ) filter, centred at 2.169 µm, were performed ( Figure 3). Cloud patterns are not observable with this filter due to the high optical depth of the gaseous CO 2 component. Subtracting k cont and Brγ images allowed us to efficiently correct the pollution due to this scattered light from the saturated dayside [18]. . From left to right, a continuum K image, an image obtained with a Bracket gamma filter, and the result of subtracting a Brγ image from the one taken with k cont . On these images, the dark regions are clouds, the bright regions are optically thinner areas between the clouds that allow thermal emission from the lower atmosphere to escape, and the outlined crescent is the saturated dayside of the planet. Images from 11 July 2012.
The correction of the scattered light with the Brγ images also helped to improve the visualization of the planetary limb and, hence, to obtain a more accurate position of the coordinates' grid. Besides, co-adding the best images and cross-correlating regions of clouds ( Figure 4) helped us to obtain an effective resolution significantly by overcoming the seeing-limited resolution. . Images A and B are produced from a selection of the best 10% taken, which were registered to a common coordinate system and finally co-added to form the present images. The images A and B were taken with approximately one hour interval between them. In this way, we can apply the cloud tracking technique in order to retrieve the wind field at the bottom layer of Venus' cloud deck (∼48 km). A and B images at the left were obtained in the first day of observations, the pair in the middle were taken at the second day and, finally, the A and B images on the right side of the figure were obtained at the last observing day.
In order to prevent image saturation, even on the nightside, when using the k cont filter, we used a Gray 5 attenuator filter (which reduces the flux between 4.8 and 5.7 magnitudes in the wavelength range from 1.03 to 2.15 µm with a 0.3 magnitude accuracy).
On the other hand, since we started observing Venus still with low elevation (airmass ∼3), we were obliged to increase the exposure time at the beginning of each day, observing the run. Accordingly, the used exposure time was 10 s until the airmass reached approximately 2.5. With this caution, we could obtain-even at the beginning of each day's run-images with enough contrast which allowed us to perform cloud tracking. With the steady decrease in the airmass, we were forced to use a 3 s exposure time to prevent nightside image saturation and maintain the optimal contrast at Venus' lower cloud deck. Orbital geometry and observational characteristics are displayed on Table 1. Table 1. Orbital geometry and circumstances of ground-based observations : (1-2) Date/UT interval; (3)(4)(5) disk aspect; (6) sub-observer longitude and latitude (planetocentric); (7-9) observing conditions and geometry. The image pixel scale in km (9) was computed at disk centre and using NICS camera resolution, at the TNG telescope, of 0.13 arcsec/pix.

Coordinated Space-Based Observations with VEx/VIRTIS-M
VEx/VIRTIS-M is a 2-channel imaging spectrometer operating in the visible and infrared (0.3-1 µm and 1-5 µm) wavelength ranges [26]. The visible channel maps the dayside at ultraviolet wavelengths of about 380 nm, probing UV features at an altitude of 66-73 km. We selected pairs of VIRTIS-M ultraviolet (UV) images at 380 nm to retrieve wind speeds from the motions of the clouds on the dayside of Venus. To improve the measurement accuracy, the time interval between images was maximized (typically ∼1 h). The motions of the UV upper cloud features were visually measured on equirectangular projections of these VIRTIS-M images, which were previously processed to enhance the visualization of the cloud patterns following the same procedure as in previous works [21,38].
Coordinated VEx observations included in this work comprised a set of UV dayside images acquired during VEx orbit numbers from 2272 to 2273, during the days 11-12 of July 2012, respectively. Since the quality of these VEx/VIRTIS-M images was variable, we discarded the images with lower signal-noise ratio (SNR) for our purpose of measuring wind speeds by manually tracking cloud tracers. The selected pairs of images were separated by a time interval of around 60 min, encompassing a latitude range from the equator to 80 • S. Planetary local hours spanned from 10 h to 17 h (see Table 2). These observations were chose to match the same observing dates as our ground-based observations with TNG/NICS (see Table 2).
The original VIRTIS-M images were navigated and processed to improve the SNR ratio using the convolution with the directional filter described by Hueso et al. [38]. As the position of the spacecraft changes along with the instant where each image was acquired, each couple of images were cylindrically or polar projected depending on whether they were sensing regions closer to the equator or to the South Pole. Figure 5 shows 2 examples of VIRTIS-M images at 380 nm after being processed and geometrically projected. In all the cases, a grid with a spatial resolution of 0.2 • was employed.

Wind Determination Methods
Tracking the motions of the clouds in images has long been regarded as a trustworthy method for estimating the wind speeds of Venus within the cloud layer [3], and the results have been repeatedly validated by comparison with in situ measurements by entry probes and balloons [40]. For winds at the cloud tops (∼70 km above the surface), these are usually evaluated using UV images of the dayside of Venus [5,8,41,42], while for the deeper winds at the middle-to-lower clouds (approximately at 47.5-56.5 km of altitude), we use images of the dayside taken at visible and near-infrared wavelengths [42,43] and infrared images at bands 1.7 and 2.3 µm for nightside clouds, seen as opacity patterns against the background emission from the deep atmosphere of Venus [3,17,22].

Cloud-Tracked Winds on the Nightside Using TNG/NICS
In the case of the ground-based k cont images taken by TNG/NICS, we used image pairs, such as those in Figure 4, and a semi-automatic technique relying on phase correlation to perform the template matching [22]. On the other hand, the navigation of the images (i.e., attributing latitude and longitude coordinates to each image's pixel) obtained with ground-based telescope benefit from the used of SPICE kernels.
Our observing strategy was similar to the one used previously by Tavenner et al. [18] at IRTF (Figure 3). SPICE Kernels: Enhancing TNG/NICS Image's Navigation Accuracy SPICE is a system developed by NASA's Navigation and Ancillary Information Facility (NAIF) to assist scientists in planning and interpreting scientific observations from spacebased instruments onboard planetary spacecrafts, such as Venus Express or Galileo, and it is also used in engineering tasks associated with space missions. This system's main focus is on solar system geometry, concerning the ephemeris of objects in the solar system; the orientations of both the object and the observer; the sizes, shapes and time conversions between systems; and the lapse rates [25,44].
To measure wind velocities by tracking the displacement of contrasting features on the observed atmospheric layer of a planet, a high level of accuracy on the geographical position of such features is paramount. With SPICE, this accuracy is achieved using a wealth of available data on the ephemeris of a great number of objects in the solar system. Ephemeris data on the trajectory of the different objects are reached by carefully calculating the light travel time between observer and target on a designated reference frame, contained in data available on hundreds of kernels (SPICE data files). The precision of such computations can reach an accuracy in the order of 4.0 × 10 −11 s [25], which in turn offer a position error of 1.2 cm at 50 A.U. Hence, the position error from SPICE image navigation on these observations will be less than a centimetre. Since SPICE is openly available to the community, its contents are also subject to a high level of scrutiny, and every piece of data is always updated with the latest results that improve these archives, which are reviewed periodically.
Using the same software as Peralta et al. [22], the images of the nightside of Venus taken by TNG/NICS were navigated with SPICE, and the position and orientation of the computed grid of coordinates (latitude, longitude, and local time) was adjusted visually to match the disk of Venus in the images using a limb-fitting technique, obtaining a precision of about 1 pixel for the position of the grid. The size of the grid was determined using the angular size of Venus predicted by SPICE and the pixel scale of the camera NICS (0.13 arcsec per pixel).

Cloud Tracking (CT) Wind Retrieval with VEx/VIRTIS-M
The long integration time needed to obtain a single VIRTIS data cube in the visible channel (∼12 min), along with the highly eccentric polar orbit of VEx with its apocenter located near south pole, limits VIRTIS-M observations (and cloud tracking) to the southern hemisphere. The spatial resolution in the images of these cubes vary between 45 km per pixel for high latitudes and about 15 km per pixel for lower latitudes. As a first step, the UV images from VIRTIS-M were processed and then projected onto equirectangular or polar coordinates (depending on the latitudinal coverage) with an angular resolution compatible with the worst spatial resolution for the images in the pair. The apparent motions of features observed in the images were measured manually, while pairs of images with a time separation of about 1 h were chosen. This time interval arose from the compromise and maximizing the precision of the speed measurements and ensuring that the shape of most of the cloud tracers remained coherent and clear for identification.

Cloud Tracking Method, Cloud Features Tracers, and Accuracy
We used the same methodology as in Machado et al. [10,11], Gonçalves et al. [12], Peralta et al. [22] to retrieve the zonal and meridional winds in the VEx/VIRTIS-M and TNG/NICS images. We obtained cloud-tracked winds from VEx/VIRTIS-M following the temporal apparent motions of UV upper cloud features on Venus' dayside (Venus/VIRTIS-M). A total of 72 cloud tracers from 2 pairs of dayside images, separated by a time interval of about 1 h (see Table 2), were identified to carry out the calculation of the wind velocities.
Concerning the near-infrared (1.74, 2.28 µm), we used the same cloud tracking technique as for the VIRTIS-M UV (380 nm) images; although, the cloud patterns to be tracked on the nightside images keep some important differences. As mentioned before, on the dayside hemisphere, at visible/UV wavelengths, features can be distinguished as moving cloud patterns in the upper clouds, particularly thanks to the unknown UV absorber, which causes high-contrast patterns to form as solar radiation is absorbed and reflected back to the observer (VEx). On the nightside images at near-infrared wavelengths, the moving features we observe are formed by opacity patterns against the bright thermal emission from the deeper and warmer atmosphere, below the cloud layer.
To perform this method in the different dynamical regimes on both cloud layers, other factors should be considered as well. Previous nightside cloud tracking winds using space-based infrared images [5,21,22] have stated that the lower clouds at the bottom cloud deck, in contrast to the upper clouds, stay coherent for longer timescales (more than 2 h) and generally suffer less variability in the zonal wind speed [22]. For these reasons, the approximately one hour time interval (between images A and B) was perfectly suitable for our cloud tracking purposes when applied to TNG/NICS observations of the bottom cloud layer.
The following expressions were used for the calculation of the wind components: where u represents the zonal wind, v is the meridional wind component, a is the Venus radius, H is the height above the surface, φ and λ are the longitude and latitude in degrees, and ∆t is the time difference between the images (in seconds). Here, zonal wind means a horizontal (velocity vector along the same altitude) system of winds parallel to equator, while the meridional component represents winds in the north-south direction. The errors associated with both components of the cloud-tracked winds were calculated in the same way as in previous works [5,[10][11][12]28]. Considering that the error in time δt is quite small, the absolute errors for both components of the winds will be given by the following expressions, taken from the general forms described in Bevington et al. [45], i.e., δu ≈ δX/∆t and δvs. ≈ δY/∆t, where δX and δY are absolute errors for the spatial displacement of the clouds. Concerning VEx/VIRTIS, the grids used for cloud tracking had a spatial resolution of 0.2 • (for both latitudes and longitudes). Hence, δX and δY were about 21 km, implying wind velocity measurements errors of ∼5 m·s −1 . Regarding ground-based TNG/NICS images, the grids used for cloud tracking purposes had a spatial resolution of 0.4 • , and the retrieved zonal and meridional winds were affected by an error of about 30 m·s −1 .

Winds at the Dayside Cloud Tops with VEx/VIRTIS-M
To perform cloud tracking with the VIRTIS-M images, we chose pairs of images separated by 60 min for the case of VEx orbit 2273, and about 45-60 min for orbit 2272. Sometimes we can observe noticeable differences between cloud patterns, these transiting from "patchy" to more "plain", from lower to higher latitudes (see Figure 5). Figure 6 shows our individual wind measurements along with the latitudinal profiles of zonally averaged speeds from Sanchez-Lavega et al. [5]. The results for the zonal component of the wind are displayed in the left graph, while the meridional component is shown on the right graph. The reference profiles from Sanchez-Lavega et al. [5] are based on cloud-tracked winds averaged along latitude bins of 2 • , and the graphs exhibit the mean value and the corresponding standard deviation. One can note that, in both zonal and meridional latitudinal wind profiles, there is a gap in the covered latitudes (between 65 and 70 • South). This was caused by the difficulty of finding reliable cloud tracers for this latitudinal region were cloud details were not easy to discern. Concerning the zonal wind (left side), we can note the presence of a nearly uniform zonal wind between the equator and about 45 • S, where the zonal velocity of the wind is approximately constant with latitude. Around 50 • S, a pronounced mid-latitude jet is apparent, a phenomenon already reported by Sanchez-Lavega et al. [5]. For higher latitudes, we observe a steady decrease in the zonal wind towards the pole. With respect to the meridional wind, results in Figure 6 (right side) exhibit no clear trend, discarding-at the time of the observations-net poleward circulation at the cloud tops. The results of a latitudinal wind profile (zonal and meridional winds), retrieved from Venus' dayside space-based VEx/VIRTIS-M data at cloud tops level (around ∼70 km altitude) are in good agreement those from Sanchez-Lavega et al. [5], though the mid-latitude jet seems more enhanced than reported in other works.

Cloud Patterns and Winds at the Nightside Lower Clouds with TNG/NICS
Due to the variable opacity of the deeper clouds, the infrared thermal radiation coming from the super-heated low atmosphere of Venus can be observed as high-contrasted markings at the lower clouds' layer in the TNG/NICS images [17,18], with dark areas being related to thicker clouds (higher local opacity) and thin clouds seen as bright regions in the images since they allow the transmission of higher infrared thermal radiation fluxes coming from the hot lower atmosphere [13,18,21,46] (see Figure 7). We also studied the motions of these markings to retrieve wind velocity field on the nightside lower clouds of Venus [21,22].

Morphology of the Nightside Lower Clouds with TNG/NICS Images
During the 3 consecutive days of our TNG/NICS observations, the morphology of the nightside lower clouds exhibited a huge diversity and fast variability (see Figure 7). Especially striking was how the morphology of the clouds changed from one day to the next, and the strong hemispherical asymmetries that were previously reported by Crisp et al. [14]. In our observations, the mid-latitude region is prone to show bright bands of clouds roughly aligned with the equator and hence correlated with the zonal wind. The limited spatial resolution of ground-based images prevents the visualization of many atmospheric phenomena identified in images from space missions, such as mesoscale waves, dark spots, mushroom-shaped clouds, vortexes, or patterns resembling shear instabilities [23,[47][48][49].
Other features with larger scales, such as the sharp dark streaks, circumequatorial belts, or the bright equatorial wall, are infrequent [23], and they seemed missing during our TNG/NICS observations. The hemispherical symmetry observed for the bright midlatitude bands during our observations and their orientation northwest-to-southeast (see Figure 7) suggests the manifestation of other large patterns, such as the bright planetary-scale streaks [50]; although, the limited spatial resolution of our images inhibits us from confirming whether these patterns are composed of thinner bright streaks [50].

Was the Planetary-Scale Cloud Discontinuity Present?
Peralta et al. [9] reported the astonishing discovery of a sharp equatorial discontinuity or "disruption" propagating at the middle and lower clouds of Venus, where no planetary wave had been identified before, except for mesoscale waves [23,51]. This discontinuity manifests as a dark band with cyclical behaviour that has remained recurrent for decades, rotating around the planet faster than the mean zonal flow. The related sharp cloud discontinuity spans along the mid-latitude region, almost at a planetary scale. Peralta et al. [9] suggests that this disruption may be the physical manifestation of a new type of Kelvin wave that propagates trapped about the equator, dissipates before arriving at the upper clouds, and significantly alters the spatial distribution of clouds' aerosols [9].
In the present work, based on ground observations, the evidence of a planetary-scale homogeneous dark clouds is clear, such as an atmospheric wall (see Figure 8). This striking atmospheric feature was coherent with all the images we obtained for the second day of observations (12 July 2012). Nevertheless, it is not clear whether this event can be directly associated with a manifestation of the cloud discontinuity reported by Peralta et al. [9], since Kelvin waves are expected to propagate faster than the zonal mean flow. During 2016, Peralta et al. [9] estimated that the zonal phase speed of the discontinuity was, on average, 23 ± 9 m·s −1 faster than the mean zonal flow-a difference hard to detect given the accuracy of our wind speeds with TNG/NICS (∼30 m·s −1 ). In fact, no faster speeds were detected during 12 July 2012 (see Figure 9).

Dynamics of the Nightside Lower Clouds with TNG/NICS Images
Following the same method as Peralta et al. [29], we performed cloud tracking using phase correlation for template matching, being able to retrieve wind speeds that were comparable with some cases measurements from space orbiters (see Figure 9). A set of cloud patterns were selected attending to their shape coherence along time and they were tracked along time to characterize their motions and, therefore, the wind speeds (see in Figure 10). For each one of the three days of observations, we retrieved the velocity based on each pair of related tracers (from image A to B). We measured both horizontal components (zonal and meridional) of the wind. For the zonal component of the wind, we built a latitudinal profile of zonally averaged velocity values in latitude bins of 5 • for each day of observations (see Figure 9). The error bars associated with each velocity were derived from the calculations described in Section 3.3. These were a direct consequence of the observing conditions, mainly from the images' resolution and the time interval between the exposures that produced images A and B. They also correspond to the standard deviation, related to the set of wind velocities retrieved within each latitude bin.
In Figure 9, one can clearly see that equatorward of 60 • , the zonal wind is around 60 m·s −1 in both hemispheres, with a variability ranging 5-10 m·s −1 . Near the 60 • latitude, in both hemispheres, there are hints of a weak jet, although this cannot confirmed out of the error bars. At higher latitudes, we observed the expected steady decrease towards the poles already reported in previous works [5,21,22,52].
In relation to the meridional flow of the wind, despite covering a wide range of latitudes and local time positions, the accuracy of our velocity measurements prevents us from discerning any clear trend in the zonally averaged profile out of the error bars (see Figure 11), consistent with previous reports from VEx [21] and Akatsuki [22].
The results concerning the zonal component of the wind velocity at the lower cloud deck of Venus is undoubtedly a major result from this analysis, given past measurements with ground-based observations [17,18,22]. We can see, in Figure 12, a consolidated latitudinal profile of the zonal wind from all cloud-tracked measurements from all days of our TNG/NICS observational run. The retrieved zonal winds have the particularity that they were measured with an improved cloud tracking method that uses phase correlation between images. In Figure 12, we can note the high consistency between previous space-based (VEx/VIRTIS) zonal wind latitudinal profile from Hueso et al. [21] and Peralta et al. [22] (Akatsuki/IR2) at lower cloud deck of Venus, and the latitudinal zonal wind velocity profile coming from ground-based observations in this study (TNG/NICS); although, in this study, the presence of an equatorial jet is not evident .  , obtained using ground-based, near-infrared observations (shown as red diamonds) from the Telescope Galileo (TNG) and the near-infrared camera (NICS). The results shown here are the averaged outcome of a 5 • binning in latitude made from all our measurements and comprising all observational days (11-13 July 2012). For comparative reasons, we also present the latitudinal profile of the zonal wind obtained using space-based data from the infrared channels of VIRTIS instrument from Venus Express space probe (green in the figure, from Hueso et al. [21]) and from the IR2 camera on the Akatsuki orbiter (blue dots, Peralta et al. [22]). We note that VEx/VIRTIS results only cover the southern hemisphere of Venus due to orbital constraints that allowed the obtention of suitable images for performing cloud tracking, only in the southern hemisphere, while the Akatsuki's equatorial orbit enable the probing of both hemispheres of Venus. For an interpretation of the references to colour in this figure caption, the reader is referred to the web version of this article .

Discussion and Conclusions
The ground-based data obtained with TNG/NICS during 11-13 of July 2012 have allowed us to monitor the wind variability on the nightside lower cloud deck of Venus. The latitudinal zonal wind profile obtained from our ground-based observations are consistent with previous space-based results. Given that, from the Earth, it is possible to observe both hemispheres of Venus simultaneously, our observations helped us to extend the monitoring of deeper winds at the northern hemisphere with a quality closer to the observations by Akatsuki/IR2 [22], and drastically improve past ground-based observations [16,17], confirming zonal winds that can be regarded as hemispherically symmetric in contrast to the cloud morphology. Our measurements during 11-13 July 2012 also discard the presence of the equatorial jets reported during 2016 with Akatsuki/IR2 data [22,48].
As already stated, our measurements of the nightside meridional winds display no clear trend at the lower cloud deck, a result in agreement with previous reports from VEx [21] and Akatsuki [22]. Nevertheless, using higher-accuracy wind measurements and the larger data set from VIRTIS-M, Gorinov et al. [52] showed that meridional winds at the nightside lower clouds exhibit abrupt changes in its sign at lower latitudes, a result to be confirmed with our further studies via more accurate measurements. In any case, results of meridional wind velocities from both ground-based observations and space missions (VEx/VIRTIS and Akatsuki/IR2) [21,22,52] seem inconsistent with a Hadley cell circulation at the lower clouds.
We detected a significant variability in the opacity of lower clouds' morphology along the different days of observations (see Figures 1 and 7), a phenomenon already reported in previous works [16,22]. Even though the cloud pattern visible in the TNG/NICS images during 12 July 2012 (see Figure 8) can be attributed to the cloud discontinuity reported by Peralta et al. [9], caution must be taken since its associated zonal speeds do not seem faster than the background zonal flow above the error bars.
The dayside cloud-tracked winds obtained based on VIRTIS-M, reported in this work, are consistent with previous results (from Sanchez-Lavega et al. [5] and with local time averaged zonal wind profiles from Horinouchi et al. [8]). The zonal wind retrieved at mid-latitudes cloud tops indicates an almost homogeneous flow of the zonal wind (around 100 m·s −1 ), with a steep decrease at the latitudes where it is expected to find the so-called cold collar (nearly 60 • S). Figure 13 shows two levels of the atmosphere of Venus sounded in the framework of the coordinated observations undertaken in this project. We can see the averaged zonal velocities of dayside cloud tops ranging around 100 m·s −1 in the mid-latitudes region with a steep decrease for higher latitudes. However, the nightside NIR observations sensing the lower clouds present an almost steady zonal wind profile of the order of 60 m·s −1 , again in the mid-latitude region with a rapid decrease at higher latitudes. The results obtained are consistent with the expected increase in the zonal wind speed from lower layers of the atmosphere to the cloud top altitude.
The coordinated observations between Venus nightside lower cloud layer with TNG/ NICS (ground-based), and VEx/VIRTIS-M dayside cloud tops observations allowed us to compare the dynamics of Venus simultaneously at the two different atmospheric altitude levels of the lower clouds and top of the upper clouds. This made it possible to also have an integrated view of the atmospheric dynamics of Venus in the dayside and nightside of the planet, since in situ wind measurements by space missions indicate that the zonal wind speeds at the deeper clouds are similar on the day and night sides: Pioneer Venus Day and Night probes exhibited zonal speeds discrepancies below 7 m·s −1 within altitudes ranging 50-60 km [40], while VEGA balloons (propagating along a similar range of altitudes from the night to the dayside) generally displayed, for the zonal speeds, Doppler residuals below 4 m·s −1 , with a maximum of 16 m·s −1 [53]. This similarity between the deeper winds at the day and night sides are also supported by numerical simulations [30]. The new modelling tool, the Venus Climate Database (VCD), from 2021 [4,54], also highly suggests in its simulations and predictions that the zonal wind at the bottom cloud level on the atmosphere of Venus is stable along all planetary local times. Therefore, we can consider that the zonally averaged profile of zonal winds at the nightside lower clouds is also representative of their dayside counterpart, and our results for the wind speeds at the upper (VIRTIS-M) and lower clouds (NICS) can be combined to estimate the vertical shear of the zonal wind on the dayside of Venus. Figure 13. Comparison between two levels of the atmosphere is latitudinal profile of the zonal wind coming from coordinated ground-based TNG/NICS Venus nightside observations (NIR at 2.28 µm) sensing the bottom of the cloud layer (∼48 km of altitude) on 11-13 July 2012, and the dayside cloud tops zonal wind velocities obtained from our coordinated space-based observations with VEx/VIRTIS-M on 11 July 2012. For the ground-based case, we used a latitudinal binning of 5 degrees, while for the space-based observations, a 1 degree binning was used. For comparison purposes, we also plot reference results from VEx/VIRTIS-M [5] at cloud tops, VEx/VIRTIS-H IR channel, which senses the bottom of the clouds [21], and also from the Akatsuki IR2 camera [22].
If we consider that the lower clouds are located within altitudes ∼48-55 km [29] and the top of the upper clouds within ∼70-74 km [42], their vertical separation to be roughly estimated to be 20 ± 4 km. Figure 14 displays our estimation of the vertical shear of the zonal wind between these two atmospheric layers, compared with the same type of profile estimated with only VIRTIS-M by Hueso et al. [21]. Note that these values of the vertical shear cannot be representative for the real vertical distribution of the shear, since the largest values (higher than mean values represented here) are concentrated between the upper and bottom of the upper clouds (60-70 km) [21]. Figure 14 shows that our estimations of (∂u/∂z) 2 during 11-13 July 2012 are consistent with those obtained by Hueso et al. [21] with VEx/VIRTIS-M data during the years 2006-2008. Some discrepancy shown at about 50 • of latitude could be a consequence of the notorious mid-latitudes jet measured by us at the top of the clouds, which appears smoother on the results from Hueso et al. [21], perhaps due to the longer time interval that produced their averaged winds and to the known temporal variability of the mid-latitudes' jets, so the observed differences may be attributed to the amount of data and time coverage used in each case. Regarding the comparison with the IPSL-VGCM (LMD) modelling profile, [55,56], it is notorious that the predicted vertical wind shear profile shows consistency with the observation-based profile between mid-latitudes; however, for higher latitudes (higher than the location of the so-called cold collar), there is some lack of agreement between observations and the prediction profile-this was already pointed out by Machado et al. [27] for the case of zonal wind comparison, due to a decrease in the winds' velocity at higher latitudes than in the observations' profile. For comparison, we also display (∂u/∂z) 2 , as estimated by Hueso et al. [21], using only VIRTIS-M data during the years 2006-2008. It is worth noting that modelling simulations predict that the zonal wind at the bottom of the cloud layer is stable (both on the dayside and on the nightside) along all the planetary local times at a given latitude Garate-Lopez and Lebonnois [30], which enables, as a first approach, an estimation of a latitudinal vertical shear wind's profile on the southern hemisphere of Venus. For comparison reasons, a prediction profile from the IPSL-VGCM is also shown in light green [55,56].
Regarding the behaviour with latitude, both profiles of (∂u/∂z) 2 exhibit a gradual increase from equator to mid-latitudes, followed by a steeper decrease up to about 60 • S. The lower values for (∂u/∂z) 2 at higher latitudes are expected since the difference between zonal speeds at distinct cloud layers becomes smaller when winds decrease towards to the poles [5,21,28,42]. Furthermore, the highest value for (∂u/∂z) 2 occurs at ∼50 • S, close to the edge of the cold collar and a region where vertical instabilities cloud be favoured [28].
In this work, we used an improved cloud tracking method based on a phase correlation between images. We also improved the accuracy of the latitudes and longitudes assignment (i.e., navigation) to the tracers marked on the ground-based images in order to optimize the determination of shifts in cloud features between images. The attributed coordinates to each tracer came from the precise navigation of telescope images by the use of SPICE kernels; in this manner, each image pixel benefited from an accurate navigation process, which increased the precision of the retrieved wind velocities.
Ground-based infrared observations allied to the evolved cloud tracking method used in this work comprise an effective tool to explore and study the lower cloud deck of Venus. Ground-based observations obtained with telescopes and infrared detectors, in Earth's atmospheric infrared transparency windows, constitute a complementary and competitive technique for studying the atmosphere dynamics of the lower cloud layer on the nightside of Venus.
In order to enhance the spatial resolution and accuracy of retrieval of wind velocity, we intend to consider using deep machine learning methods applied to ground-based observations of extended sources. This method has recently been used with success in the context of solar investigations by Baso et al. [57], Xiangchun et al. [58], Dash et al. [59].
However, adaptive optics is also a prospect to consider in future observations in order to better compensate for atmospheric turbulence. A technique, also coming from solar physics, was successfully used in obtaining a homogeneous seeing improvement over a wide field by Schmidt et al. [60]. The adaptation of this recent method also opens a window to future improvements in studies of the atmosphere of Venus based upon ground-based observations. Author Contributions: Conceptualization, P.M. and J.P.; data curation, P.M., J.P., J.E.S., F.B. and R.G.; formal analysis, P.M., J.P., J.E.S. and F.B.; investigation, P.M. and J.P.; methodology, P.M. and J.P.; resources, F.B. and M.S.; software, J.P.; supervision, P.M. and J.P.; validation, P.M. and J.P.; visualization, J.E.S., F.B., M.S. and R.G.; writing-original draft, P.M. and J.P.; writing-review and editing, J.P., P.M., J.E.S., F.B. and R.G. All authors have read and agreed to the published version of the manuscript.