High nitrous oxide emissions from temporary flooded depressions within croplands

Nitrous oxide is an important greenhouse gas and emissions from managed ecosystems are directly correlated to anthropogenic nitrogen input. Here we have measured nitrous oxide emissions from ﬂ ooded depressions within croplands and from incubated soil samples. We scaled emissions to >20,000 comparable ﬂ ooded depressions across Zealand in Denmark using a deep-learning approach based on aerial photos and satellite images. We show that ﬂ ooded depressions within cultivated ﬁ elds, representing less than 1% of the total cultivated area, can release 80 times more nitrous oxide compared to rest of the ﬁ elds. Fluxes can remain high for more than two months after fertilisation and can account for 30 ± 1% of the nitrous oxide budget during that period. This highlights the urgent need for assessment of nitrous oxide hotspots, as managing these hotspots appear to represent an important part of the overall greenhouse gas emissions from managed croplands and an ef ﬁ cient mitigation action.


N
itrous oxide (N 2 O) is considered the third most important greenhouse gas after carbon dioxide (CO 2 ) and methane with a global warming potential 298 times that of CO 2 over a 100-year timespan 1 .Agricultural N 2 O emissions resulting from manure and nitrogen-based fertilisers have dominated global anthropogenic N 2 O emissions for decades [2][3][4] .Emissions of N 2 O are a consequence of nitrogen (N) inputs to soil and occur via a range of microbial transformation pathways including nitrification, nitrifier-denitrification and denitrification 4 .Oxygen availability is a key determinant of these N 2 O producing processes 5 .Under oxic conditions, nitrification converts ammonia to hydroxylamine, nitric oxide, nitrite and nitrate, and N 2 O is produced as a result of abiotic or biotic transformations of the intermediates 6,7 .The process of NO 3 -reduction via denitrification occurs under anoxic conditions, where N 2 O is an obligate intermediate in the denitrification pathway.
Temporal and spatial variability challenge our understanding of N 2 O emissions, particularly in relation to hotspots (or short periods known as hot moments) to explain high emissions 8,9 .We consider croplands that becomes temporarily flooded as a potential hotspot for fast denitrification process, resulting in a relatively large fraction of N being converted quickly to N 2 O. "Hot moments" have been shown to be the case in relation to drought followed by rewetting for managed grassland 6 and in depressions within the first four weeks of spring thaw 7 .We consider the period immediately following fertilization as a potential hot moment, particularly if nutrient-rich near-surface water quickly can reach flooded depressions, which therefore may become hotspots.The next steps to quantify the importance of flooded fields include: (1) to quantify the order of magnitude in N 2 O emissions from flooded depressions within contrasting fields, (2) to quantify the duration of high N 2 O fluxes from such depressions and finally (3) to quantify the area-contribution of such depressions within a larger region.These three targets have all been addressed in this study.Here, we focus on the distribution and role of field hotspots for N 2 O emission associated with spring fertilization.
Denmark has relatively flat croplands with rolling hills shaped by past glaciers and water and dominated by intensive farming 10 .It has enough topography to form gradients for surface or nearsurface water to flow towards depressions within individual fields 9 .Despite the spatial variation in elevation, fields receive the same order of magnitude of N from fertilization.That means that water with high concentrations of inorganic N (mainly nitrate) can accumulate in poorly drained areas 11,12 , which typically become flooded annually during winter and spring.Small-scale depressions and flooded fields have not been included in current national greenhouse gas budgets 11 despite cultivated soils generally have been identified as potential source areas for high N 2 O fluxes particularly related to early spring fertilisation prior to the main growing season [12][13][14] .Recently, N 2 O emission hotspots at plot and landscape scale have been associated with organic soils under warm well-drained conditions 15 .This study aims to assess this knowledge gap by quantifying the area and duration of flooded mineral soils within intensively managed croplands combined with detailed in-situ N 2 O flux measurements of both flooded and non-flooded areas within fields as well as N 2 O flux measurements under controlled laboratory incubations.

