Gradients in nutrient distributions and the effects of eutrophication are common features in most coastal marine areas. These structures occur in aquatic systems due to spatial differences in hydrography, nutrient loading and key biogeochemical processes. Two statistical methods, cluster analysis and probability mapping, have been used in the present study to determine characteristics and patterns in water quality variables. Combined, these two methods provide a useful tool to statistically determine spatial homogeneity and representativity of areas and stations. A case study is presented here in which water quality variables (salinity, dissolved inorganic nitrogen, dissolved inorganic phosphorus and chlorophyll) in surface waters of the Kattegat are analysed for the time period 1993-1997. It was found that morphology, the proximity to sources of nutrient loading from land, nutrient uptake and the infrequent water exchange between the Baltic Sea in the south and the Skagerrak in the north all contribute to create distinct regions of water quality. Regions with concentrations significantly different from the overall mean are identified. In addition to identification of regions with similar characteristics, representative stations (as well as not representative stations) for the respective regions were made. This type of information can be used to design new or revise old monitoring programs.