In this paper a 3-D mesoscale PBL dynamic prognostic model for a coastal region with complex terrain was developed to simulate the 3-D flow field under a local sea-land breeze circulation. The output from the PBL model was used as an input to an Eulerian numerical model which can be used to simulate finely the temporal and spatial distributions of an air pollutant (SO,) during a sea-land breeze developing in a coastal region. With respect to its use as a diagnostic model, only a few data would be required to simulate the background winds controlled by a larger scale synoptic system, and then provide initial winds for the PBL model. Having linked the three models and defined the coefficients of turbulent diffusion in a simple form, an integrated 3-D numerical air quality model suitable for the coastal environment was designed. The period of May 29, 1986 was selected for simulating and analysing the distribution of air pollutants over the coastal area of Bohai-sea in Northern China. The results indicated that the calculated concentrations corresponded with the observed ones on the whole. Thus this linked model has been shown to be feasible and useful in practice.