Introduction

The German Bight in the south-eastern North Sea, including the Wadden Sea that is listed as a UNESCO World Heritage Site, is a coastal region of enormous ecological and economical importance, extensively used for transportation, fishing, tourism and, more recently, energy generation by offshore wind parks. Environmental conditions depend strongly on the exchange between land and sea, in particular on the balance between particle deposition and erosion, and dissolved land–sea fluxes of freshwater, nutrients and contaminants. Such fluxes can be difficult to constrain, however, as the tracking of diffuse sources is often complicated (Gallardo and Marui 2006; Burnett et al. 2006). Naturally occurring radiotracers could conveniently be used for this purpose, as their loss due to radioactive decay with increasing distance from the source can be used to constrain the rates of supply, as well as the transport away from the source by mixing and advection.

Hydrological restoration represents a major current challenge along the south-eastern North Sea coast (Reise 2005; van der Hagen et al. 2008). Moreover, weather-dependent fluctuations in the residence time of Wadden Sea tidal water masses have recently been hypothesised as pivotal in controlling regional eutrophication (van Beusekom et al. 2008), and is has been shown that pore waters from tidal flats represent an important part of the nutrient and carbon cycle in the North Sea (Billerbeck et al. 2006a, b). Therefore, additional tracers that could improve evaluations of groundwater and pore water discharge, water mass residence times or exchange rates would be a useful tool in this environment.

The suite of naturally occurring radon (222Rn) and radium isotopes (223Ra, 224Ra, 226Ra and 228Ra) have been widely used to quantify estuarine and ocean mixing processes, as well as submarine groundwater discharge to the coastal zone (e.g. Rama and Moore 1996; Moore 2000a, b; Charette et al. 2003; Burnett et al. 2006). All four natural Ra isotopes are produced by the decay of different isotopes of thorium (Th), a highly particle-reactive element. The decay of Th sorbed to particles represents a continuous source of the more soluble Ra to seawater. In turn, 222Rn is produced by the decay of 226Ra. It has been shown to be a powerful tracer of submarine groundwater discharge, especially if losses to the atmosphere can be quantified (Cable et al. 1996; Burnett and Dulaiova 2003; Stieglitz et al. 2010).

The distribution, enrichment and depletion of Ra and Rn isotopes in the water column depend on their production and radioactive decay, as well as the ambient hydrodynamic and physicochemical conditions (e.g. Key et al. 1985; Burnett et al. 2006; Ku and Luo 2008). In freshwater, Ra tends to sorb to particles, whereas in seawater it is found primarily in the dissolved form (Li et al. 1977). Therefore, desorption occurs during estuarine mixing, due mainly to the effect of increasing ionic strength (I) at higher salinity—that is, the partition coefficient of Ra varies inversely with I (Li et al. 1977; Moore et al. 1995; Webster et al. 1995; Krest et al. 1999; Charette et al. 2008). A similar effect takes place within freshwater aquifers, where radium adsorbed onto the surface of particles (Porcelli 2008) becomes mobilised when seawater infiltrates and replaces freshwater in the pore spaces. This process takes place frequently in aquifers within the tidal zone, as found in our study area.

Due to their short half-lives, both 224Ra (t1/2 = 3.7 days) and 223Ra (t1/2 = 11.4 days) can serve as tracers of hydrological processes over relatively short timescales (days, weeks), characterising exchange fluxes along the river–estuary–coastal zone continuum, or the movement of water masses over continental shelves. On such short timescales, radioactive decay does not significantly affect the activities of 226Ra (t1/2 = 1,600 years) and 228Ra (t1/2 = 5.8 years). Therefore, the combined activities of all four radium isotopes, the so-called radium quartet, can be used to study the exchange of water masses between different compartments of the coastal ocean on various temporal and spatial scales (Rama and Moore 1996).

The German Bight unifies several key processes that would be expected to be mirrored in the distribution of Ra and Rn isotopes, including sediment–water exchange in estuaries and on extensive tidal flats, some influenced by freshwater seepage, as well as mixing between nearshore and offshore water masses. Despite a wealth of information dealing with this subject matter for this region, to date there exist only very localised data based on the Ra approach (Schnetger et al. 2001; Shaw 2003) and, to our knowledge, no Ra data from the North Sea are available in the peer-reviewed literature.

