A weakly nonhydrostatic, two-layer numerical model based on the Boussinesq equations is presented which is capable of describing, among others, the generation and propagation of nonlinear weakly dispersive internal waves in the Strait of Gibraltar. The model depends on one space coordinate only, but it retains several features of a fully three-dimensional model by including a realistic bottom profile, a variable channel width, and a trapezoidal channel cross section. The nonlinear primitive Boussinesq equations include horizontal diffusion, bottom friction, and friction between the two water layers. The model is driven by a height difference of the mean interface depth between the Atlantic and the Mediterranean boundaries and by semidiurnal tidal oscillations of the barotropic transport. The model presented in this paper describes (1) the mean and tidal flow in the Strait of Gibraltar, (2) the variation of the depth of the interface during a tidal cycle, (3) the generation of strong depressions of the interface at the western sides of the Spartel Sill and the Camarinal Sill, (4) the generation of strong eastward propagating internal bores, and (5) their disintegration into trains of internal solitary waves. The surface convergence patterns associated with depressions of the interface at the Camarinal Sill, internal bores, and internal solitary waves are calculated and compared with roughness patterns visible on synthetic aperture radar (SAR) images of the first European Remote Sensing Satellite ERS 1. In total, 155 ERS 1 SAR scenes from 94 satellite overflights over the Strait of Gibraltar, which were acquired in the period from January 1992 to March 1995, have been analyzed. It is shown that the proposed model is capable of explaining the observed temporal and spatial evolution of surface roughness patterns associated with eastward propagating internal waves inside the Strait of Gibraltar as well as the observed east-west asymmetry of the internal wave field. ¿ American Geophysical Union 1996