Information on the coupled liquid water, water vapor, and heat transport under arable field conditions is still limited, particularly for unsaturated soils of semi-arid and arid regions such as New Mexico. HYDRUS-1D model was applied to evaluate various transport mechanisms associated with temporal variations in water content and soil temperature in the unsaturated zone of a sandy loam furrow-irrigated field located at Leyendecker Plant Science Research Center, Las Cruces, New Mexico. The model was calibrated using measured soil water content and soil temperature at 5-, 10-, 20-, and 50-cm depths during a 19-day period (day of the year [DOY] 85 to DOY 103, 2009) and validated for a 31-day period (DOY 104 to DOY 134, 2009). Measured and optimized soil hydraulic and thermal properties and hourly meteorological data were used in model simulations. HYDRUS-1D simulated water contents and soil temperatures correlated well with the measured data at each depth. The total liquid water flux, composed of isothermal and thermal liquid water flux, dominated the soil water movement during early periods after irrigation, whereas the contribution of total water vapor flux, composed of primarily thermal and much smaller isothermal water vapor flux, increased with increasing soil drying. During the soil drying process, the total liquid flux within 15-cm depth eventually changed to water vapor flux near the surface. The upward total liquid and vapor fluxes decreased from 5 cm, indicating that vapor flux was much higher in the layer near the soil surface. The total vapor flux in this unsaturated soil layer was approximately 10.4% of the total liquid and vapor fluxes during the simulation period.