Many sources of uncertainty exist when emu. are modeled for a gas turbine combustion system. They originate from uncertain inputs, boundary conditions, calibration, or lack of sufficient fidelity in a model. In this paper, a nonintrusive polynomial chaos expansion (NIPCE) method is coupled with a chemical reactor network (CRN) model using Python to quantify uncertainties of NOx emission in, a premixed burner. The first objective of uncertainty quantification (UQ) in this study is development of a global sensitivity analysis method based on the NIPCE method to capture aleatory uncertainty on NOx emission due to variation of operating conditions. The second objective is uncertainty analysis (UA) of NOx emission due to uncertain Arrhenius parameters in a chemical kinetic mechanism to study epistemic uncertainty in emission modeling. A two-reactor CRN consisting of a perfectly stirred reactor (PSR) and a plug flow reactor (PFR) is constructed in this study using Cantera to model NOx emission in a benchmark premixed burner under gas turbine operating conditions. The results of uncertainty and sensitivity analysis (SA) using NIPCE based on point collocation method (PCM) are then compared with the results of advanced Monte Carlo simulation (MCS). A set of surrogate models is also developed based on the NIPCE approach and compared with the forward model in Cantera to predict NOx emissions. The results show the capability of NIPCE approach for UQ using a limited number of evaluations to develop a UQ-enabled emission prediction tool for gas turbine combustion systems.