Samples were taken at seven sites (Figure 1) that represent salt marshes inundated by the Guadalquivir estuary during the 16 field campaigns conducted between March 2016 and March 2018. Sites were selected to account for the salinity gradient present in the estuary, which can range from 34 psu in the river mouth down to near 0 psu in the upper limit of the land strip adjacent to the riverside considered in our study (Huertas et al., 2018). Vegetation in the marshes is characterized by pioneer halophilous species like Sarcocornia fruticosa, Halimione portulacoides (Chenopodiaceae), and Limoniastrum monopetalum (Plumbaginaceae) (García-Murillo et al., 2014). During the monitoring period, we also identified Arthrocnemum macrostachyum, Plantago coronopus, Cotula coronopifolia, Taraxacum spp, Centaurium pulchellum, Polypogon maritimus, and Juncus sp., although plant biodiversity exhibited a marked spatio-temporal variability, which is higher during spring periods (Romagnoli, 2018).

All samplings were carried out when the estuary was under the influence of spring and rising tides. Some studies in marshes and other coastal wetlands, such as mangroves, have shown that there is a strong tidal signal in CH₄ and water chemistry (Call et al., 2015). Our sampling design might therefore have affected gas dynamics, as sample collection occurred along with the net water flux into the system. However, those tidal conditions ensured overflow in the mudflat and the presence of a significant water level (a minimum water height of 0.1 m) allowing sampling, which was conducted on foot. This strategy may have also had implications for gas distribution due to sediment dynamics disturbance. Nonetheless, it was the only way to access the flooded marshes in certain sites. Considering this caveat, we tried to minimize sediment disturbance as much as possible.

The floodplain is highly variable depending on annual rainfall (Huertas et al., 2017), and is also presently affected by the estuarine intrusion, that in turn is controlled by the tide amplitude. In fact, water levels were too low for sampling during the summer months and early fall. Therefore, data are not available for all sites during the dry season, and temporal variability of aquatic parameters has been examined between periods in which measurements could be performed.

Even though surface water coverage in the marshland changes temporarily, the mudflat is a relatively shallow system. Considering past and current records of wind speed in the area and the average depth of the water column, the Richardson number was calculated, which indicated that at each site, the flooded area was a well-mixed water body (not shown). Due to mixing conditions, methane levels measured at each sampling station were assumed to be representative of each inundated zone.

In all sites, conductivity, temperature and pH (National Bureau of Standards, NBS scale) were obtained with a Yellow Spring (YSI Incorporate) portable multiparameter probe YS6820v2. The pH probe was calibrated before samplings using the United States National Bureau of Standards buffer solutions (4 and 7). Discrete surface water samples were taken to determine dissolved CH₄ concentrations and the rest of the aquatic biogeochemical variables considered in our study such as chlorophyll (Chl a), dissolved oxygen (DO), total alkalinity (TA), inorganic nutrients (NH4+, NO2−, NO3-, PO43-, and SiO44-), dissolved organic carbon (DOC), total dissolved nitrogen (TDN), and suspended particulate matter (TSM).

The air-water CH₄ flux (F_CH4, μmol m⁻² d⁻¹) was computed as:

Where Cw and Ca are the dissolved gas concentration and the equilibrium concentration in water, based on the molar atmospheric ratio, and k (cm h⁻¹) is the gas transfer rate as a function of wind speed at 10 m height, that was calculated from k and normalized to a Schmidt number (S_c) of 600 (k₆₀₀) according to:

S_c was obtained using the formulations given by Wanninkhof (1992) interpolated at the in-situ salinity from the freshwater and seawater equations and k₆₀₀ was calculated from U₁₀ according to Cole and Caraco (1998).

where U₁₀ was calculated from uz measured at the meteorological station according to Smith (1988).