In order to close this gap and to explore the potential use of Ra and Ra isotopes as tracers of land–ocean exchange in the North Sea, we studied the concentrations of all four naturally occurring radium isotopes in potential nearshore mixing end-members, and in the open German Bight. Our study aims at finding possibly distinctive isotopic signatures of the potential end-members and, based on the results, we explore the possibilities to disentangle contributions of individual hydrological compartments of the Wadden Sea to the German Bight.

Materials and methods

Study areas

Water samples were collected in three main study areas in the German Bight (Fig. 1a, b), these being a corridor extending from the barrier islands of Langeoog and Spiekeroog (East Frisian Wadden Sea) offshore to the island of Helgoland (Fig. 1c), an open tidal flat site influenced by freshwater seepage near the settlement of Sahlenburg (Fig. 1d), and the outer Weser Estuary (Fig. 1e).

Fig. 1
figure 1

Location maps of the North Sea (a) and the German Bight (b), as well as sampling locations (squares) in the open North Sea (c), the Sahlenburg intertidal flats (d) and the Weser Estuary (e). The bathymetry is shown in metres. NS North Sea transect sites, tf tidal flat site, B Blexen, Bhv Bremerhaven, E Elbe Estuary, W Weser Estuary; JB Jade Bay, L Langeoog, S Spiekeroog, H Helgoland, Sb Sahlenburg

The Wadden Sea covers approximately 9,300 km2 of the Southern North Sea (Grunwald et al. 2009). The tides are semidiurnal, and the region can be divided into two mesotidal sections (tidal range <3.5 m), characterised by barrier islands and back-barrier tidal flats, and a macrotidal section (tidal range >3.5 m) consisting of open, unprotected tidal flats lining the inner German Bight. The back-barrier tidal flats, including estuaries, have been estimated to cover 3,364 km2 (Grunwald et al. 2009). Century-long diking activities and land reclamation have resulted in an overall paucity of true mudflats in the region (sediments with mud contents exceeding 95%; mud = <63 μm fraction), and a dominance of sand flats (mud contents of <5%) and so-called mixed flats (for various sediment classification schemes, see Flemming 2000). The effects of this Wadden Sea squeeze (Delafontaine et al. 2000; Mai and Bartholomä 2000) are less pronounced in protected embayments (e.g. Jade Bay) that can accommodate fine-grained sediment deposition, and in tidal flats nourished by estuarine silts and clays (e.g. Sahlenburg tidal flats). Sediment reworking and destabilisation are caused by burrowing organisms and deposit feeders (bioturbation), as well as physical disturbance by tidal currents and storm events (e.g. Widdows et al. 2000; Volkenborn et al. 2009).

The main freshwater input to the German Bight originates from the Elbe and Weser rivers and a number of smaller streams, draining through sluice gates in the dikes. Very localised, fresh groundwater discharge occurs, for example, near Sahlenburg (cf. Gabriel et al. 2003) or near the island of Sylt (cf. Zipperle and Reise 2005). At the site near Sahlenburg (Fig. 1b), a mainland aquifer extends ca. 0.5 km onto the tidal flats, discharging fresh groundwater at scattered springs (“sand boils”, according to Li et al. 1996). Those at the Sahlenburg site are 10−2 to 3*10−1 m in diameter. During flooding tide, when the pressure gradient increases due to the overlying body of water, the punctual groundwater discharge at the sample site decreases to a minimum value of 21 mL minute−1 (30 L day−1), but the hydraulic gradient never allows the seawater to penetrate the freshwater aquifer. At low tide, the discharge increases to about 35 mL minute−1 (50 L day−1; Kurtz 2004).

To the south of the Sahlenburg tidal flats, the Weser River estuary is associated with a mean discharge of 323 m3 s−1 (Schuchardt et al. 1993), representing the second-largest freshwater source for the German Bight, after the Elbe River. The outer estuary is characterised by a funnel-like morphology with extensive open tidal flats. The mean tidal range exceeds 3.5 m, but water levels are also influenced by strong winds. Due to high current velocities in the estuary (Schuchardt et al. 1993), the water is generally well mixed and highly turbid, with a weakly stratified salt wedge. Highly mobile fluid mud layers >1 m in thickness have been observed, particularly at slack water (Schrottke et al. 2006). Average salinities near the city of Bremerhaven (Fig. 1c, own data) were found to range from 3–15 psu, decreasing to 1 psu at the tidal freshwater reach further upstream.

