The Vimmerby Moraine is the only significant ice-marginal moraine on the eastern side of southern Sweden, but no detailed studies exist on its formation during the final deglaciation of the Fennoscandian Ice Sheet. Through ground-penetrating radar surveys and detailed sediment logging, we provide evidence for an active, oscillating ice margin during the formation of the Vimmerby Moraine, suggesting that the deglaciation of the South Swedish Uplands was, at least in some regions, dynamic. Ground-penetrating radar surveys enabled imaging of internal sediment and delineation of the bedrock surface. These were complemented by common mid-point surveys and sediment logging, as well as lithofacies analysis at three exposures in agravel pit. This approach revealed multiple subglacial till units partially separated by intercalated glacifluvial deposits. The glacifluvial sediments exhibit evidence of glaciotectonism, suggesting active overriding by the last ice sheet. Further evidence of an active ice margin is provided by the ground-penetrating radar profiles collected perpendicular to the moraine crest. These contain a series of northerly dipping reflectors, which we interpret as evidence of repeated basal freeze-on and melt-out of sediment slabs during ice margin oscillations, as has been observed at contemporary glacier margins in Iceland, Norway, and the Alps. The data presented here demonstrate that the last Fennoscandian Ice Sheet remained active around the time of the Vimmerby Moraine formation. This work highlights the benefits of including detailed sediment logging and near-surface geophysical surveys in the interpretation of deglaciation dynamics.

The North Sea Basin has been covered by ice sheets originating from both the British Isles and Scandinavia at multiple times during the Pleistocene. The Witch Ground Basin (WGB) in the central northern North Sea is a critical location in terms of interpreting Late Pleistocene glacial to glacimarine history of the North Sea since it was the location of the Witch Ground Ice Stream that was active on multiple occasions during the Mid to Late Pleistocene. We map five mega-scale glacial lineation flowsets corresponding to the changing ice flow direction of the Witch Ground Ice Stream and investigate the sedimentological fingerprint and corresponding subglacial depositional processes of this palaeo-ice stream. We show that sorted sand layers within a subglacial traction till represent periodic hydraulic jacking and ice–bed decoupling at the base of the Witch Ground Ice Stream. In contrast to previous studies that have described glacitectonites deposited below the most recent grounded ice in the WGB, we present analysis of sediment cores that recovered primarily massive diamictons without any obvious deformation structures. The most recent ice cover in the WGB (~18–16 ka) was thought to have been sourced from a localized ice cap over Orkney and Shetland. The presence of chalk clasts sourced from NW of the WGB described in this study from the stratigraphically youngest till confirms this interpretation. The transition from subglacial to glacimarine deposition, while acoustically well defined (from opaque to laminated acoustic units), appears surprisingly uniform in the recovered sediment cores, but can be differentiated based on a change in colour including mottling and banding, presence of whole intact shells, and the increased number of silt and sand lenses. 14^C dating of glacimarine muds indicate high sedimentation rates of between 80 and 260 cm ka⁻¹. The transition from glacimarine to marine deposition is represented by a comparative decrease in sedimentation rate and deposition of Holocene age sandy mud. This study demonstrates a highly dynamic Witch Ground Ice Stream in the northern North Sea during the Late Pleistocene with evolving subglacial hydrology and depositional processes at the ice stream bed that left a distinct geomorphological and sedimentological fingerprint within the WGB.

Lake Victoria is the world’s largest tropical lake and the third-largest water body, providingsignificant water resources for surrounding environments including the cultural, societal, andlivelihood needs of people in its basin and along the White Nile. The aim of this study was to usedecade-long time series of measured lake flow in the lake system and phosphorus deposition todevelop a suitable numerical model based on shallow water equations (SWE) for assessing waterquality in Lake Victoria, an increasingly important tool under climate variation. Different tech-niques were combined to identify a numerical model that included: i) a high-resolution SWEmodel to establish raindrop diffusion to trace pollutants; ii) a two-dimensional (2D) verticallyintegrated SWE model to establish lake surface flow and vertically transported wind speed flowacting on lake surface water by wind stress; and iii) a site-specific phosphorus deposition sub-model to calculate atmospheric deposition in the lake. A smooth (non-oscillatory) solution wasobtained by applying a high-resolution scheme for a raindrop diffusion model. Analysis with thevertically integrated SWE model generated depth averages for flow velocity and associatedchanges in water level profile in the lake system and showed unidirectional whole lake windblowing from the southwest to northeast. The atmospheric phosphorous deposition modelenabled water value assessment for mass balances with different magnitudes of both inflows andoutflows demonstrating annual total phosphorus at 13, 500 tons concentrating at mid-lakewestern and eastern parts. The model developed here is simple and suitable for use in assess-ing flow changes and lake level changes and can serve as a tool in studies of lake bathymetry andnutrient and pollution transport processes. Our study opens towards refining models of complexshallow-water systems