We acknowledge that the empirical relationship chosen to adjust k₆₀₀ values may not be ideal for representing wind–enhancement effects in the flooded marshes (shallow water bodies with substantial variations in extent) but it was proven to be suitable to compute air-water CO₂ fluxes in relation to other parameterizations (Morris et al., 2013). Nevertheless, to give an indication of the uncertainty this choice introduces, air-water CH₄ fluxes were also predicted considering a 0.67 exponent for the Schmidt number. Despite minor differences obtained at high wind speeds (<10%), the mean flux remained unchanged when annual values were averaged.

The program language MATLAB was used to perform statistics analysis. Probability distributions of variables were examined through a Shapiro-Wilk test. Pearson's product-moment correlation (PPMC) was used to test for significant correlations between variables. Significance levels were set at p < 0.05.

All data contained in this work are available for download from Digital.CSIC, the Institutional Repository of the Spanish National Research Council (CSIC), http://hdl.handle.net/10261/173204, http://dx.doi.org/10.20350/digitalCSIC/8588).

Monthly averaged air temperature during the samplings ranged between a minimum of about 8 and a maximum of 24.5°C in January 2017 and July 2017, respectively (Figure 2A). Rainfall was quite variable during the study period (Figure 2B), typical of the Mediterranean climatic region that is characterized by a rainy season (usually extending from October to March) and a dry season (from June to September) when rainfall is almost absent. The inter-annual precipitation pattern in this geographical area also fluctuates often, which affects the flooding regime in the wetlands (Huertas et al., 2017). Particularly, between March 2016 and March 2018 which corresponds to 3 hydrological years (starting 1 September 2015), total annual rainfall was 493 mm from September 2015 to August 2016, 539 mm from September 2016 to August 2017 and 500 mm from September 2017 to March 2018. These values are very close to the mean annual precipitation in the region, equivalent to 550 mm (Serrano Martín et al., 2008), thereby representing normal hydrological cycles. Monthly averaged precipitations during the sampling period varied from a minimum of 0 mm, registered during the months of June and September 2017, to maximum of about 150 mm in spring of 2016 and 2018 and autumn 2016 (Figure 2B). The tidal coefficient at the mouth of the estuary during the sampling dates and defined as the amplitude of the tide forecast, oscillated between 41 and 111, as observed in January 2017 and November 2016, respectively (Figure 2B). Monthly averaged horizontal wind velocity at 10 m above the surface ranged between 1.8 in January 2017 and 4.5 m s⁻¹ in March 2018 with a median value of 2.5 m s⁻¹ (Figure 2C).

Water temperature during samplings oscillated between 11.4 and 32.3°C (Figure 2A), with a clear seasonal pattern that was similar at all sites (not shown for sites individually). The range of variation of other aquatic biogeochemical variables is summarized in Table 1.

Because of meteorology and tide variations, salinity markedly changed across the marshes during the study period. It is worth noting that the Guadalquivir estuary is a mesotidal system, in which tidal influence can be noticeable up to 100 km upstream of the river mouth. Nevertheless, high freshwater discharges from an upstream dam during the rainy season interferes with the tidal effect and causes drastic changes in salinity (Díez-Minguito et al., 2013). Therefore, under heavy and prolonged rainfall events, salinity along the estuary experiences large decreases, although during the dry period (70% of the hydrological year) low discharges regularly occur to compensate for evaporation losses and the estuary is tidally dominated.

Consequently, salt marshes located closer to the river mouth (S1:S4) generally had higher salinities compared to those located on the land strip upstream of the estuary (S5:S7). With the exception of a single value of 33.6 psu in March 2016, salinity always remained below 5 psu in Site 7. Based on salinity, S1:S4 could be considered polyhaline tidal marshes whereas S5:S7 formed a group of oligo-mesohaline marshes. Clear variations for most other variables were also observed in all sites during the different samplings (Table 1).

According to Chl a, the salt marshes may be categorized as meso to eutrophic systems, with the higher pigment concentrations mostly found in the mesohaline sites (Table 1). This was also the case for DOC and nutrients (PO43- and DIN), which appeared at higher levels in S5:S7. Possibly related to the presence of high phytoplankton biomass, the flooded marshes were always well-oxygenated, and hypoxia was not detected during the samplings, with DO levels keeping values above 5 mg l⁻¹ (Table 1). Results of the Pearson correlations performed with the complete dataset of the variables measured (Table 2), confirmed that the chlorophyll concentration was negatively and significantly (p < 0.05) correlated with AOU and was therefore, directly related to the dissolved oxygen concentration.

