Benjamin,
Interesting read. I have to agree with Crandles, I think this is due to the ice thickness/growth feedback. Thinner ice or open water generates ice more efficiently.
Furthermore I've not read of PIOMAS taking into account influx from glaciers. The PIOMAS model models ocean and sea ice response to NCEP/NCAR forcings. There is no mention in the literature I've read that glacier flow is included, nor is there mention in the gridded PIOMAS data.
However PIOMAS volume is below CryoSat 2, maybe glacier discharge is a part of this. Without figures of net discharge of ice from glaciers in N Greenland and the Canadian Arctic Archipelago I can't hazard a guess as to how much of a factor this might be.
If you need more to read about PIOMAS I recommend:
"Modeling Global Sea Ice with a Thickness and Enthalpy Distribution Model in Generalized Curvilinear Coordinates." Zhang & Rothrock 2003.
http://psc.apl.washington.edu/zhang/Pubs/POIM.pdfWhich is the base paper describing the details of the model. Frankly the maths is beyond me to delve into deeply, but sect 2a describes the model and by going through the text and the general form of the equations you can get a general feel for what goes into the model.
Schweiger at al 2011 states:
Daily mean NCEP/NCAR reanalysis data are used as atmospheric forcing, i.e., 10-m surface winds, 2-m surface air temperature (SAT), specific humidity, precipitation, evaporation, downwelling longwave radiation, sea level pressure, and cloud fraction. Cloud fraction is used to calculate downwelling shortwave radiation following Parkinson and Washington [1979].
Which runs down what the current version of the model is forced with.
It's good to see someone else using Excel!