Stability of ion pairs from gypsum solubility degree of ion pair formation between the major constituents of seawater

The stability constants of the ion pairs NaSO₄^?, KSO₄^?, MgSO₄^?, CaSO₄, MgCl⁺ and CaCl⁺ were determined at 25°C and 0.7 M formal ionic strength, by measuring the solubility of gypsum (CaSO₄ · 2H₂O) in different media. The media used contained one or two of the following electrolytes: NaCl, KCl, MgCl₂, NaClO₄, Mg(ClO₄)₂, Na₂SO₄. Values for the stability constants are 1.22, 1.84, 12.3, 30.6, 0.48 and 1.20 M^?1, respectively, and the solubility product for gypsum is 2.87 · 10^?4M². The distribution of the main constituents of seawater was calculated using these results and the values of the carbonate and bicarbonate constants given by Dyrssen and Hansson (1972–1973). The solubility of gypsum in seawater as calculated and determined experimentally was 21.43 mM and 21.10 mM, respectively.     

Stability of the calcium sulphate ion pair at the ionic strength of seawater by potentiometry

The stability of the ion pair CaSO₄ was determined from measurement of the change in calcium ion activity with medium composition at constant ionic strength. A calcium selective PVC-matrix liquid membrane electrode was used to monitor the calcium ion activity. All measurements were performed at 1 atm, 25 ± 0.1°C and 0.7 M formal ionic strength. The evaluation of the stability constant depends on the degree of complexation between calcium and chloride and between sodium and sulphate. The dependence of K_CaSO4, on K_NaSO4 and K_CaCl can be described by the following relations:K_CaSO4=17.7 K_CaCl+16.5 (K_NaSO4=1.8)K_CaSO4=18.8 K_CaCl+17.3 (K_NaSO4=2) for K_CaCl=0–1A value of K_CaSO4=25.4 is suggested.     

Characterization of multiple functional groups on kaolinite

Multiple proton functional group conditional binding constants (K′_i) and their concentrations (C_i) are determined from detailed acid/base titration data. The C-K′ information is obtained for kaolinite by assuming that the distribution can be approximated by cumulative independent mono protic groups. Linear programming optimization techniques are used to Fit the data. In addition, the electrode calibration is optimized in the data fit. This adjustment is important for high and low pHs. Discrete concentration dependent pXs of about 3.4, 4.5, 6.7 and 9.8 are determined for a reference kaolinite. These correspond reasonably well to designations made by Wehrli et al. (1990, Aquatic Sciences, Vol. 52, pp. 1–31 ) to A10H2 and AlOH proton reactions at edge and surface sites and to silanol exchange. There is an ionic strength effect for one site. Long and short reaction times and reversibility affect the results.     

Kinetic and equilibrium studies of copper-dissolved organic matter complexation in water column of the stratified Krka River estuary (Croatia)

An interaction of dissolved natural organic matter (DNOM) with copper ions in the water column of the stratified Krka River estuary (Croatia) was studied. The experimental methodology was based on the differential pulse anodic stripping voltammetric (DPASV) determination of labile copper species by titrating the sample using increments of copper additions uniformly distributed on the logarithmic scale. A classical at-equilibrium approach (determination of copper complexing capacity, CuCC) and a kinetic approach (tracing of equilibrium reconstitution) of copper complexation were considered and compared. A model of discrete distribution of organic ligands forming inert copper complexes was applied. For both approaches, a home-written fitting program was used for the determination of apparent stability constants (K_i^equ), total ligands concentration (L_iT) and association/dissociation rate constants (k_i¹,k_i^- 1).A non-conservative behaviour of dissolved organic matter (DOC) and total copper concentration in a water column was registered. An enhanced biological activity at the freshwater–seawater interface (FSI) triggered an increase of total copper concentration and total ligand concentration in this water layer. The copper complexation in fresh water of Krka River was characterised by one type of binding ligands, while in most of the estuarine and marine samples two classes of ligands were identified. The distribution of apparent stability constants (log K₁^equ: 11.2–13.0, log K₂^equ:8.8–10.0) showed increasing trend towards higher salinities, indicating stronger copper complexation by autochthonous seawater organic matter.Copper complexation parameters (ligand concentrations and apparent stability constants) obtained by at-equilibrium model are in very good accordance with those of kinetic model. Calculated association rate constants (k₁¹:6.1–20 × 10³ (M s)^− 1, k₂¹: 1.3–6.3 × 10³ (M s)^− 1) indicate that copper complexation by DNOM takes place relatively slowly. The time needed to achieve a new pseudo-equilibrium induced by an increase of copper concentration (which is common for Krka River estuary during summer period due to the nautical traffic), is estimated to be from 2 to 4 h.It is found that in such oligotrophic environment (dissolved organic carbon content under 83 µM_C, i.e. 1 mg_CL^− 1) an increase of the total copper concentration above 12 nM could enhance a free copper concentration exceeding the level considered as potentially toxic for microorganisms (10 pM).     

