Long-term monitoring of Water Constituent Concentrations (WCCs) is essential for water quality assessment in coastal waters, a priority for various government agencies and environmental organizations. However, utilizing satellite observations to track the extended spatial and temporal variations of WCCs, including Chlorophyll-a (Chla), Suspended Particulate Matters (SPMs), and Coloured Dissolved Organic Matters (CDOMs), remains challenging in coastal waters due to their optical complexity and the need for accurate atmospheric correction. In this study, we examined two decades of spatial and temporal variations in Chla [mg m-3], SPM [g m-3] concentrations, and CDOM absorption at 440nm [m-1] using time series data from the Medium Resolution Imaging Spectrometer (MERIS) and the Ocean and Land Colour Instrument (OLCI) spanning from 2003 to 2023. Our research focused on the complex, shallow, and highly turbid waters of the Dutch Wadden Sea, the Netherlands. To achieve this, we employed a coupled atmospheric and water optical model known as MOD2SEA, enabling simultaneous atmospheric correction and WCC retrievals from MERIS images (2003-2012) and OLCI images (2018-2023) in the study area. Our findings reveal that SPM and CDOM variations followed a consistent seasonal pattern over the two-decade period, with SPM values ranging from 10 [g m-3] to 60 [g m-3] and CDOM absorption varying from 0 [m-1] to 0.8 [m-1]. Chla estimates, on the other hand, displayed a different trend, ranging from 5 [mg m-3] to 25 [mg m-3] until 2013. Subsequently, Chla concentrations increased, reaching 45 [mg m-3] by 2023.