We developed a numerical model, RFLUX, which uses the heat tracer method for vertical groundwater flux estimation, and applied it to the Leizhou Peninsula, South China, to provide information to inform local groundwater resource utilization and management. The temperature-depth (TD) profiles of 24 boreholes, along with the observed ground surface temperature (GST) and surface air temperature (SAT) series in recent decades, were collected in this area. Underground TD data demonstrated the capacity to identify groundwater flow patterns, and local GST and SAT data demonstrated a strong correlation with each other over monthly, seasonal, and annual scales. In the RFLUX model, the average GST and SAT series were applied as an upper boundary condition, and a nonlinear initial condition was set using an analytical solution from the literature. The model results of selected TD profiles demonstrated that the annual vertical groundwater flux was about 0.15 m a −1 , which tended to be overestimated if a linear initial condition was used. This model can be easily applied with minor modifications, considering its clear purpose and simplicity.