Gas quantification based on the response of an array of metal oxide (MOX) gas sensors in an open sampling system is a complex problem due to the highly dynamic characteristic of turbulent airflow and the slow dynamics of the MOX sensors. However, many gas related applications require to determine the gas concentration the sensors are being exposed to. Due to the chaotic nature that dominates gas dispersal, in most cases it is desirable to provide, together with an estimate of the mean concentration, an estimate of the uncertainty of the prediction. This work presents a probabilistic approach for gas quantification with an array of MOX gas sensors based on Gaussian Processes, estimating for every measurement of the sensors a posterior distribution of the concentration, from which confidence intervals can be obtained. The proposed approach has been tested with an experimental setup where an array of MOX sensors and a Photo Ionization Detector (PID), used to obtain ground truth concentration, are placed downwind with respect to the gas source. Our approach has been implemented and compared with standard gas quantification methods, demonstrating the advantages when estimating gas concentrations.