<p>The 65 km<sup>2</sup>&#160;B&#225;rdarbunga caldera is located in the NW part of the Vatnaj&#246;kull glacier in central Iceland.&#160;&#160;The caldera floor lies under 500-800 m thick ice and the rims are fully subglacial as well.&#160;&#160;The caldera subsided by 65 m during the B&#225;rdarbunga-Holuhraun eruption in 2014-2015, when about 2 km<sup>3</sup>&#160;of magma drained out from a magma reservoir at ~10 km depth leading to the largest eruption in Iceland since Laki in 1783.&#160;&#160;Deformation surveys outside the caldera have indicated inflation since soon after the end of the eruption in February 2015 and seismicity has been elevated.&#160;&#160;The extensive ice cover precludes conventional microgravity surveys or detailed surveys of caldera floor elevation.&#160;&#160;However, we have studied gravity changes by comparing results of repeated Bouguer anomaly surveys.&#160;&#160;We perform a full Bouguer correction using detailed DEMs of both the ice surface and the ice-radar-derived bedrock.&#160;&#160;Ice surface changes are also mapped, allowing the removal of effects on gravity by ice mass changes.&#160;&#160;Possible sources of significant anomalies are either changes in bedrock elevation between surveys, other more deep-seated mass changes beneath the volcano, or changes in the water table and pore pressure.&#160;&#160;Surveys were carried out using a Scintrex CG-5 in 2015, 2016, 2018 and 2019, with measurements done at 25-50 locations each time.&#160;&#160;As no benchmarks exist on the ice the spatial difference in station location of 10-20 m exists between survey years.&#160;&#160;&#160;However, post-processing provides kinematic GPS position and elevation accuracy better than 0.1 m. Analysis of the data and error sources indicate an accuracy in estimates of changes of 50-100 &#181;Gal. The results obtained indicate change with an amplitude of a few hundred &#181;Gals; over the four years between 2015-2019 a clear Bouguer anomaly increase is recorded over the caldera relative to the surrounding area.&#160;Sharp gradients in the gravity difference near the caldera boundary point to a shallow source, consistent with the gravity signal arising from or near the ice-bedrock boundary.&#160;&#160;This indicates fast resurgence at B&#225;rdarbunga since 2015. The elevation of bed reflections delineated from radio echo sounding profiles (~2 MHz), measured within the caldera in June 2015 and accurately repeated in June 2019, further supports this.&#160; The suggested deformation mechanisms can be compared to geodetic observations outside the caldera for further evaluation. If all the signal is interpreted in terms of magma movements, a rise of the caldera floor by several meters and the inflow of 0.2-0.3 km<sup>3</sup>&#160;of new magma is inferred.</p>