Severe extratropical winter storms are a recurrent feature of the European climate and cause widespread socioeconomic losses. Due to insufficient long-term data, it remains unclear whether storminess has shown a notable response to changes in external forcing over the past millennia, which impacts our ability to project future storminess in a changing climate. Reconstructing past storm variability is essential to improving our understanding of storms on these longer, missing timescales. Peat sequences from coastal ombrotrophic bogs are increasingly used for this purpose, where greater quantities of coarser grained beach sand are deposited by strong winds during storm events. Moving inland however, storm intensity decreases, as does sand availability, muting potential paleostorm signals in bogs. We circumvent these issues by taking the innovative approach of using mid-infrared (MIR) spectral data, supported by elemental information, from the inorganic fraction of Store Mosse Dune South (SMDS), a 5000-year-old sequence from a large peatland located in southern Sweden. We infer past changes in mineral composition and thereby, the grain size of the deposited material. The record is dominated by quartz, whose coarse nature was confirmed through analyses of potential local source sediments. This was supported by further mineralogical and elemental proxies of atmospheric input. Comparison of SMDS with within-bog and regionally relevant records showed that there is a difference in proxy and site response to what should be similar timing in shifts in storminess over the-100 km transect considered. We suggest the construction of regional storm stacks, built here by applying changepoint modelling to four transect sites jointly. This modelling approach has the effect of reinforcing signals in common while reducing the influence of random noise. The resulting Southern Sweden-Storm Stack dates stormier periods to 4495-4290, 3880-3790, 2885-2855, 2300-2005, 1175-1065 and 715-425 cal yr BP. By comparing with a newly constructed Western Scotland-Storm Stack and proximal dune records, we argue that regional storm stacks allow us to better compare past storminess over wider areas, gauge storm track movements and by extension, increase our understanding of the drivers of storminess on centennial to millennial timescales.

Seismic refraction models should routinely be reported with their associated uncertainty. Tomographic solutions are widespread, but estimating uncertainties in these via Monte Carlo simulation places great demands on computer resource, hence this task is often omitted. By considering the Plus-Minus method of seismic refraction interpretation, we use Monte Carlo simulations to evaluate the uncertainty in seismic refraction results and determine the sources of uncertainty that are most impactful on the reliability of the output model. Our analysis considers the impact of survey mislocation (i.e., geophones misplaced from a planned position) and interpretational problems (i.e., misidentification of first-break picks and uncertainty in identifying crossover distances) on the overall uncertainty in inferred unit thicknesses and seismic velocities. These are considered for synthetic data with varying subsurface velocity structure, and for field data collected at a shallow (< 50 m) bedrock site in north Wales (UK). Analysis of synthetic data shows that the impact of the aforementioned errors on thickness estimates is similar to 1000 times that on velocity estimates. Of all permutations tested, the most significant impact on thickness uncertainty was the accuracy of first-break picks, with the variance in target thickness estimates increasing roughly exponentially with first-break pick uncertainty. It is therefore prudent to minimise such uncertainty through appropriate survey practice (e.g., maximising source energy, taking multiple shots for stacking) and to properly define the resultant uncertainty in unit thickness and velocity estimates. The simplicity of the Plus-Minus method makes it an effective tool for highlighting the errors that would impact more sophisticated interpretation approaches, such as tomography or Full Waveform Inversion. The results from such analysis can be directly applied in straightforward environmental or engineering investigations and can be used to inform more advanced refraction methods. As such, the practice we highlight should be considered for any refraction interpretation.