Sampling

Open North Sea transect

During cruise 234 of the RV Heincke (18–22 July 2005), water samples were taken at six locations (NS0 to NS5) along an open North Sea transect of ~60 km (Fig. 1c). Locations NS0 and NS1 were east and west respectively of the Otzumer Balje tidal inlet, on the shoreface of the islands of Spiekeroog and Langeoog. Locations NS2 to NS5 extended N–NE further offshore up to Helgoland. The water samples were analyzed for 223Ra, 224Ra, 226Ra and 228Ra isotopes (albeit long-lived isotopes not on all samples). It should be noted that the wind situation changed drastically during the sampling period, with calm conditions during the sampling of NS0 and NS1, and a storm event during the following transect to Helgoland.

Tidal flat

On 18 June 2004 and again 18 June 2006, surface seawater samples were collected from the Sahlenburg tidal flat near the permanent coastline (rectangle d in Fig. 1b; label tf in Fig. 1d) immediately after low and high tide (samples tf in Table 1, taken at flooding and ebbing tide). Groundwater samples were taken from nearby scattered freshwater springs at low tide on 18 June 2004. Here, all the samples were analyzed for 223Ra, 224Ra, 226Ra and 228Ra isotopes, as well as 222Rn.

Table 1 Ra isotopes in groundwater (Sahlenburg), and waters of the southern North Sea (NS North Sea transect), a German Wadden Sea tidal flat at Sahlenburg and the Weser Estuary (B Blexen, Bhv Bremerhaven) at different tidal stages. Flooding and ebbing refer to the inundation stage directly after low and high tide respectively. n.d. Not determined, n sample size

Weser Estuary

In the Weser Estuary, water samples were collected near the city of Blexen from aboard the RV Heincke on 21 July 2005 (Fig. 1e). An earlier batch of samples had been taken on 3 July 2002 near the city of Bremerhaven, directly from the shore (Handelshafen). In both cases, sampling was every 1–2 h over one tidal cycle, i.e. including high, low, flooding and ebbing tide phases.

Activities of 223Ra and 224Ra were determined in all samples from the two locations. 226Ra and 228Ra activities were determined only in the Blexen samples taken at high tide, as the high suspended particle loads from low tide samples hampered filtration, and high particle loads interfere with subsequent analysis by gamma spectroscopy.

226Ra and 228Ra measurements

From each sample, 40 L of water was filtered through a 1 μm pore size filter (Hytrex filter cartridges, 12 inch length). Radium was co-precipitated as Ba(Ra)SO4 and stored sealed for 3–4 weeks to enable the short-lived daughters 214Pb, 214Bi and 228Ac to equilibrate with their parent radionuclides 226Ra and 228Ra (Moore 1984). The activities of 226Ra and 228Ra were then determined by gamma-ray spectrometry using a high-purity germanium well-type detector. Standard aliquots of known activities analyzed under conditions identical to those of the water samples were used as reference material.

223Ra and 224Ra measurements

Radium was extracted usually from 40 to 90 L of sample water in each case, these volumes having been shown to be in the optimum range for minimising counting errors (Garcia-Solsona et al. 2008). Occasionally, lower volumes were processed because of high particle loads or limited availability (fresh groundwater). The samples were passed through a manganese-oxide-impregnated acrylic fibre at flow rates <0.5 L minute−1, which results in a good quantitative adsorption of Ra (Reid et al. 1979). Activities of 223Ra and 224Ra were then determined using a radium delayed coincidence counter (RaDeCC), a nuclear counter system suitable for low-activity measurements in natural waters (Moore and Arnold 1996). Analyses of the Bremerhaven (2002) estuarine samples were performed in air as carrier gas. The decay of these samples was monitored by repeated (every 2–3 days) measurements over 2 weeks to verify the concentrations of individual radionuclides via their half-life, and to exclude artefacts of counts from the 224Ra-decay chain that may occur in the 219Rn channel (Garcia-Solsona et al. 2008). All other samples were analyzed in He as carrier gas. Due to generally low absolute 224Ra activities in these samples, no significant spillover of counts from the 224Ra channel to the 223Ra channel was expected, and no second counting for 223Ra was performed in this case.

