# Regulation of nitrogen fixation from free-living organisms in soil and leaf litter of two tropical forests of the Guiana shield

Leandro van Langenhove, Thomas Depaepe, Sara Vicca, Joke van den Berge, Clément Stahl, Elodie A Courtois, James Weedon, Ifigenia Urbina, Oriol Grau, Dolores Asensio, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Leandro van Langenhove, Thomas Depaepe, Sara Vicca, Joke van den Berge, Clément Stahl, et al. Regulation of nitrogen fixation from free-living organisms in soil and leaf litter of two tropical forests of the Guiana shield. Plant and Soil, 2019, 10.1007/s11104-019-04012-1. http-02282288

HAL Id: hal-02282288

https://hal.science/hal-02282288v1

Submitted on 21 Aug 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Regulation of nitrogen fixation from free-living organisms in soil and leaf litter of two tropical forests of the Guiana shield

Leandro Van Langenhove · Thomas Depaepe · Sara Vicca · Joke van den Berge · Clement Stahl · Elodie Courtois · James Weedon · Ifigenia Urbina · Oriol Grau · Dolores Asensio · Josep Peñuelas · Pascal Boeckx · Andreas Richter · Dominique Van Der Straeten · Ivan A. Janssens

Received: 29 June 2018 / Accepted: 25 February 2019 / Published online: 1 April 2019
© The Author(s) 2019

## Abstract

*Background and aims* Biological fixation of atmospheric nitrogen (N₂) is the main pathway for introducing N into unmanaged ecosystems. While recent estimates suggest that free-living N fixation (FLNF) accounts for the majority of N fixed in mature tropical forests, the controls governing this process are not completely

understood. The aim of this study was to quantify FLNF rates and determine its drivers in two tropical pristine forests of French Guiana.

*Methods* We used the acetylene reduction assay to measure FLNF rates at two sites, in two seasons and along three topographical positions, and used regression analyses to identify which edaphic explanatory variables,

Responsible Editor: Bernardo M Flores.

**Electronic supplementary material** The online version of this article (https://doi.org/10.1007/s11104-019-04012-1) contains supplementary material, which is available to authorized users.

L. Van Langenhove (✉) · S. Vicca · J. van den Berge ·

E. Courtois · I. A. Janssens

Centre of Excellence PLECO (Plants and Ecosystems),
Department of Biology, University of Antwerp, Wilrijk, Belgium
e-mail: leandro.vanlangenhove@uantwerpen.be

T. Depaepe · D. Van Der Straeten

Laboratory of Functional Plant Biology, Department of Biology,
Ghent University, Ledeganckstraat 35, B-9000 Ghent, KL,
Belgium

C. Stahl

INRA, UMR Ecology of Guiana Forests (Ecofog),
AgroParisTech, Cirad, CNRS, Université des Antilles, Université
de Guyane, 97387 Kourou, French Guiana

E. Courtois

Laboratoire Ecologie, évolution, interactions des systèmes
amazoniens (LEEISA), Université de Guyane, CNRS, IFREMER,
French Guiana, 97300 Cayenne, France

J. Weedon

Department of Ecological Science, Vrije Universiteit Amsterdam,
De Boelelaan 1085, 1081 HVAmsterdam, The Netherlands

I. Urbina · O. Grau · D. Asensio · J. Peñuelas

CSIC, Global Ecology Unit CREAF-CSIC-UAB,
08193 Bellaterra, Catalonia, Spain

I. Urbina · O. Grau · D. Asensio · J. Peñuelas

CREAF, 08193 Cerdanyola del Vallès, Catalonia, Spain

P. Boeckx

Department of Applied Analytical and Physical Chemistry,
Faculty of Bioscience Engineering, Isotope Bioscience Laboratory
– ISOFYS, Ghent University, Coupure Links 653, 9000 Ghent,
Belgium

A. Richter

Department of Microbiology and Ecosystem Science, University
of Vienna, Althanstr. 14, 1090 Vienna, Austria
including carbon (C), nitrogen (N), phosphorus (P) and molybdenum (Mo) content, pH, water and available N and P, explained most of the variation in FLNF rates. *Results* Overall, FLNF rates were lower than measured in tropical systems elsewhere. In soils seasonal variability was small and FLNF rates differed among topographies at only one site. Water, P and pH explained 24% of the variation. In leaf litter, FLNF rates differed seasonally, without site or topographical differences. Water, C, N and P explained 46% of the observed variation. We found no regulatory role of Mo at our sites.

*Conclusions* Rates of FLNF were low in primary rainforest on poor soils on the Guiana shield. Water was the most important rate-regulating factor and FLNF increased with increasing P, but decreased with increasing N. Our results support the general assumption that N fixation in tropical lowland forests is limited by P availability.

**Keywords** Free-living nitrogen fixation · Tropical forest · French Guiana · Nutrients · Phosphorus · Molybdenum

Introduction

Nitrogen (N) availability is a limiting factor for plant
growth in a wide range of terrestrial ecosystems world-
wide (LeBauer and Treseder 2008) and restricts the
amount of carbon (C) that can be assimilated and stored
in the terrestrial biosphere (Hungate et al. 2003;
Penuelas et al. 2013; Zaehle et al. 2015). Globally,
biological N fixation is the most important natural path-
way for introducing previously inert N - namely N₂ gas
-, into unmanaged ecosystems (Galloway et al. 1995;
Vitousek et al. 2013). Nitrogen is fixed by microorgan-
isms, known as diazotrophs, that can reduce gaseous N
(N₂) into ammonia. Diazotrophs are divided in two
groups: symbiotic and free-living. Symbiotic N fixers
are generally found in root nodules and live in a mutu-
alistic relationship with higher plants that allocate C to
the N fixers in exchange for N. Free-living N fixers are
hetero- or autotrophic bacteria or archaea inhabiting
water, soil, rocks and leaf litter (Dynarski and Houlton
2018). Global terrestrial N inputs from biological N
fixation have been estimated at 60 Tg y⁻¹ (Vitousek
et al. 2013), and biome-level comparisons suggest that
tropical rain forests may fix more N than any other
unmanaged ecosystem (Galloway et al. 2004; Reed
et al. 2011). Until recently, a large proportion of tropical

N fixation was attributed to symbiotic organisms due to
the high abundance of leguminous trees (Losos and
Leigh 2004), typically associated with symbiotic N
fixers. The relative importance of symbiotic N fixation,
however, has been questioned because it is facultative
(Menge et al. 2009) and may decline to near zero in
mature tropical rainforests (Barron et al. 2011;
Batterman et al. 2013; Sullivan et al. 2014). In this
context, N fixation by free-living organisms is thus
increasingly considered to be more important, with re-
cent estimates from multiple rain forests suggesting that
a substantial amount of N is indeed fixed by free-living
diazotrophs (e. g. Reed et al. 2007, Cusack et al. 2009,
Černá et al. 2009, Barron et al. 2009, Wurzburger et al.
2012, Matson et al. 2014).

Because of its relevance for ecosystem functioning, N
fixation and its rate-controlling factors have been the focus
of previous research. Nitrogen fixation is energetically
expensive, requiring much ATP for both the reaction itself
(Gutschick 1981) and for maintaining an oxygen-poor
intracellular environment for nitrogenase, the enzyme re-
sponsible for N fixation (Robson and Postgate 1980). The
oxygen-poor environment can be created, either by an
upregulated respiration rate, or developing cellular struc-
tures that can limit the entry of oxygen (Robson and
Postgate 1980). This regulation is needed because oxygen
binds to nitrogenase and inactivates the enzyme. Free-
living N fixation (FLNF) is thought to decrease as N
availability increases in the environment (Hartwig 1998),
because uptake of inorganic N is less costly than N fixation
(Gutschick 1981; Menge et al. 2009)).

Low phosphorus (P) availability has been reported to
limit FLNF in several tropical environments (e.g.,
Pearson and Vitousek 2002; Reed et al. 2007), an ob-
servation often attributed to the diazotrophs's high P
requirement. Nonetheless, studies of P addition to trop-
ical forests have reported contradictory responses. Some
reported that FLNF in soils and plant litter increased in
response to P additions (Benner et al. 2007; Reed et al.
2007), while others found no effect (Pérez et al. 2008;
Barron et al. 2009). Recent studies suggest that the
stimulation of N fixation by P addition may be due to
molybdenum (Mo) contamination of the commonly
used superphosphate fertiliser (Barron et al. 2009;
Wurzburger et al. 2012). Molybdenum is a rock-
derived trace element required to produce the FeMo
cofactor necessary for the functioning of the most com-
mon group of nitrogenases (Igarashi and Seefeldt 2003).
Molybdenum is thus associated more with fundamental
enzymatic requirements than with the high energy con-
sumption of N fixation, but limitation of both Mo and P
has been documented in some forests (Barron et al.
2009; Wurzburger et al. 2012; Reed et al. 2013).

Energy in the form of external organic carbon (C) is
another likely factor regulating the activity of free-living
N fixers. Like other heterotrophic microbes, diazotrophs
rely on extracellular organic C for respiration, and C
supply may be more limiting than nutrient availability in
free-living fixers (Hofmockel and Schlesinger 2007),
even though the C cost of FLNF is not well-quantified
(Dynarski and Houlton 2018). Lastly, seasonal variation
in soil moisture can play a large role in the regulation of
N fixation rates, because the rates are positively corre-
lated with moisture content (Roskoski 1980). Seasonal
differences in N fixation have indeed been observed, but
the direction of these seasonal changes differed among
studies. Studies conducted at different tropical sites have
reported both higher (cf Reed et al. 2007) and lower (cf
Matson et al. 2014) rates in the wet season compared to
the dry season with likely factors besides moisture, such
as changes in oxygen supply, causing the discrepancy in
results.