Dissociation constants of protonated cysteine species in seawater media

The dissociation constants (pK₁, pK₂ and pK₃) for cysteine have been measured in seawater as a function of temperature (5 to 45 °C) and salinity (S = 5 to 35). The seawater values were lower than the values in NaCl at the same ionic strength. In an attempt to understand these differences, we have made measurements of the constants in Na–Mg–Cl solutions at 25 °C. The measured values have been compared to those calculated from the Pitzer ionic interaction model. The lower values of pK₃ in the Na–Mg–Cl solutions have been attributed to the formation of Mg²⁺ complexes with Cys²⁻ anions

Mg²⁺ + Cys²⁻ = MgCys

The stability constants have been fitted to

after corrections are made for the interaction of Mg²⁺ with H⁺.The pK₁ seawater measurements indicate that H₃Cys⁺ interacts with SO₄²⁻. The Pitzer parameters β⁰(H₃CysSO₄), β¹(H₃CysSO₄) and C(H₃CysSO₄) have been determined for this interaction. The formation of CaCys as well as MgCys are needed to account for the values of pK₂ and pK₃ in seawater.The consideration of the formation of MgCys and CaCys in seawater yields model calculated values of pK₁, pK₂ and pK₃ that agree with the measured values to within the experimental error of the measurements. This study shows that it is important to consider all of the ionic interactions in natural waters when examining the dissociation of organic acids.     

Aggregation of riverine colloidal iron in estuaries: A new kinetic study using stopped-flow mixing

The kinetics of aggregation of riverine colloidal iron have been examined using a stopped-flow technique which probes the first few seconds of mixing between river and sea water end members. A significant fraction, up to 80%, of the colloidal iron is aggregated during the first 1–2 s of mixing, indicating that the aggregation process is much faster than previously thought. Most of the aggregation induced by seawater results from the divalent cations Mg²⁺ and Ca²⁺, with the overall ionic strength having little influence. At equal concentrations of 27 mM, the rate of aggregation by alkaline earth cations increased with ionic size Mg²⁺ < Ca²⁺ < Sr²⁺ < Ba²⁺. The aggregation rates were indifferent to the anion (Cl⁻ or SO₄²⁻) present. Very high aggregation rates were also induced by the common water treatment coagulant Al₂(SO₄)₃ at concentrations in the range 20–30 µM Al(III), several orders of magnitude lower that those used for the seawater cations. Our results support the view that specific chemical interactions between the cations and the colloid particle surface, rather than simple electrical effects, control the colloid stability.     

Dissociation constants of protonated methionine species in seawater media

The dissociation constants (pK₁ and pK₂) for methionine have been measured in artificial seawater as a function of salinity (S = 5 to 35) and temperature (5 to 45 °C). The seawater pK₂ values were lower than the values in NaCl at the same ionic strength while the pK₁ values in seawater were lower only at S = 35. In an attempt to understand these differences, we have made measurements of the constants in Na–Mg–Cl solutions at 25 °C. The measured values have been used to determine the formation of Mg²⁺ complexes and Pitzer interaction parameters for Mg²⁺ with methionine. The seawater model with the interaction parameters accounts for the differences between the value of pK₁ and pK₂ between NaCl and seawater. This study demonstrates that it is important to consider all of the ionic interactions in natural waters when examining the dissociation of organic acids.     

The speciation of Fe(II) and Fe(III) in natural waters

