Understanding the spatial heterogeneity within environmental biofilms can provide an insight into compartmentalization of different functions in biofilm communities. We used a non-destructive and label-free method by combining Synchrotron Radiation-based Fourier Transform Infrared Microspectroscopy (SR-FTIR) with Confocal Laser Scanning Microscopy (CLSM) to distinguish the spatial chemical changes within multispecies biofilms grown from natural storm waters in flow cells. Among the different surfaces tested for biofilm growth and optimal imaging, mylar membranes were most suited and it enabled successful spatial infrared imaging of natural biofilms for obtaining reliable and interpretable FTIR spectra. Time series analysis of biofilm growth showed that influx of water during biofilm growth, results in significant changes in biofilm formation. Early biofilms showed active nutrient acquisition and desiccation tolerance mechanisms corresponding with accumulation of secreted proteins. Statistical approach used for the evaluation of chemical spectra allowed for clustering and classification of various regions of the biofilm. Microheterogeneity was observed in the polymeric components of the biofilm matrix, including cellulose, glycocalyx and dextran-like molecules. Fructan and glycan-rich regions were distinguishable and glycocalyx was abundant in the strongly adhering peripheral regions of biofilms. Inner core showed coexistence of oxygen dimers and ferrihydrite that will likely support growth of Fe (II)-oxidising bacteria. The combined SR-FTIR microspectroscopy and CSLM approach for complex natural biofilms described here will be useful both in understanding heterogeneity of matrix components and in correlating functions of juxtaposed microbial species in complex natural biofilms with physicochemical microenvironment to which they are exposed.