Ocean waves have long been a research topic, and numerous formulas of the ocean wave spectrum have been widely developed to provide new prospects for ocean experiments and the advancement of radar probing. Nonetheless, the wave spectra developed by researchers fall short of the standards set by remote sensing specialists, mainly due to the limitation of capturing the surface roughness influenced by both short and long waves, prompting ongoing efforts to develop a model covering a diverse scale of wavenumbers in the absence of a generally recognized reference formula. In response, a standard two-scale formulation of wave frequency spectra (WFS) was introduced, where wave height and spectral peak period are determined by sea state. The proposed composite WFS holds the potential to incorporate directional spreading, contributing to the angular distribution of ocean wave energy in the form of directional WFS, making it applicable for ocean modeling. In an effort to investigate directionality effects, an array of well-established spreading functions, including cosine-squared, half-cosine 2s, parameterized half-cosine 2s, hyperbolic secant-squared, and composite structured functions, has been developed here for numerical modeling of random ocean media (ROM) surface roughness and synthesized for their spectral scattering distribution mode, encompassing scattering pattern, scattering orientation, and fractal roughness properties. Nonetheless, the generated ROM models, each varying due to inherent limitations in directional WFS formulation and numerical approximations, demonstrate indeterminacy and unpredictability in surface features, posing challenges for accurately synthesizing ROM roughness patterns. These challenges intensify under varying sea states and different directional WFS formulas, leading to a situation where no single synthesized composite ROM model consistently outperforms the others, rendering them imprecise frameworks for analyzing roughness patterns and investigating texture electromagnetic interactions within the realm of remote sensing. As an approach, a pattern-sensitive fusion method is proposed, employing a multi-scale transform domain (MTD) fusion scheme that leverages the learning potential of a deep super resolution network. The objective is to fuse the reconstructed ROM roughness models, generating an optimal roughness while maintaining their scattering pattern, scattering orientation, and dominant directionality, pivotal for texture consistency and, consequently, the backscattering properties from the synthetic aperture radar (SAR) viewpoint. To validate the reliability of ROM modeling and its roughness synthesis, including the texture fusion and raw data generation, a comprehensive objective quality assessment technique is utilized. These assessments demonstrate the complete consistency of the simulation results with the underlying spectral theory, highlighting their potential contribution to projects related to ocean radar probing and remote sensing.INDEX TERMS Directional wave frequency spe...