Besides direct precipitation, the only water source to the Doñana marshland is the adjacent estuary and thus, the increased salinity in the flooded marshes reflected the effective entry of estuarine waters through the tides. Salinity was significantly (p < 0.05) and positively correlated with TSM, DOC and TDN, indicating that the dissolved and suspended matter measured in the marshes could partially have a fluvial origin.

Dissolved methane ranged between 5.7 and 12,290 nmol l⁻¹ and generally increased in the mesohaline marshes (S5:S7, Table 1). The spatio-temporal variability of CH₄ in all sites is shown in Figure 3 where salinity has also been plotted. Methane values at site S1 were below 18 nmol l⁻¹ during the entire sampling period, except for a single value of 22 nmol l⁻¹ observed in March 2018, coinciding with a drop in salinity (~4 psu) due to heavy monthly precipitations and a moderate estuarine intrusion (marked by the tidal coefficient) during the sampling date (Figure 2B). The highest methane levels in this site were observed when salinity exhibited values below 15 psu, as it occurred in May and December of 2016, March 2017 and March 2018 (Figure 3).

In S2, CH₄ was invariably higher than in S1 with a median value of 278 nmol l⁻¹ (Table 1 and Figure 3). A minimum of around 50 nmol l⁻¹ was measured in March, April and December 2016 when salinity was above 14 psu. A second minimum of 100 nmol l⁻¹ was observed in March 2018 despite the low salinity present in the flooded marshes, which may be attributable to a dilution effect caused by the high monthly rainfall (Figure 2B). The maximum concentration of dissolved methane in this site (617 nmol l⁻¹) was detected in the summer of 2017 at a salinity of 19 psu. Therefore, even though salinity was likely to affect methane production, other processes were involved in CH₄ dynamics in this particular sector of the marshes.

The range of variation of dissolved methane in S3 decreased with respect to that in S2, oscillating between 13.1 and 319.3 nmol l⁻¹ (Table 1) and with an average concentration of 104 nmol l⁻¹. Apparent relationship with salinity was not evidenced in this site, as maximum levels of methane were not associated with decreased salinity and vice versa (Figure 3).

In S4, the average concentration of methane during the sampling period was 65 nmol l⁻¹, with a maximum value of 114 nmol l⁻¹ found in December 2016 at a salinity of 4 psu (Figure 3). A second maximum was also observed in March 2018 at a similar salinity. However, high CH₄ concentrations (~90) were still measured under significant estuarine intrusions, such as those occurring in June 2016 and May 2017 which can be evidenced by the salinity rises (~35 and ~24 psu, respectively, Figure 3).

In the mesohaline marshes (S5:S7) that were characterized by salinities below 10, except in March 2016, methane concentration clearly increased in relation to the polyhaline sites (Table 1 and Figure 3). Average concentrations of CH₄ were 338.5, 261.2, and 431.3 nmol l⁻¹ in S5, S6, and S7, respectively. However, it must be noted that an exceptionally high methane level, equivalent to nearly 12,300 nmol l⁻¹, was found in S5 during July 2017, which has been excluded of the average local computation. As a general trend, methane concentration rose in these marshes during the summer of 2017 when water temperatures reached 32°C (Figure 2A).

Overall, the influence of temperature on dissolved CH₄ was observed in all sites, as the gas levels increased during the spring and summer months, compared to the levels measured over the winters (Figure 3). The temperature dependence was confirmed by the significant (p < 0.05) and positive correlation obtained between this variable and CH₄ (Table 2). When methane concentration was plotted against temperature (Figure 4A), the thermal effect on CH₄ levels was especially obvious in the mesohaline marshes, where sulfate inhibition is likely to be suppressed.