The interactions of Fe(II) and Fe(III) with the inorganic anions of natural waters have been examined using the specific interaction and ion pairing models. The specific interaction model as formulated by Pitzer is used to examine the interactions of the major components (Na⁺, Mg²⁺, Ca²⁺, K⁺, Sr²⁺, Cl⁻, SO₄⁻, HCO₃⁻, Br⁻, CO₃²⁻, B(OH)₄⁻, B(OH)₃ and CO₂) of seawater and the ion pairing model is used to account for the strong interaction of Fe(II) and Fe(III) with major and minor ligands (Cl⁻, SO₄²⁻, OH⁻, HCO₃⁻, CO₃²⁻ and HS⁻) in the waters. The model can be used to estimate the activity and speciation of iron in natural waters as a function of composition (major sea salts) and ionic strength (0 to 3 M). The measured stability constants (K_FeX^*) of Fe(II) and Fe(III) have been used to estimate the thermodynamic constants (K_FeX) and the activity coefficient of iron complexes (γ_FeX) with a number of inorganic ligands in NaClO₄ medium at various ionic strengths: In(K_FeX/γ_Feγ_X) = InK_FeX − In(γ_FeX) The activity coefficients for free ions (γ_Fe, γ_x) needed for this extrapolation have been estimated from the Pitzer equations. The activity coefficients of the ion pairs have been used to determine Pitzer parameters (B_FeX, B_FeX⁰, C_FeX^φ) for the iron complexes. These results make it possible to estimate the stability constants for the formation of Fe(II) and Fe(III) complexes over a wide range of ionic strengths and in different media. The model has been used to determine the solubility of Fe(III) in seawater as a function of pH. The results are in good agreement with the measurements of Byrne and Kester and Kuma et al. When the formation of Fe organic complexes is considered, the solubility of Fe(III) in seawater is increased by about 25%.     

Ionic activity coefficients in aqueous mixtures of NaCl and MgCl2

Methods developed earlier, based on hydration numbers for individual ionic species, have been extended to the calculation of ionic activity coefficients in aqueous systems of two electrolytes MX and NX₂ with a common unhydrated anion (X⁻). The data required include the mean activity coefficients of MX and NX₂ in the mixtures, together with the osmotic coefficient. The procedure is illustrated by a calculation of γ_Na, γ_Mg, and γ_Cl in a mixture of NaCl and MgCl₂ closely approximating the composition of seawater with salinity of 35‰.     

Depth profiles of volatile iodine and bromine-containing halocarbons in coastal Antarctic waters

Measurements of bromoform (CHBr₃), diiodomethane (CH₂I₂), chloroiodomethane (CH₂ICl) and bromoiodomethane (CH₂IBr) were made in the water column (5–100 m depth) of the Southern Ocean within 0–40 km of the Antarctic sea ice during the ANTXX1/2 transect of the German R/V Polarstern, at five locations between 70–72°S and 9–11°W in the Antarctic spring/summer of 2003–2004. Some of the profiles exhibited a very pronounced layer of surface sea-ice meltwater, as evidenced by salinity minima and temperature maxima, along with surface maxima in concentrations of CHBr₃, CH₂I₂, CH₂ICl and CH₂IBr. These results are consistent with in situ surface halocarbon production by ice algae liberated from the sea ice, although production within the sea ice followed by transport cannot be entirely ruled out. Additional sub-surface maxima in halocarbons occurred between 20 and 80 m. At a station further from shore and not affected by surface sea-ice meltwater, surface concentrations of CH₂I₂ were decreased whereas CH₂ICl concentrations were increased compared to the stations influenced by meltwater, consistent with photochemical conversion of CH₂I₂ to CH₂ICl, perhaps during upward mixing from a layer at  70 m enhanced in iodocarbons. Mean surface (5–10 m) water concentrations of halocarbons in these coastal Antarctic waters were 57 pmol l^− 1 CHBr₃ (range 44–78 pmol l^− 1), 4.2 pmol l^− 1 CH₂I₂ (range 1.7–8.2 pmol l^− 1), 0.8 pmol l^− 1 CH₂IBr (range 0.2–1.4 pmol l^− 1), and 0.7 pmol l^− 1 CH₂ICl (range 0.2–2.4 pmol l^− 1). Concurrent measurements in air suggested a sea-air flux of bromoform near the Antarctic coast of between 1 and 100 (mean 32.3, median 10.4) nmol m^− 2 day^− 1 and saturation anomalies of 557–1082% (mean 783%, median 733%), similar in magnitude to global shelf values. In surface samples affected by meltwater, CH₂I₂ fluxes ranged from 0.02 to 6.1 nmol m^− 2 day^− 1, with mean and median values of 1.9 and 1.1 nmol m^− 2 day^− 1, respectively.     

Stability constants of NaSO4−, MgSO4, MgF+, MgCl+ ion pairs at the ionic strength of seawater by potentiometry

