Retrieving water quality variables from remotely sensed reflectance spectra (Rrs) with global-scale, cross-waterbody applicability remains a grand challenge in environmental remote sensing. Here we introduce Hierarchical Spectral Ensemble (HSE), a machine learning pipeline for concurrent retrieval of four ecologically relevant water quality parameters (chlorophyll-a (Chl-a), total suspended solids (TSS), coloured dissolved organic matter absorption coefficient at 440 nm (aCDOM(440)), and Secchi depth (Zsd)), applied to the GLORIA 2022 global dataset of 7,572 co-located hyperspectral in situ ground truth measurements spanning six continents. HSE leverages a combination of five diverse base learners trained with 10-fold out-of-fold stacking, and integrated with a diversity-constrained non-negative least squares (NNLS) meta-learner, supported by a 291-dimensional feature suite spanning spectral, spatial, and temporal representations. We propose Interpolation-based Spectral Data Augmentation (ISDA), an ecologically inspired, class-balance oversampling methodology applied to the training set after the initial stratified split to prevent information leakage. Evaluated on the held-out global test set, HSE attains R2 = 0.822, 0.704, 0.841, and 0.962 for Chl-a, TSS, aCDOM(440), and Zsd respectively (mean R2 = 0.832; all in original physical units following back-transformation). Stratified analysis on a per-quartile basis reveals negative values of R2 in low-concentration regimes for Chl-a, TSS, and aCDOM(440), illustrating how global metrics can mask failures in retrieval across a majority of the distribution, specifically in the oligotrophic regime. This constitutes the primary known limitation of the method: reliable prediction in low-concentration, low-optical-signal regimes remains elusive and is demonstrated only in high-concentration, eutrophic samples where a measurable optical signal exists.