WiNDS: An H$\alpha$ kinematics survey of nearby spiral galaxies -- Vertical perturbations in nearby disk-type galaxies

We present the Waves in Nearby Disk galaxies Survey (WiNDS) consisting of 40 nearby low inclination disk galaxies observed through H$\alpha$ high-resolution Fabry Perot interferometry. WiNDS consists of 12 new galaxy observations and 28 data archived observations obtained from different galaxy surveys. We derive two-dimensional line-of-sight velocity fields that are analyzed to identify the possible presence of vertical velocity flows in the galactic disks of these low-inclination late-type galaxies using velocity residual maps, derived from the subtraction of an axisymmetric rotation model to rotational velocity map. Large and globally coherent flows in the line-of-sight velocity of nearly face-on galaxies can be associated with large vertical displacement of the disk with respect to its mid-plane. Our goal is to characterize how frequent vertical perturbations, such as those observed in the Milky Way, arise in the Local Universe. Our currently available data have allowed us to identify 20$\%$ of WiNDS galaxies with strong velocity perturbations that are consistent with vertically perturbed galactic disks.


INTRODUCTION
Over the last ten years, oscillatory perturbations towards the outskirts of the Solar Neighborhood in the Galactic disk have been reported several times (e.g. López-Corredoira et al. 2002;Momany et al. 2006;Widrow et al. 2012;Slater et al. 2014;Price-Whelan et al. 2015;Xu et al. 2015;Antoja et al. 2018;Gaia Collaboration et al. 2021). In particular, using 11k main-sequence stars from Sloan Digital Sky Survey (SDSS, York et al. 2000), Widrow et al. (2012) detected asymmetries in both, i) the bulk velocity of Solar Neighborhood stars associated to a breathing pattern, i.e. compression and rarefaction motions; and ii) in the vertical stellar number density distribution. The latter was related to a bending pattern, i.e. local displacements of the disk from the midplane (Widrow et al. 2014). This study was followed up by Xu et al. (2015) who analyzed measurements of stellar number counts of main-sequence stars located at Galactic latitudes 110 • < l < 229 • , as a function of Galactocentric radius. They showed that the amplitude of the perturbations in the vertical stellar number counts, i.e. the displacement of the disk with respect to its midplane, increases towards the outskirts of the Galaxy. This type of oscillatory perturbation on the stellar and/or gaseous component of the disk is known as a corrugation pattern. Thanks to the analysis of the full 6D phasespace information for more than six million stars from Gaia Collaboration et al. (2018), Antoja et al. (2018) showed that, indeed, our own Galactic disk is undergoing phase mixing of a non-equilibrium configuration, a perturbation likely asso-ciated with the interaction of the Milky Way with Sagittarius dwarf galaxy in the past (e.g. Gómez et al. 2013;Laporte et al. 2018bLaporte et al. , 2019Laporte et al. , 2020. More recently Gaia Collaboration et al. (2021) confirmed that the Milky Way disk has very complex dynamics, with vertical velocity perturbations that can be partially described as a bending wave.
Interestingly, evidence of such complex vertical oscillatory patterns on external galaxies, such as those observed in the Milky Way, are still extremely limited. Vertically perturbed disks have been extensively observed, especially in edge-on galaxies. Studies of 21 cm neutral hydrogen line observations show that most galaxies that have extended HI disks, with respect to the optical, are warped. Using a sample of 26 edge-on galaxies located in different environments with inclinations i > 75 • and blue diameters larger than 1.5 , from Westerbork HI Survey of Spiral and Irregular Galaxies (WHISP, van der Hulst et al. 2001), García-Ruiz et al. (2002) showed that all galaxies in their sample with an extended HI disk (∼ 76%) are warped with respect to the optical disk. They also found that warped disk galaxies are more often found in less dense environments, indicating that tidal interactions are not the only mechanism producing Galactic warps. However, in denser environments, the disk vertical perturbations typically show larger amplitudes than those in less dense environments. This is likely due to the fact that, in denser environments, galaxies are more likely to undergo violent tidal interactions with companion galaxies.
For observations in optical bands, Ann & Park (2006) analyzed a sample of 325 galaxies from the SDSS, which consists of a majority of late-type spiral galaxies with inclinations i > 84 • and isophotal major axis length D 25 > 1 at µ B = 25 mag arcsec −2 . The study showed that 73% of the sample contained a warped disk: 51% corresponds to Sshaped and 22% to U-shaped perturbations. For the analysis of the environment-dependence on the frequency of warps, Ann & Park (2006) analyzed 75 galaxies with redshift data from the SDSS DR3, of which 56 have companion galaxies. For the subset of galaxies that have companions, 64% show warped disks. For galaxies without companions, 81% show warped disks. This result is consistent with those presented in García-Ruiz et al. (2002), i.e, in poor environments a higher frequency of warped disks is found than in rich environments. Vertical perturbations have also been studied in edge-on galaxies using dust lines. Narayan et al. (2020) analyzed a sample of five nearby objects with inclinations i > 80 • and also showed that such perturbations are common.
S-shaped warps are expected to be a very common feature of both interacting and isolated late-type galaxies. Gómez et al. (2016Gómez et al. ( , 2017 studied in detail the origin and evolution of vertical perturbation in the disk of late-type galaxies using high-resolution cosmological magneto-hydrodynamical simulations of the Auriga Project . By means of the construction of mass-weighted mean height, Z , and mean vertical velocities maps, V Z , for both the stellar and cold gas components, these studies found that 70% of the Milky Way-like models showed strongly vertically perturbed disks. Interestingly, while only half of the vertically perturbed disks showed integral-sign (S-shaped) warps, the remaining half showed a more complex vertical structure, i.e. corrugation patterns, indicating that such structures are expected to be common.
A possible reason behind the lack of detected vertical patterns relates to the complexity behind their observation. On edge-on galaxies, such pattern can be easily confused with Sshaped warps due to projection effects. However, as discussed in Gómez et al. (2017), such a corrugation pattern should be detectable along line-of-sight velocity (V los ) of nearly face-on galaxies. Due to the oscillating nature of a corrugation, global patterns on the (V los ) field can be directly linked to global patterns in a Z map. Until recently, previous efforts to detect corrugations in the cold gas disk component of nearby galaxies were based on H α emission through long-slit spectroscopy. Alfaro et al. (2001) analyzed NGC 5427, a galaxy with an inclination angle of 30 • , and found wavy-like variations in the velocity profile of the ionized gas. Sánchez-Gil et al. (2015) analyzed a sample of four spiral galaxies and showed that two presented corrugation patterns. However, due to the limited coverage associated with the long-lit spectroscopic observations, such patterns could be confused with the effects of local perturbations such as fountain flows. In a recent work, Gómez et al. (2021) presented a full 2D velocity map of the low inclination galaxy, VV304a. Their study was based on Fabry-Perot interferometer H α observations. This technique has the great advantage of providing high spectral resolution in a narrow frequency range, along with good spatial resolution (in this work ∼ 2 ). As such, it allows resolving local velocity perturbations with amplitudes of the order of 10 km s −1 . The study showed, for the first time, that VV304a presents global and coherent velocity perturbations that are consistent with a corrugation pattern. They also showed that these velocity perturbations cannot be described by the effects associated with the presence of axisymmetric features on its disk, such as the bar and spiral structure. Thus, they conclude that these perturbations are a direct consequence of the gravitational interaction with its similar mass companion galaxy, VV304b.
The characterization of vertical perturbations in nearby disk galaxies can provide very valuable information about their recent interaction with their environment. Several mechanisms can be behind the formation of warps and corrugation patterns (Sellwood 2013;Gómez et al. 2016Gómez et al. , 2017. One of them is the tidal distortion of a pre-existing disk by an exter-nal perturber, such as the cases of our own Milky Way (Antoja et al. 2018;Laporte et al. 2019) and VV304a (Gómez et al. 2021). The strong tidal torques being exerted on a preexisting disk as relatively massive satellites fly by or merge can induce strong vertical perturbations (Ostriker & Binney 1989;Quinn et al. 1993;Velazquez & White 1999;Bailin 2003;Kazantzidis et al. 2009;Gómez et al. 2013;D'Onghia et al. 2016;Gómez et al. 2017;Laporte et al. 2018a,b). In addition, the torques associated with a non-spherical mass distribution of dark matter can also trigger the formation of vertical patterns (Debattista & Sellwood 1999;Jiang & Binney 1999;Shen & Sellwood 2006;DeBuhr et al. 2012;Yurin & Springel 2015;Gómez et al. 2016;Laporte et al. 2018a). Other possibilities may include misaligned accretion of cold gas due to the cooling of hot gas halo, infalling from the cosmic web, or being left by gas-rich host-satellite interactions (Sancisi et al. 2008;Jiang & Binney 1999;Roškar et al. 2010;Aumer et al. 2013;Radburn-Smith et al. 2014;Gómez et al. 2017) and ram-pressure of the surrounding intergalactic material (Haan & Braun 2014). The objective of this work is to search for velocity perturbations in the disks of nearby latetype galaxies. To achieve this goal we use high-resolution observations of near face-on galaxies obtained with a Fabry-Perot interferometer. For this purpose, we present the Waves in Nearby Disk galaxies Survey (WiNDS), which consists of a data set with high spectral resolving power (R ∼ 10000) and large spatial coverage, making them ideal for investigating the kinematics of the ionized gas in galactic disks. This study is the first step in an effort that aims to characterize the history of interactions of nearby galaxies, and estimate the frequency with which vertical disturbances such as those observed in the Milky Way arise in the Local Universe.
The paper is organized as follows. In Section 2 we introduce the WiNDS data sample and we describe our new observations. In Section 3, we discuss the data reduction procedure. Section 4 and Section 5 describes the data analysis and selection and quantification criteria of bending modes. Finally, we present our results in Section 6 and in Section 7 we provide a discussion and our conclusions. In Appendix A the comments for each individual galaxy and maps of the new observations made in this work are given. The residual velocity maps are shown in Appendix B. The rotation curves are displayed in Appendix C. In Appendix D the image processing method is described.

WINDS DATA SAMPLE
In this Section, we present our survey WiNDS, an ongoing observational campaign currently comprising 40 nearby nearly-face-on spiral galaxies. WiNDS consist of 12 new observational data and we complement it with additional archival data of 28 late-type galaxies from different surveys. The entire data sample contains 3D data cubes obtained using Fabry-Perot interferometer with a resolving power at the H α rest wavelength of R ≈ 10000. The galaxies in the WiNDS sample were selected according to three main criteria: • Distance: nearby spiral galaxies were selected to have systemic velocities v ≤ 3000 km s −1 , corresponding to less than 45 Mpc.
• Inclination: all galaxies have previously estimated low inclination angles i ≤ 40 degrees. Both, morphologically and kinematically estimated inclination were considered in this step.
• Size: all galaxies have projected diameters, a p , between 2 arcmin ≤ a p ≤ 4 arcmin.
The goal of this selection criteria is to allow us to resolve velocity perturbations in the observed disks with amplitudes as small as 10 km s −1 . In what follows we discuss WiNDS new observations in detail. We also briefly describe the additional surveys utilized to complement WiNDS. In Table 1 we list some of the main properties of the WiNDS sample. The blue histograms in Figure 1 show the distribution of the most relevant parameters for the final WINDS sample. The solid curves depict the smooth continuous approximation of the underlying histograms. These smooth histograms were obtained using a Gaussian Kernel Density Estimator. The green histograms in the same figure show the parameter distribution of the vertically perturbed galaxy candidates within the WiNDS sample (see Section 5). Although the WiNDS sample is not complete in terms of, e.g., mass and magnitude, it allows performing a first systematic search of vertically perturbed disks in the nearby Universe. The 40 selected galaxies span different environments such as clusters, groups, and field galaxies, and a wide variety of morphologies (Sa to Sc). Note however that this is an evolving project, and more objects will be added in the future to complement the available sample. Figure 2 shows the optical images of the WiNDS sample, obtained from SDSS DR9 using the g, r and z bands.

WiNDS: New data cubes
The new observations for WiNDS were obtained using the Fabry-Perot interferometer at the Observatoire de Haute Provence Observatory (OHP, France). In addition to the selection criteria mentioned in Section 2, we focus on galaxies with previously reported H α observations (James et al. 2004;Sheth et al. 2010). The new data cubes include galaxies with inclination angles ranging from 6 • to 33 • , except for NGC 2500 which has a i ≈ 40 • . The sample contains galaxies with a range of morphological types (between Sa to Sc) and a wide absolute B-band magnitude range (−17.9 ≤ M B ≤ −21.9 mag). It also contains both isolated and interacting galaxies belonging to different environments. The 12 new galaxies observed as part of WiNDS are mainly located in the Northern Hemisphere, the list of 12 galaxies is reported in Table 1. In particular, NGC 2763 is the only galaxy in our sample observed from the Southern Hemisphere using the SOAR Adaptive Module Fabry-Perot (SAM-FP, Mendes de Oliveira et al. 2017). This galaxy was added to the sample at a later time, based on the availability of the SAM-FP.

WiNDS: Observation of new data cubes
Most WiNDS observations for the new galaxies were performed in February 2019 using the 1.93-m telescope at the OHP. The observations were made with a mean integration time of 2 hours per galaxy. We note that due to poor weather conditions three galaxies in the sample were observed with a seeing slightly over 3 arcsec.
The new datacubes were obtained through the Gassendi HAlpha survey of Spirals (GHASP) instrument, a focal reducer containing a Fabry-Perot scanning interferometer which has a large field-of-view (FoV), of 5.9 × 5.9 arcmin 2 , an aperture ratio of f /3.9, and a pixel scale 0.68 arcsec pixel −1 . The Free Spectral Range (FSR) for the two used interferometers with interference orders of 798 and 2600 are 375.9 km s −1 and 115.4 km s −1 , respectively. The resolving power reached in our survey was R = 10000 − 28600, which translates into a velocity resolution of 11.6 − 3.6 km s −1 , respectively. As discussed in Gómez-López et al. (2019) the detector, Image Photon Counting System (IPCS), has the advantage of a zeroreadout noise, which allows a very fast scan of the interferometer through the entire FSR. The use of an IPCS makes it possible to carry out short exposures (10 seconds) and to make a large number of cycles (typically 20) to average the variations in observation conditions. Within new observation of WiNDS, the H α line can be found within the range going from 6563 Å to 6627 Å. A total of seven filters, each with a Full Width at Half Maximum (FWHM) of 15 Å, are available to isolate the H α emission line.
For the wavelength calibration of the data, a cube is obtained before and after the observation of each galaxy to account for possible changes in observing conditions. Each calibration cube is obtained using the Neon emission line (6598.95 Å), which is isolated using a narrow-band interference filter (∼ 15 Å) centered on this wavelength. This line is used because it is intrinsically narrow and close to the wavelength of the redshifted H α line.
An additional galaxy, NGC 2763, was observed in February 2016 using the 4.1-meter telescope at the Southern Astrophysical Research (SOAR), Cerro Pachón, in Chile. The instrument used is called SOAR Adaptive Module Fabry-Perot (SAM-FP). It is a restricted-use instrument providing a new mode of operation of SAM for spectroscopy using a Fabry-Perot (FP) and SAM-Imager (SAMI). The SAM-FP provides a FoV of 3 × 3 arcmin 2 , a CCD pixel scale of 0.045 arcsec pixel −1 (after electronic binning on the detector) and 48 wavelength channels, in the third dimension of the datacube. For the SAM-FP observations, the interference order was 609 with a FSR of 492.6 km s −1 and resolving power of R = 10000 which translates into velocity sampling of 10.3 km s −1 per channel (for a resolution of ∼ 20.6 km s −1 ), see Mendes de Oliveira et al. (2017) for more details.

Additional Data Archive
The new data cubes described in Section 2.1 were complemented with 28 galaxies previously observed, available at the Fabry-Perot interferometer database 7 and Herschel database 8 at the Laboratoire d'Astrophisique de Marseille. The galaxies were selected from the GHASP (G), SINGS-H α (S), VIRGO-H α (V) and HRS-H α (H) surveys, each of them briefly discussed below. Some basic parameters, describing the corresponding observations, are listed in Table 1.
The data cubes available from the previously mentioned surveys were reprocessed in order to homogenize the galaxy analysis. This was performed considering the same signalto-noise and smoothing method for all data cubes. Data cubes from the GHASP, VIRGO-H α , and SINGS-H α surveys are available in the Fabry-Perot database website. On the other hand, HRS-H α data cubes are available in the Herschel Database.
2.2.1. GHASP GHASP (Epinat et al. 2008a,b) survey consists of 3D H α data cubes of 203 late-type galaxies (spiral and irregular). This sample contains a wide range of morphological types, from Sa to Irr, with absolute magnitudes in the range of −16 ≤ M B ≤ −22 mag, stellar masses ranging from 10 9 M to 5 ×10 11 M and are located in nearby low-density environments. From the total sample of 203 galaxies, 83 are strongly barred galaxies (SB or IB) and 53 moderately barred galaxies (SAB or IAB). According to our selection criteria described in Section 2, thirteen galaxies were selected and are reported in Table 1.

SINGS-Hα
The Legacy survey Spitzer Infrared Nearby Galaxies Survey (SINGS, Kennicutt et al. 2003) consists of 75 nearby galaxies selected with a wide coverage of morphological type (E to Im), luminosity types (IR quiescent to luminous IR galaxies) and CO/HI ratio covering over 3 orders of magnitudes. These galaxies are located in different environments such as galaxy groups, clusters, and low-density fields. The SINGS survey was created to characterize the infrared emission in a wide range of galaxy properties and star formation environments. The SINGS-H α data is a H α kinematic followup survey of SINGS (Daigle et al. 2006a;Dicaire et al. 2008) and consists of 65 late-type galaxies which present HII regions. According to the criteria described in Section 2, the selected galaxies are reported in Table 1.

VIRGO-Hα
The 3D H α data cubes from the Virgo survey, VIRGO-H α , (Chemin et al. 2006) consists of 30 spiral and irregular galaxies. This survey is a subsample of the Virgo Cluster Catalog (Binggeli et al. 1985(Binggeli et al. , 1993, which is the nearest cluster to the Milky Way and counts ∼ 1400 members, mainly dwarftype galaxies. The galaxies considered in VIRGO-H α have an apparent magnitude greater than B 0 t = 12 mag, wide morphological types (S0/a to Im), and inclination angles between 25 • and 89 • . All galaxies are located in the cluster's core and its extension towards M49. According to our criteria selection described in Section 2, the three selected galaxies are reported in Table 1.

HRS-Hα
The Herschel Reference Survey, HRS-H α , (Gómez-López et al. 2019) consists of 152 star-forming galaxies observed using the 1.93-m telescope at the OHP. This survey aimed to study the relationship between the baryonic and dynamical mass of galaxies. The HRS-H α galaxies sample spans a wide range of morphologies (from Sa to Sm-types, including Blue Compact Dwarfs) and stellar masses (10 8 M ≤ M star ≤ 10 11 M ).
From the HRS-H α sample, and excluding overlap with previous surveys (GHASP, SINGS-H α and VIRGO-H α ), we found seven galaxies that follow our selection criterion. The names of the six targets are listed in Table 1.
Note that not all of the galaxies come from the same source. Table 2 shows our observations broken-down by instrument, telescope, survey of origin, and number of galaxies observed. The most relevant observations parameters are listed in Table  3. For completeness, in the same table we also list the properties of the additional observations extracted from the data archive.     The WiNDS data sample was reduced using the pipeline based program COMPUTEEVERYTHING 9 and the RE-DUCWIZARD 10 interface, following the steps described in Daigle et al. (2006b) and Epinat et al. (2008b).
The data reduction process includes a) the wavelength calibration, b) the night sky lines subtraction, c) the astrometry process, d) the adaptive spatial binning using the 2D Voronoi tessellation, whose implementation finally produces a smoothing specially adapted to a given SNR, e) the generation of the 2D momentum maps, f) a semi-automatic cleaning of the velocity fields and g) velocity dispersion correction. In the following we briefly describe these processes: a) The calibration process consists of applying a phase map to provide the same wavelength/velocity origin to each profile, for each pixel of the observation datacube.
The 2D phase map is computed from the two calibration datacubes. During the phase correction process, the individual exposures are re-centered with respect to each other, using field stars, to minimize telescope drifts and instrument bend.
b) The OH emission lines are the main foreground sky contamination. The sky subtraction is done by dividing the collapsed datacube into a galaxy-dominant and a sky-dominant region. The night-sky emission is interpolated from this galaxy-free region using a two-degree polynomial, and subsequently removed from the whole datacube.
c) The next step is astrometry which is done through the KOORDS task in the KARMA 11 package (Gooch 1996), comparing star fields between XDSS 12 R-band images and our continuum images of each galaxy.

