Processing data
1. Fluorescence data were measured on a Horiba AquaLog
- Fluorescence was measured as signal/reference beam.
- Internal excitation- and emission-correction factors were applied
- Blanks (ultrapure water) with identical instrument settings were subtracted
2. Inner filter effects were eliminated with the absorbance-based method 
3. Fluorescence signals were normalized using the Raman peak area of ultrapure water at 350nm.
4. First-order Rayleigh scatter bands were removed and replaced with missing numbers and Raman scatter bands and the second-order Rayleigh scatter band were replaced with interpolated data. 
5. Fluorescence emission wavelengths were restricted to the range 300-620 nm, and excitation was restricted to the range 250-620 nm.
Model fitting
1. Start: Processed fluorescence data above.
2. Model were obtained with the following options
- Starts: 100
- Convergence: 1e-6
- Constraints: nonnegativity in all modes
- Initialization: random orthogonolized numbers

About the data
1. Each *.csv file contains the wavelength information in nm as column and row headers (columns=excitation, row=emission) and the fluorescence observations processed as described above.
2. Each file is named with its sample identifier.
3. The files "metadata_wastewater_for model.csv" contain information on the sample type and all other useful information for the description of the samples detailed in the main text.
4. The *.ods (OpenDocument Spreadsheet) file contains two spreadsheets with the loadings of the PARAFAC model.
5. The zip-file contains  *.mat (MATLAB) files that hold data equivalent to the *.ods-files in a format that is compatible to the drEEM toolbox (dreem.openfluor.org)