Results and discussion
In total 102 partly flooded fields (from less than 0.4 to 2.5 ha) have been identified within a 9000 km 2 area in the central part of Zealand in East Denmark in 2018 and 2019 (see area marked in red in Fig. 1).Subsequently, the entire island of Zealand was surveyed to identify the number and area of similar flooded fields for the period March-April 2019.For that we used a deeplearning approach to identify and segment single ephemeral water bodies in remotely sensed images at 3-m spatial resolution (Fig. 1).
Due to the dispersed emergence of ephemeral water bodies on cultivated land, their appearances/shapes could be distinguished based on a deep learning segmentation and a combination of optical and topographic remotely sensed information, due to the markedly different appearance from their surrounding cropland areas that are largely uniform in spectral properties and appearance (see "Methods" for detailed description).In total 20,063 individual flooded areas within fields were identified with an average size (± one std) of 802 ± 150 m 2 and in total equal to 23.1 ± 1.3 km 2 .That means that flooded areas within fields represent 0.50% of the total farmed area of Zealand equal to 4,973 km 2 .Of these flooded areas, the rim area (5 m rim of each flooded area) represented 0.26% of the farmed area (and equals 13.2 ± 0.8 km 2 ), which in the following is considered a key hotspot for N 2 O emissions.See Table S1 for accuracy metrics.
The duration of flooding during the two-month study period of March/April 2019 was assessed by mapping the near surface wetness index (see methods and supplementary Figs. 1, 2) and indicated that individual flooded depressions within fields varied in size but that the area-integrated area being flooded over the two months did not change significantly.This is consistent with the general pattern observed in the field.
Spring N 2 O bursts.During spring 2018/2019 after fertilisation, fluxes of N 2 O were measured once at 102 different flooded fields including measurements in the central part of the flooded area, at the rim of the flooded area and at drained adjacent locations representing a background value for N 2 O characterising the remaining part of the field (each location in 3 replicates).This measurement campaign is considered a snapshot of the spatial variability of N 2 O fluxes measured within two months after fertilisation and prior to the main growing season (regardless of main crop and time since fertilization).Flux measurements were made using static chambers and by measuring the N 2 O concentrations for 12 min in each chamber with a photo-acoustic gas monitor (INNOVA 1312, Lumasense Inc., Ballerup, Denmark) 16 .Flux measurements were accompanied by soil temperature and soil moisture measurements as well as soil samples taken for soil pH measurements and nitrate analysis.See details under "Methods".
N 2 O emissions in the rim area of the flooded area (the rim defined as being completely saturated and without standing water of more than 0-1 cm) were significantly higher compared to any other part of the fields (Fig. 2).Emission rates in the centre of the flooded depression remained generally lower but emitted highly variable amounts of N 2 O compared to the upland parts of the fields.On average, N 2 O fluxes from the rim area of the depressions were 652 ± 118 μg N m −2 h −1 , corresponding to more than 80 times higher average N 2 O flux than at adjacent locations (8 ± 2 μg N m −2 h −1 ).The latter value was consistent with previous numbers reported from Denmark 11,[17][18][19] .The N 2 O fluxes corresponded to soil temperatures between 2-7 °C and a mean soil pH of 6.4 ± 0.3 without significant differences between flooded and non-flooded parts of the fields.Soil moisture measurements for the centre and the rim of flooded fields revealed complete water saturation and water contents by volume varied between 45 and 52% corresponding to the soil porosity.For adjacent locations, the soil water content was in average 32 ± 8.3% by volume.Corresponding nitrate concentrations revealed a significant and positive correlation of high fluxes and high concentrations of nitrate in the rim area (see Supplementary Fig. 4).
Laboratory fertilisation experiment.During the autumn 2022, a sensitivity test was made based on intact soil cores collected from a typical loamy soil from one of the adjacent locations to a flooded field.A water and nitrate addition experiment was made simulating three water content levels (freely drained, ambient and total flooded conditions) with nitrate addition in three treatments: low (10 kg ha −1 ), normal fertilisation level (50 kg ha −1 ) 20 and triple fertilisation level (150 kg ha −1 ).The treatments corresponded to an addition of 3.8, 19.2 and 57.7 mg N as nitrate per core.The latter level is made to mimic split fertilisation practises 20 , which is commonly applied in the area and/or the fact that fertilisation near flooded fields is high as tracks in many cases (based on aerial photos, supplementary Figure 3) show that machineries are driving around rather than through flooded areas of fields.All treatments were made in three replicates; in total 27 cores were pre-incubated for 2 weeks and subsequently incubated at 7 °C for 75 days.N 2 O emissions were measured using an Isotopic N 2 O Analyzer (914-0027 from Los Gatos Research, Quebec, Canada).
Mean laboratory fluxes (Table 1) show that under normal fertilisation rates of nitrate, field and laboratory observations under completely flooded conditions are in the same order of magnitude; around 400 μg N 2 O-N m -2 h -1 .Field measurements are higher, which is in line with the fact that soil cores in the incubation represent emissions from N 2 O produced and released from the top 12 cm of the ploughing layer only.In contrast, measurements in the fields represent N 2 O produced and released from the entire soil profile; presumably from the top 30 cm mainly, representing the ploughing layer.
The results (Table 1) illustrate that the very same soil can release from 0.7 to 1072 ug N as N 2 O per m 2 as a mean over 8 weeks only reflecting the nitrate input and the environmental conditions in terms of drainage/flooding.For a "normal" nitrate addition, the flooded cores released on average 370 times more N 2 O as compared to ambient conditions, while under "high" nitrate addition the flooded cores released 29 times more N 2 O as compared to ambient conditions and 536 times as compared to dry conditions.These incubation-based values can be compared to the factor of 80 observed under field measurements, which indeed is in the low end of the increased rate due to flooding.Other soil characteristics measured (± st.dev) show that the bulk density was 1.4 ± 0.1 g cm -3 , the soil organic carbon content was 1.59 ( ± 0.09)% and the nitrogen content was 0.18 (±0.01)%.The average C:N ratio was 8.85 and soil pH was 6.25 (±0.03).These soil characteristics indicate that it is a common cultivated soil in Zealand 21 , which can be considered a Luvisol representing a soil type covering more than 95% of the area in Zealand 21 .
Landscape integrated N 2 O fluxes.The importance of N 2 O hotspots in regards to greenhouse gas mitigation is on the agenda in most countries with intensive farming [12][13][14][15] .Thus, scaling the above results to a larger area as Zealand is important for assessing the wider significance of this phenomenon, but is also related to sources of uncertainty.Consequently, we here present a conservative estimate based on the following discussion and criteria's: 1.The study period (2018/2019) was relatively normal in terms of precipitation.However, both years reflect an increasing tendency of warmer summers (particularly the summer 2018) and less precipitation during the main growing season and more episodic precipitation events during autumn 22 .
2. We assume that the 102 fields represent an average in relation to land use practice, although they cover a range of fields depending on crop types, crop rotation, specific farmer fertilisation strategies etc.Many of these variables will vary from year to year but will typically do so within a four-year cycle.Thus, snapshot measurements of 102 flooded fields show that a typical agricultural soil in Zealand in Denmark can release 80 (during incubation up to 140) times more N 2 O solely depending on nitrate availability if kept under completely water-saturated conditions.3. Field measurements were completed within a two-month period which equal the period where flooding was observed.The duration of flooding was verified based on an assessment of a near-surface moisture index as derived from (see supplementary figure 1).From this we consider the snapshot measurements of 102 lakes to represent a mean flux across two months.This is further justified by the incubation experiment revealing high fluxes for at least 8 weeks after nitrogen addition (Fig. 3).4. We have identified 20,063 similar partly flooded fields across Zealand.We assume that N 2 O emissions from the 102 flooded areas within fields are representative for the 20,063 flooded areas across croplands in Zealand.This is justified by the overall variation and mean value of soil development and soil characteristics (including nitrate  Table 1 Mean fluxes in μg N 2 O-N m −2 h −1 (standard deviations given in brackets) for the eight weeks of the laboratory experiment for three nitrogen levels (Low = 10 kg ha -1 , Normal = 50 kg ha -1 and High = 150 kg ha -1 ).Scaling these results to a larger area as Zealand; considering (1) flooded fields during a two-month period, (2) only the 5-m rim of flooded fields and (3) an average rate of N 2 O emissions being 80 times higher as compared to the rest of the field, we conclude that these hotspots represent 30 ± 1% of the total N 2 O emission during that period.The uncertainty range refer solely to the uncertainties related to the area estimates of flooded fields (see Fig. S2).We consider the 30% being a minimum value, which is most likely markedly higher.Factoring in the entire flooded area (and not only the 5-m rim area) will increase the percentage of N 2 O emission from flooded fields, however, results presented also reveal a larger variability in fluxes if the entire flooded area would be included.
The results also highlight that the total area of flooded area within fields is not a suitable measure to estimate the total N 2 O loss, as a strong gradient in emissions may appear along the topography near flooded fields.The extremely high N 2 O fluxes are measured at the rim of the depressions (presumably partly due to receiving water from uphill soils) whereas lower fluxes are measured in the central part of the depressions.These lower fluxes measured in the central part may point at that nitrate reduction at the rim is fast and that water is depleted with respect to nitrate before reaching the central parts of the depressions and that nitrate from deeper sources (groundwater) is less important for the observed high N 2 O fluxes.The depression size-dependency and scaling of N 2 O fluxes across fields with e.g., different soil types, crops and management strategies are some of the most important unknowns in order to estimate the environmental and economic loss of N 2 O (and nitrogen) in depressions.
In general, main factors controlling the formation of hot moments of N losses have been identified to be temperature, precipitation and fertilisation.Roy et al. ( 2021) 13 concluded, based on a US country-level assessment on nitrogen budgets, that a variety of hot moments and hotspots collectively may account for about 63% of the total surplus nitrogen balance.The nitrogen loss as N 2 O was part of the budget, but not specifically quantified due to lack of hotspot measurements 13 2022) 24 used an ecosystem model to quantify how land management practices and weather can control nitrous oxide emissions.They concluded that spring-thaw-related N 2 O pulses in northern agroecosystem regions of the United States occur for a short period, but increase the emissions from croplands and grasslands by up to 16% at the national level.
The above mentioned references clearly demonstrate a growing awareness of hotspots (and hot moments) being critical and important for robust N budgets for croplands, but also that hotspots, as flooded depressions within fields, remain a major knowledge gap which need to better understood and quantified in order to improve mitigation actions and optimise the use of fertiliser.We conclude that the extreme importance of a relatively small area within field for area-integrated N 2 O emissions is completely new in literature and suggests an urgent need to address management of these hotspots.Our study focuses on the end of non-growing season/start of the growing season.However, the non-growing season has recently been shown to account for 38-43% of the annual N 2 O emissions depending on crop type based on a study of Canadian croplands 25 .Taking the conclusions from the present study may markedly increase the importance of the non-growing season and the transition to the growing season.The same study 25 also showed that the percentage of added N released as N 2 O (the emission factor) varied from 0.29-7.32%for yearly emissions but only 0.12-2.02%for the growing season.Our experimental incubation suggests that more than 5-8% of added nitrate is being released as N 2 O to the atmosphere (see Supplementary Table 1).Such high emission factors stand in contrast to the commonly accepted factor of 1% used in national greenhouse gas inventories (NGGI) 11 .However, due the limited area of flooding (less than 1%), results presented here cannot be used to argue that the currently used emission factor based on an average for individual fields needs to be adjusted.In addition, despite the focus on the non-growing season as crops start to grow, plants will compete with microorganisms for N and thereby reduce the N 2 O production.Finally, incubation experiments are excellent for keeping environmental factors constant during measurement campaigns, but more sites representing different soil types, climate and management practice are important to reduce uncertainty before scaling our numbers to national or larger region N 2 O budgets.This is relevant in the light of only one soil type being included in the incubation experiment.Although a representative site and soil type in terms of soil characteristics has been chosen for the incubation experiment, variations in soil type are expected to have consequences for observed N 2 O fluxes and the overall N 2 O budget.Similarly, we have in this study focused on an actual time period (an actual year and spring time only), but year to year variations can be critical for both the area within field being flooded, but also the time period being flooded, and also the environmental conditions while flooded.Particularly, soil temperature, precipitation events and nitrate availability while flooded are expected to impact the net N 2 O budget for flooded fields.
Implications and recommendations.Due to gravity, even minor depressions are the final destination for water runoff from surrounding upland fields, with reduced infiltration rates and longer residence times of water compared to the upland areas.This means that they are more prone to being waterlogged and becoming anoxic, particularly in the non-growing season.Hydrological models and continuous measurements of nitrate are needed to identify: (1) the upland areas, (2) runoff patterns towards depressions (runoff versus groundwater) and (3) the residence time of water and nitrate in depressions.
Large precipitation events (producing waterlogged areas or floods) occur mainly during the non-growing season.They are a result of cyclic patterns and are considered to increase in frequency and intensity in the future 26 .In relation to N 2 O hot moments, early spring fertilisation events are considered very critical.A similar hot moment is likely after the "first precipitation event" following a prolonged drought 27 .This is a typical autumn situation, where the timing and intensity of rain will determine how much nitrate is being flushed out from the part of the field surrounding the depressions 28,29 .Currently, there are no estimates of the overall proportion of nitrate in fields that can be mobilised and potentially transported to depressions and subsequently converted to N 2 O and free nitrogen (N 2 ) by denitrifying microorganisms.It is important to note, that in contrast to N 2 O, the release of N 2 is not an issue for the climate.But any loss of nitrogen represents a loss of resources seen from a farming perspective.This is also the case in paddy soils where high fluxes have been measured after soil drainage 30 but here nitrification is considered the main process as a source of N 2 O. Without a loss, the nitrogen pool would otherwise have been kept in the soil system and increased growth and crop yields in the following growing season.
Thus, N 2 O hotspots are an essential focus point for planning actions to reduce the overall greenhouse gas footprint of the agricultural sector.Currently, N 2 O hotspots are not included in national N 2 O budgets 11 .The lack of measurements is likely the main reason for this avoidance, probably reflecting the logistical problems of obtaining solid N 2 O flux measurements with flux chambers in muddy/flooded fields.The results shown here suggest that nitrate pulses and associated N 2 O bursts can have a duration of days to several weeks, which highlights that hightemporal measurements are necessary to capture the N 2 O bursts.This has also been proven in wetland ecosystems 5 and the technology for implementing similar setups in cultivated fields is available.
Policies for taking cultivated land out of production in order to reduce net greenhouse gas emissions, as a means for securing a successful green transition, are on the political agenda right now.In many regions, farmers and environmentalists agree that a promising option to reduce net greenhouse gas emissions is to convert existing croplands to nature.This agreement between farmers and environmentalists is based on a willingness of spending public money to cover some of the cost for such land use changes.Therefore, the critical question is which part of the cropland should be converted in order to obtain the highest benefits for the climate, nature and farmers.From a farmer's perspective, the periodically flooded fields are often considered the least cost-effective to operate, least accessible (part of the time), most in need of re-sowing, and, despite the hassle, usually supplying no or low yields.Consequently, flooded areas within fields appear to be giant N 2 O emitters and represent a low hanging fruit in relation to greenhouse gas mitigation actions.
We conclude that fluxes of N 2 O in and near flooded depressions within fields are extremely high when nitrate levels are high (particularly after fertilisation) and fluxes in the rim area of flooded fields are 80-140 times higher compared to similar soil types being effectively drained.Fertilisation (which is often repeated during springtime) can be considered a hot moment and flooded fields an example of extreme hotspot, showing very high nitrate values and corresponding extremely high N 2 O emission for at least 8 weeks.During such conditions, laboratory results suggest that more than 5-8% of added nitrate is being released as N 2 O to the atmosphere.Flooded depressions within fields are found widespread in Zealand but represent only an area of less than 1%.Scaling these field-measured fluxes to Zealand in Denmark show that during a two-month period, the rim area of flooded fields representing 0.26% of the farmed land can release 30% of the total N 2 O emission.Consequently, land use changes taking flooding of fields into account are very likely a strong winwin opportunity for climate, nature and farmers.

