An objective definition of seasons for the Mediterranean region is performed using the mean intra-annual variations of 12 climatological parameters. The aim is to achieve a climatologically appropriate determination of the onset and cessation dates, as well as the duration, of each season, better than the conventional seasons (winter, spring, summer and autumn). Moreover, possible changes of the defined seasons during the 70-year study period 1949-2018 are investigated over the climatologically sensitive Mediterranean region. The data used are daily NCEP/NCAR Reanalysis grid point values of precipitation rate, convective precipitation rate, 2 m temperature, total cloud cover, 2 m zonal and meridional wind, 500 and 1,000 hPa geopotential height, 500 and 850 hPa temperature, 850 hPa specific humidity and precipitable water over the Mediterranean region, for the period 1949-2018. Firstly, Principal Component Analysis is applied to the mean intra-annual variations of the above parameters in order to reduce the dimensionality and then k-means Cluster Analysis is applied to the resultant Principal Components in order to group dates with similar regimes of the above parameters. This methodology is applied for the 70-year period 1949-2018 as well as for the five overlapping 30-year subperiods