The HYDRUS-1D and HYDRUS (2D/3D) computer software packages are widely used finite-element models for simulating the one-and two-or three-dimensional movement of water, heat, and multiple solutes in variably saturated media, respectively. In 2008, Šimůnek et al. (2008b) described the entire history of the development of the various HYDRUS programs and related models and tools such as STANMOD, RETC, ROSETTA, UNSODA, UNSATCHEM, HP1, and others. The objective of this manuscript is to review selected capabilities of HYDRUS that have been implemented since 2008. Our review is not limited to listing additional processes that were implemented in the standard computational modules, but also describes many new standard and nonstandard specialized add-on modules that significantly expanded the capabilities of the two software packages. We also review additional capabilities that have been incorporated into the graphical user interface (GUI) that supports the use of HYDRUS (2D/3D). Another objective of this manuscript is to review selected applications of the HYDRUS models such as evaluation of various irrigation schemes, evaluation of the effects of plant water uptake on groundwater recharge, assessing the transport of particle-like substances in the subsurface, and using the models in conjunction with various geophysical methods.Abbreviations: CRS, cosmic-ray sensing; EC, electrical conductivity; ERT, electrical resistivity tomography; FEM, finite element mesh; GPR, ground-penetrating radar; GUI, graphical user interface; TDT, time-domain transmissometry.The HYDRUS-1D and HYDRUS (2D/3D) software packages (Šimůnek et al., 2008b) are finite-element models for simulating the one-and two-or three-dimensional movement of water, heat, and multiple solutes in variably saturated media, respectively. The standard versions, as well as various specialized add-on modules, of the HYDRUS programs numerically solve the Richards equation for saturated-unsaturated water flow and convection-dispersion type equations for heat and solute transport. The flow equation incorporates a sink term to account for water uptake by plant roots as a function of water and salinity stress. Both compensated and uncompensated water uptake by roots can be considered. The heat transport equation considers movement by both conduction and convection with flowing water. The governing convection-dispersion solute transport equations are written in a relatively general form by including provisions for nonlinear, nonequilibrium reactions between the solid and liquid phases and linear equilibrium reactions between the liquid and gaseous phases. The transport models also account for convection and dispersion in the liquid phase as well as diffusion in the gas phase, thus permitting the models to simulate solute transport simultaneously in both the liquid and gaseous phases. Hence, both adsorbed and volatile solutes, such as pesticides and fumigants, can be considered.The solute transport equations further incorporate the effects of zero-order production, first-o...