The stability of the ion pair NaSO₄ was determined by measuring the change in sodium activity with medium composition at constant ionic strength, using a sodium-sensitive glass electrode. The stability constants of MgSO₄ and MgCl⁺ were determined indirectly from measurements of the stability of MgF⁺ in different media. All measurements were performed at 1 atm pressure, 25 ± 0.1 °C and 0.7 M formal ionic strength. The stability constants for NaSO₄^?, MgSO₄, MgF⁺ and MgCl⁺ are 1.8 ± 0.1, 6.3 ± 0.1, 22.9 ± 0.1 and 0.34 ± 0.02 M^?1, respectively.     

An improved model for the calculation of CO2 solubility in aqueous solutions containing Na, K, Ca, Mg, Cl, and SO4

An improved model is presented for the calculation of the solubility of carbon dioxide in aqueous solutions containing Na⁺, K⁺, Ca²⁺, Mg²⁺, Cl⁻, and SO₄²⁻ in a wide temperature–pressure–ionic strength range (from 273 to 533 K, from 0 to 2000 bar, and from 0 to 4.5 molality of salts) with experimental accuracy. The improvements over the previous model [Duan, Z. and Sun, R., 2003. An improved model calculating CO₂ solubility in pure water and aqueous NaCl solutions from 273 to 533K and from 0 to 2000 bar. Chemical Geology, 193: 257–271] include: (1) By developing a non-iterative equation to replace the original equation of state in the calculation of CO₂ fugacity coefficients, the new model is at least twenty times computationally faster and can be easily adapted to numerical reaction-flow simulator for such applications as CO₂ sequestration and (2) By fitting to the new solubility data, the new model improved the accuracy below 288 K from 6% to about 3% of uncertainty but still retains the high accuracy of the original model above 288 K. We comprehensively evaluate all experimental CO₂ solubility data. Compared with these data, this model not only reproduces all the reliable data used for the parameterization but also predicts the data that were not used in the parameterization. In order to facilitate the application to CO₂ sequestration, we also predicted CO₂ solubility in seawater at two-phase coexistence (vapor–liquid or liquid–liquid) and at three-phase coexistence (CO₂ hydrate–liquid water–vapor CO₂ [or liquid CO₂]). The improved model is programmed and can be downloaded from the website http://www.geochem-model.org/programs.htm.     

Accurate measurement of O2, N2, and Ar gases in water and the solubility of N2

The degree of atmospheric saturation for O₂, Ar, and N₂ gases in water can be determined to accuracies of ±0.1–0.3% using mass spectrometry to determine the gas ratios and Winkler titrations for oxygen analysis. We describe methods used to obtain this level of accuracy and precision. Oxygen accuracy of ±0.1% can be obtained by careful attention to standardization using KIO₃ standards that have been corrected for impurities. Accurate O₂/Ar and O₂/N₂ gas ratios (±0.1–0.2%) are obtainable by measuring the mass ratios against the atmosphere if the effect of different gas concentrations on the performance of the mass spectrometer are taken into account. Oxygen and argon saturation values have been determined previously to accuracies of less than or equal to ±0.1%, but published estimates of the saturation value for nitrogen differ by more than 1%. We have redetermined the N₂ saturation value at 19°C and zero salinity to be 0.92% greater than the results reported in the work of Weiss (1970).     

Methane-dominated thermochemical sulphate reduction in the Triassic Feixianguan Formation East Sichuan Basin, China: towards prediction of fatal H2S concentrations