In contrast, no statistical correlation between methane and salinity was found (Table 2). Nonetheless, when the relationship between the gas saturation levels in the salt marshes and salinity was more closely examined, separately considering the polyhaline (S1:S4) and mesohaline sites (S5:S7), two different behaviors of CH₄ dynamics and their relation to this variable were distinguished (Figure 4B). Polyhaline marshes were characterized by a CH₄ saturation range of 252–36,735% (average 5,170%) and a scattered distribution with salinity (Figure 4B). On the contrary, mesohaline marshes exhibited a larger CH₄ saturation range (374–620,007%, average 31,541%) and a positive relationship with salinity was noticeable (Figure 4B).

In order to better elucidate the effect of ecosystem metabolism on the dynamics of CH₄ in the salt marshes, the relationship between CH₄ saturation and the oxygen availability was also considered. As shown in Figure 4C, log₁₀%CH₄ vs. %O₂ displayed a scattered distribution but a positive relationship, indicating that high levels of CH₄ were present in well-oxygenated inundated marshes.

Moreover, dissolved CH₄ concentration was shown to be strong and significantly (p < 0.05) correlated with Chl a (Table 2). The increase of methane with the phytoplankton biomass was evident in both the polyhaline and mesotrophic sites (Figure 5A) and in the mesohaline and eutrophic sites (Figure 5B). The trend found between DOC and Chl a in the salt marshes (Figures 5C,D) suggests that (Figures 5C,D) organic matter in the water column may have originated partially from phytoplankton productivity. Nonetheless, the fact that high values of DOC (>10 mg l⁻¹) were still present under a low Chl a concentration (<1.5 μg l⁻¹) in the two groups of sites, indicate the contribution of external sources to the organic matter content in the marshes. Increases in organic and suspended matter were accompanied by rises in dissolved CH₄ (Figures 5E–H). Significant and positive correlations with salinity were found for both variables (Table 2) and between salinity and TDN, supporting the hypothesis that estuarine inputs also contribute to organic matter accumulation and suspended matter entry in salt marshes.

