Basel cell carcinoma (BCC) is the most common human skin malignancy. Its incidence has increased significantly in Australia, Europe and North America over the past decade. A number of modalities are currently used for treatment of BCC, including cryosurgery which offers a potential for high cure rate, low cost, minimal bleeding and good cosmetic effect. However, cryosurgery is not used frequently for BCC because no current method exists to design adequate treatment parameters. We present a numerical analysis on the thermal history of the target tissue during cryosurgery of a nodular BCC using liquid nitrogen (LN2) spray. The model uses Pennes equation to describe the heat transfer within the target tissue. A convective thermal boundary is used to describe the heat interaction between the tissue and LN2, and the apparent heat capacity method is applied to address the tissue phase change process. A parametric study is conducted on the convective heat transfer coefficient (hs: 104~106 W/m2K), cooling site area (rs/R0: 0.5~1.0) and spray time (t: 0~30 sec.), with the objective to understand the thermal history during tissue freezing, including lethal temperature (-50 C) and cooling rate (CR). Results demonstrate that propagation of the lethal isotherm is sensitive to the convective heat transfer coefficient, hs, with a range of 104~5104 W/m2K. Increasing the cooling site area can significantly enhance cooling efficiency, producing dramatic increase in the amount of tissue encompassed by the lethal isotherm. The cooling rate (CR) shows a highly dynamic distribution during the cooling process: the highest CR drops quickly from 140 C/sec. (t=0.5 sec.) to 20 C/sec. (t=5 sec.). The highest CR is initially located close to the cooling site but moves toward the inside of the tissue as treatment proceeds. The model presented herein provides a simulation tool for treatment planning of cryosurgery using LN2 spray, in which the protocol parameters, e.g. cooling site area and spray time, can be determined for an optimal outcome. The quantitative predictions on the propagation of lethal isotherm and the distribution of CR should help to optimize cryosurgery efficacy.