It arises a puzzle in \NGC\, how to secularly maintain the counter-rotating disc from $0.2$ to $7\,$pc unambiguously detected by recent ALMA observations of molecular gas. Upon further analysis of disc dynamics, we find that the Kelvin-Helmholtz (KH) instability (KHI) results in an unavoidable catastrophe of the disc developed at the interface between the reversely rotating parts, and demonstrate that a close binary of supermassive black holes provides tidal torques as the unique external sources to prevent the disc from the KH catastrophe. We are led to the inescapable conclusion that there must be a binary black hole at the center of NGC 1068, to prevent it from the KH catastrophe. The binary is composed of black holes with a separation of $0.1\,$pc from GRAVITY/VLTI observations, a total mass of $1.3\times 10^{7}\:M_{\odot}$ and a mass ratio of $\sim 0.3$ estimated from the angular momentum budge of the global system. The KHI gives rise to forming a gap without cold gas at the velocity interface which overlaps with the observed gap of hot and cold dust regions. Releases of kinematic energies from the KHI of the disc are in agreement with observed emissions in radio and $\gamma$-rays. Such a binary is shrinking with a timescale much longer than the local Hubble time via gravitational waves, however, the KHI leads to an efficient annihilation of the orbital angular momentum and speed up merge of the binary, providing a new paradigm of solving the long term issue of "final parsec problem". Future observations of GRAVITY+/VLTI are expected to be able to spatially resolve the CB-SMBHs suggested in this paper.

Close binaries of supermassive black holes (CB-SMBHs) with separations of $\lesssim 0.1$pc as the final stage of galaxy mergers are sources of low frequency gravitational waves (GW), however, they are still elusive observationally because they are not spatially resolved. Fortunately, reverberation as echoes of broad emission lines to ionizing continuum conveys invaluable information of the dynamics of broad-line regions (BLRs) governed by supermassive black holes in the central regions of active galactic nuclei (AGNs). In this paper, we demonstrate how to composite the hybrid 2-dimensional transfer functions of binary BLRs around the CB-SMBHs in AGNs, providing an opportunity of identifying them from reverberation mapping (RM) data. It is found that there are variation-coupling effects in the transfer functions, arising from the coupling of CB-SMBH light curves in the Fourier space. We provide semi-analytical formulations of the transfer functions for kinematic maps of the gas. For cases with the simplest variation-coupling effects, we make calculations for several BLR models and reveal significant distinctions from those of single active black holes. In principle, the difference is caused by the orbital motion of the CB-SMBH systems. In order to search for CB-SMBHs in time-domain space, selection of target candidates should focus on local AGNs with H$\beta$ double-peaked profiles and weaker near-infrared emission. High-fidelity RM-campaigns of monitoring the targets in future will provide opportunities to reveal these kinematic signatures of the CB-SMBHs and hence for measurements of their orbital parameters.