The pCO₂ levels in the marshes were also well-correlated to the log of %CH₄ and to %O₂ (Figures 6A,B). pCO₂ values spanned three orders of magnitude, ranging between 0.5 and 4,427 μatm in S3 and S5, respectively (Table 1). Extremely low pCO₂ levels, could be attributed to the high pH values measured in the marshes (~10). The inundated mudflat remained undersaturated with respect to atmospheric CO₂ (average concentration of ~405.4 μatm during the sampling period estimated with the Mauna Loa trends, https://www.esrl.noaa.gov/gmd/ccgg/trends/) half of the time. In particular, S3 always behaved as a sink for atmospheric CO₂ (Table 2).

CH₄ emissions exhibited a marked variability during the sampling period in all sites, varying from ~3 μmol m⁻² d⁻¹ in site 1 in March 2016, to a maximum of ~12,700 μmol m⁻² d⁻¹ in site 5 in July 2017 (Figure 7). Overall, a tendency of higher methane effluxes was observed during the summer months, although emissions increased in the mesohaline marshes (from 5.6 to 12,715 μmol m⁻²d⁻¹ ,average 637 μmol m⁻² d⁻¹) in relation to the polyhaline marshes (from 2.6 to 720 μmol m⁻²d⁻¹, average 104 μmol m⁻² d⁻¹).

Coastal wetlands, and particularly estuarine systems are intense sources of CH₄ to the atmosphere (Borges and Abril, 2011; Upstill-Goddard and Barnes, 2016) where emissions are sustained by fluvial inputs and sedimentary methanogenesis, fueled by a high organic matter deposition (Borges et al., 2017; Rosentreter et al., 2018a).

We recently reported a data set of dissolved CH₄ concentrations in the Guadalquivir estuary transect that flooded the Doñana salt marshes from March 2016 to March 2017 (Huertas et al., 2018). For this period, that partially encompasses the one used in this study, the CH₄ concentration ranged from 14 nmol l⁻¹ in the river mouth (average salinity during the study >25) to 750 nmol l⁻¹ upstream to the limit of site 5 (average salinity ~5). These levels of dissolved methane corresponded to a saturation range of 520–30,800% (average 2,285%) (Huertas et al., 2018), which is comparable to that measured in the Doñana polyhaline marshes (S1:S4) corresponding to 252–36,735% (average 5,170%) but much lower than that observed in the mesohaline marshes (374–620,007%, average 31,541%). This indicates that even though estuarine inputs of CH₄ in the salt marshes may occur by lateral transport, through riverine overflow, a local additional source of CH₄ must contribute to the observed high CH₄ values. The methane over-saturations and concentrations measured in Doñana fall within the range reported in temperate estuarine systems characterized by a salinity gradient (Upstill-Goddard and Barnes, 2016; Borges et al., 2018b) and in mangroves (Rosentreter et al., 2018b and references therein).

The most likely additional source of CH₄ to the Doñana tidal flat would be sedimentary methanogenesis. Warm and waterlogged soils provide ideal conditions for methanogenesis (Roehm, 2005). Deep standing water enriched in dissolved organic matter and covered by plant litter creates an anaerobic environment that triggers CH₄ production (Ding and Cai, 2007), with this process being favored by temperature increases. Our data show a direct and strong relationship between the levels of dissolved methane, organic matter content and water temperature, as it has been observed in other brackish coastal wetlands (Welti et al., 2017). Methanogenesis is substantially driven by temperature (Martin and Moseman-Valtierra, 2017) and it is therefore the most dominant environmental driver of CH₄ emissions in wetlands (Yvon-Durocher et al., 2014) specifically in coastal salt marshes (Abdul-Aziz et al., 2018). The overall pattern found here therefore suggests a methanogenesis control by temperature in the sediment, which has been described as muddy with numerous stems and roots and scarce micro-fauna (Ruiz et al., 2004). In addition, vegetation in the marshes may have supported a methanogens population that could contribiute to an increase of CH₄ concentrations during warmer conditions, which has been found in Spartina dominated wetlands (Burke et al., 2002; Abdul-Aziz et al., 2018). The organic-matter rich sediment in conjunction with the organic inputs from the heterotrophic Guadalquivir estuary (Flecha et al., 2015) and those generated in situ by the phytoplankton productivity, as indicated by the association found between Chl a and DOC, provide optimal conditions for sedimentary methanogenesis to proceed, with temperature acting as a major regulating driver. The higher temperatures during the late spring and summer months possibly resulted in increased microbial activity and greater primary productivity (Chl a and temperature were significantly and positively correlated), leading to high organic inputs, sediment carbon contents, and subsequently higher substrate availability. Large suspended matter loads that are characteristic in the turbid Guadalquivir River (Flecha et al., 2015) would also provide higher silt deposition and enhance methane generation, in agreement with what has been observed in shallow coastal areas (Borges et al., 2017). The pattern between TSM and dissolved CH₄ in the salt marshes supports this hypothesis, particularly in the mesohaline marshes where very high levels of dissolved methane were measured at a high Chl a concentration, DOC content and TSM.

The couplings between pCO₂, CH₄, DOC, and %O₂ also indicate that the dynamics of dissolved gases in the tidal flat were driven by the oxidation of the organic matter produced by either phytoplankton, as described elsewhere (Borges et al., 2017; Welti et al., 2017) or from allochthonous sources (Upstill-Goddard et al., 2017). Here, the relationship observed between pCO₂ and CH₄ is quite frequent (Borges et al., 2015, 2018a; Rosentreter et al., 2018b) and denotes a common source and similar production dynamics, as both gases are produced through organic matter respiration.

In the polyhaline marshes, large CH₄ concentrations were measured when DOC levels were high, but this association was not as clear as in the mesohaline sites. This was probably due to the general inhibitory effects of abundant sulfate electron acceptors and reducing bacteria generally found in saline environments (Alongi et al., 2001; vanDijk et al., 2015). In fact, even though we did not observe a linear relationship between salinity and CH₄ levels, coinciding with the trend described in other brackish systems (Chuang et al., 2017; Welti et al., 2017), the lower CH₄ levels in the Doñana salt marshes were mainly associated to high salinity values, consistent with existing literature on coastal wetlands (Bartlett et al., 1987; Poffenbarger et al., 2011). Nevertheless, large CH₄ levels can still occur in wetlands even when sulfate reduction is relevant (Lee et al., 2008; Alongi and Brinkman, 2011).

Similarly, distinct associations between methane and oxygen have been reported in aquatic systems. It is certain that methanogens are obligate anaerobes, and thus a negative relationship %CH₄ vs. %O₂ would be expected, as previously shown in other coastal wetlands (Livesley and Andrusiak, 2012), including the Guadalquivir estuary (Huertas et al., 2018). But alternatively, positive associations have been described in Amazon floodplain lakes (Devol et al., 1990), in freshwater lakes (Bogard et al., 2014), African savannah rivers (Upstill-Goddard et al., 2017) and shallow coastal zones (Borges et al., 2017). This feature has been attributed to several causes. First, the diffusion of CH₄ produced in the underlying sediment into the aerated water column under low levels of turbulence since the CH₄ diffusion term in these conditions, exceeds combined CH₄ losses via oxidation, and water-to-air exchange (Upstill-Goddard et al., 2017). The contribution of high macrophyte-related productivity has also been suggested as a cause for positive relationships with dissolved oxygen due to direct CH₄ production (Stanley et al., 2016). Additionally, methanogens have a higher sensitivity to temperature than methanotrophs do (Dunfield et al., 1993), which leads to higher rates of CH₄ production in relation to consumption at elevated temperatures within a certain range of salinity values. Furthermore, methanogenesis in aerobic conditions by methanogen archaea fixed at the surface of phytoplankton cells has also been reported (Grossart et al., 2011). Therefore, peaks of dissolved CH₄ that correlate positively with oxygen, phytoplankton biomass, and productivity are a recurrent feature in aquatic ecosystems (Bogard et al., 2014 and references therein).

Considering the shallowness of the tidal flat in the Doñana, it is plausible to assume that benthic productivity could have been decoupled from CH₄ production (Borges et al., 2015). Moreover, methane could also have been transported from the creek onto the marshes during the flood tide as it occurs in estuarine systems (Rosentreter et al., 2018b), whereas the DO signal would be related to the primary productivity occurring in the mudflat. Unfortunately, tidal pumping cannot be evaluated due to the lack of measurements in the creek, but it could be an important mechanism for CH₄ import to the marshes (Rosentreter et al., 2018b). The correlation of methane with a Chl a concentration may support a strong benthic CH₄ source, fueled by sinking phytoplankton material (Borges et al., 2017) rather than an aerobic CH₄ production.

Therefore, according to the patterns found in our study, temperature, and phytoplankton activity seem to be the main drivers of CH₄ formation when salinity does not exceed 17 psu (see for instance temporal changes in methane levels in site 5 and the overall relationship %sat CH₄ vs. salinity). Above this salinity threshold, sulfate inhibition would constrain methane formation in the marshes.

Nevertheless, our findings also indicate that besides sedimentary methanogenesis, inflowing estuarine inputs would provide an additional source of methane to the salt marshes. Offsetting these sources, losses by microbial oxidation, lateral advection towards the estuary and dilution by heavy rain events, must also be incorporated as mechanisms affecting temporal variability of methane in the system. Unfortunately, the relative importance of all sources and sinks of dissolved methane cannot presently be estimated.

The mean air-water CH₄ flux during the sampling period was 324 μmol m⁻² d⁻¹, although a high variation in the emissions was observed (±1,394 μmol m⁻² d⁻¹). This was a consequence of the different ranges of emissions found between the polyhaline and mesohaline marshes, as the methane fluxes were generally lower in the former compared to the latter, in agreement with what has been observed in systems characterized by marked salinity gradients (Poffenbarger et al., 2011; Welti et al., 2017). Nevertheless, the mean flux provided here falls within the ranges of CH₄ emissions reported in European estuarine zones (Upstill-Goddard and Barnes, 2016), in temperate fluvial catchments (Upstill-Goddard et al., 2017) and in mid-latitude mangroves (Rosentreter et al., 2018a). Few studies have quantified and investigated air-water CH₄ exchange and dynamics in tidal flats (Dean et al., 2018). The fluxes measured in our study are slightly lower than those observed in subtropical tidal wetlands (Tong et al., 2013; Welti et al., 2017), coherent with a temperature control on methanogenesis (Abdul-Aziz et al., 2018). Lower emissions can therefore be expected from temperate coastal wetlands compared to tropical and subtropical coastal wetlands. In fact, the higher methane fluxes in the Doñana were found during the summer months in all sites, as increasing temperature increases microbial activity. Nevertheless, methane emissions from the Doñana salt marshes are still of the same order of magnitude than those measured in Indian and Australian mangroves (Linto et al., 2014; Call et al., 2015; Rosentreter et al., 2018b).

The average air-water CH₄ flux estimated here is almost 5-fold higher than the mean value computed in the Guadalquivir estuary transect that floods Doñana and equivalent to 66.2 μmol m⁻²d⁻¹ (Huertas et al., 2018). The biogeochemistry in the estuary is totally different than that in the marshes where large loads of inorganic suspended particulate matter cause a considerable turbidity that limits phytoplankton growth (Ruiz et al., 2013). The low primary productivity in the estuarine waters would constrain labile organic matter supply for sedimentary methanogenesis, whereas in the shallow and lower turbulent salt marshes organic inputs from in situ photosynthesis, or from external sources, would fuel the process. This hypothesis however, remains to be tested.

The Doñana salt marshes behaved as a moderate CH₄ emitter under the environmental conditions present during the sampling period, which corresponds to a normal hydrological cycle for this geographical area. Cycles characterized by heavy precipitations, which are not so uncommon (Díaz-Delgado et al., 2016; Huertas et al., 2017), are expected to affect methane dynamics and increase the emissions due to the overall reduction in salinity over the inundated marshes. This effect has been observed in tropical mangrove-dominated estuaries (Rosentreter et al., 2018b). Moreover, it is important to acknowledge that the approach used here to compute methane fluxes, does not quantify CH₄ ebullition fluxes, whose contribution to total CH₄ emissions may be significant (Baulch et al., 2011), particularly during low water periods. Furthermore, our sampling regime did not account for the potential importance of lateral tidal pumping as a source of CH₄ and higher concentrations (and fluxes) during low tide, as it has been reported in other intertidal systems (Call et al., 2015; Maher et al., 2015). In fact, tidal differences may account 5-fold in the estimated CH₄ diffusive water to air fluxes (Maher et al., 2015). CH₄ may also be transferred directly from sediments to the atmosphere (Borges and Abril, 2011). At low tide, this can take place through diffusion or ebullition and in vegetated tidal flats, such as Doñana, plants may act as both active and passive conduits of the CH₄ (Foster-Martinez and Variano, 2016). All these pathways were not determined in our study and total CH₄ emissions must be seen as a lower limit.

Our study is, however, foundational because it underlines the drivers controlling CH₄ dynamics in the Doñana marshes. These wetlands are threatened by climate change and many local stressors associated to anthropogenic activities (Green et al., 2017). In the Mediterranean climatic region, higher minimum temperatures, more extreme high temperature events in summer and less precipitation are projected to occur by the end of the 21st century (Giorgi and Lionello, 2008). How the expected warming will affect methane dynamics in the ecosystem complex formed by the estuary and the marshes cannot be anticipated. Increased temperatures can cause methane to be produced at a faster rate but a rise in the saline intrusion associated to the projected sea level rise may offset methane release. Sustained observations would then be required to constrain the gas budget in this iconic coastal wetland.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.