Land-ocean interface has been a challenge in determining the boundary conditions for earth system modelling. Due to the uncertainty of material fluxes, the localization effect on counter gradient convergence, and the existence of self-adapted structures, land-ocean interface is unable to be monitored ergodically under the framework of classic linear theory. There is a great difficulty in quantifying real parameters for this interface, which is the key connection between linear and nonlinear systems under harmonic conditions, to better improve the simulation and prediction capability of earth system modelling. In a previous work, based on a two-dimensional multi-constituent tidal model and its adjoint model, together with assimilation using satellite altimeter data and tidal gauge array data, the spatial and temporal variation in bottom friction coefficient was estimated in an improved way. Inspired by the abovementioned study, this contribution proposed a suggestion of complex parameterization for land-ocean interface, in combination with remote sensing observation and numerical modelling.