“…Recently, this problem has been applied for irrigation systems design, see (Haghverdi et al, 2015) where linear programming is used as one of the methods for delineating management zones among others as Kmeans and Isodata, methods that are still being used. These other methods are classified as clustering methods, see (Ortega et al,2002), (Jaynes et al,2005) and (Jiang et al (2011), but their major drawback is the resulting fragmentation of the zones, because these methods generate oval shaped and disjoint zones. Although the problem of defining management zones in presence of site specific variability has been studied in previous works, to the best of our knowledge, an important characteristic that has not been considered yet is variability in time of the chosen soil property.…”