The strain energy release rates of adhesively-bonded pultruded GFRP joints were determined experimentally. The crack propagated in the adherend along paths outside the symmetry plane accompanied by fiber bridging. A new method, designated the "extended global method", was introduced to facilitate mode partitioning in the mixed-mode experiments. Non-linear finite element models were developed in order to quantify the effect of the observed fiber bridging on crack propagation. An exponential traction-separation cohesive law was used to model the fiber bridging zone and calculate the energy release rate due to the fiber bridging, while the virtual crack closure technique was used for calculation of the fracture components at the crack tip. Experimental, analytical and numerical analyses were used to establish quasi-static mixed-mode failure criteria for crack initiation and propagation. The derived mixed-mode failure criteria can be used for simulating progressive crack propagation in other joint configurations comprising the same adhesive and adherends. (C) 2014 Elsevier Ltd. All rights reserved.