Thermal Response Testing is now a well-known and widely-used method allowing the determination of the local thermal or geometrical properties of aBorehole Heat Exchanger (BHE), those properties being critical in the design of GSHP systems. The analysis of TRTs is an inverse problem that hascommonly been solved using an approximation of the ILS solution. To do this, however, the heat rate during a TRT must be kept constant, or least be nontime-correlated, during the test, which is a challenging constraint. Applying temporal superposition to the ILS model is a way to account for varying power,although it requires the use of an optimization algorithm to minimize the error between a parametrized model and experimental values.In this paper, the Newton-Raphson method is applied to the time-superposed ILS for parameter estimation in TRTs. The parameter estimation is limitedto the effective thermal conductivity and the effective borehole resistance. Analytical expressions of the first and second derivatives of the objective function,chosen as the sum of quadratic differences, are proposed, allowing to readily inverse of the Hessian matrix and speed the convergence process.The method is tried for 9 different TRTs, 2 of which are reference datasets used for validation of the method (Beier et al., 2010). Differences betweenestimated and reference thermal conductivities are of 3.4% and 0.4% for the first and second reference TRTs, respectively. The method is shown to be stableand consistent: for each of the 9 TRTs, 11 realizations are performed with different initial values. Convergence is reached in all cases and all realizationslead to the same final values for a given TRT.The proposed convergence method is about 70% to 90% faster than the Nelder-Mead simplex and require about 8 times less iterations in average. Theconvergence speed varies between 0.3 to 13.6 s with an average of 3.7 s for all TRTs.