Geophysical surveys provide an efficient and non-invasive means of studying subsurface conditions in numerous sedimentary settings. In this study, we explore the application of three geophysical methods to a proglacial environment, namely ground penetrating radar (GPR), seismic refraction and multi-channel analysis of surface waves (MASW). We apply these geophysical methods to three glacial landforms with contrasting morphologies and sedimentary characteristics, and we use the various responses to assess the applicability and limitations of each method for these proglacial targets. Our analysis shows that GPR and seismic (refraction and MASW) techniques can provide spatially extensive information on the internal architecture and composition of moraines, but careful survey designs are required to optimise data quality in these geologically complex environments. Based on our findings, we define a number of recommendations and a potential workflow to guide future geophysical investigations in analogous settings. We recommend the initial use of GPR in future studies of proglacial environments to inform (a) seismic survey design and (b) the selection of seismic interpretation techniques. We show the benefits of using multiple GPR antenna frequencies (e.g., 25 and 100 MHz) to provide decimetre scale imaging in the near surface (e.g., < 15 m) while also enabling signal penetration to targets at up to ∼40 m depth (e.g., bedrock). This strategy helps to circumvent changes in radar signal penetration resulting from variations in substrate conductivity or abundant scatterers. Our study also demonstrates the importance of combining multiple geophysical methods together with ground-truthing through sedimentological observations to reduce ambiguity in interpretations. Implementing our recommendations will improve geophysical survey practice in the field of glacial geology and allow geophysical methods to play an increasing role in the interpretation of glacial landforms and sediments.

Sediments deposited during glacial-interglacial cycles through the Early to Mid-Pleistocene in the North Sea are chronologically poorly constrained. To contribute to the chronology of these units, amino acid racemization (AAR) and strontium (Sr) isotope analyses have been performed on samples from four shallow borings and one oil well along a transect in the northern North Sea. D/L Asp (aspartic acid) values obtained through reverse-phase liquid chromatography in the benthic foraminiferal species Elphidium excavatum is focused on because of consistent results and a good stratigraphic distribution of this benthic species. For the Early Pleistocene, an age model for the well 16/1–8, from the central part of the northern North Sea based on Sr ages allows for dating of the prograding wedges filling the pre-Quaternary central basin. A regional calibration curve for the racemization of Asp in Elphidium excavatum is developed using published ages of radiocarbon-dated samples and samples associated with the previously identified Bruhnes/Matuyama (B/M) paleomagnetic boundary and a Sr age from this study. Based on all the available geochronological evidence, samples were assigned to marine oxygen isotope stages (MIS) with uncertainties on the order of 10–70 ka. Sr ages suggest a hiatus of <2 million years (Ma) possibly due to non-deposition or low sedimentation between the Utsira Formation (Pliocene) and the Early Pleistocene. An increase in sedimentation rates around 1.5 ± 0.07 Ma (∼MIS 51) may partly be due to sediment supply from rivers from the south-east and partly due to the extension of ice sheet around 1.36 ± 0.07 Ma from the Norwegian coast to the central North Sea. A possible basin-wide glaciation occurred around 1.1 Ma (∼ MIS 32) (upper regional unconformity/top of unit Q4 in this study), resulting in erosion and regional unconformity. Two interglacials in the Norwegian Channel have been dated: the Radøy Interglacial to 1.07 ± 0.01 Ma (possibly MIS 31, the ‘super interglacial’), and the Norwegian Trench Interglacial to 0.50 ± 0.02 Ma (possibly MIS 13). A massive till unit identified at the same stratigraphic level in all shallow borings may partly represent an extensive MIS 12 glaciation. This study shows that the combined use of amino acid racemization data and Sr isotope chronology can refine the chronological ambiguities of Quaternary North Sea sediments related partly to the impact of glacial processes.

In 1964, exploration drilling in the German Sector of the North Sea hit a gas pocket at ∼2900 m depth below the seafloor and triggered a blowout, which formed a 550 m-wide and up to 38 m deep seafloor crater now known as Figge Maar. Although seafloor craters formed by fluid flow are very common structures, little is known about their formation dynamics. Here, we present 2D reflection seismic, sediment echosounder, and multibeam echosounder data from three geoscientific surveys of the Figge Maar blowout crater, which are used to reconstruct its formation. Reflection seismic data support a scenario in which overpressured gas ascended first through the lower part of the borehole and then migrated along steeply inclined strata and faults towards the seafloor. The focused discharge of gas at the seafloor removed up to 4.8 Mt of sediments in the following weeks of vigorous venting. Eyewitness accounts document that the initial phase of crater formation was characterized by the eruptive expulsion of fluids and sediments cutting deep into the substrate. This was followed by a prolonged phase of sediment fluidization and redistribution widening the crater. After fluid discharge ceased, the Figge Maar acted as a sediment trap reducing the crater depth to ∼12 m relative to the surrounding seafloor in 2018, which corresponds to an average sedimentation rate of ∼22,000 m3/yr between 1995 and 2018. Hydroacoustic and geochemical data indicate that the Figge Maar nowadays emits primarily biogenic methane, predominantly during low tide. The formation of Figge Maar illustrates hazards related to the formation of secondary fluid pathways, which can bypass safety measures at the wellhead and are thus difficult to control.