Tea color is a key factor in consumer purchasing decisions and reflects product quality. Thus, modeling and controlling the soursop slices drying process are essential to enhance sensory value and reduce investigation costs. In this study, mathematical techniques were applied to evaluate the changes in color values (L*, a*, b*, triangle E*, and BI) during heat pump drying using a high-capacity pilot drying equipment. Experimental data were collected from assessing color changes at drying temperatures of 20 degrees C, 30 degrees C, 40 degrees C, and 50 degrees C over time. The data were fitted into regression equations to assess the kinetics of the process. Results revealed that each kinetic model could appropriately describe the color value changes for one or more drying processes at different temperatures. Notably, both fractional conversion and Zero-order models were chosen to describe the change in L* value during drying at 20 degrees C and 30 degrees C, while the Power-law model was used for L* value changes at 40 degrees C and 50 degrees C. At the same temperature ranges, another finding showed that the experimental data for a* and b* values followed the Combined kinetic model (R-2 > 0.92). Kinetic models describing the changes in BI and Delta E* during drying at various temperatures were also identified based on statistical parameters (R-2, chi(2), and RMSE). The activation energy required for changes in color values was calculated, revealing that the a* value was the most sensitive to temperature, whereas the sensitivities of L*, b*, triangle E*, and BI were comparable. The study discussed the nature, modeling methods, and optimization of drying processes, offering insights into both quantity and quality.