The amounts of short-lived 223Ra and 224Ra radionuclides were derived from the emanation of their daughter radionuclides 220Rn and 219Rn respectively. In order to limit device-specific errors, each measurement was repeated two or three times using multiple RaDeCC channels. Blank determination was performed using a manganese-oxide-impregnated fibre that had been flushed with de-ionised water. A final measurement of the samples was conducted after several weeks in order to determine the supported 223Ra and 224Ra. The parent of 223Ra, 227Ac, has been shown to be somewhat elevated in coastal areas, potentially causing a small but regionally variable amount of supported 223Ra (Geibert et al. 2008).

222Rn measurements

In each case, a 10 mL water sample was collected in a syringe pre-filled with 10 mL of a non-evaporating, water-immiscible scintillation cocktail (di-isopropylnaphthalene, Hidex Maxilight) to prevent outgassing from the aqueous phase. 222Rn activities were determined in the home laboratory, using a portable single tube liquid scintillation counter (Triathler® LSC). As 222Rn has a much higher affinity for the organic solvent, it was extracted from the aqueous phase by short, vigorous shaking of the sample. The mixture was then allowed to stand for 3 h in order to separate the phases and to equilibrate radon with its α-emitting daughters 218Po and 214Po (Purkl and Eisenhauer 2004). The excess Rn concentration, 222Rnex, was determined in a second extraction step after 4 weeks, when 222Rn had equilibrated with its parent radionuclide 226Ra (Purkl and Eisenhauer 2004). Blank levels for this method are extremely low (a few counts per day), making the LSC particularly suitable for Rn measurements on very small volumes.

Salinity

Salinity data for the open North Sea transect were obtained from the physical oceanography group onboard the RV Heincke. Other salinity data are own conductivity-based measurements.

Error handling

The uncertainties for the radium isotope analyzes were determined by applying the principle of error propagation, including decay, uncertainty in the efficiency determination, and chance coincidence corrections (for the RaDeCC), following the strategy presented in Giffin et al. (1963) and Moore and Arnold (1996). The error of the blank correction, considered to be small compared to the effects of efficiency determination, was not propagated for the short-lived Ra isotopes. The detector efficiency for the long-lived Ra isotope measurements was 12 and 10% for 226Ra and 228Ra respectively. The efficiency of the RaDeCC varied between 27% (224Ra) and 21% (223Ra) in helium, and it was only ~25% (224Ra) and ~17% (223Ra) in air, consistent with the findings of Moore (2008). The moisture content of the fibres and, consequently, efficiency were found to be more variable when air was used as carrier gas. The uncertainties for 222Rn were determined by simple counting statistics.

Results

228Ra and 226Ra

The offshore activities of the long-lived Ra isotopes were in general lower than in nearshore waters (see Fig. 2, Table 1). Activities of 228Ra (in dpm 100 kg−1) were 35.1 ± 3.5 during low tide at nearshore site NS0 (presumably a signal from the back-barrier tidal flat), decreasing to 28.2 ± 2.8 at high tide, and only 23.0 ± 2.3 at the offshore site NS4 near Helgoland. 226Ra activities were 13.8 ± 1.4 dpm 100 kg−1 during low tide at site NS0, 11.8 ± 1.2 during high tide at site NS0, and 8.2 ± 0.8 at offshore site NS4.

Fig. 2
figure 2

Maximum activities of natural Ra isotopes (dpm 100 kg−1) and corresponding salinity in the sampling area (cf. Table 1 for sampling dates). In all cases, maximum activities were recorded during ebbing or at low tide. a Salinity (psu), b 223Ra, c 224Ra, d 226Ra, e 228Ra. For orientation, see geographic features in Fig. 1b

Irrespective of the tidal phase, the values that were recorded on the tidal flat (tf) were similar to the signal originating from the nearshore site NS0 at low tide along the transect (cf. differences between ebbing and flooding tide were within the measurement uncertainties for both isotopes at tf). 228Ra was found to vary from 39.3 ± 3.8 to 46 ± 4.5, and 226Ra activity from 15.5 ± 1.5 to 16.5 ± 1.7 dpm 100 kg−1. Furthermore, the long-lived Ra activities in the fresh groundwater discharging on the tidal flat were quite similar, with 228Ra and 226Ra values of 42.1 ± 4.1 and 21.3 ± 2.2 dpm 100 kg−1 respectively. The 228Ra and 226Ra activities in the Weser Estuary at high tide were 46.3 ± 4.6 and 17.1 ± 1.7 dpm 100 kg−1 respectively (Table 1), again very similar to the other nearshore values, but significantly different from the offshore values.

