Finite element electrothermal modeling is employed to study shading-induced failure in monolithically integrated thin-film photovoltaic modules. A key element is spatial nonuniformity in current-voltage characteristics, which causes inhomogeneous current flow when part of a module is under reverse bias due to shading. Time-dependent calculations show that spots with lower reverse breakdown voltage experience greater current density and localized thermal runaway that can cause permanent damage, resulting in ohmic shunts. Such failure events lead to performance loss, especially when the module is in normal operating conditions and shunted currents lead to abnormally high module temperatures. Our simulations of temperature distributions, current-voltage characteristics, and electroluminescence are compared to data from the literature for copper indium gallium diselenide devices. (C) 2016 Elsevier Ltd. All rights reserved.