New sour pools have recently found in the Lower Triassic Feixianguan Fm carbonate reservoirs in the East Sichuan Basin in China with H₂S up to 17.4% by volume. A recent blowout from a well drilled into this formation killed hundreds of people as a result of the percentage concentrations of H₂S. In order to assess the origin of fatal H₂S as well as the cause of petroleum alteration, H₂S concentrations and the isotopes, δ³⁴S and δ¹³C have been collected and measured in gas samples from reservoirs. Anhydrite, pyrite and elemental sulphur δ³⁴S values have been measured for comparison. The high concentrations of H₂S gas are found to occur at depths >3000 m (temperature now at 100 °C) in evaporated platform facies oolitic dolomite or limestone that contains anhydrite nodule occurrence within the reservoirs. Where H₂S concentrations are greater than 10% its δ³⁴S values lie between +12.0 and +13.2‰ CDT. This is within the range of anhydrite δ³⁴S values found within the Feixianguan Fm (+11.0 to +21.7‰; average 15.5±3.5‰ CDT). Thus H₂S must have been generated by thermochemical sulphate reduction (TSR) locally within the reservoirs. Burial history analysis and fluid inclusion data reveal that the temperature at which TSR occurred was greater than about 130–140 °C, suggesting that the present depth-temperature minimum is an artifact of post-TSR uplift. Both methane and ethane were actively involved in TSR since the petroleum became almost totally dry (no alkanes except methane) and methane δ¹³C values become significantly heavier as TSR proceeded. Methane δ¹³C difference thus reflects the extent of TSR. While it is tempting to use a present-day depth control (>3000 m) to predict the distribution of H₂S in the Feixianguan Fm, this is an invalid approach since TSR occurred when the formation was buried some 1000–2000 m deeper than it is at present. The likelihood of differential uplift across the basin means that it is important to develop a basinal understanding of the thermal history of the Feixianguan Fm so that it is possible to determine which parts of the basin have been hotter than 130–140 °C.     

On the relevance of iron adsorption to container materials in small-volume experiments on iron marine chemistry: Fe-aided assessment of capacity, affinity and kinetics