224Ra and 223Ra

North Sea transect

The highest 224Ra and 223Ra activities along the transect were measured in the outflow of the back-barrier system (site NS0, Table 1) at low tide, with 20.3 ± 6.1 for 224Ra and 1.2 ± 0.1 dpm 100 kg−1 for 223Ra. Towards the open North Sea, activities decreased to 2.0 ± 0.6 for 224Ra and 0.10 ± 0.02 dpm 100 kg−1 for 223Ra (Fig. 2b, c), slightly rising again near Helgoland (site NS5). The salinity increased slightly from 32 psu at site NS1 nearshore to 33 psu at site NS4 further offshore.

Tidal flat and fresh groundwater

The highest activities of the short-lived Ra isotopes were measured on water samples collected from the tidal flat near Sahlenburg (up to 85.3 ± 6.3 for 224Ra and 8.0 ± 1.2 dpm 100 kg−1 for 223Ra; see Table 1). The 224Ra activities from the tidal flat, even during flooding tide when only 34.3 ± 2.2 dpm 100 kg−1 was recorded, were much higher than in the fresh groundwater discharging from nearby sand boils (224Ra: 5.1 ± 0.9 dpm 100 kg−1). The difference was less pronounced for the 223Ra activities (3.6 ± 0.5 to 8.0 ± 1.2 on the tidal flat vs. 2.6 ± 1.3 in the fresh groundwater), associated with a very low 224Ra/223Ra ratio (Table 1). Distinctly higher 224Ra and 223Ra activities were detected at ebbing tide, with values rising to 85.3 ± 6.3 and 8.0 ± 1.2 dpm 100 kg−1 respectively. The salinity of seawater samples varied from 28 psu (flooding tide) to 26 psu (ebbing tide), whereas the salinity of the groundwater was found to be 0.7 psu, indicating its origin from the freshwater aquifer that extends onto the intertidal flat.

Weser Estuary

At the Bremerhaven (Bhv) and Blexen (B) sites (Fig. 1e), temporal variations of activities were found to be related to the tidal phase (Table 1, Fig. 3). The 224Ra activities were 26.1 ± 8.2 to 36.5 ± 6.1 for B and 31.6 ± 1.1 to 36.1 ± 1.0 dpm 100 kg−1 for Bhv. The 223Ra activities ranged from 1.8 ± 0.1 to 2.5 ± 0.2 for B and 3.0 ± 0.5 to 4.0 ± 0.4 dpm 100 kg−1 for Bhv. A peak in 224Ra activity was observed with decreasing salinity during ebbing tide at the Bhv site (Fig. 3d). A similar trend would be expected for 223Ra, but errors for low 223Ra activity are too large to allow a clear statement. The salinity was in the range of 15–18 psu at high tide and 6–12 psu at low tide.

Fig. 3
figure 3

Temporal variations of short-lived Ra isotopes (black) and salinity (white) in the Weser Estuary at Blexen (a 223Ra, b 224Ra; 21 July 2005) and Bremerhaven (c 223Ra, d 224Ra; 3 July 2002). HT High tide, LT low tide

222Rn

Preliminary 222Rn measurements carried out on the intertidal flat revealed that the activity of 222Rn in seawater was significantly lower than that in fresh groundwater (234.4 ± 40.4 dpm kg−1). This was the case at both flooding (46.9 ± 8.8 dpm kg−1) and ebbing tide (66.7 ± 11.2 dpm kg−1).

Discussion

226Ra and 228Ra

Extensive continental shelf areas have been reported to be an important source of 228Ra to the ocean (Moore 1981), prominent examples being the Arctic and the Yellow Sea (Rutgers van der Loeff et al. 1995, 2003; Kim et al. 2005). The North Sea as a whole resembles these larger examples, while the Wadden Sea, with its extensive intertidal flats, constitutes a distinct situation within the North Sea, with exceptionally strong water and sediment dynamics. To our knowledge, we report here the first values for long-lived radium isotopes from the North Sea, the Wadden Sea and the Weser Estuary. As these are coastal environments influenced by a large shelf area, we would expect to find comparatively high values determined largely by the residence time of Atlantic waters in this shallower North Sea sector. Indeed, open ocean 228Ra concentrations are typically much lower than coastal concentrations as a consequence of decay and mixing (Sarmiento et al. 1982).

