On the other hand, geophysical fluid mechanics (meteorology, oceanography) is oriented more towards the comprehension and prediction of the relevant processes at broader temporal and spatial scales. Environmental fluid mechanics targets its major concern somewhere between these two extremes with the aim of assessing the potential environmental hazard impact and helping in decision-making processes for proposed project solutions (Cushman-Roisin et al. 2008). Dominant forcing and its intensities in the mixing processes affecting the effluent plume on its path from the diffuser

orifice to the arbitrary downstream profile are highly variable. Therefore, the concept of separating the far-field and near-field zones with different dominant forcing is widely adopted (Fischer et al. 1979). In the case of a submarine public sewage outfall, the near-field domain in the vicinity of the outfall see more diffuser ranges from the inflow point up to the sea AC220 order surface or the neutral buoyancy layer, where further effluent plume rise is interrupted, after which plume dynamics is mainly in the horizontal direction (Akar & Jirka 1994a,b). Therefore, the integral solution of the problem is usually obtained through the combination of two structurally different numerical

models. Plume propagation in the far field is modelled with a 2D or 3D oceanographic numerical model using initial concentration fields calculated from the near-field model (Wood et al. 1993, Akar & Jirka 1994a,b, Pun & Davidson 1999). Using this approach medroxyprogesterone one can avoid a high-resolution numerical grid within the far-field model required for resolving the near-field mixing process. In this study we have slightly modified the previously described methodology in order to assess the influence of bora-induced density changes on effluent

plume dynamics. This approach consists of two steps: a) temporal changes in the vertical density distribution along the water column at the positions of the analysed submarine outfall diffusers are obtained from 3D numerical model simulations; b) mixing processes in the near-field are resolved using a numerical model constructed according to Featherstone (1984), with previously calculated vertical density profiles. More details on the nearfield numerical model used are given in section 3. In the Rijeka Bay area (Figure 1) tourist activities and the bathing season are at their height in the summer, at the end of June and during the first half of July. Many projects and construction activities related to the implementation of municipal sewage systems as well their improvements are currently in progress. Three new submarine outfalls L, O, MNJ (Figure 1) are envisaged for construction, and plans are in hand to extend the one already in existence R (Figure 1). The basic hydraulic characteristics of these four outfalls are given in Table 1 (Andročec et al. 2009). The hydraulics of the diffusers were taken into account according to the methodology of Fischer et al.