d)
For the estimation of The Signal-to-Noise ratio (S/N), which is related to the flux in the line, to the spectral resolution and the r.m.s. in the continuum. We use adapted spatial binning through Voronoi tessellation, as described in Daigle et al. (2006b). This allows us to obtain the highest spatial resolution for a given S/N, which is the main advantage of this binning technique. The Voronoi bins are constructed iteratively from a given pixel, by accreting adjacent pixels until the resulting spectrum reaches the desired S/N. Pixels where the S/N is greater than the chosen threshold value will not be affected by the binning process, 9 /https://projets.lam.fr/projects/computeeverything 10 https://projets.lam.fr/projects/fpreducwizard 11 KARMA tools package is available on website https://www.atnf.csiro.au/computing/software/karma/ 12 ESO Online Digitized Sky Survey https://archive.eso.org/dss/dss while regions of low S/N will emerge from noise without being contaminated by adjacent regions of higher S/N. For the WiNDS sample, the objective is to obtain a S/N ≥ 7 per spatial bin, where S/N is considered as the square root of the flux in the line. We implemented a Hanning smoothing of a one-spectral-element of the spectrum which preserves the flux.
e) From this Voronoi binned datacube, the different momentum maps are computed as explained in Daigle et al. (2006b): after identifying the H α line boundaries, the continuum emission around the line is estimated and subtracted from the spectrum. The moments of the line are then estimated. Namely, the line or monochromatic flux is the zeroth-order moment (intensity integrated within the line boundaries), the line-of-sight velocity (V los ) is the first-order moment (intensity-weighted velocity sum within line boundary), and the LoS velocity squared velocity dispersion (variance) as the second-order moment (difference between the intensity-weighted squared velocity and squared of the intensity-weighted velocity within the line boundaries).
f) In some cases, the lowest S/N regions are strongly affected by night sky line contamination. Thus a semiautomatic cleaning is performed. Indeed, to achieve the desired signal-to-noise ratio in the outer-most galactic regions very large Voronoi bins are produced. These are mainly associated with sky subtraction residuals and background emission. Those bins are first semiautomatically erased and then manually deleted when the semi-automatic process is not sufficient.
g) Finally, velocity dispersion maps were corrected by instrumental width. This correction was done by subtracting the mean dispersion of the instrument contribution to the observed velocity dispersion map, as follows: where σ obs corresponds to the observed velocity dispersion and σ inst is the mean dispersion of the instrument contribution. Moreover, σ inst is a function of the FSR, the interference order p, and the resolving power R: The results of the data reduction process are presented below for a WINDS subsample where panels a), b) and c) show the XDSS B-band image, the H α line-of-sight velocity map, and the H α velocity dispersion map, respectively. In addition, the results of the data reduction for the remaining new observations made in this work, where XDSS blue image ( (2)  . coherent kinematic perturbations on our sample of late-type, low-inclination galaxies. In particular, we seek for perturbations that are consistent with bending modes produced by warps and corrugation patterns. In this Section, we discuss the procedure followed to analyze our resolved velocity fields.
As discussed in Gómez et al. (2021), the observed V los field of each galaxy is bound to contain contributions from three different velocity components: the distributions within the disk plane, radial V R and rotational, V rot ; and the perpendicular velocity to the disk plane, V Z . More precisely, each V los can be described as: where V sys is the systemic velocity of the galaxy and it is considered as a fixed value. The polar coordinates (r, θ) in the plane of the galaxy are measured from the position angle (PA), inclination (i) and sky position (α, δ) of the rotation center. The first step to identify global perturbations in the data is to subtract from the V los map an axisymmetric model of V rot . To obtain such an axy-symmetric velocity model we assume that, at first order, the contributions from the radial and the vertical velocities to V los are negligible. Therefore, the observed velocity Equation 3 reduces to: To obtain V rot (R) we use the tilted-ring method (Begeman 1987), which assumes that the galaxy can be analyzed using concentric rings along the major axis, which are described by the parameters V r rot , i r , PA r . Here the supra index r indicates that we are referring to the rings. In this work, we use an improved tilted-ring method, described in detail in Epinat et al. (2008a), and derive V rot (R) for each galaxy bin is modeled using the modified Zhao function (Epinat et al. 2008a), This four-parameter model was especially chosen for its versatility to adjust to a very diverse set of different rotation curve shapes. The four parameters involved in the V rot model are v t and r t which correspond to the velocity and radius when the rotation curve changes from an increasing velocity to a flat regime, respectively; a and g which are related with sharpness of the turnover. As discussed in Epinat et al. (2008a), for each galaxy the fitting procedure requires a set of initial values for the Zhao function parameters as well as for (α, δ), i and PA. In our work, the values of (α, δ) are set to the peak of the emission in the continuum image and set as fixed values. For the systemic velocity, PA, and i, initial values were extracted from the literature. The resulting values of i and PA, obtained as a result of the fitting procedure, are referred to as kinematically inferred values. Note that, in those cases where disks are projected nearly perfectly face-on (i < 10 • ), the value of i is not allowed to vary and, thus, is kept fixed to their estimated morphological value. Instead, v t , r t , a, g parameters are always allowed to vary. Finally, V rot for each Voronoi binning is estimated through the minimization of χ 2 based on the Levernberg-Marquardt method (Press et al. 1992), computing an iterative 3.5σ clipping on the observed bin-centroid velocity field. All kinematical parameters (i k , PA k ,V sys ) estimated for WiNDS galaxies are listed in Table 4. The estimation of errors in the determination of kinematic parameters for the WiNDS sample was calculated using the power spectrum of the residual velocity field and the application of a Monte Carlo method, as described in more detail in Epinat et al. (2008a). It is worth noting that the degeneracy between V max and i when fitting our models (Begeman 1989), especially for very low inclined disk, does not affect our results as we are focusing our analysis on residual velocity fields, V res . The H α V max values derived in this work are listed in Table 4 Once an unperturbed axisymmetric model of V los is obtained for each galaxy, we generate their corresponding V res fields. This is done by subtracting the properly V los model from the observed V los maps. It is worth recalling that global and coherent features in a V res field of a late-type galaxy could be the result of kinematic perturbations induced by features such as a bar or spiral structures, but also from improper model parameters (such as center, systemic velocity, position angle of major axis, inclination or rotation curve model). The multi-polar signatures expected for ill-defined parameters (Warner et al. 1973) are not observed in our data.
It is also worth noticing that small variations in the magnitude of the V res could be attributed to uncertainties in the   In order to verify the sensitivity of the method used to the initially chosen and fixed inclination value, the model is rerun considering the extreme values of the inclination that are within the error range shown in Table 4.
The contributions from in-plane flows are expected to be small in low inclination galaxies such as those studied in this work. Using a suite of numerical simulations, Gómez et al. (2021) characterized the contribution from in plane flows to the V los driven by non-axisymmetric structures. Their models consisted of a disk galaxy, projected into an inclination of 35 • , considering a bar with similar characteristics to that of the Milky Way and spiral overdensities that ranged from 100% to 1000% density contrast (∆ρ/ρ) with respect to the disk background density. Even for ∆ρ/ρ = 1000 the resulting velocity perturbations in the corresponding V res field were 10 km s −1 .

MODES
In this Section we describe the selection criteria applied to select those galaxies within the WiNDS sample that show velocity perturbations consistent with a bending mode, such as a warp or corrugation pattern. To avoid selecting discrete and local vertical perturbations, such as those associated, e.g., fountain flows we have established the main criteria to select strong candidates are the following: • We focus on galaxies that present an extended H α coverage, extending for at least 0.7 R opt . This criterion allows us to globally explore the kinematics of the disks, especially on the outer regions where bending modes typically show their stronger amplitudes.
• We focus on galaxies where global perturbations show amplitudes that are > 10 km s −1 . This enhances the chances that the observed perturbations in this low inclination sample of disks are not mainly driven by the axisymmetric components of the galaxies.
• We applied a Gaussian low-pass filter to the V res images in Fourier space with cutoff frequency equivalent to a spatial distance of ≈ 500 pc to highlight relevant and large perturbations (see Appendix D.1). This allows us to smooth out local and discrete perturbations, and focus on large-scale and coherent features. Based on these images, using a circular grid centered at the center of the galaxy, evenly spaced by 0.3 R opt as shown in Fig  3, we select as potential candidates all those galaxies where the perturbations cover at least an azimuthal extension of 60 degrees. We consider perturbations that at a given azimuthal angle cover a radial extension 0.3 R opt as corrugation patterns, whereas the remaining as typical warps.
Examples of the resulting filtered images are shown in Figure 3. The upper panel shows an example of well-defined global and coherent velocity perturbations. The bottom panel shows the case where the map was consistent with an unperturbed velocity field. According to the these criteria eight galaxies show V res maps consistent with a vertical perturbed disk.
We note that previous studies have analyzed velocity dispersion maps, V disp , to identify vertically perturbed disk galaxies. For example, Jiménez-Vicente & Battaner (2000) used V disp to detect plausible shells or chimneys in the galactic disk of NGC 5668. This was done by comparing the geometry of the V res and V disp maps and associating local perturbations in both maps to vertical motions in the galaxy. In our work, V disp maps are obtained as a sub product of the data reduction process (see Section 3) and, thus, we will present them together with the V res maps. However, as opposed to Jiménez-Vicente & Battaner (2000), our goal here is to identify global rather than localized and discrete velocity perturbations. As a result, V disp maps were not explored in detail. We will further analyze these maps in follow-up work.
We emphasize that our selection criteria cannot confirm or rule out the presence of vertical perturbations in our disks. Although the objective is to detect global and coherent perturbations, it should be noted that our method does not exclude the presence of fountain flows, which are shown as local and discrete perturbations but rather, both types of perturbations could coexist. Nonetheless, note that we have carefully checked the H α intensity profiles in all galaxies that were selected as potential vertically-perturbed objects. We found that all candidates within the WiNDS sample show, in general, well-behaved Gaussian profiles. That is, there is no presence of strong and intense multiple components in the H α profiles due to the presence of other velocity components (such as those observed in merging systems such as HCG31, for example, Amram et al. 2007).

RESULTS
In this Section, we analyze and discuss the resulting V res fields of the eight WiNDS galaxies that fulfilled our quantification criteria described in Section 5. However, we discuss the remaining new observed galaxies in the Appendix A.1 and we present in Appendix B.1 the V res maps of galaxies that do not accomplish the criteria described in Section 5.
In Figure 4, we show the results obtained for the galaxy NGC 3344. This galaxy has a very low morphologycal inclination angle (i ∼ 18 • , Paturel et al. 1997) and, thus, represents an ideal candidate to search for vertical perturbations on its galactic disk.  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel.
show the resulting V los map, obtained following the procedure outlined in Section 3. The map shows a good H α coverage of the overall disk (1 R opt ), despite being limited by the FoV of the GHASP instrument. The amplitude of V los reaches ≈ 100 km/s. Panel c) shows the H α velocity dispersion map, normalized to σ 95th = 44.2 km s −1 13 Panel d) shows the residual velocity field V res , obtained after subtracting the axisymmetric rotational velocity model from V los . Interestingly the V res map reveals a global, strong, and coherent oscillating-like pattern throughout the entire disk of the galaxy, with an amplitude of the order of ±30 km s −1 . To compare the disk morphology to the structure observed on its residual velocity field, on panel e) we show the galaxy B-band image with the contours obtained from the V res field. Notice that no clear correlation between axisymmetric patterns and velocity perturbations can be observed. The extended black box or slit, placed across the 13 A normalization of the velocity dispersion to the 95th percentile was used for all galaxies in order to avoid outlying dispersion measurements and better show possible correlations between residual velocity fields and velocity dispersion. disk on panel e), highlights the multiple transitions from positive to negative values experienced by V res as a function of galactocentric distance. This is better shown on panel g) where we highlight the behaviour of V res across this particular slit (black crosses). Here the black line represents the corresponding smoothed data using a moving average function with a 7 data points window, while the dashed vertical line indicates the galaxy centre. The direction of the slit was chosen to highlight transitions of the residual velocity values. NGC 3344 is a galaxy that is cataloged as isolated, with no nearby galaxies within a projected radial distance of 150 kpc and radial velocity difference ∆ ≤ 500 km s −1 . We recall that NGC 3344 has a very mildly inclined disk, i ∼ 18 • . Thus, the contribution from in-plane flows to the resulting V res , associated with e.g. bars and/or spiral structure, is expected to be small. As previously discussed, a plausible interpretation of these V res perturbations is associated with the contribution from a vertical corrugation pattern. Such corrugation pattern could be the result of several mechanisms, such as close interactions with satellites (Sancisi et al. 2008;Jiang & Binney 1999;Roškar et al. 2010;Aumer et al. 2013;Radburn-Smith et al. 2014;Gómez et al. 2017), torques associated with either misaligned triaxial DM halos or even DM overdensity wakes (Ostriker & Binney 1989;Quinn et al. 1993;Velazquez & White 1999;Bailin 2003;Kazantzidis et al. 2009;Gómez et al. 2013;D'Onghia et al. 2016). Vertical perturbations on this galaxy have been previously reported by other authors. Indeed Briggs (1990) reports the presence of a warp in the outer regions of NGC 3344 after analysing observations of the HI 21 cm line. Verdes-Montenegro et al. (2000) presented a detailed study of this galaxy, also using HI 21 cm lines. The authors suggest a possible relationship between the inner and outer star-forming rings with a strong spiral structure and the warp in the HI layer. However, they were not able to clearly link the warp structure with the perturbation in the more internal regions. Our analysis of NGC 3344, with a spectral resolution ∼ 10 km s −1 , reveals a very complex structure of NGC 3344. NGC 3344 is a very interesting candidate of the vertically perturbed galaxy that will be analyzed in greater detail on a follow-up study.
In Figure 5 we now show the results obtained for NGC 628 (M74). This is a nearly face-on galaxy with an inclination angle of 6.5 • ; thus, an ideal candidate to search for possible vertical perturbations. On panel a) we display a B-band image of the galaxy, obtained from the XDSS. NGC 628 is a nonbarred galaxy with two main spiral arms emerging from its bulge. The H α observation coverage of NGC 628 corresponds to a mean radius of 12.5 kpc, which represents 0.8 R opt . The H α velocity map (panel b), shows an amplitude of ∼ 60 km/s, as expected from a low inclination galaxy. For comparison we also show on panel f) the HI velocity map from The HI Nearby Galaxy Survey (THINGS) (Walter et al. 2008, THINGS), observed with Very Large Array (VLA), which is consistent with the H α velocity map. Interestingly, studies based on neutral hydrogen (Kamphuis & Briggs 1992) have detected the presence of an elongated warp structure at around 12 projected distance from the nucleus. On panel c) we show the derived H α velocity dispersion map, normalized to σ 95th = 14 km s −1 . The V res field of NGC 628 is shown on panel d). Interestingly, global and coherent velocity perturbations, of the order of 20 km s −1 , can be clearly observed. Panel e) allows us to correlate NGC 628-disk morphology with its residual velocity map. The velocity perturbations on this galaxy are more clearly highlighted on panel g) where we follow the mean V res along the highlighted slit, suggesting a warp-like structure.
As shown on its B-band image (panel a), an extended tail is observed in the southwest direction of the disk. This substructure has been previously reported in Kamphuis & Briggs (1992). Although NGC 628 has currently no close companions, a dwarf galaxy system (UGC 1176 and UGC 1171) is located at 140 kpc. The overall mass of the satellite system has been reported to lie on the order of magnitude of the total NGC 628 mass (Michałowski et al. 2020). This system may have played an important role in the evolution of NGC 628. As discussed by Briggs (1986), if the dwarf system is at rest at its current location concerning for to NGC 628, it would fall into its host in about 1 -2 Gyr. However, a more recent study by Michałowski et al. (2020) suggests that a tidal origin for the asymmetric HI tail located on the south-western outskirt is an unlikely scenario. They argue that the tail does not resemble recent tidal features, whereas an older feature would wind almost symmetrically around the galaxy. In addition, recently induced tidal tails typically show m = 2 patterns, with two nearly symmetrical arms. Instead, NGC 628 shows a strong m = 1 like spiral structure. Nonetheless, they do not rule out a possible previous interaction with the neighboring dwarf pair. It is interesting to compare with the results presented in Gómez et al. (2017), based on fully cosmological simulations. In particular, their model S13 shows similar morphology to NGC 628, with a strong m = 1 arm extending from the inner galactic regions. This model displays a vertically perturbed stellar and cold gas disk, but no recent interaction with massive satellites (M sat /M host > 0.03). Instead, S13 experienced a massive gas-rich merger 5 Gyr in the past. As a result, its pre-existing disk is destroyed but a new thin disk quickly forms thereafter due to the re-accretion of misaligned gas. This gas accretion gives rises to both the vertical perturbations observed in this modeled galaxy and its strong m = 1 morphology. This suggests that the velocity perturbations observed in NGC 628, as well as its m = 1 spiral morphology, could be linked to the recent accretion of misaligned cold gas.
In Figure 6 we focus on NGC 6946. This is a spiral galaxy classified as SABcd with a kinematical inclination of i ∼ 17 • , an optical radius of ≈ 10 kpc, and an estimated total mass of 1.2 × 10 11 M (Pisano & Wilcots 2000). Optical and neutral hydrogen observations of NGC 6946 have revealed its strong spiral arms, a particularly high star formation rate, a nuclear starburst, and a weak bar (Pisano & Wilcots 2000;Schinnerer et al. 2006). Panel a) shows its XDSS B-band image, where we can clearly see its spiral nature. NGC 6946 is considered to be isolated (Tully 1988). Using HI observations, previous studies have detected a dozen of low mass irregular dwarf galaxies companions (Pisano & Wilcots 2000;Karachentsev et al. 2000Karachentsev et al. , 2005. However, at the present day NGC 6946 does not show any sign of undergoing strong direct gravitational interactions with its satellite population. Its two most massive companions, UGC 11583 and L149 have estimated total masses of 1.5 × 10 9 and 3.0 × 10 8 M and projected distances of 83 and 75 kpc, respectively (Pisano & Wilcots 2000). The H α V los maps of NGC 6946, shown in panel b), has an amplitude of ∼ 150 km s −1 . Following our selection criteria, the map reveals a significant H α coverage, reaching at least a radius of 7 kpc, corresponding to 0.7 R opt . As in the case of NGC 3344, the coverage for NGC 6946 is limited by the FoV of the GHASP instrument. For comparison, we show (on panel f) the HI velocity map, extracted from Walter et al. (2008). Both maps are consistent, even though they are considering different components of the ISM and come from different instruments and techniques. Note, however, that our observations have a better spatial resolution, of 3 arcsec against the 6 -12 arcsec resolution in the HI observations. The H α velocity dispersion map is shown on panel c), normalized to σ 95th = 55.2 km s −1 . On panels d) and e) we show the galaxy V res map, which reveals global and coherent perturbations, reaching amplitudes 35 km s −1 . The V res velocity contours, overlaid on the B-band image (panel e), allow us to directly contrast this velocity structure with the morphology of NGC 6946. The oscillating nature of this velocity field is better highlighted on panel g), where we show the mean V res along the highlighted slit on panel e). It is worth recalling NGC 6946 spiral nature, and its low but not negligible inclination. As a result, a contribution from in-plane flow to the V res field can be expected. However, as shown by Gómez et al. (2021), for disks as inclined as NGC 6946 (i ∼ 40 • ), even an spiral structure 1000 times denser than its mean background disk density cannot generate velocity perturbations with amplitudes 10 km/s. This indicates that the observed global  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel. velocity perturbation could be partially linked to vertical velocity flows.
It is interesting to consider what mechanisms could be driving the perturbations in this galaxy. As previously discussed, NGC 6946 does not have nearby massive satellites that could be directly tidally interacting with its disk. However, even low mass satellites at relatively large galactocentric distances could significantly perturb an embedded disk through the excitation of dark matter overdensity wakes (Weinberg 1998;Vesperini & Weinberg 2000). For example, using cosmological simulations of Milky Way-mass galaxies, Gómez et al. (2016) studied the onset and evolution of a strong vertical pattern in the disk. The vertical pattern in this model, with an amplitude of ∼ 50 km s −1 , is the result of a satellite-host halodisk interaction. Interestingly, the satellite had a total mass of ∼ 5 per cent of the host and a pericentre distance of 80 kpc. The satellite was not massive enough to directly perturb the galactic disk but the density field of the host dark matter halo responded to the satellite passage and strongly amplified its perturbative effects. A similar scenario could be taking place in the case of NGC 6946, considering its low mass com-panions at a relatively large projected galactocentric distance. Other possible mechanisms, in addition to the previously discussed contribution from in-plane flows, could be related to smooth accretion of misaligned cold gas.
In Figure 7 we show NGC 4321 (M100), a grand-design nearby spiral galaxy belonging to the Virgo Cluster, classified as SABbc (RC3). The galaxy has a small bulge and two welldefined, prominent and, symmetric spiral arms, as shown on its B-band image (panel a). The nucleus of the galaxy is compact and bright. NGC 4321 is a very low-inclination galaxy, with a kinematical inclination of i ≈ 25 • and an optical radius of R opt = 16.3 kpc from RC3. Panel b) shows that the coverage of the V los distribution, obtained from the H α observations, is extensive and reaches nearly to its optical radius; i.e. 13 kpc ∼ 0.8R opt . The V los map shows an amplitude of ≈ 150 km s −1 . For completeness, we show on panel c) the H α velocity dispersion map, normalized to σ 95th = 33.1 km s −1 .
The derived residual velocity map, obtained after subtracting the axisymmetric velocity model is shown on panels d) and e). As in previous examples, this map reveals global and coherent velocity flows with amplitudes that can reach  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel.
∼ 30 km s −1 . Panel f) more clearly shows the radial variations of the residual velocity along the highlighted slit in panel e). Note the oscillatory behavior of V res across the disk. These perturbations, observed in the ionized gas of NGC 4321, could be the result of a recent interaction with one of its satellite galaxies. Indeed, NGC 4321 has two dwarf companion galaxies, NGC 4328 and NGC 4323, located at projected distances of 24 kpc and 28 kpc from their host, respectively (Hernandez et al. 2005). Knapen et al. (1993) showed that the HI component of NGC 4321 is mostly contained within the optical disk. However, they also identify a large HI extension that could be the result of the tidal interaction with NGC 4323. Indeed, Knapen et al. (1993) find important differences in the behaviour of its rotation between the approaching and the receding sides. They suggest these differences could be caused by deviations from circular motions in the outer disk probably due to a close passage of its companion galaxy. They also suggest this interaction could be the cause of the observed asymmetry in the total HI distribution. NGC 4321 also shows diffuse stellar extension in two directions, one of them towards NGC 4323, likely the result of ongoing interaction with this satellite (Hernandez et al. 2005).
NGC 5194 (M51a) is a well-known nearby grand design spiral galaxy classified as non-barred by RC3, with a very low inclination angle i = 20 • . NGC 5194 is tidally interacting with its companion NGC 5195 (M51b), an early-type SB0 galaxy. NGC 5194 shows an intense star formation activity at its center and along its spiral arms. Panel a), on Figure  8, clearly shows its strong spiral structure with signs of being lopsided, likely due to the interaction with its companion. Panel b) shows that the H α observations provide wide coverage of the disk, which allows us to obtain its velocity field up to distances of 12 kpc ∼ 0.8 R opt . The resulting V los map shows a velocity range of ≈ 150 km s −1 . It is interesting to compare our H α velocity map with the one derived from HI observations, obtained from THINGS using VLA (Walter et al. 2008), shown in panel f). As seen in NGC 628 and NGC 6946, the HI velocity map of NGC 5194, constrained to the same FoV as the H α data, shows a very good agreement with the ionaized gas velocity map, despite differences in both the spatial and spectral resolutions. In panel c), we show the velocity dispersion map normalized to σ 95th = 15 km s −1 .  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel.
The residual velocity map is shown in panels d) and e). It is interesting to note the large velocity perturbations throughout the disk, reaching amplitudes of 40 km s −1 . We recall the very low inclination angle of this galaxy, thus rendering it unlikely that these velocity perturbations are mainly due to in-plane velocity flow. Panel e) allows comparing the substructure of the V res map with the galaxy morphology as seen in the B-band image. In panel g) we follow the smoothed residual velocity profile along the highlighted slit. As before, this allows appreciating the radial oscillatory behavior of V res , consistent with what is expected for a corrugated disk. Shetty et al. (2007) presented a deep kinematic study of NGC 5194 using full 2D velocity distributions from interferometric CO and Fabry-Perot H α observations. They report a complex velocity field, with variations of the disk's PA and inclination along the galactic radius that strongly suggest a vertically perturbed disk and very significant out-of-plane motions. An HI warp has also been previously reported in NGC 5194 (Rots et al. 1990;Oikawa & Sofue 2014;Haan & Braun 2014), which has been attributed to the tidal interaction with NGC 5195. The presence of vertical velocity flows in this very in-teresting galaxy will be further analyzed in a follow-up study, following the procedure described in Gómez et al. (2021).
In Figure 9 we present NGC 5668, a nearly face-on latetype spiral galaxy with an estimated inclination angle 18 • (Schulman et al. 1996) and R opt of ∼ 14.3 kpc, located at 27.6 Mpc. The galaxy is classified as a SA(s)d by RC3 and, on its optical image (panel a), it shows a weak bar or oval structure on its inner 12 region. Observations of NGC 5668 were extracted from the GHASP sample (Epinat et al. 2008a). The H α coverage reaches out to 20 kpc corresponding to ∼ 1.4 R opt . On panel b) we show the resulting V los field, with a velocity amplitude reaching ≈ 70 km s −1 . The velocity dispersion map, normalized to σ 95th = 57.5 km s −1 , is shown on panel c).
The V res map, shown on panel d), reveals global perturbations with amplitudes > 20 km s −1 km/s which are also consistent with a warped H α disk. In panel e) we present the optical image with the overlay of the residual velocity map. The warp is better shown in panel f), where we can see the mean velocity profile along the highlighted slit on panel e). Note the large amplitude of this kinematical perturbation reaching  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel.
velocities > 20 km s −1 at the disk outskirts. High-resolution observations in the 21-cm line of the neutral hydrogen velocity field, from the Very Large Array (VLA), have detected a kinematic warp that begins, at least, at the end of the optical radius (Schulman et al. 1996). Note that the H α warp is very well aligned with the previously reported HI warp, indicating that the perturbation extends even within the optical radius of the galaxy. In addition to the HI kinematic warp, Schulman et al. (1996) detected high-velocity wings beyond the double-horned 21 cm profile, related to High-Velocity Clouds (HVCs) in the disk and halo of the galaxy. As discussed by Schulman et al. (1996), the high-velocity wings located outside the optical disk may be the result of infalling material, comparable to the Magellanic Stream observed in the Milky Way.
The origin of the kinematic warp, and the infalling material, could be related to a past interaction with its neighbor galaxy UGC 9380, currently located at a projected distance of ∼ 200 kpc of NGC 5668. It is worth mentioning that Jiménez-Vicente & Battaner (2000) also analyzed the 2D velocity field of this galaxy through its H α emission using Fabry-Perot in-terferometer. However, they did not cover the region where we detect the kinematic warp. Instead, they reported the detection of HVCs and localized high residual velocity dispersion regions, associating them to shell/chimney regions. Even though we present a velocity dispersion map in panel c), in this work we have not attempted to recover these regions, typically associated to vertical motions of ionized gas due to starforming processes.
In Figure 10 we show NGC 2763. This low inclination disk galaxy (i ≈ 29.5 • ) is the only object in WiNDS observed in the Southern Hemisphere, and one of the first galaxies observed with the SAM-FP instrument (Mendes de Oliveira et al. 2017). The B-band image on panel a) reveals two main spiral arms, each with multiple sub-arms, and a small bar. In addition, it shows a significant lopsided structure. The galaxy shows an extended H α emission across the disk, reaching up to 0.9 R opt , with R opt ≈ 9.9 kpc. Panel b) shows its V los maps, with an amplitude of ≈ 100 km s −1 . The velocity dispersion map, normalized to σ 95th = 37.2 km s −1 is shown on panel c).
The resulting V res maps, shown on panel d), reveals a very complex structure with global residual velocity flows reach-  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel.
ing amplitudes 20 km s −1 . In panel e) we can observe the B-band image of NGC 2763 with the residual velocity map overlapped, highlighting the region where complex residual velocity perturbations are observed. Interestingly, no detailed kinematical studies of NGC 2763 velocity field, nor companion galaxies within 150 kpc, have been reported. We have confirmed the lack of massive nearby companions by performing a systematic search within a radial projected distance of 250 kpc and radial velocity difference ∆ ≤ 1000 km s −1 using NASA/IPAC Extragalactic Data base. Thus, NGC 2763 is considered as an isolated galaxy. As such, it is not clear whether its significantly perturbed morphology and kinematics are the results of a previous minor merger event or due to significant misaligned smooth gas accretion.
NGC 5713 is oriented nearly face-on, with inclination i ≈ 10 • , and shows a significant lopsided morphology. This can be seen in panel a) of Figure 11, where we show its XDSS Bband image. NGC5713 has an R opt ∼ 13.5 kpc, it is a barred and multiple arm spiral galaxy (SABb-type de Vaucouleurs et al. 1991), and is located at 33.4 Mpc. The H α V los map, shown on panel b) covers approximately a region of 0.9 R opt , and shows an amplitude of ∼ 100 km s −1 . The velocity dispersion map, normalized to σ 95th = 54.2 km s −1 , is shown on panel c).
The V res (panel d) reveals a very peculiar structure, with several global and coherent perturbations reaching peak velocities ≥ 30 km s −1 . The galaxy is currently undergoing a very strong tidal interaction with its similar mass Sab-type companion, NGC 5719 (Vergani et al. 2007), and shows two HI tidal tails. Both galaxies (NGC 5713 and NGC 5719) show an optical radius at a similar galactocentric distance and are connected by two HI tidal bridges. The projected distance between these two galaxies is ∼ 77 kpc. NGC 5719 shows an ionized gas disk that counter rotates with respect to its neutral gas counterpart. According to Vergani et al. (2007), this counter-rotating disk is the result of the accretion of NGC 5713 material onto NGC 5719 during their closest approach. Interestingly NGC 5719 shows a nearly edge-on configuration with a clearly visible warped disk. In addition, NGC 5713 is interacting with the Sm-type dwarf galaxy, PGC 135857, which is located at a projected distance of ∼ 64 kpc (Vergani et al. 2007). These interactions are likely behind the observed  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel.
NGC 5713 strongly perturbed velocity field. The interacting system will be compared against numerical models in a follow-up study.