Our activities for 228Ra in the open North Sea (site NS4, 23.0 ± 2.3 dpm 100 kg−1) and in the Weser Estuary (23.0 ± 2.3 dpm 100 kg−1) are comparable to observations for the Siberian shelves incorporating the Ob, Yenisey and Lena estuaries, with values ranging from 3.3 dpm 100 kg−1 towards the open Arctic Ocean up to 40.9 dpm 100 kg−1 in the estuary of the river Ob (Rutgers van der Loeff et al. 2003). Activities reported from another large shelf region, the Yellow Sea, range from 3.4 ± 0.8 to 68.9 ± 1.1 dpm 100 kg−1, with values often at the higher end of the range (Kim et al. 2005). From a slightly smaller shelf region in the northeast Gulf of Mexico, values ranging from 13.6 to 137.3 dpm 100 kg−1 have been reported by Moore (2003). This suggests that the North Sea is no exception to the general trend of high 228Ra activities on extended continental shelves due to release from sediments. Indeed, similar activities have been reported from the South Atlantic Bight (Moore 2000b), the Amazon shelf (Moore et al. 1995), and the Mississippi and Atchafalaya rivers (Krest et al. 1999).

In contrast to 228Ra, 226Ra with a half-life of 1,600 years is typically less variable in surface waters, but it has previously been found to show some variation with water mass (Ku and Luo 1994, 2008). Typical values of about 7 dpm 100 kg−1 have been reported from the North Atlantic (Broecker et al. 1976), with relatively little variation as a result of the long half-life of the isotope. Our data for the German Bight of the North Sea (site NS4: 8.2 ± 0.8 dpm 100 kg−1) reveal a similar concentration for 226Ra, confirming the North Atlantic origin of North Sea waters, with little subsequent modification of the signal.

Summarising our findings for long-lived Ra isotopes at a larger scale, we have demonstrated that 226Ra concentrations in the open North Sea are very similar to those in the North Atlantic, the former being only marginally more elevated. 228Ra values indicate enrichment over the large shelf area of the North Sea, consistent with findings from other shelf regions.

Somewhat unrelated to the large-scale trends are the long-lived Ra isotope activities in our fresh groundwater, which originates from an aquifer related to a buried valley from a past glacial period (Gabriel et al. 2003). With 42.1 ± 4.1 and 21.3 ± 2.2 dpm 100 kg−1 for 228Ra and 226Ra respectively, they deviate considerably from values reported for other coastal aquifers. For example, Rama and Moore (1996) found 228Ra activities as high as 176 dpm 100 kg−1 in a seep of the North Inlet salt marsh near Georgetown, South Carolina. More recently, Breier and Edmonds (2007) measured high 226Ra (95–446 dpm 100 kg−1) and 228Ra (74–576 dpm 100 kg−1) activities in the freshwater end-member of a spring discharging in the Corpus Christi Bay groundwater system of Texas. We explain the low radium activity in our groundwater by the very low salinity (0.7 psu), leading to a stronger adsorption onto particles. The lack of pronounced enrichment of long-lived Ra isotopes in fresh groundwater discharging in the intertidal zone limits the use of 228Ra and 226Ra isotope signals as tracers of fresh localised groundwater inputs to the German Bight. However, recirculated saline pore water or groundwater would contribute to radium enrichment in the Wadden Sea and, eventually, the open North Sea.

223Ra and 224Ra

