Two types of signals are clearly visible in continuous GPS (cGPS) time-series in Iceland, in particular in the vertical component. The first one is a yearly seasonal cycle, usually sinusoid-like with a minimum in the spring and a maximum in the fall. The second one is a trend of uplift, with higher values the closer the cGPS stations are to the centre of Iceland and ice caps. Here, we study the seasonal cycle signal by deriving its average at 71 GPS sites in Iceland. We estimate the annual and semi-annual components of the cycle in their horizontal and vertical components using a least-squares adjustment. The peak-to-peak amplitude of the cycle of the vertical component at the studied sites ranges from 4 mm near the coastline up to 27 mm at the centre of the Vatnajökull, the largest ice cap in Iceland. The minimum of the seasonal cycle occurs earlier in low lying areas than in the central part of Iceland, consistent with snow load having a large influence on seasonal deformation. Modelling shows that the seasonal cycle is well explained by accounting for elastically induced surface displacements due to snow, atmosphere, reservoir lake and ocean variations. Model displacement fields are derived considering surface loads on a multilayered isotropic spherical Earth. Through forward and inverse modelling, we were able to reproduce a priori information on the average seasonal cycle of known loads (atmosphere, snow in non-glaciated areas and lake reservoir) and get an estimation of other loads (glacier mass balance and ocean). The seasonal glacier mass balance cycle in glaciated areas and snow load in non-glaciated areas are the main contributions to the seasonal deformation. For these loads, induced seasonal vertical displacements range from a few millimetres far from the loads in Iceland, to more than 20 mm at their centres. Lake reservoir load also has to be taken into account on local scale as it can generate up to 20 mm of vertical deformation. Atmosphere load and ocean load are observable and generate vertical displacements in the order of a few millimetres. Inversion results also shows that the Iceland crust is less rigid than the world average. Interannual deviation from the GPS seasonal cycle can occur and are caused by unusual weather conditions over extended period of time.