Iron chemistry in seawater has been extensively studied in the laboratory, mostly in small-volume sample bottles. However, little has been reported about iron wall sorption in these bottles. In this paper, radio-iron ⁵⁵Fe was used to assess iron wall adsorption, both in terms of capacity, affinity and kinetics. Various bottle materials were tested. Iron sorption increased from polyethylene/polycarbonate to polymethylmetacrylate (PMMA)/high-density polyethylene/polytetrafluoroethylene to glass/quartz, reaching equilibrium in a 25–70 h period. PMMA was studied in more detail: ferric iron (Fe(III)) adsorbed on the walls of the bottles, whereas ferrous iron (Fe(II)) did not. Considering that in seawater the inorganic iron pool mostly consists of ferric iron, the wall will be a factor that needs to be considered in bottle experiments.The present data indicate that for PMMA with specific surface (S)-to-volume (V) ratio S/V, both iron capacity (42 ± 16 × 10^− 9 mol/m² or 1.7 × 10^− 9 mol/L recalculated for the S/V-specific PMMA bottles used) and affinity (log K_Fe'W = 11.0 ± 0.3 m²/mol or 12.4 ± 0.3 L/mol, recalculated for the S/V-specific PMMA bottles used) are of similar magnitude as the iron capacity and -affinity of the natural ligands in the presently used seawater and thus cannot be ignored.Calculation of rate constants for association and dissociation of both Fe'L (iron bound to natural occurring organic ligands) and Fe'W (iron adsorbed on the wall of vessels) suggests that the two iron complexes are also of rather similar kinetics, with rate constants for dissociation in the order of 10 ⁻⁴–10^− 5 L/s and rate constants for association in the order of 10⁸ L/(mol s). This makes that iron wall sorption should be seriously considered in small-volume experiments, both in assessments of shorter-term dynamics and in end-point observations in equilibrium conditions. Therefore, the present data strongly advocate making use of iron mass balances throughout in experiments in smaller volume set-ups on marine iron (bio) chemistry.     

A high accuracy method for determining nitrogen, argon and oxygen in seawater

An improved gas chromatographic system was constructed to analyze oceanic dissolved N₂, Ar and O₂ with a higher accuracy and shorter analytical time. To obtain a higher accuracy of N₂, Ar and O₂ measurements, the following was added to the system: (I) an air trapping system; (II) a N₂–CO₂ trapping system after the operation of the air trapping system; (III) an active carbon column system for separating N₂ and CO₂ completely and (IV) the introduction of automatic valves controlling most of the system. Compared to previous studies, the precision of the measurements of N₂, Ar and O₂ concentrations was higher at 0.04%, 0.05% and 0.02%, respectively, and our analytical time was shorter at 600 s. Using the improved analytical technique, concentrations of N₂ (C_N2, 561.69–611.81 μmol/kg) and Ar (C_Ar, 15.126–16.238 μmol/kg), saturation states of N₂ (ΔN₂, − 5.1–0.9%) and Ar (ΔAr, − 7.0 to − 1.1%) from 0 m to 3000 m depth in the western North Pacific were observed during March 2005. Based on these data, we propose a new concept for estimating the amount of bubble injection (B). The total error in calculating B was estimated to be about 20%. We estimated B from 12 to 43 μmol/kg in this region using the observational values of N₂ and Ar. As each water mass had a significantly different value of B even with an error of 20%, it is possible to use it as an index of sea surface state for when each water mass is produced in the sea surface mixed layer. Moreover, based on our values of B, we estimated preformed dissolved oxygen (DO) (C_preDO, 309–332 μmol/kg) and the saturation state of C_preDO (ΔpreDO, − 7.0 to − 1.2%) in this region. Thus, the difference between C_preDO and DO content in the ocean interior may be a more useful index for biogenic organic decomposition in the ocean field compared to Apparent Oxygen Utilization (AOU). Until now, the estimation of oceanic uptake of anthropogenic CO₂ has used AOU as a major parameter. Therefore, it may be necessary to re-evaluate the oceanic uptake of anthropogenic CO₂ based on our new concept of B.     

Diffusion of seawater ions. Part II. The role of activity coefficients and ion pairing

A theoretical evaluation of basic thermodynamic relationships reveals that variation of activity coefficients, ion pairing and electrical interactions must be considered when modelling ionic diffusion in seawater. The contributions of ion-pair formation and change in activity coefficient along the diffusion path were studied experimentally by conducting diffusion experiments in which solutions of KCl, NaCl, MgCl₂, Li₂SO₄, K₂SO₄, Na₂SO₄ and MgSO₄, at an ionic strength of 0.7, were allowed to diffuse into distilled water. The study reveals that the thermodynamic factor, required to correct for changes in the activity coefficient along the diffusion path, is significant for all the salts studied. Agreement between a simple diffusion model, which does not include ion pairing, and observed data was good for completely dissociated salts, but poor for salts which are known to form ion pairs at the concentration levels studied. The diffusion of MgSO₄, 0.425 of which is associated at I = 0.7, was successfully modelled by assuming that the diffusion coefficient of the MgSO₄⁰ ion pair is different from the diffusion coefficient of the dissociated salt. The diffusion coefficient of this ion pair is estimated to be 1.9 × 10⁻⁵ cm² s⁻¹ at 30°C, as compared to 0.49 × 10⁻⁵ cm² s⁻¹ for the dissociated salt. It is suggested that the high mobility of this ion pair could cause magnesium enrichment in pore water of sulfate depleted sediments.     

Rapid sediment accumulation and microbial mineralization in forests of the mangrove Kandelia candel in the Jiulongjiang Estuary, China

Rates of sediment accumulation and microbial mineralization were examined at three Kandelia candel forests spanning the intertidal zone along the south coastline of the heavily urbanized Jiulongljiang Estuary, Fujian Province, China. Mass sediment accumulation rates were rapid (range: 10–62 kg m⁻² y⁻¹) but decreased from the low- to the high-intertidal zone. High levels of radionuclides suggest that these sediments originate from erosion of agricultural soils within the catchment. Mineralization of sediment carbon and nitrogen was correspondingly rapid, with total rate of mineralization ranging from 135 to 191 mol C m⁻² y⁻¹ and 9 to 11 mol N m⁻² y⁻¹; rates were faster in summer than in autumn/winter. Rates of mineralization efficiency (70–93% for C; 69–92% for N) increased, as burial efficiency (7–30% for C; 8–31% for N) decreased, from the low-to the high-intertidal mangroves. Sulphate reduction was the dominant metabolic pathway to a depth of 1 m, with rates (19–281 mmol S m⁻² d⁻¹) exceeding those measured in other intertidal deposits. There is some evidence that Fe and Mn reduction-oxidation cycles are coupled to the activities of live roots within the 0–40 cm depth horizon. Oxic respiration accounted for 5–12% of total carbon mineralization. Methane flux was slow and highly variable when detectable (range: 5–66 μmol CH₄ m⁻² d⁻¹). Nitrous oxide flux was also highly variable, but within the range (1.6–106.5 μmol N₂O m⁻² d⁻¹) measured in other intertidal sediments. Rates of denitrification were rapid, ranging from 1106 to 3780 μmol N₂ m⁻² d⁻¹, and equating to 11–20% of total sediment nitrogen inputs. Denitrification was supported by rapid NH₄ release within surface deposits (range: 3.6–6.1 mmol m⁻² d⁻¹). Our results support the notion that mangrove forests are net accumulation sites for sediment and associated elements within estuaries, especially Kandelia candel forests receiving significant inputs as a direct result of intense human activity along the south China coast.