This paper deals with thermoelectric problems including the Peltier and Seebeck effects. The coupled elliptic and doubly quasilinear parabolic equations for the electric and heat currents are stated, respectively, accomplished with powertype boundary conditions that describe the thermal radiative effects. To verify the existence of weak solutions to this coupled problem (Theorem 1), analytical investigations for abstract multi-quasilinear elliptic-parabolic systems with nonsmooth data are presented (Theorem 2 and 3). They are essentially approximated solutions based on the Rothe method. It consists on introducing time discretized problems, establishing their existence, and then passing to the limit as the time step goes to zero. The proof of the existence of time discretized solutions relies on fixed point and compactness arguments. In this study, we establish quantitative estimates to clarify the smallness conditions. 2010 Mathematics Subject Classification. 35K51, 35R05, 35J62, 35Q79.