Lowland tropical forests found on highly weathered old soils are typically considered to be N-rich, because they can accumulate, recycle and export large quantities of N (Vitousek and Sanford 1986; Hedin et al. 2009; Brookshire et al. 2012). This open N cycle is also corroborated by enriched δ¹⁵N values of soils and plant tissues, due to high gaseous and leaching losses of depleted ¹⁵N sources (Amundson et al. 2003). Simultaneously, these lowland tropical forests have usually been described as P limited because of biogeochemical theory predicting that P limitation should be prevalent in old, strongly weathered soils (Walker and Syers 1976; Wardle et al. 2004). There is also a wealth of indirect evidence including high N availability (Brookshire et al. 2012), high N:P ratios in leaves (Vitousek 1984) and correlations between forest properties and soil fertility at continental scale (Quesada et al. 2012; ter Steege et al. 2006. Locally, however, nutrient availability in tropical landscapes can vary with topography, but the magnitude and direction of this variation is variable and influenced by, e.g., terrain steepness and rainfall, leading to differences in physical denudation rates and solute transportation (Porder et al. 2005; Weintraub et al. 2015). Differences in soil type and redox status along topographical gradients may also affect nutrient availability (Tiessen et al. 1994).

Here we present results from a study carried out in mature tropical rainforests of French Guiana, where we measured FLNF in both soil and leaf litter. The rolling hills of French Guianese tropical forests, part of the Guiana Shield and perched upon parent material that is amongst the oldest and most weathered in the world (Hammond 2005), are characterised as very poor in mineral nutrients (van Kekem et al. 1996), with topography inducing spatially heterogeneous nutrient availabilities. We studied three distinct topographical positions at two different forest sites in the wet and dry season, with the aim of maximizing the range in soil nutrients and moisture. Our goals were (I) to determine and compare rates of FLNF in soils and leaf litter to other tropical forest studies, (II) to evaluate whether or not N fixation rates differed between sites or between seasons, and (III) to identify which environmental factors (or combination thereof) best explained the spatial and seasonal variation in N fixation rates.

**Materials and methods**

Study sites

The study was conducted at two primary forest sites in French Guiana, in the research stations of Paracou and Nouragues. Paracou is situated 15 km from the coast (5°15′N, 52°55′W), while Nouragues is located inland (4°05′N, 52°41′W). Annual rainfall quantities at both sites were similar, with Paracou receiving an average of 3100 mm year⁻¹ for the period 2004–2015 (Aguilos et al. 2018) and Nouragues receiving an average of 2990 mm year⁻¹ (Bongers et al. 2001). Mean annual air temperature is near 26 °C for both sites (Gourlet-Fleury et al. 2004; Bongers et al. 2001). The French Guianese climate is characterized by a wet and a dry season due to the north/south movement of the Inter-Tropical Convergence Zone. The region receives heavy rains from December to July and a dry period, typically characterized by less than 100 mm rainfall month⁻¹, from August to November (Aguilos et al. 2018).

Topography at both sites is undulating with maxi-
mum slopes of approximately 15°. The elevational dif-
ference between hill summits and valleys is 20–50 m
over horizontal distances of 200–400 m. Soils at the
Paracou site are schist soils with veins of pegmatite
along the Bonidoro series, a Precambrian metamorphic
formation (Epron et al. 2006). Soils at the Nouragues
site are also derived from the same Bonidoro series, and consist of mainly Caraib gneiss (Bongers et al. 2001). This Precambrian geological substrate is particularly low in P content compared to the generally younger, nutrient-richer soils of western Amazonia (Hammond 2005; Grau et al. 2017) and therefore soils at both sites are classified as nutrient-poor Acrisols (FAO-ISRIC-ISSS 1998). Soils at Paracou range from loamy sand to sandy loam, while soils at Nouragues contain much more clay and span the range of sandy loam to silty clay according to the USDA texture classification chart (Fig. S1).

## Experimental design

To maximize the natural variation in soil nutrient availability the experimental plots encompassed a topographical gradient. At each forest site twelve plots were installed, distributed over four hillslopes, with each hillslope having three plots located at three topographical levels: (1) bottom, i.e. just above the creek running through the valley, (2) slope, i.e. the intermediate section of the elevation, and (3) top, i.e., where the slope evens out and becomes the hilltop. Each plot was of 20 × 20 m in size. In total, 24 plots spread over two sites and three topographical categories per site were studied. Distances between the plots were 10–100 m and in each plot five soil and litter samples were collected. Four samples were collected at 2 m distance from each corner and a fifth sample was taken in the middle of the plot. A total of 120 sampling points (2 study sites, 3 topographical categories per site, 4 plots per category and 5 samples per plot) were sampled for both soil and leaf litter in both the wet and dry season (240 samples in total). Courtois et al. (2018) have provided more detailed information on the experimental design.

## N fixation

Nitrogen fixation rates were determined using the acetylene reduction assay (Hardy et al. 1968). Samples were collected in May and September 2016, in the wet and dry season, respectively. Samples of leaf litter were collected manually from the soil surface and soil samples were collected with a 2-cm diameter corer to a depth of 5 cm after removing all litter from the surface.

All samples were placed in clear 100 ml borosilicate jars. The jars were sealed with rubber septa and 10 ml of air was replaced with 10 ml of acetylene gas (welding

grade, Air Liquide) to create a 10% headspace concentration by volume. The samples were then incubated in situ at ambient forest light (no direct sunlight) and temperature for 18 h. Sample moisture was not manipulated in any way, but was determined after the incubation as the weight loss after oven drying at 70 °C during 48 h.

After incubation, a subsample from the sample headspace was injected into a pre-evacuated 12 ml borosilicate vial (Labco Limited, Ceredigion, UK) and shipped to Ghent, Belgium for analysis. Ethylene concentrations were measured using laser-based photo-acoustic spectroscopy (ETD-300, Sensor Sense, The Netherlands). Parallel acetylene blanks (no leaf litter or soil) were created to assess background levels of ethylene in the acetylene gas (1.5 – 0.4 nl ethylene ml⁻¹ acetylene gas), which were subtracted from the sample ethylene concentrations. Controls for ethylene production in the soil or litter in the absence of acetylene gas were also assayed and were consistently below the detection limit of 0.01 nl ethylene ml⁻¹ air. Soil and leaf litter samples that over the incubation time produced ethylene concentrations below the detection limit (0.01 nl ethylene ml⁻¹ air) were recorded as half of this value.

We converted the rate of ethylene production, expressed as nmol ethylene g⁻¹ sample h⁻¹, into N fixation rates, expressed as kg N fixed ha⁻¹ y⁻¹ using the densities of the soil (to a depth of 5 cm) or leaf litter, and the theoretical ratio of 3 mol ethylene produced per mole N fixed (Benner et al. 2007; Cusack et al. 2009; Matson et al. 2014; Pearson and Vitousek 2001; Reed et al. 2007). The latter being based upon the conclusion that reducing 3 mol of acetylene to ethylene is equivalent to the 6 electron transfer involved in the reduction of one mole of N₂ to ammonium (Seitzinger and Garber 1987). We attempted to measure uptake of ¹⁵N through the pool dilution assay (Furnkranz et al. 2008) in a subset of soil samples to gauge the actual ratio of moles ethylene produced per mole N fixed, but we failed to do so due to a combination of low soil FLNF rates and high background levels of N in these environmental samples. Other authors have reported encountering similar issues and could not measure ¹⁵N uptake in soil (Matson et al. 2014) or leaf litter (Menge et al. 2009) samples. Soil samples for determining bulk density to a depth of 5 cm were taken with Kopeck rings. The samples were oven dried at 105 °C for 48 h and bulk density was calculated by dividing soil weight by Kopeck ring volume. Leaf litter was collected in a 0.5 m² wooden frame at five
locations per plot and once per season for determining
litter density. The litter was dried at 105 °C for 48 h and
the density was calculated by dividing weight by area
(kg ha⁻¹).

Chemical analyses

Total C, N, P and Mo contents in the soil and litter were determined on the same samples that were used for acetylene reduction. Ethylene production was first measured (see N fixation), and afterwards samples were oven dried at 70 °C for two days and then ground in a ball mill (Retsch, Germany). Total C and N contents were determined by dry combustion with an elemental analyser (Flash 2000, Thermo Fisher Scientific, Germany). Total P (mg kg⁻¹) and Mo (mg kg⁻¹) contents were determined by the sequential digestion of the ground soil and litter samples in heated strong acid (69% HNO₃ and 30% H₂O₂), followed by analysis on an iCAP 6300 Duo ICP optical emission spectrometer (Thermo Fisher Scientific, Germany).

Soil texture and nutrient availability were determined on a composite sample made of three soil cores per sampling point, each core to a depth of 15 cm. Texture was determined on fresh soil using the hydrometric method (Gee and Bauder 1986). Soil particles were dispersed with sodium metaphosphate and the amounts of sand, silt and clay were determined using a hydrometer. Soil samples for measuring nutrient availability were collected in May and September 2015, sieved (2 mm) and split into three subsamples. The first subsample was extracted with 1 M KCl in a 1:2.5 w:v ratio for pH measurement and determination of available N. On this extract pH<sub>KCl</sub> was measured with a pH meter (HI 2210–01, Hanna Instruments, USA), after which the extract was filtered through a 42 µm filter and the filtrate's concentration of NH<sub>4</sub><sup>+</sup> and NO<sub>3</sub><sup>-</sup> was determined colorimetrically (SAN++ continuous flow analyser, Skalar Inc., The Netherlands). Available N (mg kg<sup>-1</sup>) was defined as the sum of the NH<sub>4</sub><sup>+</sup> and NO<sub>3</sub><sup>-</sup> concentrations. The second subsample was used to determine the gravimetric water content by measuring weight loss after drying at 105 °C during 48 h. The third subsample was used to determine extractable P and Mo. Soils were oven dried at 60 °C for 48 h after which extractable P was determined with Bray-P acid fluoride extraction (Bray and Kurtz 1945). Available Mo was determined through resin extraction (Wurzburger et al. 2012). The samples were mixed with water in a 1:6 ratio

and five 2 cm² strips of anion-exchange membrane (VWR Chemicals, USA) were added. This mixture was stirred for 24 h and the strips were then rinsed and eluted with 10% HNO3. Available P and Mo contents were determined with a iCAP 6300 Duo ICP optical emission spectrometer (Thermo Fisher Scientific, Germany).

Literature comparison

To compare the FLNF rate results we found to those of other authors we performed a database search similar to the search carried out by Dynarski and Houlton (2018). We searched Web of Science using the keywords nitrogen, free-living, asymptotic, fixation, soil, leaf litter and tropical forest. From the resulting studies we selected those that were performed in natural terrestrial tropical ecosystems. For studies that presented results from multiple time points or seasons, we averaged the data over the course of the study period. Studies that did not report any measure of variance were assigned a standard error of ¼ of the mean. Reported rates of FLNF were converted to nmol ethylene produced g⁻¹ dry substrate h⁻¹ in order to compare N fixation rates between studies. When N fixation rates were presented on a per area basis we used the bulk density of the N fixation substrate (soil or leaf litter) and the ethylene: N₂ conversion factor indicated in the study to convert to nmol ethylene produced g⁻¹ dry substrate h⁻¹. When no conversion factor was indicated, we assumed the standard 3: 1 conversion factor (Hardy et al. 1968). Results of this database search are summarised in supplementary information Table S4.

Data analysis

To assess the differences in N fixation rates and soil and
leaf litter variables between sites and seasons we used
linear mixed effects regression models (LMER),
analysing soil and leaf litter data separately. We used
site (Paracou or Nouragues), season (Wet or Dry) and
topographical position (Bottom, Slope or Top) as fixed
factors and plot identity as a random factor. The validity
of the linear models’ assumptions (linearity, normality
of residuals, no influential outliers, homoscedasticity)
were evaluated with standard functions of R (R core
team 2017, version 3.4.3), including diagnostic plots.
Prior to analysis, data were log transformed if their
distribution was right-skewed to improve normality of
model residuals. Multiple comparisons within a factor were analysed using Tukey post hoc tests. We performed principal component analyses to visualize the correlations of previously standardised soil and leaf litter variables according to site, season and, if necessary, topography. We observed that a large proportion of our measured samples yielded ethylene production rates that fell below the detection limit (25% of samples in soil). To investigate if there were site-specific or seasonal patterns in the occurrence of values below detection limit, we transformed our ethylene production rates into binomial data (1 for measured rate and 0 for below detection limit rate) and analysed the resulting data with binomial generalized linear model (GLM) with season and site as factors.

To identify which set of physico-chemical variables significantly contributed to the observed variation in N fixation rates we performed a forward stepwise model selection, i.e. starting from a null model and retaining the predictor variable that led to the largest decrease in the Akaike information criterion (AIC), corrected for sample size (AICc). This process was iterated until no additional predictor reduced the model AICc by at least two units. This procedure was performed for the dataset as a whole, as well as for each combination of site and season that was shown to be significantly different in the mixed-effects models (see above) in soil and in leaf litter. For these analyses the measurements of FLNF that were below the detection limit were excluded to assess which predictor variables participate in regulating the FLNF that were detectable and thus participate to the ecosystem scale N fixation. The potential predictor variables for both leaf litter and soil were standardised to a mean = 0 and standard deviation = 1 prior to the model fitting procedure to avoid potential issues in interpretation and numerical stability due to differences in magnitude between variables. Potential predictor variables were gravimetric water content, C:N ratio, N:P ratio, total carbon, total nitrogen, total phosphorus and total molybdenum. For soil we additionally included available nitrogen, phosphorus and molybdenum, soil pH and soil texture. We present the best-fit model for each data subset, based on this forward stepwise procedure.

As an additional check of robustness we used an Akaike weight approach to assess the importance of predictor variables. We summed Akaike weights computed for all possible first-order models containing a given predictor, to obtain a measure of the relative variable importance for this predictor (Burnham and

Anderson 2002). We did this for all abovementioned predictors. This approach yielded very similar results to the forward stepwise model selection, confirming the robustness of our analysis, but for readability reasons are not presented nor discussed in this paper. These results are, however, shown in supplements (Figs. S2 and S3). All analyses were carried out with the software package R (version 3.4.3) using packages lme4 (Bates et al. 2015), MASS (Venables and Ripley 2002), ggfortify (Tang et al. 2016) and AICcmodavg (Mazerolle 2017).

## Results

Rates of FLNF were on average $0.015 \pm 0.003$ (standard error) nmol ethylene $g^{-1} h^{-1}$ or $0.57 \pm 0.10$ kg N $ha^{-1}$ $y^{-1}$ in soil and $0.25 \pm 0.04$ nmol ethylene $g^{-1} h^{-1}$ or $0.09 \pm 0.01$ kg N $ha^{-1}$ $y^{-1}$ in leaf litter. Per unit mass, FLNF rates were thus, on average, tenfold higher in leaf litter than in the top 5 cm of the soil (Table 1). However, when reported per unit area, FLNF rates were lower in litter than soil, due to the huge difference in density between the sample types (Table 1).

*Soil - N fixation* Overall, soil FLNF rates did not differ between seasons ($P = 0.27$). The effect of topography on soil FLNF rates differed between Paracou and Nouragues (Site x Topography interactive effect; $P = 0.021$) (Fig. 1a, b). In Paracou, soil FLNF rates were 20 times higher (+/-SE = 6 to 65 times higher) in the bottom and slope plots than in the top plots ($P = 0.047$ and $P = 0.009$, respectively), but did not differ significantly between the bottom and slope plots. In Nouragues no differences were observed in the FLNF rates of the bottom, slope or top plots. A considerable proportion of the soil FLNF rates were below the detection limit (25%, out of 230 samples). Although the LMER did not reveal a significant effect of season on the soil FLNF rates, season did affect the number of samples whose FLNF rates were below the detection limit (Fig. S4 A & B), albeit differently in Paracou than in Nouragues (Site x Season interactive effect $P = 0.004$). Whereas in Paracou the number of soil FLNF rates that were below the detection limit did not differ between wet and dry season, in Nouragues values below the detection limit occurred primarily in the dry season (18 out of 52 samples versus 2 out of 56 samples, in the dry and wet season, respectively, $P < 0.001$). When conducting the
<table><caption>**Table 1** Range, median, mean and standard error (SE) of FLNF rates at forest sites Paracou and Nouragues</caption><thead><tr><th rowspan="3"></th><th rowspan="3"></th><th rowspan="3"></th><th colspan="9">Paracou</th><th colspan="9">Nouragues</th></tr><tr><th colspan="3">Wet</th><th colspan="3">Dry</th><th colspan="3">Wet</th><th colspan="3">Dry</th></tr><tr><th>Bottom</th><th>Slope</th><th>Top</th><th>Bottom</th><th>Slope</th><th>Top</th><th>Bottom</th><th>Slope</th><th>Top</th><th>Bottom</th><th>Slope</th><th>Top</th></tr></thead><tbody><tr><td rowspan="5">Soil</td><td rowspan="4">Ethylene<br>production<br>(nmol g<sup>-1</sup> h<sup>-1</sup>)</td><td>Minimum</td><td>BDL*</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td></tr><tr><td>Median</td><td>0.006</td><td>0.005</td><td>BDL</td><td>0.015</td><td>0.008</td><td>0.001</td><td>0.002</td><td>0.003</td><td>0.004</td><td>0.014</td><td>0.006</td><td>0.012</td><td>0.012</td><td>0.012</td><td>0.012</td><td>0.012</td><td>0.012</td><td>0.012</td><td>0.012</td></tr><tr><td>Mean</td><td>0.012</td><td>0.017</td><td>0.002</td><td>0.014</td><td>0.013</td><td>0.003</td><td>0.003</td><td>0.005</td><td>0.028</td><td>0.055</td><td>0.010</td><td>0.037</td><td>0.037</td><td>0.037</td><td>0.037</td><td>0.037</td><td>0.037</td><td>0.037</td><td>0.037</td></tr><tr><td>Maximum</td><td>0.044</td><td>0.191</td><td>0.022</td><td>0.029</td><td>0.056</td><td>0.018</td><td>0.006</td><td>0.034</td><td>0.404</td><td>0.298</td><td>0.086</td><td>0.460</td><td>0.460</td><td>0.460</td><td>0.460</td><td>0.460</td><td>0.460</td><td>0.460</td><td>0.460</td></tr><tr><td>SE</td><td>0.003</td><td>0.009</td><td>0.001</td><td>0.002</td><td>0.002</td><td>0.001</td><td>0.000</td><td>0.002</td><td>0.020</td><td>0.030</td><td>0.004</td><td>0.021</td><td>0.021</td><td>0.021</td><td>0.021</td><td>0.021</td><td>0.021</td><td>0.021</td><td>0.021</td></tr><tr><td rowspan="5">N fixation rate<br>(kg N ha<sup>-1</sup> y<sup>-1</sup>)</td><td>Minimum</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td></tr><tr><td>Median</td><td>0.296</td><td>0.183</td><td>0.008</td><td>0.738</td><td>0.350</td><td>0.060</td><td>0.084</td><td>0.119</td><td>0.129</td><td>0.537</td><td>0.228</td><td>0.395</td><td>0.395</td><td>0.395</td><td>0.395</td><td>0.395</td><td>0.395</td><td>0.395</td><td>0.395</td></tr><tr><td>Mean</td><td>0.565</td><td>0.724</td><td>0.091</td><td>0.663</td><td>0.528</td><td>0.149</td><td>0.100</td><td>0.176</td><td>0.742</td><td>1.301</td><td>0.405</td><td>1.279</td><td>1.279</td><td>1.279</td><td>1.279</td><td>1.279</td><td>1.279</td><td>1.279</td><td>1.279</td></tr><tr><td>Maximum</td><td>2.082</td><td>7.836</td><td>1.215</td><td>1.449</td><td>2.168</td><td>0.858</td><td>0.181</td><td>1.157</td><td>8.944</td><td>10.513</td><td>3.635</td><td>15.850</td><td>15.850</td><td>15.850</td><td>15.850</td><td>15.850</td><td>15.850</td><td>15.850</td><td>15.850</td></tr><tr><td>SE</td><td>0.168</td><td>0.379</td><td>0.060</td><td>0.077</td><td>0.099</td><td>0.054</td><td>0.010</td><td>0.055</td><td>0.460</td><td>1.073</td><td>0.171</td><td>0.741</td><td>0.741</td><td>0.741</td><td>0.741</td><td>0.741</td><td>0.741</td><td>0.741</td><td>0.741</td></tr><tr><td rowspan="5">Litter</td><td rowspan="4">Ethylene<br>production<br>(nmol g<sup>-1</sup> h<sup>-1</sup>)</td><td>Sign.</td><td>a</td><td>a</td><td>b</td><td>a</td><td>a</td><td>b</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td><td>a</td></tr><tr><td>Minimum</td><td>0.026</td><td>0.025</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td></tr><tr><td>Median</td><td>0.186</td><td>0.143</td><td>0.112</td><td>0.111</td><td>0.053</td><td>0.000</td><td>0.109</td><td>0.059</td><td>0.030</td><td>0.087</td><td>0.045</td><td>0.043</td><td>0.043</td><td>0.043</td><td>0.043</td><td>0.043</td><td>0.043</td><td>0.043</td><td>0.043</td></tr><tr><td>Mean</td><td>0.524</td><td>0.433</td><td>0.403</td><td>0.376</td><td>0.290</td><td>0.001</td><td>0.346</td><td>0.300</td><td>0.150</td><td>0.207</td><td>0.046</td><td>0.139</td><td>0.139</td><td>0.139</td><td>0.139</td><td>0.139</td><td>0.139</td><td>0.139</td><td>0.139</td></tr><tr><td>Maximum</td><td>3.018</td><td>2.500</td><td>3.976</td><td>2.981</td><td>4.554</td><td>0.007</td><td>2.898</td><td>2.407</td><td>1.252</td><td>0.930</td><td>0.120</td><td>2.112</td><td>2.112</td><td>2.112</td><td>2.112</td><td>2.112</td><td>2.112</td><td>2.112</td><td>2.112</td><td>2.112</td></tr><tr><td rowspan="5">N fixation rate<br>(kg N ha<sup>-1</sup> y<sup>-1</sup>)</td><td>Minimum</td><td>0.011</td><td>0.005</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td><td>BDL</td></tr><tr><td>Median</td><td>0.052</td><td>0.045</td><td>0.044</td><td>0.028</td><td>0.015</td><td>0.000</td><td>0.042</td><td>0.017</td><td>0.008</td><td>0.043</td><td>0.015</td><td>0.019</td><td>0.019</td><td>0.019</td><td>0.019</td><td>0.019</td><td>0.019</td><td>0.019</td><td>0.019</td><td>0.019</td></tr><tr><td>Mean</td><td>0.166</td><td>0.127</td><td>0.142</td><td>0.109</td><td>0.102</td><td>0.000</td><td>0.126</td><td>0.112</td><td>0.045</td><td>0.111</td><td>0.017</td><td>0.065</td><td>0.065</td><td>0.065</td><td>0.065</td><td>0.065</td><td>0.065</td><td>0.065</td><td>0.065</td><td>0.065</td></tr><tr><td>Maximum</td><td>0.889</td><td>0.580</td><td>1.332</td><td>0.898</td><td>1.707</td><td>0.002</td><td>0.615</td><td>0.979</td><td>0.418</td><td>0.800</td><td>0.053</td><td>0.867</td><td>0.867</td><td>0.867</td><td>0.867</td><td>0.867</td><td>0.867</td><td>0.867</td><td>0.867</td><td>0.867</td></tr><tr><td>SE</td><td>0.054</td><td>0.034</td><td>0.062</td><td>0.046</td><td>0.083</td><td>0.000</td><td>0.044</td><td>0.048</td><td>0.021</td><td>0.041</td><td>0.002</td><td>0.042</td><td>0.042</td><td>0.042</td><td>0.042</td><td>0.042</td><td>0.042</td><td>0.042</td><td>0.042</td><td>0.042</td></tr><tr><td></td><td></td><td>Sign.</td><td>a</td><td>a</td><td>ab</td><td>b</td><td>b</td><td>c</td><td>a</td><td>abd</td><td>ac</td><td>bcd</td><td>ce</td><td>be</td><td>be</td><td>be</td><td>be</td><td>be</td><td>be</td><td>be</td><td>be</td></tr></tbody></table>

Rates measured in sample types soil and leaf litter are given for the wet and dry season separately and split up by landscape unit. Rates are expressed both on a mass basis as ethylene production rates and on an area basis as N fixation rates. Significant differences (*P* < 0.05) within a single site and sample type are indicated by differing letters. Calculation of median, mean and standard error and statistical analyses were performed from data including BDL* measurements. *BDL, below detection limit of 0.01 nl ethylene produced per ml air. This equates to a detection limit of 1*10⁻⁴ nmol g⁻¹ h⁻¹ and 7*10⁻⁴ kg N ha⁻¹

LMER analysis on a subset of the Nouragues soil FLNF rates, excluding samples where FLNF rates were below the detection limit, we identified an effect of season (*P* = 0.008), with the highest rates occurring in the dry season. (Fig. S5 B).

**Soil – Environmental variables** As each site had different soil FLNF rates relating to their topography, we performed a principal component analysis (PCA) for each site separately, allowing us to visualize site specific relationships among soil parameters. In Paracou, PC1 and PC2 explained 63.3% of the variation, with PC1 explaining 36.3% of the variation and containing most of the topographically induced variation (Fig. 2a). Clay content

correlated with total P and available N and was highest on the slopes. PC2, explaining 27% of the variation, contained the seasonally-induced variation (Fig. 2a). Total C correlated with total N and was higher in the dry season, while total Mo correlated well with pH and was higher in the wet season. Available Mo correlated well with moisture and was also higher in the wet season. Available P was not heavily affected by season, but was higher in the bottom plots than on either the slopes or tops. In Nouragues PC1 and PC2 together explained 68.6% of the variation (Fig. 2b). The first principal component (PC1) explained 55.2% of the variation and, just as in Paracou, mostly contained topographically-induced variation. In Nouragues soil properties for top plots were clearly
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Area-based N fixation rates for soil at Paracou (a) and Nouragues (b) and for leaf litter at Paracou (c) and Nouragues (d). Rates represent the means (±1 SE) of the samples collected from the bottom, slope and top landscape positions in the wet (May) and dry (September) seasons (N = 20). Letters denote significant</figcaption></figure>

different from bottom and slope plots. Clay content corre-
lated with total P and available N and all were highest in
the top plots. In contrast to Paracou, however, total Mo
correlated with available P and both were higher on the
bottom and slope plots than on the top plots. Moisture,
total C and total N correlated well with each other and
clustered together at a ±45° angle from PC1. Available
Mo and C:N were correlated along PC2, which explained
13.4% of the variation. Lastly, pH also varied topograph-
ically and was highest on the bottom and slope plots, yet
the variation was small. For all Pearson correlations, see
Table S3 A and B.

*Leaf litter – N fixation* Leaf litter N fixation rates were not different between sites. On average, litter FLNF rates in the wet season were 7.5 to 13.5 times (± SE) higher (*P* < 0.001) than in the dry season. This

differences amongst different seasons and landscape units within
a single site for soil or leaf litter at the P < 0.05 level. Significance
testing was performed by mixed effects regression models using
the log transformed values as measured values were non-normal.
Data include measurements that were below the detection limit

difference between both seasons was further empha-
sized by the distinct seasonal difference in the propor-
tion of FLNF measurements that were below the detec-
tion limit (Fig. S4 C & D): the vast majority (97%) of the
16% (33 out of 235) of litter FLNF measurements that
were below the detection limit were from the dry season
(binomial GLM, P < 0.001). Nonetheless, the impact of
season interacted with that of topographical position
(Season x Topography interactive effect P < 0.001)
(Fig. 1c, d). In the wet season FLNF rates were similar
among topographic positions, but in the dry season,
rates on the top plots were 2 to 7.5 times (± SE) lower
(P < 0.05) than rates on the bottom or slope plots, which
did not differ from each other.

*Leaf litter – Environmental variables* Because the LMER stated that FLNF rates were similar between sites
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Principal Component Analysis of soil variables in Paracou (a) and Nouragues (b). Points on the graph are plot averages. C = total C, N = total N, P = total P, Mo = total Mo, C:N = C:N ratio, P = available P, N in = available N, Mo in = available Mo, Moisture = water content, pH = pH, Clay = % clay, Sand = % sand and BD = bulk density</figcaption></figure>

we analysed the leaf litter stoichiometry of both sites together. Together, the first two principal components of the PCA accounted for 73.3% of the variation (Fig. 3). PC1 explained 49.8% of the variation and distinguished between the wet and dry season. Along this component, moisture, N and P content all correlated positively with each other and negatively with C:N and N:P ratio, indicating that moisture, N and P content was higher in the wet season, while C:N and N:P ratios were higher in the dry season. This was confirmed by linear effects regression analysis (Table S2). PC2, accounting for 23.5% of the variation, correlated positively with C content. The absence of a site effect on PC1 and PC2 indicates that leaf litter chemistry was similar for both sites. For all Pearson correlations, see Table S3 C.

*Drivers of FLNF in soils and litter* In soils, stepwise regression analyses indicated that soil water content, available P and pH were the primary drivers of N fixation rates (Table 2a), explaining 24% of the variation in soil FLNF rates. Because soil FLNF rates differed significantly between both sites, we ran the stepwise regression analysis again for each site separately. In Paracou, P exerted a strong effect, as both total P and available P had a positive influence on FLNF rates. Available N was negatively

related to FLNF, and together with P, explained 36% of the variation. For Nouragues a model containing only total N, total C or water content explained 40% of the variation in FLNF rates. These three parameters were strongly correlated and the regression analysis deemed the models containing either one of them singularly equally valid (Table 2a).

In litter, across both sites and seasons, about 46% of the variation in FLNF rates was explained by water and N:P ratio (Table 2b). While water had a positive effect on FLNF rates, N:P ratio showed a negative effect. As FLNF rates differed between seasons we ran the stepwise analysis again for each season separately. In the wet season 33% of the variation could be explained by water, C content, P content and C:N ratio. The model for the dry season explained about 68% of the variation and was dependent on the positive influence of water and the negative influence of N content.

## Discussion

Overall, the sum of FLNF rates of soils and leaf litter measured in this study fall into the lower end of the 0.1–
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Principal Component Analysis of leaf litter variables in both Paracou and Nouragues and for both wet and dry season. Dots on the graph are averaged for plot. C = total C, N = total N,</figcaption></figure>

P = total P, Mo = total Mo, C:N = C:N ratio, N:P = N:P ratio and Moisture = water content

<table><caption>**Table 2** Results of the forward stepwise selection analyses for N fixation rates in soil (A) and leaf litter (B). For these analyses we excluded the measurements that were below the detection limit</caption><thead><tr><th>Site</th><th>Season</th><th>Parameters</th><th></th><th></th><th>R<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="2">Par & Nou</td><td rowspan="2">Wet & dry</td><td>+0.57 Water*</td><td>+0.20 Available P*</td><td>+0.17 pH*</td><td rowspan="2">0.240</td></tr><tr><td>F<sub>1,169</sub> = 15.6</td><td>F<sub>1,179</sub> = 4.8</td><td>F<sub>1,180</sub> = 4.5</td></tr><tr><td rowspan="2">Par</td><td rowspan="2">Wet & dry</td><td>+1.1 Total P***</td><td>+0.34 Available P**</td><td>-0.51 Available N**</td><td rowspan="2">0.362</td></tr><tr><td>F<sub>1,88</sub> = 44.7</td><td>F<sub>1,91</sub> = 7.2</td><td>F<sub>1,68</sub> = 9.3</td></tr><tr><td rowspan="3">Nou</td><td rowspan="3">Wet & dry</td><td>+0.25 Total N*</td><td></td><td></td><td rowspan="3">0.398</td></tr><tr><td>F<sub>1,88</sub> = 5.0</td><td></td><td></td><td></td></tr><tr><td>+0.25 Total C*</td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>F<sub>1,88</sub> = 4.7</td><td></td><td></td><td rowspan="2">0.396</td></tr><tr><td></td><td></td><td>+0.25 Water*</td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>F<sub>1,88</sub> = 5.0</td><td></td><td></td><td>0.401</td></tr></tbody></table><table><thead><tr><th>Site</th><th>Season</th><th>Parameters</th><th></th><th></th><th>R<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="2">Par & Nou</td><td rowspan="2">Wet & Dry</td><td>+0.76 Water***</td><td>-0.32 N:P**</td><td></td><td rowspan="2">0.457</td></tr><tr><td>F<sub>1,40</sub> = 24.67</td><td>F<sub>1,195</sub> = 8.34</td><td></td></tr><tr><td rowspan="2">Par & Nou</td><td rowspan="2">Wet</td><td>+1.71 Water**</td><td>+0.65 Total P**</td><td>+1.38 Total C**</td><td>+ 0.81 CN*</td><td rowspan="2">0.328</td></tr><tr><td>F<sub>1,65</sub> = 9.10</td><td>F<sub>1,117</sub> = 10.69</td><td>F<sub>1,117</sub> = 8.64</td><td>F<sub>1,115</sub> = 5.87</td></tr><tr><td rowspan="2">Par & Nou</td><td rowspan="2">Dry</td><td>+2.63 Water***</td><td>-0.79 Total N*</td><td></td><td></td><td rowspan="2">0.676</td></tr><tr><td>F<sub>1,82</sub> = 26.26</td><td>F<sub>1,73</sub> = 6.70</td><td></td><td></td><td></td></tr></tbody></table>

Each row lists the parameters that were included into the model that provided the best fit, based on AICc values. Columns one and two give information regarding the dataset that was used. For soils (A) a model including data from both sites (Par, Paracou; Nou, Nouragues) and seasons was made along with models for each site separately because N fixation rates were different between sites. For leaf litter (B) a model including data from both sites and seasons was made along with models for each season separately because seasonality affected N fixation. Parameters that had a significant effect (* = p < 0.05, ** = p < 0.01 and *** = p < 0.001) on N fixation rates are given with their F values. Plus or minus numbers represent beta values. R² was calculated for the fit of the modelled and measured fixation rates

<footnote n="1">¹ The stepwise selection identified total N, total C and water content as equally important for determining FLNF rates in Nouragues soils. This is because these three variables are highly correlated and behave similarly in the model. Explanations as to why this occurs are offered in the text</footnote>
60 kg N ha⁻¹ y⁻¹ range reported for FLNF on tropical forest floors worldwide (Reed et al. 2011), and much below the more recent estimate of 15–36 kg N ha⁻¹ y⁻¹ fixed in tropical forests (Pajares and Bohannan 2016) and the average ethylene production rate of 5.32 nmol ethylene g⁻¹ h⁻¹ reported in a recent meta-analysis on N fixation rates in tropical forest ecosystems (Dynarski and Houlton 2018). The FLNF rates found in the present study are much lower than those found in Costa Rica (Reed et al. 2007; Reed et al. 2010; Reed et al. 2013), Panama (Barron et al. 2009; Wurzburger et al. 2012), Puerto Rico (Cusack et al. 2009) or in the younger sites along a Hawaiian chronosequence (Crews et al. 2000). An overview of FLNF rates reported by these authors can be found in Table S4. Sullivan et al. (2014) measured N fixation rates in Costa Rica, about 35 km away from where Reed et al. (2007) conducted their study. These authors reported soil and leaf litter FLNF rates that were lower than those previously published by Reed et al. (2007), but while their reported soil FLNF rates were similar to those found in our study, their reported leaf litter FLNF rates were still twice as high as the rates we found. Other studies of FLNF rates conducted in less fertile tropical forests, e.g. on the higher altitudes of an altitudinal transect in Ecuador (Matson et al. 2014) or in the older sites of a Hawaiian chronosequence (Crews et al. 2000), reported FLNF rates comparable to those found in the present study. In Hawaii, Crews et al. (2000) measured decreasing rates of leaf litter FLNF on increasingly older soils and attributed this decrease to a combination of lower concentrations of geologically cycled nutrients, such as P, and high N pools at their oldest sites. This, coupled with increases in FLNF rates after P fertilization, led them to one of their main conclusions; that low P availability limited FLNF rates at their older sites. It would, however, be inaccurate to conclude that soil P alone determines N fixation since, for example, in Panama (Wurzburger et al. 2012) two of the studied sites along a total soil P gradient containing high (AVA) and medium (Gigante) levels of total soil P displayed the lowest rates of FLNF in their study.

## Soil free-living N fixation

At both sites moisture was important in the regulation of FLNF (Table 2a). Water is essential for all microbes, but for diazotrophs in particular water plays an important role in regulating activity. Not only does their

metabolism require sufficient amounts of water, but nitrogenase activity is hindered in the presence of oxygen (O₂) (Hill 1988). Oxygen binds to nitrogenase and inhibits its function (Hartmann and Burris 1987) and because increased soil moisture decreases soil O₂ concentrations, water content can have a large impact on soil FLNF rates. As diazotrophs are mainly heterotrophic the thickness of the soil water film, which is important for the diffusion rates of extracellular enzymes and soluble organic-C substrates and is directly affected by soil water content (Davidson and Janssens 2006), will also play a role.

In Paracou, we observed no seasonal effect on soil FLNF rates and although rates were typically reported to be higher in the wet season than in the dry season (e. g. Hofmockel and Schlesinger 2007; Reed et al. 2007), Matson et al. (2014) have shown that the opposite can also occur. They postulated that it is likely that moisture fluctuations were not directly responsible for their observed seasonal changes in FLNF rates, just as in our study it might not have contributed to unchanging rates. Their reasoning was based on the fact that cyanobacteria can fix N at moisture levels as low as 6% and in one study reached maximum N fixation rates at soil moistures between 22 and 42% (Jones 1977), though it did not specify whether this was in sandy or clayey soil. As soil moisture in Paracou always remained between 15 and 44% (Table S1) it is likely that diazotrophs had sufficient moisture to keep oxygen levels low and continue N fixation.

At Nouragues, the effect of soil moisture on FLNF rates was high (40% of variability explained by soil moisture alone, Table 2a) and seasonality caused an interesting pattern in the distribution of very low FLNF rate measurements (below detection limit). The vast majority of below detection limit measurements occurred during the dry season (Fig. S4 B) and when analysing only data for which FLNF rates were above the detection limit, we found a significant season effect in Nouragues with quantifiable rates (i.e. in samples above the detection limit) actually being higher in the dry season than in the wet season (Fig. S5 B). This meant that the range of FLNF rates was broader during the dry season than during the wet season (Fig. S5 A), indicating that N fixation hotspots were of increasing importance. These are typically found in tropical ecosystems (Pajares and Bohannan 2016) and reflect the very small-scale spatial heterogeneity of abiotic factors affecting the dynamics of the diazotroph community (Reed et al. 2010). As the soil dried out at the onset of
the dry season the heterogeneity of water-filled pore space increased, leading to the creation of aerobic and anaerobic, as well as dry and mesic microsites that co-occured on a small spatial scale (Sexstone et al. 1985).

Besides water content, available P and pH were the most important predictors for soil FLNF at our sites. The pH in our soils is very low, ranging between 3.8 and 4.2 (Table S1). This is a relatively small range, yet pH has been identified to affect bacterial community composition and diversity (Rousk et al. 2010; Tripathi et al. 2014) or soil CO₂ effluxes (Courtois et al. 2018) at the local scale even when changes are small. We were unable to find studies assessing the effect of pH on FLNF in tropical ecosystems, but nitrogenase activity of soil diazotrophs has been shown to increase with increased pH in a German pine forest situated in north-east Bavaria (Limmer and Drake 1996). Roper and Smith (1991) found that nitrogenase activity of bacteria extracted from clayey Australian soils reached its peak around pH 7.5, matching the nitrogenase pH optimum found by Pham and Burgess (1993), and decreased in more acidic soils. Their study, however, was carried out on soils that at the start of the experiment were only slightly acidic after which the pH was decreased during the course of the study. Their result is therefore not necessarily representative for N fixing microbial communities that developed in acidic soils, such as in our study. However, taking into account an enzymatic pH optimum of 7.5 (Pham and Burgess 1993) and decreased nitrogenase activity in more acidic temperate soils (Zahran 1999), a positive relationship between soil pH and FLNF rates in our acidic tropical forest sites is plausible.

Phosphorus can play a pivotal role in regulating N fixation rates and is often limiting the rate of this process in highly weathered lowland tropical soils (e. g. Camenzind et al. 2018; Reed et al. 2011). In Paracou, both P predictors were positively correlated with the FLNF rates, suggesting a higher activity of N fixing microbes when more P is present in the soil. Additions of P to the forest floor have shown to increase both diversity and abundance of diazotrophs in tropical soils (Reed et al. 2010), as well as increases in FLNF rates (e. g. Benner et al. 2007; Reed et al. 2013). The positive correlation between available P and total P on the one hand, and FLNF rates on the other, is in line with the longstanding idea that increased P benefits FLNF rates in tropical lowland forests. In addition to P, the model for Paracou also identified available N as a predictor, indicating that higher available N is linked to decreased

FLNF rates. This is in line with what was observed in previous tropical forest studies, where additions of N led to the decrease of N fixation rates (e. g. Barron et al. 2009; Crews et al. 2000; Cusack et al. 2009). This observation is generally attributed to the theory that many heterotrophic N₂ fixers are facultative fixers and able to down-regulate their fixation pathway when other sources of N are available and more cost-efficient to acquire (Menge et al. 2009). The topographic effect on FLNF rates in Paracou, effectively resulting in lower rates on the top than on bottom and slope plots is likely caused by the interplay of P and N in the Paracou soils as they also varied with topographic position. Total P was highest on the slope and available P highest on the bottom, while both were lowest on the top. Because of the higher clay content of the Paracou slope soils there are likely more aluminium and iron oxides (Tiessen et al. 1994) that are able to occlude more P, resulting in higher soil total P contents as seen on the slope. On the more sandy bottom landscape position metal oxides that can occlude P are likely scarcer and the higher water content, especially during rain events when runoff causes disproportionate changes in water status, leads to reduction of iron oxides (Colombo et al. 2014), liberating occluded P and provoking higher concentrations of available P.

In Nouragues, we found that water content, C and N performed equally well, each individually explaining 40% of the observed variation in FLNF rates, suggesting that the absolute concentrations of C and N in the soil explain a substantial part of the variation in FLNF rates. It is surprising that the model identified a positive effect of N content on FLNF rates, given that most studies associate increased N with decreased or unchanged fixation rates (Camenzind et al. 2018). Additionally, we would expect that available N, which was also included as potential variable but was not selected by the model, would be identified as a variable affecting FLNF rates instead of N content because no decomposition is needed before assimilation. Likely, the identification of N by the model is purely mathematical and caused by its near perfect correlation with C content (Table S3 B). The positive relationship of C with FLNF rates, which is predominantly carried out by heterotrophic diazotrophs (Sprent and Sprent 1990), is the result of decomposition of organic material and subsequent assimilation of additional C providing the energy necessary to carry out fixation. In Nouragues, the high correlation between C content and moisture (Table S3
B) partly explains their equal importance in the model. Their correlation might be due to the influence of soil moisture on decomposition rates, especially when it is very wet (negative relation) or very dry (positive relation). At high moisture levels, soil organic carbon and nitrogen stocks increase because of the slower decomposition in water-saturated soil (Van Sundert et al. 2018). Moreover, soils in Nouragues are clay-rich, exacerbating this effect as soils containing more clay stabilize and store more soil organic matter than sandy soils (Reis et al. 2014), such as those in Paracou. Both increased water and C can be beneficial for free-living diazotrophs as many species possess heterotrophic anaerobic metabolisms (Dixon and Kahn 2004) and proliferate under oxygen poor and carbon rich conditions.

In Nouragues the topographical patterns of N and P were different from Paracou, with much higher total P concentrations that occurred on the top landscape position instead of on the slopes, and much smaller differences in available P concentrations among the topographic positions. This likely played a role in the regulation of soil FLNF rate, causing them to remain high on the top landscape position. Because of this strong difference in topography effects between Nouragues and Paracou, we cannot draw general conclusions about landscape-scale variations in soil FLNF rates across topographies in French Guianese tropical forests. Instead, our results support the idea that soil FLNF rates at our sites varied in function of water and nutrient availability, similar to what was reviewed by Dynarski et al. (2018).

Lastly, it is interesting to note that total or available Mo was never selected as explanatory variable for soil FLNF rates, in spite of playing a regulatory role in FLNF at other sites (Barron et al. 2009; Reed et al. 2013; Wurzburger et al. 2012). Molybdenum concentrations in our study sites were slightly higher than those found elsewhere (Gupta and Lipsett 1981), and available Mo was ten times higher than in Panama (Wurzburger et al. 2012), which was the only study reporting the effect of available molybdenum on FLNF. Likely, the concentration of available Mo at our sites is high enough to preclude any regulatory role of Mo in this P-limited environment.

## Leaf litter free-living N₂ fixation

Overall, leaf litter FLNF rates did not vary with topography and were instead driven mainly by water content and the N:P ratio of the leaf litter, as shown by the overall

model explaining 46% of the variation (Table 2b). The negative influence of N:P ratio is in line with what was postulated by Reed et al. (2007) and likely illustrates a stoichiometric and energetic balance shift. In a decomposing leaf a decreasing N:P ratio leads to P becoming comparatively more abundant, shifting the stoichiometric balance to move away from P limitation to N limitation, favouring diazotrophs. Additionally, increased P in the environment may alleviate energetic constraints (see above) and enable diazotrophs to invest the required energy into fixing N. As the general model included data from both seasons, a large influence of water content, as evidenced by the largest beta value (Table 2b), was to be expected because litter FLNF rates were significantly affected by seasonality; wet season FLNF rates were up to nearly two times higher than dry season FLNF rates (Table 1) and the number of samples for which FLNF rates were below the detection limit was much smaller in the wet than in the dry season (Fig. S4 C & D). As expected, litter dried out severely during the dry season and water content decreased from an average ~64 to ~40% across both sites. Given that diazotrophs grow and proliferate in aqueous environments (Scott 1957), it is likely that leaf litter water content in the dry season dropped below a threshold of minimum water required for diazotroph proliferation, resulting in a net drop of FLNF rate. This response of diazotroph activity to leaf litter moisture was already observed in a sub-tropical karst forest, where researchers found that decreases in leaf litter moisture resulted in decreased FLNF rates (Li et al. 2018). The large beta value (2.63, Table 2b) associated with water content in the dry season model hints towards the disruptive effect of water shortage. This model explained 67% of the variation in litter FLNF rates and included only water content and litter N, which, just as the N:P ratio in the general model, had a negative effect on the FLNF rates. As mentioned earlier, N assimilation is cheaper than N fixation from an energetic point of view and when more N is present fixation will likely be down regulated (Menge et al. 2009). On the top landscape position in Paracou, also during the dry season, leaf litter water content was lowest while N:P ratio was highest. This likely led to a complete collapse of N fixation (Fig. 1c), because N was in ample supply compared to P and water was scarce. During the wet season, litter stoichiometry explained 33% of the variation in FLNF rates which likely reflects, at least partially, the interactions between seasonal changes in labile C availability, P content and N demand (Reed et al. 2007). Litterfall peaks at the onset of the dry season
(Chave et al. 2010; Wagner et al. 2013) and once the wet season starts the daily rainfall provides a vehicle for the movement of readily decomposable, dissolved organic C (DOC) within the litter layer (Courtois et al. 2018). The N₂-fixing microbial community in the litter layer is dominated by heterotrophic microorganisms (Sprent and Sprent 1990) and the influx of litter DOC provides an easily accessible energy source for these diazotrophs. Additionally, in our study, we found that the concentration of both P and N in leaf litter was higher in the wet season than in the dry season, similar to what was found in another study in Costa Rica (Wood et al. 2005). However, relative to N content, P content increased more towards the wet season, resulting in a lower N:P ratio during the wet season. In combination with more labile C input, this could stimulate N fixation through the relief of energetic constraints and the added advantage of being able to fix N compared to assimilating N from an environment where it is, comparatively, less abundant than in the dry season. This process was already observed in tallgrass prairie soils (Eisele et al. 1989).

In contrast to litter P content, which changed significantly across seasons yet showed no significant change along the topographical gradient, litter Mo content changed significantly between the two seasons and along the topographical gradient. Just as in soils, however, and in spite of its importance for N fixation (Kaiser et al. 2005; Seefeldt et al. 2009), we found no evidence to support a regulatory role for Mo content. Similar as in soils, litter Mo concentrations were ten-fold higher than those reported from other tropical sites (Barron et al. 2009; Bowell and Ansah 1993; Reed et al. 2013; Wurzburger et al. 2012) and likely too high to render a regulatory role to Mo in litter FLNF rates.

Lastly, in both soils and leaf litter our models were unable to explain more than 67% of the observed variation, and in most cases only around 40%. This means that often the majority of variation in FLNF rates could not be explained by the variables we selected and thus it is very likely that other factors not measured in this study contribute to the regulation of FLNF in our tropical lowland forests. Knowing the diazotroph community composition, in both soil and leaf litter, could enhance our understanding of their nutritional and environmental needs and help us estimate at what point parameters such as moisture, pH and P or N availability are beneficial or, inversely, detrimental. The likelihood that within a single diazotroph community both aerobic and anaerobic lifestyles can occur (Dixon and Kahn

2004) and that N fixing Archaea within the community may possess different nutritional requirements than bacterial diazotrophs (Leigh 2000) are additional reasons to study the diazotroph community composition. Also, we did not assess the iron (Fe) or the vanadium (V) availabilities at our sites while they have the potential of participating in FLNF rate regulation (Zhang et al. 2016). While Mo is a necessary co-factor of most nitrogenase enzymes (Igarashi and Seefeldt 2003), Fe is found in all known nitrogenases and the occurrence of ‘iron-only’ (Fe-Fe) nitrogenases has been widely documented (e. g. Yang et al. 2014; Zheng et al. 2018). The role of V in regulation of FLNF is understudied, but it is certain that the occurrence of an alternative enzymatic co-factor, namely the vanadium-iron (V-Fe) cofactor, is widespread and V availability may therefore also play a regulatory role in nitrogenase biosynthesis (Hu et al. 2012). A follow-up study investigating the community composition of the diazotrophs, specifically looking into the prevalence of aerobic and anaerobic organisms, combined with research into the occurrence of V-Fe nitrogenases and V and Fe availabilities in the soil and leaf litter could enable us to explain more of the observed variation than we presently could.

## Conclusion

Our study has shed light on the drivers behind FLNF in tropical soil and leaf litter on the Guiana Shield and has shown that the rates of FLNF are much lower than those estimated for most tropical forests elsewhere. Water, N and P played the main roles in determining FLNF rates in both sample types, while pH only regulated in the soil. The effect of seasonality differed between sample type and differences in FLNF regulation between sites could be observed in soils, but not in litter. Despite having been shown to influence N fixation rates in other mature tropical forests, the micro nutrient Mo played no role in the regulation of FLNF rates at our sites in French Guiana.

In the sandy site, Paracou, the stimulating effect of P and the inhibiting effect of N were the main drivers behind soil FLNF, but in the clayey site, Nouragues, soil FLNF was mainly stimulated by water or C content. Our models for soil FLNF could not explain more than 40% of the observed variation, illustrating the complexity of predicting
fixation rates upon measured variables in an envi-
ronment that is highly heterogeneous on a regional,
local and even micro scale. In leaf litter we also
identified water, N and P as main drivers, but the
underlying mechanisms that caused variation may
have been different compared to the soil. In the leaf
litter we observed no differences in FLNF rates
between sites, but during the dry season litter rates
exhibited a drastic decline that was mainly related
to water insufficiency and the inhibiting effect of
N. During the wet season water was still of impor-
tance, but now stimulating effects of C and P also
came into play. It is important to note that our litter
wet season model explained only about 30% of the
observed variation.

It is likely that in both soil and leaf litter
diazotroph community composition and iron or
vanadium availability had an influence on FLNF
rates making them interesting to measure in future
studies. In the larger framework of global change,
where N deposition is expected to increase
(Penuelas et al. 2013), P deposition to tropical
forests may change (Gross et al. 2016) and the
possibility of a drier Amazon basin (IPCC 2013)
may cause disruptive changes to the FLNF rates in
soil and litter. Nutrient addition studies may offer
clues as to the response of FLNF to changes in N
or P supply, but testing the influence of climatic
changes in situ calls for a very specific type of
studies.

**Acknowledgements** This research was supported by the Euro-
pean Research Council Synergy grant ERC-2013-SyG.

610028-IMBALANCE-P. We thank the staff of the Nouragues station, managed by USR mixte LEEISA.

(CNRS; Cayenne), and the Paracou station, managed by UMR
Ecofog (CIRAD, INRA; Kourou).

Both research stations received support from “Investissement d’Avenir” grants managed by Agence Nationale de la Recherche (CEBA: ANR-10-LABX-25-01, ANAEE France: ANR-11-INBS-0001). We thank Tom Van der Spiet, Anne Cools and Miguel Portillo Estrada for their help with analyses and Philippe Gaucher, Nadine Calluy, Lore Verryckt and Kevin Van Sundert for their help with logistics and statistics.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

References

Aguilos M, Hérault B, Burban B, Wagner F, Bonal D (2018) What drives long-term variations in carbon flux and balance in a tropical rainforest in French Guiana? Agric For Meteorol 253-254:114-123. https://doi.org/10.1016/j.agrformet.2018.02.009

Amundson R, Austin AT, Schuur EAG, Yoo K, Matzek V, Kendall C, Uebersax A, Brenner D, Baisden WT (2003) Global patterns of the isotopic composition of soil and plant nitrogen. Global Biogeochem Cy 17. https://doi.org/10.1029/2002gb001903

Barron AR, Wurzburger N, Bellenger JP, Wright SJ, Kraepiel AML, Hedin LO (2009) Molybdenum limitation of asymptotic nitrogen fixation in tropical forest soils. Nat Geosci 2:42–45. https://doi.org/10.1038/Ngeo366

Barron AR, Purves DW, Hedin LO (2011) Facultative nitrogen fixation by canopy legumes in a lowland tropical forest. Oecologia 165:511–520. https://doi.org/10.1007/s00442-010-1838-3

Bates D, Mächler M, Bolker B, Walker S (2015) Fitting linear mixed-effects models using lme4. 67:48. https://doi.org/10.18637/jss.v067.i01

Batterman SA, Hedin LO, van Breugel M, Ransijn J, Craven DJ, Hall JS (2013) Key role of symbiotic dinitrogen fixation in tropical forest secondary succession. Nature 502:224–227. https://doi.org/10.1038/nature12525

Benner JW, Conroy S, Lunch CK, Toyoda N, Vitousek PM (2007) Phosphorus fertilization increases the abundance and Nitrogenase activity of the Cyanolichen Pseudocryptillaria crocata in Hawaiian montane forests. Biotropica 39:400–405. https://doi.org/10.1111/j.1744-7429.2007.00267.x

Bongers F, Charles-Dominique P, Forget P, Thèry M (2001) Nouragues dynamics and plant-animal interactions in a Neotropical rainforest. Kluwer Academic Publishers, Dordrecht, The Netherlands https://doi.org/10.1007/978-0-015-9821-7

Bowell RJ, Ansah RK (1993) Trace element budget in an African savannah ecosystem. Biogeochemistry 20:103–126. https://doi.org/10.1007/bf00004137

Bray RH, Kurtz LT (1945) Determination of Total, organic, and available forms of phosphorus in soils. Soil Sci 59:39–46

Brookshire ENJ, Gerber S, Menge DNL, Hedin LO (2012) Large losses of inorganic nitrogen from tropical rainforests suggest a lack of nitrogen limitation. Ecology letters 15:9–16. https://doi.org/10.1111/j.1461-0248.2011.01701.x

Burnham K, Anderson D (2002) Model selection and multimodel inference: a practical information-theoretic approach. Springer-Verlag, New York

Černá B, Rejmánková E, Snyder JM, Šantrůčková H (2009) Heterotrophic nitrogen fixation in oligotrophic tropical marshes: changes after phosphorus addition. Hydrobiologia 627:55–65. https://doi.org/10.1007/s10750-009-9715-y

Chave J, Navarrete D, Almeida S, Álvarez E, Aragão LEOC, Bonal D, Châtelet P, Silva-Espejo JE, Goret JY, von Hildebrand P, Jiménez E, Patiño S, Peñuela MC, Phillips OL, Stevenson P, Malhi Y (2010) Regional and seasonal patterns of litterfall in tropical South America. Biogeosciences 7:43–55. https://doi.org/10.5194/bg-7-43-2010
Colombo C, Palumbo G, He J-Z, Pinton R, Cesco S (2014) Review on iron availability in soil: interaction of Fe minerals, plants, and microbes. J Soils Sediments 14:538–548. https://doi.org/10.1007/s11368-013-0814-z

Courtois EA, Stahl C, van den Berge J, Bréchet L, van Langenhove L, Richter A, Urbina I, Soong JL, Peñuelas J, Janssens IA (2018) Spatial variation of soil CO2, CH4 and N2O fluxes across topographical positions in tropical forests of the Guiana shield. Ecosystems. 21:1445–1458. https://doi.org/10.1007/s10021-018-0232-6

Crews TE, Farrington H, Vitousek PM (2000) Changes in Asymbiotic, heterotrophic nitrogen fixation on leaf litter of Metrosideros polymorpha with long-term ecosystem development in Hawaii. Ecosystems 3:386–395. https://doi.org/10.1007/s100210000034

Cusack DF, Silver W, McDowell WH (2009) Biological nitrogen fixation in two tropical forests: ecosystem-level patterns and effects of nitrogen fertilization. Ecosystems 12:1299–1315. https://doi.org/10.1007/s10021-009-9290-0

Davidson EA, Janssens IA (2006) Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440:165–173. https://doi.org/10.1038/nature04514

Dixon R, Kahn D (2004) Genetic regulation of biological nitrogen fixation. Nat Rev Microbiol 2:621–631. https://doi.org/10.1038/nrmicro954

Dynarski KA, Houlton BZ (2018) Nutrient limitation of terrestrial free-living nitrogen fixation. The New phytologist 217: 1050–1061. https://doi.org/10.1111/nph.14905

Eisele KA, Schimel DS, Kapustka LA, Parton WJ (1989) Effects of available P and N:P ratios on non-symbiotic dinitrogen fixation in tallgrass prairie soils. Oecologia 79:471–474. https://doi.org/10.1007/bf00378663

Epron D, Bosc A, Bonal D, Freycon V (2006) Spatial variation of soil respiration across a topographic gradient in a tropical rain forest in French Guiana. J Trop Ecol 22:565–574. https://doi.org/10.1017/s0266467406003415

FAO (1998) World Reference Base for Soil Resources, by ISSS-ISRIC-FAO. World Soil Resources Report No. 84., Rome, p 88

Furnkranz M, Wanek W, Richter A, Abell G, Rasche F, Sessitsch A (2008) Nitrogen fixation by phyllosphere bacteria associated with higher plants and their colonizing epiphytes of a tropical lowland rainforest of Costa Rica. The ISME journal 2:561–570. https://doi.org/10.1038/ismej.2008.14

Galloway JN, Schlesinger WH, Levy H, Michaels A, Schnoor JL (1995) Nitrogen fixation: Anthropogenic enhancement-environmental response. Global Biogeochem Cy 9:235–252. https://doi.org/10.1029/95GB00158

Galloway JN, Dentener FJ, Capone DG, Boyer EW, Howarth RW, Seitzinger SP, Asner GP, Cleveland CC, Green PA, Holland EA, Karl DM, Michaels AF, Porter JH, Townsend AR, Vörösmarty CJ (2004) Nitrogen cycles: past, present, and future. Biogeochemistry 70:153–226

Gee GW, Bauder JW (1986) Particle-size analysis. In: methods of soil analysis part 1: physical and mineralogical methods 2nd edition edn. American Society of Agronomy & soil science Society of America, Madison, Wisconsin, USA

Gourlet-Fleury S, Guehl JM, Laroussinie O (2004) Ecology and management of a neotropical rainforest. Lessons drawn from Paracou, a long-term experimental research site in French Guiana. Elsevier, Meppel, the Netherlands

Grau O, Peñuelas J, Ferry B, Freycon V, Blanc L, Desprez M, Baraloto C, Chave J, Descroix L, Dourdain A, Guitet S, Janssens IA, Sardans J, Hérault B (2017) Nutrient-cycling mechanisms other than the direct absorption from soil may control forest structure and dynamics in poor Amazonian soils. Sci Rep 7:45017. https://doi.org/10.1038/srep45017

Gross A, Turner BL, Goren T, Berry A, Angert A (2016) Tracing the sources of atmospheric phosphorus deposition to a tropical rain Forest in Panama using stable oxygen isotopes. Environ Sci Technol 50:1147–1156. https://doi.org/10.1021/acs.est.5b04936

Gupta UC, Lipsett J (1981) Molybdenum in soils, plants, and animals. In: Brady NC (ed) advances in agronomy, vol 34. Academic press, pp 73-115. https://doi.org/10.1016/S0065-2113(08)60885-8

Gutschick VP (1981) Evolved strategies in nitrogen acquisition by plants. American naturalist:607–637

Hammond D (2005) Tropical forests of the Guiana shield: ancient forests in a modern world. CABI Publishing, Cambridge

Hardy RW, Holsten RD, Jackson EK, Burns RC (1968) The acetylene-ethylene assay for N(2) fixation: laboratory and field evaluation. Plant Physiol 43:1185–1207

Hartmann A, Burris RH (1987) Regulation of nitrogenase activity by oxygen in Azospirillum brasilense and Azospirillum lipoferum. J Bacteriol 169:944–948

Hartwig UA (1998) The regulation of symbiotic N2 fixation: a conceptual model of N feedback from the ecosystem to the gene expression level. Perspectives in Plant Ecology, Evolution and Systematics 1:92–120. https://doi.org/10.1078/1433-8319-00054

Hedin LO, Brookshire ENJ, Menge DNL, Barron AR (2009) The nitrogen paradox in tropical Forest ecosystems. Annu Rev Ecol Evol Syst 40:613–635. https://doi.org/10.1146/annurev.ecolsys.37.091305.110246

Hill S (1988) How is nitrogenase regulated by oxygen? FEMS Microbiol Rev 4:111–129. https://doi.org/10.1111/j.1574-6968.1988.tb02738.x

Hofmockel KS, Schlesinger WH (2007) Carbon dioxide effects on heterotrophic dinitrogen fixation in a temperate pine Forest abbreviations: FACE, free air carbon dioxide enrichment; NPP, net primary production. Soil Sci Soc Am J 71:140–144. https://doi.org/10.2136/sssaj2006.110

Hu Y, Lee CC, Ribbe MW (2012) Vanadium Nitrogenase: a two-hit wonder? Dalton transactions (Cambridge, England: 2003) 41:10 https://doi.org/10.1039/c1dt11535a

Hungate BA, Dukes JS, Shaw MR, Luo Y, Field CB (2003) Atmospheric science. Nitrogen and climate change. Science (New York, NY) 302:1512-1513. https://doi.org/10.1126/science.1091390

Igarashi RY, Seefeldt LC (2003) Nitrogen fixation: the mechanism of the Mo-dependent nitrogenase. Critical reviews in biochemistry and molecular biology 38:351–384. https://doi.org/10.1080/10409230391036766

IPCC (2013) Climate change 2013: the physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change [stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex and P.M. Midgley (eds.)]. Cambridge, United Kingdom and New York, NY, USA. https://doi.org/10.1017/CBO9781107415324
Jones K (1977) The effects of moisture on acetylene reduction by Mats of blue-green algae in sub-tropical grassland. Ann Bot 41:801–806. https://doi.org/10.1093/oxfordjournals.aob. a085355

Kaiser BN, Gridley KL, Brady JN, Phillips T, Tyerman SD (2005) The role of molybdenum in agricultural plant production. Ann Bot 96:745–754. https://doi.org/10.1093/aob/mci226

van Kekem AJ, Pulles JHM, Khan Z (1996) Soils of the rainforest in Central Guyana. Series 2, Georgetown, Guyana

LeBauer DS, Treseder KK (2008) Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology 89:371–379. https://doi.org/10.1890 /06-2057.1

Leigh JA (2000) Nitrogen fixation in methanogens: the archaeal perspective. Current issues in molecular biology 2:125–131

Limmer C, Drake HL (1996) Non-symbiotic N2-fixation in acidic and pH-neutral forest soils: aerobic and anaerobic differen- tials. Soil Biol Biochem 28:177–183. https://doi.org/10.1016 /0038-0717(95)00118-2

Losos E, Leigh EG (2004) Tropical Forest diversity and dynam- ism: findings from a large-scale plot network. Univ, Chicago Press, Chicago

Matson AL, Corre MD, Burneo JI, Veldkamp E (2014) Free-living nitrogen fixation responds to elevated nutrient inputs in trop- ical montane forest floor and canopy soils of southern Ecuador. Biogeochemistry 122:281–294. https://doi. org/10.1007/s10533-014-0041-8

Mazerolle MJ (2017) AICcmodavg: model selection and multimodel inference based on (Q)AIC(c)

Menge DN, Levin SA, Hedin LO (2009) Facultative versus obli- gate nitrogen fixation strategies and their ecosystem conse- quences. Am Nat 174:465–477. https://doi.org/10.1086 /605377

Pajares S, Bohannan BJ (2016) Ecology of nitrogen fixing, nitri- fying, and denitrifying microorganisms in tropical Forest soils. Front Microbiol 7:1045. https://doi.org/10.3389 /fmicb.2016.01045

Pearson HL, Vitousek PM (2001) Stand dynamics, nitrogen accu- mulation, and symbiotic nitrogen fixation in regenerating stands of Acacia Koa. Ecol Appl 11:1381–1394. https:// doi.org/10.1890/1051-0761(2001)011[1381: SDNAAS]2.0.CO;2

Pearson HL, Vitousek PM (2002) Soil phosphorus fractions and symbiotic nitrogen fixation across a substrate-age gradient in Hawaii. Ecosystems 5:587–596. https://doi.org/10.1007 /s10021-002-0172-y

Penuelas J et al (2013) Human-induced nitrogen-phosphorus im- balances alter natural and managed ecosystems across the globe. Nat Commun 4:2934. https://doi.org/10.1038 /ncomms3934

Pérez S, Pérez C, Carmona M, Farina J, Armesto J (2008) Effects of labile phosphorus and carbon on non-symbiotic N2 fixa- tion in logged and unlogged evergreen forests in Chiloé Island, Chile. Rev Chil Hist Nat 81:12

Pham DN, Burgess BK (1993) Nitrogenase reactivity: effects of pH on substrate reduction and carbon monoxide inhibition. Biochemistry 32:13725–13731. https://doi.org/10.1021 /bi00212a043

Porder S, Asner GP, Vitousek PM (2005) Ground-based and remotely sensed nutrient availability across a tropical

landscape. Proc Natl Acad Sci U S A 102:10909–10912. https://doi.org/10.1073/pnas.0504929102

Quesada CA, Phillips OL, Schwarz M, Czimczik CI, Baker TR, Patiño S, Fyllas NM, Hodnett MG, Herrera R, Almeida S, Alvarez Dávila E, Arneth A, Arroyo L, Chao KJ, Dezzeo N, Erwin T, di Fiore A, Higuchi N, Honorio Coronado E, Jimenez EM, Killeen T, Lezama AT, Lloyd G, López-González G, Luizão FJ, Malhi Y, Monteagudo A, Neill DA, Núñez Vargas P, Paiva R, Peacock J, Peñuela MC, Peña Cruz A, Pitman N, Priante Filho N, Prieto A, Ramírez H, Rudas A, Salomão R, Santos AJB, Schmerler J, Silva N, Silveira M, Vásquez R, Vieira I, Terborgh J, Lloyd J (2012) Basin-wide variations in Amazon forest structure and function are medi- ated by both soils and climate. Biogeosciences 9:2203–2246. https://doi.org/10.5194/bg-9-2203-2012

R CT (2017) R: A language and environment for statistical com- puting., 3.4.3 edn., Vienna, Austria

ter Steege H, Pitman NCA, Phillips OL, Chave J, Sabatier D, Duque A, Molino JF, Prévost MF, Spichiger R, Castellanos H, von Hildebrand P, Vásquez R (2006) Continental-scale patterns of canopy tree composition and function across Amazonia. Nature 443:444–447. https://doi.org/10.1038 /nature05134

Reed SC, Cleveland CC, Townsend AR (2007) Controls over leaf litter and soil nitrogen fixation in two lowland tropical rain forests. Biotropica 39:585–592. https://doi.org/10.1111 /j.1744-7429.2007.00310.x

Reed SC, Townsend AR, Cleveland CC, Nemergut DR (2010) Microbial community shifts influence patterns in tropical forest nitrogen fixation. Oecologia 164:521–531. https://doi. org/10.1007/s00442-010-1649-6

Reed SC, Cleveland CC, Townsend AR (2011) Functional ecolo- gy of free-living nitrogen fixation: a contemporary perspec- tive. Annu Rev Ecol Evol S 42:489–512. https://doi. org/10.1146/annurev-ecolsys-102710-145034

Reed SC, Cleveland CC, Townsend AR (2013) Relationships among phosphorus, molybdenum and free-living nitrogen fixation in tropical rain forests: results from observational and experimental analyses. Biogeochemistry 114:135–147. https://doi.org/10.1007/s10533-013-9835-3

Reis CES, Dick DP, Caldas JS, Bayer C (2014) Carbon seques- tration in clay and silt fractions of Brazilian soils under conventional and no-tillage systems. Sci Agric 71:292–301. https://doi.org/10.1590/0103-9016-2013-0234

Robson RL, Postgate JR (1980) Oxygen and hydrogen in biolog- ical nitrogen fixation. Annu Rev Microbiol 34:183–207. https://doi.org/10.1146/annurev.mi.34.100180.001151

Roper MM, Smith NA (1991) Straw decomposition and nitroge- nase activity (C2H2 reduction) by free-living microorgan- isms from soil: effects of pH and clay content. Soil Biol Biochem 23:275–283. https://doi.org/10.1016/0038-0717 (91)90064-Q

Roskoski JP (1980) Nitrogen fixation in hardwood forests of the northeastern United States. Plant Soil 54:33–44. https://doi. org/10.1007/bf02181997

Rousk J, Bááth E, Brookes PC, Lauber CL, Lozupone C, Caporaso JG, Knight R, Fierer N (2010) Soil bacterial and fungal communities across a pH gradient in an arable soil. The ISME journal 4:1340–1351. https://doi.org/10.1038 /ismej.2010.58
Scott WJ (1957) Water relations of food spoilage microorganisms. In: Mrak EM, Stewart GF (eds) advances in food research, vol 7. Academic press, pp 83-127. https://doi.org/10.1016/S0065-2628(08)60247-5

Seefeldt LC, Hoffman BM, Dean DR (2009) Mechanism of Mo-dependent Nitrogenase. Annu Rev Biochem 78:701–722. https://doi.org/10.1146/annurev.biochem.78.070907.103812

Seitzinger SP, Garber JH (1987) Nitrogen-fixation and N-15(2) calibration of the acetylene-reduction assay in coastal marine sediments. Mar Ecol Progr Ser 37:8–73. https://doi.org/10.3354/Meps037065

Sexstone AJ, Revsbech NP, Parkin TB, Tiedje JM (1985) Direct measurement of oxygen profiles and denitrification rates in soil Aggregates1. Soil Sci Soc Am J 49:645–651. https://doi.org/10.2136/sssaj1985.03615995004900030024x

Sprent JI, Sprent P (1990) Nitrogen fixing organisms. Chapman and Hall, London, UK

Sullivan BW, Smith WK, Townsend AR, Nasto MK, Reed SC, Chazdon RL, Cleveland CC (2014) Spatially robust estimates of biological nitrogen (N) fixation imply substantial human alteration of the tropical N cycle. Proc Natl Acad Sci U S A 111:8101–8106. https://doi.org/10.1073/pnas.1320646111

Tang Y, Horikoshi M, Li WtUItVSRoPRP (2016) ggfortify: Unified Interface to Visualize Statistical Result fo Popular R Packages. The R Journal 8:12

Tiessen H, Chacon P, Cuevas E (1994) Phosphorus and nitrogen status in soils and vegetation along a toposquence of catastrophic rainforests on the upper Rio Negro. Oecologia 99: 145–150. https://doi.org/10.1007/bf00317095

Tripathi BM, Lee-Cruz L, Kim M, Singh D, Go R, Shukor NAA, Husni MHA, Chun J, Adams JM (2014) Spatial scaling effects on soil bacterial communities in Malaysian tropical forests. Microb Ecol 68:247–258. https://doi.org/10.1007/s00248-014-0404-7

Van Sundert K, Horemans JA, Stendahl J, Vicca S (2018) The influence of soil properties and nutrients on conifer forest growth in Sweden, and the first steps in developing a nutrient availability metric. Biogeosciences 15:3475–3496. https://doi.org/10.5194/bg-15-3475-2018

Venables WN, Ripley BD (2002) Modern applied statistics with S, Fourth edn. Springer, New York

Vitousek PM (1984) Litterfall, nutrient cycling, and nutrient limitation in tropical forests. Ecology 65:285–298. https://doi.org/10.2307/1939481

Vitousek PM, Sanford RL (1986) Nutrient cycling in moist tropical Forest. Annu Rev Ecol Syst 17:137–167. https://doi.org/10.1146/annurev.es.17.110186.001033

Vitousek PM, Menge DN, Reed SC, Cleveland CC (2013) Biological nitrogen fixation: rates, patterns and ecological controls in terrestrial ecosystems. Philos Trans R Soc Lond Ser B Biol Sci 368:20130119. https://doi.org/10.1098/rstb.2013.0119

Wagner F, Rossi V, Stahl C, Bonal D, Hérault B (2013) Asynchronism in leaf and wood production in tropical forests: a study combining satellite and ground-based measurements. Biogeosciences 10:7307–7321. https://doi.org/10.5194/bg-10-7307-2013

Walker TW, Syers JK (1976) The fate of phosphorus during pedogenesis. Geoderma 15:1–19. https://doi.org/10.1016/0016-7061(76)90066-5

Wardle DA, Walker LR, Bardgett RD (2004) Ecosystem properties and Forest decline in contrasting long-term Chronosequences. Science (New York, NY) 305:509-513. https://doi.org/10.1126/science.1098778

Weintraub SR, Taylor PG, Porder S, Cleveland CC, Asner GP, Townsend AR (2015) Topographic controls on soil nitrogen availability in a lowland tropical forest. Ecology 96:1561–1574. https://doi.org/10.1890/14-0834.1

Wood TE, Lawrence D, Clark DA (2005) Variation in leaf litter nutrients of a Costa Rican rain forest is related to precipitation. Biogeochemistry 73:417–437. https://doi.org/10.1007/s10533-004-0563-6

Wurzburger N, Bellenger JP, Kraepiel AM, Hedin LO (2012) Molybdenum and phosphorus interact to constrain asymbiotic nitrogen fixation in tropical forests. PLoS One 7:e33710. https://doi.org/10.1371/journal.pone.0033710

Yang J, Xie X, Wang X, Dixon R, Wang Y-P (2014) Reconstruction and minimal gene requirements for the alternative iron-only nitrogenase in *Escherichia coli*. Proc Natl Acad Sci 111:E3718–E3725. https://doi.org/10.1073/pnas.1411185111

Zaehle S, Jones CD, Houlton B, Lamarque J-F, Robertson E (2015) Nitrogen availability reduces CMIP5 projections of twenty-first-century land carbon uptake*. J Clim 28:2494–2511. https://doi.org/10.1175/jcli-d-13-00776.1

Zahran HH (1999) Rhizobium-legume Symbiosis and nitrogen fixation under severe conditions and in an arid climate. Microbiol Mol Biol Rev 63:968–989

Zhang X, McRose DL, Darnajoux R, Bellenger JP, Morel FMM, Kraepiel AML (2016) Alternative nitrogenase activity in the environment and nitrogen cycle implications. Biogeochemistry 127:189–198. https://doi.org/10.1007/s10533-016-0188-6

Zheng Y, Harris DF, Yu Z, Fu Y, Poudel S, Ledbetter RN, Fixen KR, Yang ZY, Boyd ES, Lidstrom ME, Seefeldt LC, Harwood CS (2018) A pathway for biological methane production using bacterial iron-only nitrogenase. Nat Microbiol 3:281–286. https://doi.org/10.1038/s41564-017-0091-5

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.