Methods
Mapping flooded fields.This study demonstrates a framework for mapping ephemeral water bodies in surface depressions of agricultural fields with satellite imagery at very high spatial and temporal resolution.This combination of rich information in both space and time, which is otherwise not feasibly achievable from common freely available satellite remote sensing systems, is obtained from the use of the constellation of 130 nano-satellites from the Planet company.We trained a deep convolutional network with the widely used U-Net architecture 31,32 for semantic segmentation to predict pixel-wise ephemeral water coverage on remotely sensed PlanetScope images at 3-m spatial resolution supported by a detailed (40 cm) LiDAR-based DEM (digital elevation model).We expected the fusion of optical and elevation data to help discriminating water bodies from their surroundings.
Therefore, we devised a strategy based on a set of rules to identify ephemeral water bodies in satellite images -, according to which we manually delineated and annotated 1824 individual objects across Zealand that were verified using very highresolution nationwide aerial photography provided by the Danish Agency for Data Supply and Infrastructure (SDFI).These othophotos were collected at 12.5 cm spatial resolution as a patchwork of 4-channel RGBNir images varying between 29 March -18 April, which despite being useful for supporting visual assessment during the training data collection, lack the temporal consistency for large-scale analysis.We ensured that the labels represent a wide range of landscape characteristics, by defining and balancing five subcategories depending on how water and vegetation are spatially distributed.The collected water bodies showed a distinct visual contrast and a clear surface depression corresponding to the object coverage, which excludes surfaces with soil saturation without notable topographic variation.In addition, we only labelled objects that showed a high temporal variability that we visually assessed in very highresolution aerial photography covering a period of seven years, as well as based on the presence of permanent vegetation, such as trees, shrubs, or emergent plant types, which we excluded from the delineation.The manually labelled water bodies were used to train and validate the deep-learning model, which was then used to predict ephemeral water bodies covering the whole study area.
The input features given to the deep learning model were PlanetScope level 3B Ortho Tile imagery at a 3-m spatial resolution and daily temporal resolution that consists of three visible (VIS: 0.40-0.75µm) and one near-infrared (NIR: 0.75-1.30µm) bands with a maximum of 20% cloud cover (Fig. 1).We selected a three-week temporal window from 29 March to 18 April to compose our median PlanetScope mosaic with an average of 36 cloudless mosaics for each 25-km x 25-km tile after cloud removal.This period corresponds to an observed peak in the water signal, which is useful to maximise the contrast between the ephemeral water bodies and their surrounding cropland, and minimise noise from occasional snow cover.
We obtained the temporal dynamics of the ephemeral water bodies using six months of PlanetScope data over 102 fieldcollected sites that showed a high intra-annual variability (Fig S1 ), with a strong peak during March and April, followed by a sudden decrease -in contrast to permanent small water bodies that displayed a more stable trajectory.We used a single band thresholding technique of the NIR band to calculate the water fraction (WF), based on the bimodal distribution function representing reflectances over water and land surfaces.The water and land limit thresholds (WL: 1300 and LL: 4500) were empirically identified to separate the pure water and land signals from the mixed pixels, as WF(%) = (LL − NIR) × 100/(LL-WL) 31 .We calculated the temporal variability from all >80% cloud-free daily data using a 5-day median between January and July 2019.
We supplemented our PlanetScope composites with data from the Danish Elevation Model (DK-DEM) at a 0.4 m resolution showing elevation in relation to the average sea level (Fig. 1), which is based on 415 billion points covering the entire country, collected by aerial LiDAR and is provided by the Danish Agency for Data Supply and Infrastructure (SDFI).We included a digital Terrain Model (DTM), a Digital Surface Model (DSM), and a DTM-derived slope layer.The DSM contains information about the height of surface objects, such as vegetation or buildings in relation to mean sea level, which are removed from the DTM.We hypothesized that information about differences in terrain and surface object height can be useful in differentiating areas of low reflectance, such as shadows that are other than ephemeral water bodies.
The deep-learning model was trained with gradient descent using a learning rate of 0.0003 and an adaptive learning rate optimization algorithm (Adam), which computes separate learning rates for each parameter 33 .We employed the crossentropy loss function during the 500 epochs of training with a batch size of 16, a patch size of 256, along with batch normalization 34 and using rectified linear units (ReLU) as the activation function 35 .To increase the variability of the dataset, we augmented the training data with artificial and transformed training images (e.g.vertical flip, horizontal flip, rotation, gaussian noise, pixel dropout) 36 .After each epoch, the model was evaluated on an independent validation set (20% of the labels) and we selected the final model with the lowest validation loss value of 0.5670.The predicted map conveys the location of each detected ephemeral water body and its area coverage in m 2 .
After the prediction, we used morphological filtering with a 6×6 structuring element to eliminate noise resulting from image artefacts and edge effects, and to smoothen the occasionally fuzzy object boundaries caused by low contrast conditions.We used the independent validation set to calculate the prediction accuracy of the filtered product (Fig. S").The total accuracy of predicting water body pixels was 72%, the balanced total accuracy of background as well as water body pixels was 86%, and the precision of predicting water body pixels was 50% (Table S1).Additionally, we used the derived error rates to proportionally adjust our rea estimates and calculate their 95% confidence intervals 37 .Ephemeral water bodies represent areas being a continuum from flooded fields to wet soils.Given the nature of ephemeral water bodies, where no exact borders can be drawn between target and background in relation to the spatial extent, the pixel-based accuracy obtained is considered satisfactory.At instance level, we report 34% precision and 34% recall (34% Fscore) when associating predicted and labelled individual water bodies with a minimum intersection over union (IoU) of 20%.The true positives have an average IoU of 46% with their corresponding labelled objects.Additionally, we report an overall water body counting error of only 3% on the validation set (216 labelled objects, 223 predicted), showing that the model has a good balance between false positives and false negatives.The post-processed results became the basis for calculating the total area coverage, as well as estimating the area under the rim-effect using a 5-m buffer from each object outline, as established by observations from the field.
The uncertainty of mapping ephemeral water bodies in general, and with a coverage of <150 m 2 in particular, is comparatively high, because of their dynamic nature, as they may emerge after a rainfall event but are often rapidly drained away afterwards.This can cause a variation within image composites using different dates, which is inherent to the use of very high-resolution satellite image analysis for regional scale analysis.An additional source of potential inaccuracy is caused by relatively large fields that are uniformly wet or show a wet, narrow swath along the field boundary, creating contrast conditions that are suboptimal for accurate detection.
N 2 O flux measurements.In situ fluxes of N 2 O were measured along transects at each flooded area within fields and after the first fertilisation event (in March).The general strategy for farmers is to fertilise first time when possible followed by two more intense fertilisation events at the end of March or start of April.The typical total load of N 150-200 kg N for these three fertilisation events.The main growing season is considered from mid April.The high number of sites in this study reflects the focus on robust mean values rather than replicates reflecting crop type, time since fertilisation or general management practices.Sites were consequently selected randomly within the study area.
At each of the 102 locations, N 2 O fluxes were measured at three positions along a transect (centre, rim and upslope) in three replicates.Flux chambers made of PVC tubes (diameter of 24 cm, height of 20 cm) were inserted 2 cm into the soil.Flux measurements were made by placing a lid on the tube and measuring the N 2 O concentration for 12 min in each chamber with a photo-acoustic gas monitor (INNOVA 1312, Lumasense Inc., Ballerup, Denmark).The closed-chamber technique is known to create a bias by altering the diffusion gradient between soil and chamber headspace (Livingston and Hutchinson, 1995).However, several studies have shown that this bias can be overcome by applying a non-linear regression method to describe the gas exchange (Forbrich et al. 2010).In our case, the flux in each chamber was estimated by fitting the increase in N 2 O concentrations to a nonlinear function for measurements starting 2 min after the lid was placed.For N 2 O flux calculations, concentrations of N 2 O were converted from ppm to μmol m -3 using the ideal gas equation and subsequently the flux was estimated using: where F is the N 2 O flux (μmol m -2 hour -1 ), C is the concentration of N 2 O (μmol m -3 ), t is time (hours), V is the volume (m 3 ) and A is the area of the soil cylinder (m 2 ).The volume was based on height measurements for each of the flux measurements.
Nitrate and pH measurements.Soil samples were collected from 0-10 cm around the flux tube (n = 3) and kept cold at 7 °C and dark until processed in the laboratory within 3-4 days.Samples (2 g) were subsequently mixed with 10 ml deionized water and shaken for 10 min and the eluent analysed for nitrate content on standard for Ion Chromatography (Dionex Aquion, Thermo Scientific, Waltham, MA, USA) and pH values measured using a VWR pHenomenal (CO3100L).
Laboratory incubations.In-situ soil samples were taken on October 6 th , 2022, on a conventional agricultural field near Sorø, Denmark, at 55°44´N, 11°59´E.Samples (n = 40) were collected in plastic tubes (12 cm in length, 7 cm in diameter), immediately transported to the laboratory and pre-incubated for almost four weeks at 7 °C.
Prior to drying and wetting all samples had soil water content of approximately 16.5 ± 1.1% by weight representing the ambient conditions.Nine random samples were selected for free drainage and drying.The samples were dried at room temperature for seven days.After drying for seven days, the water content was 11.3 ± 1.3%.Nine other samples were wetted with deionized water 24 h prior to gas measurement to ensure completely water-filled pores and at least 2 mm standing water on top of the soil surface.The remaining nine samples represent ambient conditions.
For N-treatment, we used laboratory grade KNO 3 dissolved in deionized water to represent three levels of N input (10 kg ha -1 , 50 kg ha -1 and 150 kg ha −1 ).For N-treatment, we used sterile BD Plastipak 1 mL syringes to inject 10 mL of KNO 3 + H 2 O in each sample.Five injections were made at 5 cm depth and 5 injections were made as the syringe was moved to the surface.From this we assume that the entire core was homogeneously treated with N. Injections were made on the 9th of November 2022.
Soil bulk density was calculated as a weight difference after complete drying at 105 °C of samples after the incubation.The content of C and N was determined on elemental analyser (CE 1110, Thermo Electron, Milan, Italy) and pH was measured in double deionized water (Silex H 2 O, 1:2.5 weight to volume) and measured using a potentiometer after 3 min of calibration.
Measurements of N 2 O production were conducted over eight weeks (starting 9th of November 2022) where soil samples were placed in 2 L hermetically sealed glass bottles.Between measurements samples were placed at 7 °C in a fridge open to a ventilated room but covered by parafilm (to avoid evaporation).During measurements samples were out of the refrigerator for 9 min.Room temperature was between 12 and 17 °C.From each bottle, two Versilon SE200 tubes of 120 cm were connected to an isotopic N 2 O-analyzer.On the inlet tube a 0.1 mm Whatman paper filter was attached.Luer lock "fittings" were attached to the end.Gas was continuously sampled using Isotopic N 2 O Analyser, 914-0027, from Los Gatos Research (Quebec, Canada), with a flow rate of 143 mL min −1 , and specific volume of 50 cm 3 .The analyser logged N 2 O in ppm at 1 Hz frequency, each measurement was 8 min.Each sample was measured daily for 17 days, then every second day for seven days and then twice a week until week 10.
The first 200 s of each measurement was discarded as the analyser needed time to mix up residue air and sample from the bottle.Final N 2 O fluxes were calculated similarly as described for field measurements.

Fig. 1
Fig. 1 Examples of flooded fields in the study area Zealand in Denmark.a Total study area Zealand in Denmark, the study area of 102 flooded depressions is marked as a red square and shown in detail in (b).c 1-3: Examples of flooded fields with various contrast and background conditions as they appear at various spatial resolutions, showing the resulting segmented boundaries.

Fig. 2
Fig. 2 The frequency of N 2 O fluxes measured at three categories of fields.a Example of a flooded field, which for the particular field represents about 0.2% of the entire field equal to the mean across Zealand during intensive fertilisation in March/April.b Histogram showing the frequency of mean N 2 O fluxes measured at three categories of fields: the central part of the flooded field with standing water, the rim area of the flooded part and the ambient field surrounding the flooded part.
concentrations, soil pH, %C and C/N) from the 102 fields and soil maps covering Zealand and the area representing the 102 flooded areas.
. Wallman et al. (2022) 15 measured N 2 O emissions related to fertiliser treatments and observed extreme N 2 O bursts associated with the actual fertilisation.However, the study did neither include measurements under anaerobic conditions or lateral movement of fertiliser across the field.Wagner-Riddle et al. (2020) 19 specifically included hotspots and hot moments (fertilisation events) as part of recent review of mitigation actions to minimise nitrogen emissions from agroecosystems.Kim et al. (2013) 23 made a metaanalysis showing a strong dependency of N 2 O emissions on soil nitrogen availability which is aligned with our field observations.They showed that N 2 O emissions increased linearly with increased N input in the range of 10-150 kg N per ha, but also that N 2 O emissions increased non-linearly (almost exponentially) with N inputs larger than 150 Kg N. Del Grosso et al. (

Fig. 3
Fig. 3 Measurements of N 2 O fluxes for 74 days during incubation experiment.Daily measurements of N 2 O fluxes for the first 17 days and less frequently afterwards for (a): drained, ambient and (b): flooded conditions (standard errors are shown as bars).Please see Supplementary Fig. 5 for details for all three soil moisture levels.