The decrease in short-lived Ra isotopes observed along the North Sea transect is consistent with the release from a main source of short-lived radium in the Wadden Sea and its estuaries, and subsequent decay towards the open North Sea. Although the decrease of 224Ra and 223Ra with distance from the coastline resembles the results of studies done elsewhere (e.g. the Gulf of Mexico, Moore and Krest 2004; the Amazon shelf, Moore et al. 1995), the absolute magnitude of Ra activities in the open North Sea is distinctly smaller. At a distance of 60 km offshore, only a negligible activity of short-lived Ra was measured, whereas in the Amazon shelf and the Gulf of Mexico, the decrease of the coastal Ra signal could be followed over tens to hundreds of kilometres. As a consequence of the fast decay of short-lived Ra isotopes, and the strong dilution of Ra derived from adjacent nearshore sources (i.e. the Weser Estuary, tidal flats, and submarine groundwater), short-lived Ra isotopes may be in the range of the detection limit beyond Helgoland. In principle, the decrease of short-lived Ra with distance from shore would allow constraining horizontal mixing according to previous studies (Moore 2000b). However, the potential importance of Ra diffusing from sediments in the offshore area would need to be assessed in order to understand the overall Ra budget in the German Bight of the North Sea, and the limited amount of data on the offshore transect of the present study is not sufficient to derive a robust mixing estimate.

On the tidal flat, distinct variations in 224Ra and 223Ra activity were measured as a function of the tidal stage. The peak activities might be related to the turbulent flushing of the sediment bed due to tidal disturbance, and possibly due to bioirrigation. As the vast tidal flats are exposed to air during low tide, which temporarily prevents a loss of newly produced short-lived Ra, the inundation during flooding tide causes an extensive flushing of the upper sediment layers. We hypothesise that this flushing is assisted by bioirrigation due to burrowing benthic organisms, causing a stripping of radionuclides from particles after ebb tide. Also, pore water from tidal flats lost to seawater during ebbing tide (Billerbeck et al. 2006a, b) is likely to represent an important source of Ra to the Wadden Sea.

In the Weser Estuary, the increasing Ra activities observed during ebbing tide most probably result from the radionuclide fraction that is adsorbed to particles under low salinity conditions, and stripped off when brought in contact with the tidal river water of higher ionic strength. Li et al. (1977) showed that Ra is released in estuaries during the mixing between freshwater and seawater. High loads of suspended particles in the Weser Estuary, especially at low tide (Schrottke et al. 2006), may cause a release of Ra within the water column, as a consequence of both increasing salinity and decreasing particle concentrations, whereas hydraulic action of the tides can deliver short-lived Ra from pore water and groundwater further upstream, in combination leading to a delayed and more complex response of Ra isotopes. Our results can be seen only as a first indication of possible processes, and more detailed sampling campaigns would be required to establish a Ra budget for the Weser Estuary.

222Rn

During both flooding and ebbing tide, the 222Rn signal was clearly lower in seawater than in fresh groundwater. Turbulent mixing during inundation would support the mixing of 222Rn in the overlying water and, consequently, outgassing at the water–atmosphere interface. We expect that the overlying water body of the intertidal flat would be too shallow to develop a vertical 222Rn activity gradient, especially as the intertidal flat is partially exposed to air during low tide. Our small dataset provides at least a preliminary indication of activity levels that should prove useful in the future design of more detailed studies. A particularly promising application is the linking of the noble gas 222Rn to investigations of methane (Grunwald et al. 2009), as this could help to constrain regional budgets of gas release.

Conclusions

This paper presents a first dataset on four radium radioisotopes (the radium quartet) in the German Bight of the North Sea. On a regional scale, the isotope distributions were found to track primarily the mixing between nearshore and offshore water masses. It was shown that localised fresh groundwater discharge is most probably not an important source of the short-lived 223Ra and 224Ra isotopes in this environment. Rather, 223Ra and 224Ra are released from the vast tidal flats during inundation, with a potential contribution of diffuse fluxes from saline or brackish groundwater. We found very similar activities of both short- and long-lived Ra isotopes in all potential water mass end-members, except for the waters of the open German Bight. The Weser Estuary, the tidal flats near Sahlenburg, and discharge originating from the back-barrier tidal flats near Spiekeroog all have virtually identical activities of 228Ra and 226Ra. These activities decrease towards the open German Bight, revealing a potential for constraining offshore fluxes of material originating in the tidal zone on a scale of tens of km, but limiting the potential discrimination of signals from sources like the Weser Estuary, as well as back-barrier and unprotected tidal flats. For the establishment of a Ra budget for the region, the signatures from nearby the Elbe River and from possible Ra fluxes from the sea floor in the open North Sea would be required, especially as the Elbe has been shown to be the main contributor of fresh water to our study area in historical times (Scheurle and Hebbeln 2003).

As an additional finding, we demonstrate for the first time that 228Ra in the open German Bight is elevated, comparable to the situation in other regions with extended submerged shelf areas.