DISCUSSION AND CONCLUSIONS
In this study, we have presented the Waves in Nearby Disk galaxies Survey (WiNDS), which currently consists of 40 nearby low-inclination spiral galaxies, including objects with diverse morphological features. The WiNDS sample contains galaxies within a distance limit of 45 Mpc and absolute B magnitude between -17.0 and -21.9 mag, all observed in H α with Fabry-Perot interferometer. These observations allow us to obtain very detailed velocity maps of the ionized gas distribution, with velocity sampling precision as low as 3 km s −1 , which translates into a resolution of 6 km s −1 . Our sample was selected such that galaxies have an inclination angle i 40 • . WiNDS consists of 12 new data observed, and includes archival data for 28 low-inclination late-type galaxies, extracted from the GHASP, SINGS-H α , VIRGO-H α , HRS-H α surveys. For each galaxy in the final sample, we derived their line-of-sight velocity field, an unperturbed axy-symmetric velocity model, kinematical parameters, the corresponding residual velocity fields, V res and rotation curves (see Appendix C.1).
Our main goal in this work was to examine the resulting V res fields to search for evidence of large, global, and coherent kinematic perturbations in our sample of low-inclined late-type galaxies. In particular, we seek perturbations that are consistent with those produced by warps and corrugation patterns. Through three criteria we identify clear signatures of perturbations. First, we are interested in objects that present a wide H α coverage (≥ 0.7R opt ) so that the kinematics of the disks could be globally explored. Second, we focused on galaxies with V res amplitude that are > 10 km s −1 , which allows us to be more confident in observing perturbations in the low-inclination disks that are not the result of the axisymmetric components of the galaxies. Finally, we also searched for perturbations in the V res fields that show global and coherent velocity flows, avoiding local and discrete perturbations that could be linked to, e.g., fountain flows.
We emphasize that our selection criteria cannot confirm nor rule out the presence of vertical perturbations in our disks. In-  (e). The red (blue) area corresponds to velocities above (below) the midplane of the galaxy. The derived kinematic inclination is indicated in the lower right of each panel. stead, our goal is to select potentially strong candidates for follow-up studies. Note however that, as shown by Monari et al. (2016) and Gómez et al. (2021), the amplitude of the observed velocity flows in these low-inclined disks are, typically, much too large to be driven either by spiral structure or by a bar. Within the subset of galaxies with vertical perturbations considered as candidates, those displaying discrete and local perturbations that could be associated to fountain flows are not excluded since such signatures can coexist with the global and coherent perturbations that we aim at detecting in this study.
A vertical perturbation such as a warp or a corrugation pattern is manifested by an extended and oscillatory vertical displacement of the disk with respect to the overall mid-plane. Its characterization can allow us to constrain the recent interaction history of a galaxy. Previous studies based on cosmological simulations (e.g. Gómez et al. 2017) have identified four main mechanisms behind such perturbations: distant fly-by massive companions, close satellite encounters, re-accretion of cold gas from progenitors of a gas-rich major merger and accretion of misaligned cold gas. Out of 40 observed low inclination galaxies, we find that approximately 20 percent of the sample shows clear signatures of global and coherent perturbations. Adding to this statistic the results of our pilot study based on VV304a (Gómez et al. 2021) brings this number to 22 percent. We note however that, as further discussed below, some of our observations were either not sufficient or not adequate to identify velocity perturbation on the disk outskirts and, thus, this percentage could be higher. In fact, out of the 40 galaxies in the WiNDS sample, only 70% show good H α coverage. If we only consider this subset, the fraction of vertically perturbed galaxy candidates rises to 32%.
From the galaxies that show signs of potential vertical perturbation, six of them are interacting with a satellite galaxy and/or belong to a group. Thus, the observed kinematic perturbations are likely directly linked with recent environmental interactions. The rest of the perturbed subsample of galaxies are considered to be isolated. Therefore, their perturbations could be the result of previous minor mergers or misaligned cold gas accretion. A deeper and more detailed study must be performed for each galaxy to constrain the origin of their perturbed velocity field. As shown in Fig. 1, the vertically per-turbed candidates within WiNDS sample (green distribution) shows no preferential distribution of the main parameters with respect to the overall WiNDS sample.
Previous studies, either in HI (García-Ruiz et al. 2002) or in the optical (Reshetnikov & Combes 1998;Ann & Park 2006, for the edge-on galaxies), that have characterized the frequency which vertically perturbed disks arise in the local Universe, have found that approximately 70% of them present a warped disk, typically displaying S-shaped configuration. However, evidence for more complex corrugation patterns in external galaxies, such as those hinted in this work, have been extremely scarce. Indeed, only a handful of previous studies have previously reported corrugations on external galaxies (e.g. Sánchez-Gil et al. 2015;Narayan et al. 2020;Gómez et al. 2021). In agreement with observations, cosmological simulations of high-resolution late-type galaxies within Milky Way-sized halos, Gómez et al. (2016Gómez et al. ( , 2017) estimated that 70% should show strongly vertically perturbed disks, and 35% should present a corrugated structure. The scarcity in the detection of corrugation patterns on external galaxies so far was not that surprising since most studies attempting to characterize stellar disk vertical structures have been focused on edge-on systems in which corrugation is difficult to detect due to projection effects. Our study shows that high resolution velocity maps, obtained with techniques such as Fabry-Perot interferometer, allows us to reveal much more complex velocity structures on external galactic disks. Furthermore, due to the low inclination angles of the disks in our sample, these are likely linked with vertical flows and consistent with corrugations patterns. Eight galaxies from WiNDS show potential vertical patterns, of which 25% are in a close interaction with a massive companion, 50% show nearby dwarf satellite galaxies, and 25% can be regarded as isolated.
However, the fraction of low inclination galaxies with detected velocity flow, i.e. ≈ 20%, is currently significantly lower than the expected fraction of vertically perturbed disks. Several factors can be playing a role here. The first and most obvious reason is that a fraction of ≈ 30% of the galaxies in our sample did not show a wide H α coverage. Such cases either presented a patchy and poor H α emission distribution, or it was concentrated within the inner ∼ 0.5 R opt . As a result, we were not able to globally explore the kinematics of their disks, especially on the outer disk regions where warps and vertical perturbations are expected to be stronger. Second, the weather conditions in ≈ 12.5% of our observations were not ideal. As a result, the observations were not sufficient to generate kinematic maps accurate enough to properly resolve their velocity fields. Considering these caveats, it is not surprising that we have been able to detect kinematic perturbation, consistent with a vertical perturbed disk, in a smaller fraction of galaxies than previously reported.
In a follow-up study deeper and better observation will be presented. Following Gómez et al. (2021) we will compare our observations against detailed analytic kinematic models that account for the axisymmetric perturbations measured on each galaxy and thus constrain the origin of their velocity perturbations. NGC 1058 is a nearly normal Sc face-on spiral galaxy with an inclination angle of 6 • . NGC 1058 does not present well-defined arms in H α emission and reveals a weaker emission in the innermost than in outer regions. The galaxy is located at a distance of 9.93 Mpc, the most recent work distance value is considered. It has an angular diameter of 3 arcmin with apparent magnitude in B band about 12 mag. NGC 1058 belongs to the group NGC 1023 being the least bright member. This galaxy was observed in H α line (Sánchez-Gil et al. 2015) to analyze if it presents vertical flows using the long-slit technique using different P.A. and finding that velocity peaks are associated with star-forming regions. Two supernovae 1961V and 1969L have been reported in the outer disk.
NGC 2500 is a nearby spiral galaxy with a short bar aligned with its minor kinematical axis, classified as SB(rs)d, located at 9.79 kpc and inclination angle 40 • . This galaxy was observed using H α filter by Epinat et al. (2008b) who find diffuse emission in their observed H α maps. NGC 2500's corrugated velocity patterns were studied using long-slits with H α observations in Sánchez-Gil et al. (2015), where the vertical displacements do not seem to be related with an H α emission peak. Like NGC 1058, this galaxy is part of the sample by Sánchez-Gil et al. (2015) but unlike in NGC 1058 where they find speed peaks, in NGC 2500 it was not conclusive.
NGC 3147 is a non-barred galaxy and it is classified as SA(rs)bc. It is considered as the best type 2 Seyfert candidates using optical and X-ray observations, simultaneously (Bianchi et al. 2012).
NGC 3184 nearby grand design spiral SABc type. Located at a distance of 10.05 Mpc, it has a nearly face-on orientation with inclination angle 16.7 • and 184.5 • deviated in this article using the methodology described previously. This galaxy belongs to the group NGC 3184 being the brightest with 9.7 mag. Other members are NGC 3198, NGC 3432, and NGC 3319.
NGC 3423 nearby spiral galaxy Sc type, absolute B-band magnitude of 12 mag and angular size 3.56 arcmin and is located at 15 Mp. No previous analysis studies in particular of this galaxy.
NGC 3485 nearby spiral galaxy located at 28.17 Mpc with absolute B-band magnitude 12.7 mag . The kinematical inclination derived in this work is 26 • . There are no individual works for this galaxy previously. In this work, the H α emission is faint and is it located in the galaxy arm.
NGC 3642 Non-barred galaxy with three rings and dusty spiral arms classified SAbc-type. The galaxy shows a warped outer disk in HI line (Verdes-Montenegro et al. 2002). NGC 3642 belongs to a group composed of 5 galaxies. NGC 3642 is the brightest member with elliptical galaxies NGC 3610 and NGC 3613 with similar magnitude. Another member is NGC 674 and NGC 3683 are fainter by more than 2 magnitudes. It is an extended galaxy and shows around internal disk from which a spinal arm appears that forms an external disk NGC 4136 nearby face-on galaxy with a bar and a ring as well as a well-developed spiral structure in its outer disk with low-luminosity, classified SBc type (Gusev et al. 2003).
NGC 4900 is a spiral barred galaxy, classified as SB(rs)c located in the Virgo cluster at a distance of 13.3 Mpc and at about 12 • south-east of M87. NGC 4900 presents HII nucleus using CO observations (Lee et al. 1998). To highlight the larger velocity amplitudes in the V res map, we use a filter on the residual field through the Fourier Transform and Gaussian low-pass filter. The image filtering process is carried out following the process described in Fig. 30.
Let be a spatial image f(x,y) with dimension N × N, the Discrete Fourier Transform (DFT) of f called F(u,v) is defined as  where D(u,v) corresponds to the Euclidean distance from (u,v) to the origin of the frequency plane and D 0 is cutoff frequency in pixel in Fourier space.