Next Article in Journal
Breeding Colony Dynamics of Southern Elephant Seals at Patelnia Point, King George Island, Antarctica
Next Article in Special Issue
From Standard Weather Stations to Virtual Micro-Meteorological Towers in Ungauged Sites: Modeling Tool for Surface Energy Fluxes, Evapotranspiration, Soil Temperature, and Soil Moisture Estimations
Previous Article in Journal
Erratum: Zhao, H.; Zhang, H.; Miao, C.; Ye, X.; Min, M. Linking Heat Source–Sink Landscape Patterns with Analysis of Urban Heat Islands: Study on the Fast-Growing Zhengzhou City in Central China. Remote Sens. 2018, 10, 1268
Previous Article in Special Issue
Substantially Greater Carbon Emissions Estimated Based on Annual Land-Use Transition Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combining High-Resolution Remote Sensing Products with a Crop Model to Estimate Carbon and Water Budget Components: Application to Sunflower

1
CESBIO, Université de Toulouse, CNES/CNRS/INRAE/IRD/UPS, 31400 Toulouse, France
2
Agence De l’Environnement et de Maîtrise de l’Energie (ADEME), CEDEX 1, 49004 Angers, France
3
INRAE, UMR 1248 AGIR, 31326 Castanet-Tolosan, France
4
INRAE, USC 1439 CESBIO, 31100 Toulouse, France
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(18), 2967; https://doi.org/10.3390/rs12182967
Submission received: 3 August 2020 / Revised: 7 September 2020 / Accepted: 9 September 2020 / Published: 11 September 2020
(This article belongs to the Special Issue Remote Sensing of Land–Atmosphere Interactions)

Abstract

:
The global increase in food demand in the context of climate change requires a clear understanding of cropland function and of its impact on biogeochemical cycles. However, although gas exchange between croplands and the atmosphere is measurable in the field, it is difficult to quantify at the plot scale over relatively large areas because of the heterogeneous character of landscapes and differences in crop management. However, assessing accurate carbon and water budgets over croplands is essential to promote sustainable agronomic practices and reduce the water demand and the climatic impacts of croplands while maintaining sufficient yields. From this perspective, we developed a crop model, SAFYE-CO2, that assimilates high spatial- and temporal-resolution (HSTR) remote sensing products to estimate daily crop biomass, water and CO2 fluxes, annual yields, and carbon budgets at the parcel level over large areas. This modeling approach was evaluated for sunflower against two in situ datasets. First, the model’s output was compared to data acquired during two cropping seasons at the Auradé integrated carbon observation system (ICOS) instrumented site in southwestern France. The model accurately simulated the daily net CO2 flux (root mean square error (RMSE) = 0.97 gC·m−2·d−1 and determination coefficient (R2) = 0.83) and water flux (RMSE = 0.68 mm·d−1 and R2 = 0.79). The model’s performance was then evaluated against biomass and yield data collected from 80 plots located in southwestern France. The model was able to satisfactorily estimate biomass dynamics and yield (RMSE = 66 and 54 g·m−2, respectively). To investigate the potential application of the proposed approach at a large scale, given that soil properties are important factors affecting the model, a sensitivity analysis of two existing soil products (GlobalSoilMap and SoilGrids) was carried out. Our results show that these products are not sufficiently accurate for inclusion as inputs to the model, which requires more accurate information on soil water retention capacity to assess water fluxes. Additionally, we argue that no water stress should be considered in the crop growth computation since this stress is already present because of remote sensing information in the proposed approach. This study should be considered a first step to fulfill the existing gap in quantifying carbon budgets at the plot scale over large areas and to accurately estimate the effects of management practices, such as the use of cover crops or specific crop rotations on cropland C and water budgets.

Graphical Abstract

1. Introduction

Sunflower (Helianthus annuus L.) was introduced by Spaniards in Europe during the 16th century as an ornamental plant species before its oil started to be used for food during the 19th century. Currently, sunflower is cultivated on more than five continents, with Ukraine and Russia being the largest producers, followed by the EU-28. In France, sunflower crops cover 590,000 ha, with an annual production of 1.36 M tonnes, resulting in an average yield of 2.3 t·ha−1 (source: FAO, http://www.fao.org/faostat/en, 10-year average between 2007 and 2016). In France, this oilseed crop species is grown mainly in the central-western (Nouvelle Aquitaine) and southwestern (Occitanie) regions. Sunflower is a summer crop; it is mainly grown under rainfed conditions and has a low water-use [1]. In southwestern France, the temperate Mediterranean-like climate (dry and warm summers) results in sunflower crops being subjected to frequent and severe water stress conditions that impact crop productivity at both field and landscape/regional scales. These variations should be quantified to better manage croplands. Moreover, in this context, accurate estimation of yield, carbon (C), and water budgets at the field level over large areas remains challenging. The heterogeneous character of landscape in terms of climatic and soil conditions as well as the diversity of agricultural practices make the upscaling of agronomic models difficult. Thus, a precise understanding of the cropland C budget is needed, as well as agronomic indicators (e.g., water-use efficiency) useful for monitoring crop species, such as sunflower.
Crop models were first used for research purposes [2,3] in the early 1970s. During the 1980s, they became more diversified and complex to meet management needs and to support field-scale decision-making ([4,5] and farming applications [6]). Today, a considerable number of crop models exist. They can be classified into three groups:
  • C-driven models, or what we call the “de Wit school” crop growth simulation models [7], are based on the assimilation of C by photosynthetic parts of the plant, such as the SUCROS [8] or GECROS [9] models. These models, on the basis of their hierarchical architecture, are complex, require numerous parameters, and are difficult to calibrate; thus, they are difficult to upscale to large areas.
  • Water-driven models are based on a predetermined relation between crop growth rate and transpiration according to water productivity parameters. Compared to that of the previous class, the approach of this class results in models whose structure is less complex, leading to fewer parameters. Among the most representative ones are the AquaCrop [10] or AqYield [11] models.
  • Radiation-driven models, which have been constructed on the basis of [12]’s approach, convert intercepted or absorbed solar radiation into aboveground or total biomass via a radiation-use efficiency parameter, significantly reducing the number of photosynthetic parameters needed. Most crop models are radiation driven, including the CERES [13], EPIC [14], STICS [15], SAFY [16], and SAFY-WB [17] models.
The model presented in the current study belongs to both the first and third groups. Indeed, it evolved from (i) the SAFY-WB model [17], which belongs to the third group and allows the simulation of biomass production, yield, and water fluxes and budget, and from (ii) the SAFY-CO2 model [18], which calculates biomass production as the difference between plant assimilation and respiration.
Regarding sunflower, crop growth models are either generic, with a set of parameters specific to sunflower (AquaCrop, CropSyst [19], EPIC, STICS, or WOFOST [20]) or specific to sunflower (OilCrop-SUN [21]; SUNFLO [22]).
To reduce the number of model parameters and to increase the robustness of the models upon upscaling, several remote sensing data assimilation approaches for crop models have been developed. The benefits of using HSTR satellite products in crop models have already been demonstrated in several studies [18,23,24,25,26]. These approaches are either based on establishing relations between remote sensing products and crop productivity [27] or on calibrating the model’s parameters [28,29,30]. Such approaches have already been applied to sunflower crops for biomass estimation [31,32]. Those studies showed that accounting for the spatial heterogeneity in sunflower development is crucial for the accurate estimation of crop dynamics; nevertheless, these models are not able to estimate cropland C budgets. In this context, we propose a newly developed agro-meteorological model, the Simple Algorithm for Yield, Evapotranspiration, and CO2 fluxes estimates (SAFYE-CO2), that assimilates HSTR products, and we validate the model for use in sunflower. Adapted from models developed by [17,18], this model allows us to simulate cropland growth, yield, and CO2 and water fluxes at a large scale because it requires little input data (e.g., little to no crop management data).
Initially, the Simple Algorithm For Yield estimates (SAFY) model was developed to estimate crop biomass and production and was first validated for use in winter wheat [16], followed by sunflower and maize [23,31]. Those studies demonstrated the benefits of assimilating the HSTR green area index (GAI) derived from optical remote sensing satellite data for crop monitoring by calibrating several phenological parameters as well as light-use efficiency. Later, this model was coupled with a water module based on the FAO-56 approach [33] to estimate cropland evaporation (Ev) and transpiration (TR) and thus water stress and water needs. This resulted in the SAFY-WB (water balance) model, which was validated for use in wheat [17] and maize [34]. More recently, Pique et al. [18] showed the benefit of using HSTR GAI products for assessing the components of the C budget over winter wheat crops by adapting the SAFY model to a C module, which led to a new model called SAFY-CO2. In line with these studies, a new model, coupling the SAFY-CO2 model with the same water module as that in the SAFY-WB model, has been constructed and is evaluated here. The new model is called the SAFYE-CO2 model (E stands for evapotranspiration (ETR)), which aims to assess daily biomass production and CO2 and water fluxes, as well as annual yields and the C and water budget over croplands. The development and validation of the SAFYE-CO2 model for sunflower was driven by the following perspectives: (i) Simulation of the water flux and of the soil water provides an opportunity to estimate the effects of water stress on crop growth, which could improve biomass and CO2 flux estimates for summer crops (e.g., sunflower); and (ii) SAFY-CO2 has already been validated for use in winter wheat [18], and validation for use in sunflower would allow us to simulate C and water budgets for sunflower-winter wheat crop rotations, which are very common in southwestern France.
In terms of upscaling the model over large areas and simulating the effects of water stress on crop production and C budgets, the objective of this study was twofold: (i) To analyze whether a coupled remote sensing/modeling approach, such as that involving SAFYE-CO2, is able to accurately reproduce sunflower function, i.e., biomass, yield, and CO2 and water fluxes, at the ground–atmosphere interface; and (ii) to evaluate the impact of current soil map products on the simulation accuracy as well as the additive value of coupling SAFY-CO2 with a water module to improve sunflower simulation.
In the following sections, the materials and methods are described. Afterward, the results section starts with a validation of the SAFYE-CO2 model against biomass, yield, and CO2 and water flux measurements at the Auradé integrated carbon observation system (ICOS) instrumented site. A sensitivity analysis of the model to the soil parameters was then performed. Finally, SAFYE-CO2 was validated via biomass and yield data collected from more than 80 fields during five years, covering a wide range of soil and climatic conditions. Maps of the components of the C budget over a cropping season are then presented. The study ends with a discussion on the benefits and limitations of our approach, raising issues concerning the assessment of accurate C and water budgets at the plot scale over large areas and discussing potential improvements for simulating those variables.

2. Materials

The proposed study used a large amount of data from several sources to calibrate, validate, and upscale our model. In addition to the meteorological inputs and soil property data needed for the model, high-resolution remote sensing data were required for calibration. The validation of the model is based on CO2 and water flux measurements, biomass production, and yield. The last two originated either from destructive measurements carried out at the FR-Aur site (see Section 2.2) or from surveys carried out in 2013, 2014, and 2015 [35]. The following sections present these data as well as the study area.

2.1. Study Area

The study area is located in southwestern France and encompasses 4 administrative departments (Haute-Garonne, Gers, Aude, and Ariège; see Figure 1). It is part of the Regional Spatial Observatory—Southwest (OSR-SO, https://osr.cesbio.cnrs.fr/) located near Toulouse and includes the instrumented ICOS agricultural site of Auradé (FR-Aur). The region is characterized by a temperate and Mediterranean-like climate, with dry and warm summers (see Table 1). The Garonne River flows across the study area from south to north and has shaped the topography into terraces and hill landscapes that resulted from the erosion of old terraces. This has led to heterogeneous soils called “boulbènes” (silty soil on terraces) and “terrefort” (clayey-limestone soils on hills), characterized by contrasting textures (especially with respect to the clay content) and with variations in depth within each type of soil.
According to the 2015 French Land Parcel Identification System (LPIS, or the Registre Parcellaire Graphique in French), agricultural activity represents 66% of the study area (focus in Figure 1), which includes forests, grasslands, orchards, etc. Sunflower crops constituted 19.1% of the agricultural area.
In the region, sunflower is usually sown between April and May and is harvested between the end of August and September. The climatic conditions during the crop cycle of sunflower are characterized by high radiation and strong evaporative demand, while low and erratic rainfall occurs in the summer (Table 1). The selected study area is very seldom irrigated, like elsewhere in France (<5% of the sunflower area [36]). Owing to its deep root system, its thick and hairy epidermis, and its adaptive capability (through a reduction of the photosynthetic rates), sunflower is moderately sensitive to water deficit. However, hydric stress occurring before pre-flowering can lead to early senescence, ultimately causing losses in yield and oil content.

2.2. Flux Site and Experimental Instrumentation

Since 2005, the FR-Aur site (Figure 1) has been highly instrumented for the continuous collection of micrometeorological, meteorological, soil, and vegetation data (see [37] for more details) in the context of the OSR-SO. FR-Aur is also part of the European Research Infrastructure Consortium ICOS (https://www.icos-cp.eu/), which is a European network documenting long-term greenhouse gas concentrations in the atmosphere and greenhouse gas exchange between the surface and the atmosphere in several ecosystems under different pedoclimatic conditions (crops, grasslands, forests, etc.). The plot area is 23.5 ha and is highly heterogeneous because of the topography (17.3 m height difference), which results in different soil depths (from 100 cm to more than 200 cm) and thus different crop yield potentials. The FR-Aur field belongs to a conventional grain farm where only mineral fertilization is applied. The plot is characterized by a rapeseed/winter wheat/sunflower/winter wheat rainfed rotation. Sunflower was grown during the 2006–2007 and 2015–2016 cropping seasons (see Table 1 for sowing and harvest dates). Hereafter, the crop years will be abbreviated by “CY” followed by the harvest year (e.g., CY-07 for the 2006–2007 crop season).

2.2.1. Meteorological Data

The SAFYE-CO2 model requires meteorological data as the main input, while vegetation remote sensing was used to constrain the modeled dynamics. Specifically, air temperature (TAIR) at 2 m, incoming global radiation (RG), rain, and reference ETR (ET0) were needed. These variables were provided either from in situ measurements acquired at the FR-Aur site for local validation or from SAFRAN (système d’analyse fournissant des renseignements adaptés à la nivologie) reanalysis data [38] for the spatial analysis.
At FR-Aur, meteorological data are recorded every 30 min and averaged on a daily basis (see [37] for more details). These data were used to calculate the ET0 using the Penman–Monteith formula [33].
The SAFRAN meteorological data were used as inputs to simulate the sunflower experiments of 2013, 2014, and 2015. SAFRAN reanalysis data are provided over a grid of 8 × 8 km² over France. In this study, the climatic data were estimated for each field using a bilinear interpolation of the 4 closest grid points.
The climatic conditions of all simulated years are summarized in Table 1. Each variable was averaged from the 15th of April to the 15th of September, corresponding to the average sunflower growing season in the region. The crop cycles of CY-07 and CY-16 were quite similar, although 2016 was slightly hotter and sunnier than 2007 (Table 1), leading to slightly higher biomass production in 2016. In 2013 and 2014, the growing seasons were colder compared to the season during 2015. The 2013 growing season (CY-13) was the wettest, and the CY-14 growing season was the driest. Note for the sake of clarity, only sowing and harvest dates of CY-07 and CY-16 are indicated in Table 1. Indeed, more than one field was simulated for the others CY (respectively 9, 21, and 49 for CY-13, CY-14, and CY-15), thus the sowing and harvest dates vary by plot, which explains the empty cases of Table 1.

2.2.2. Biomass and Yield Data

Destructive measurements were made at FR-Aur to characterize the sunflower dry aboveground mass (DAM) dynamics during CY-07 and CY-16 (see Table 2). In CY-07, 30 plants were collected monthly inside the flux tower footprint area distributed on either side of the eddy covariance tower towards the prevailing winds. The DAM was obtained by drying samples in an oven at 65 °C for 48 h. The total DAM produced (before harvest) was measured by cutting plants from a 24-m² area and measuring them. The grain yield (YLD) was provided by the farmers for CY-07, and the grains were separated from the heads and weighed to obtain yield for CY-16.
Extensive field campaigns were conducted in 2013, 2014, and 2015 to collect biomass and yield data across the study area (Figure 1). In 2013, 9 elementary sampling units (ESUs, 20 × 20 m2) corresponding to 8 different fields were monitored during the growing season. For each ESU and at 4 times during crop development (the last sampling was performed just before harvest), 5 samples of aboveground vegetation consisting of 10 adjacent plants were collected (except for the first date, when only 3 plants were collected), oven dried at 50 °C for 48 h, and then weighed. DAM was calculated by dividing the total biomass by the area sampled, calculated as the product of the interrow width and the length of the sampling area. On the last sampling date, the grains were separated from the heads and weighed to calculate the yield.
In 2014 and 2015, yield data were collected for more than 70 and 116 fields, respectively. They were provided by the farmers, and because some of the yield values were averaged over several fields, only yield data corresponding to one field were considered (see Table 2). Two or three visits were made to each field in the 2014 and 2015 campaigns to characterize the crop growth status (e.g., canopy heterogeneity, presence of weeds, pest and disease pressure). From those observations, we defined quality scores, and low-quality data were discarded. These data corresponded to fields where (i) vegetation-free areas larger than 2 m2 were observed, and (ii) weeds covered more than 50% of the field. After discarding all these data (115), we kept data for 70 fields for the model evaluation (Table 2). None of the fields monitored in 2013, 2014, and 2015 were irrigated.

2.2.3. Flux Data

Turbulent fluxes of CO2, water vapor (ETR and latent heat), sensible heat, and momentum were measured continuously using the eddy covariance method [39,40,41] by the combination of an open-path fast infrared gas analyzer (LI-7500, LI-COR) and a 3-D sonic anemometer (CSAT3, Campbell Scientific). The flux tower was installed near the center of the field (see Figure 1). EdiRe software (Clement, 1999) was used to calculate the ETR and net CO2 fluxes (see [18] and [37] for details). Flux filtering, quality controls, and gap filling procedures were performed in accordance with the CarboEurope-IP recommendations (http://www.carboeurope.org/).
The net ecosystem exchange (NEE) of CO2 was partitioned into gross primary productivity (GPP) and ecosystem respiration (RECO) according to the methods proposed by [42] and adapted for croplands by [37]. During the period when bare soil was exposed (deduced from field observations), the GPP was set to 0, while the RECO was set equal to the NEE.

2.2.4. Soil Data

To run the SAFYE-CO2 model, data on the soil moisture at field capacity (θFC) and the permanent wilting point (θWP) are needed for each considered field to estimate the daily soil available water content (AWC). This information can be derived either from in situ sampling and laboratory analysis or can be inferred from soil texture products by applying pedotransfer functions (PTFs). In this study, we used both methodologies: data from in situ soil sampling at FR-Aur and data from the global soil databases GlobalSoilMap (GSM [43]) and SoilGrids (SG [44]) at both FR-Aur and at the other 78 fields. These soil products were used to estimate the θFC, θWP, and soil depth for all the considered fields in this study, including those at FR-Aur. The global soil databases provide the fraction of sand, clay, and silt at a 90- (GSM) or 250-m (SG) resolution as well as the soil depth. SG is based on machine learning and remote sensing data, while GSM uses a data-mining technique and a bottom-up approach [45].
We used the PTFs proposed by [46] (Equations (9) and (10)) to calculate the volumetric soil moisture at field capacity (θFC) and at permanent wilting point (θWP). The estimations of the θFC, θWP, and soil depth for FR-Aur using in situ samples, SG, and GSM products are summarized in Table 3 for the 3 soil layers considered in the proposed modeling approach (see Section 3.1).
The soil types of all the fields monitored in 2013, 2014, and 2015 were characterized based on the textures derived from the GSM and SG products (for the top 30 cm of the soil) and considering the United States Department of Agriculture (USDA) classification system (see Figure 2).
Considering all monitored fields, the GSM products suggest that 64% of the soils are clay loams, 23% are loams, 11% are silty clay loams, and the remaining soils vary between clay, silty clay, and silty loams. Compared with the GSM products, the SG products suggest less heterogeneity among the soil types and classify 87% of the soils as clay loams, 12% as loams, and less than 1% as silty clay loams.

2.3. Satellite Data and Products

2.3.1. Multisatellite Optical Images

The authors of [47,48] highlighted the usefulness of assimilating HSTR satellite products into crop models, as it increased the accuracy in reproducing crop dynamics and/or reduced the number of model parameters. Thus, the accuracy of our approach is partly determined by the availability of remote sensing images. Indeed, to successfully calibrate some of the model parameters (see Section 3.3), a certain number of images at key dates had to be available. An extensive set of HSTR products derived from several sensors was thus used in this study. The main characteristics in terms of the spatial resolution, temporal sampling, and spectral bandwidth of each sensor are detailed in Table 4, and the number of images (all in all, 147) from each sensor is shown in Figure 3. For the purpose of the study, a temporal series of remote sensing products, averaged at the plot scale, were needed. Products of different spatial (from 8 to 30 m) and temporal (from 1 to 16 days) scales were thus used in order to get the most complete time series (Figure 3).
The authors of [23] showed that the advantages in terms of time sampling of using GAI estimates from multiple sensors outweigh the drawbacks in terms of sensor measurement mismatch. In their study, the authors compared the normalized difference vegetation index (NDVI) of different sensors and showed that, during the growing season (NDVI > 0.2) and over various crops (sunflower, maize, wheat, etc.), linear relationships (R values between 0.97 and 0.99 and relative root mean square error (rRMSE) values between 5.3% and 10.9%) exist between pairs of sensors (including DEIMOS-1, Formosat-2, Landsat-8, and SPOT4 (Take5) [49]). The KALIDEOS processing chain (http://kalideos.cnes.fr) for atmospheric, radiometric, and geometric corrections was used to derive surface reflectance values from the satellite data. The mean geometric correction accuracy was close to 0.2 pixels [50], which is satisfactory considering the area of the studied fields.

2.3.2. From Image Reflectance to GAI Estimates

The remotely sensed GAI estimates were derived from surface reflectance values using the Biophysical Variables Neural NETwork (BV-NNET) tool [51]. BV-NNET is a trained artificial neural network based on the inversion of the radiative transfer model PROSAIL [52]. The GAI estimates take into account the spectral and directional information of each sensor, which are of paramount importance when working with multiple satellites. The GAI value estimated by BV-NNET is an “effective GAI”, which does not consider the clumping effect and could lead to underestimation compared to that of the “true GAI”, mainly for high GAI values [53]. For each simulated field, the mean GAI was calculated taking into account all the pixels of the considered plot after the application of an offset of 10 m to avoid edge effects.

3. Methods

3.1. The SAFYE-CO2 Model

The SAFYE-CO2 model was constructed by coupling the SAFY-CO2 model [18] with the water module of the SAFY-WB model [17], which was adapted from the FAO-56 method [33]. SAFY-CO2 is a parsimonious agro-meteorological model running on a daily basis. It is driven by TAIR and RG as meteorological inputs. It simulates crop dynamics (GAI), crop production (DAM, YLD), CO2 fluxes (GPP, RECO, and NEE), and the annual net ecosystem C budget (NECB). Coupling with the water module, which can easily be added or removed, is needed to estimate crop TR and soil Ev and thus soil water content (SWC). From those data, soil water stress and plant water needs can be calculated. Adding this module requires two additional meteorological inputs, rain and ET0.
Since all equations involved in the SAFYE-CO2 model have already been described and detailed in other studies (SAFY [16]; SAFY-WB [17]; SAFY-CO2 [18]), we will summarize here the main methods and the interactions between the water module and the SAFY-CO2 model. All related parameters discussed in this section are summarized in Table 5.
In the SAFYE-CO2 model, GPP is first estimated as a function of the RG, a climatic efficiency (εC), the absorbed photosynthetically active radiation (APAR) calculated from the GAI, and the light-interception coefficient (KEXT), a function of the light-use efficiency taking into account the fraction of the diffuse radiation (ELUEA and ELUEB; see [18]), and a function accounting for the decrease in canopy photosynthesis capacity during senescence (sR10). The daily estimates of GPP are modulated by (i) a thermal stress function defined by three cardinal temperatures limiting crop growth (TMIN, TMAX, TOPT), and (ii) a water stress function (WtSV, see below). Net primary productivity (NPP) is estimated as the difference between the GPP and autotrophic respiration (RA). RA is estimated by the methods proposed by [54], which separates RA into two components: maintenance respiration (RM) and growth respiration (RGR). RM is computed as a function of the NPP of the previous day, a “Q10-type” RM equation [55] (with two crop-specific parameters: Q10 and R10), and the sR10 coefficient. RGR is estimated considering the daily amount of biomass produced and using a RGR coefficient parameter (YG) as proposed by [56] and later improved by [57]. Heterotrophic respiration (RH) is estimated from a simple Q10 first-order exponential equation (Q10_H) and a reference respiration parameter (RhREF).
CO2 fixed by the vegetation is converted into dry biomass via a plant C parameter (CVEG) and partitioned between aboveground DAM and belowground parts according to the methods proposed by [58]. This process depends on three parameters (f0, f, and c) that were adapted for sunflower crops based on data from [59].
Like for SAFY, the daily increase in the GAI is estimated from the daily increase in DAM, the specific leaf area (SLA), and partition-to-leaf function (defined by parameters PlA and PlB). The decrease in the GAI starting with senescence is estimated based on a threshold from the growing degree days (starting from crop emergence as detected by the remote sensing data), from the senescence parameter (SENA), and from the rate of senescence (SENB). The final yield is simply deduced from the maximum value of DAM multiplied by a harvest index (HI). The set of equations used to describe the previous methods are detailed in the study by [18].
In SAFY-CO2, the soil is a single layer whose characteristics (parameters) are used to calculate RH. In the water module of SAFYE-CO2, the soil is divided into 3 layers: A surface layer of constant depth (DEVAP), subject to Ev and TR; a root layer, subject to TR only, that increases in depth during the growing season at a rate (VR) that tracks crop growth and limited by a maximum root depth (DROOTS); and a constant depth (D) deep layer. The last one is used to initialize the water content of the root layer during root growth and to control for the effects of percolation and drainage. The soil module is detailed in the paper of [17]. The soil AWC is essentially updated daily, taking into account rain (inputs), Ev, and TR simulated based on the FAO-56 methods [33]. The AWC is defined by the water fluxes (input and output) and the total water capacity (TWC) (i.e., TWC = θFC − θWP). TR depends on a maximal TR coefficient (KcbMAX), a reduction coefficient (KTRP), a critical relative humidity parameter (Dfe), the GAI, and the maximum relative filling rate between the surface and root layer.
The water module also computes two water stress functions, one that drives Ev (WtSE, Equation (1), considering the AWC of the surface layer) and another that drives the growth of vegetation (WtSV, Equation (2), considering the surface and root layers) through GPP estimation. These functions represent the main link between the vegetation and the soil module and are defined as follows:
WtS E = ( SM θ WP θ FC θ WP ) ,
WtS E = ( SM θ WP θ FC θ WP ) ,
RH TOT = max { RH evap , RH root } ,
where SM (m3·m−3) is the SWC of the surface layer and where θFC (m3·m−3) and θWP (m3·m−3) are the water content at field capacity and at the permanent wilting point, respectively. Dfe is the critical relative filling rate for water stress, and RHevap and RHroot are the relative humidity of the surface and root layers, respectively.
In the present study, the calculation of Ev was modified from that of [17] according to the methods of [60] because differences were observed between simulated and observed ETR during bare soil periods and because of the dependence of Ev on ET0. The reason is that the FAO-56 method for estimating Ev was developed to reproduce homogeneous grass stands, whereas the soil is often completely bare during the fallow period. Additionally, [60] asserted that, in addition to being ET0 dependent, Ev depends on the hydraulic state of the soil surface layer. When ET0 is high, the first several centimeters of the soil rapidly become dry, slowing the capillary rise of the remaining water. A new equation for estimating Ev was thus adopted:
Ev = K BS ( Ev A ET 0 + Ev B ) ET 0 ( 1 F cover ) WtS E ,
where ET0 (mm·d−1) is the reference evapotranspiration, KBS (−) is a correction coefficient for bare soil, EvA (1/(mm·d−1)) and EvB (−) are the evapotranspiration coefficients, and FCOVER is the vegetation cover fraction. FCOVER is estimated from the maximal vegetation fraction cover (KCOV), an exponential coefficient (ECOV), and the GAI. FCOVER can be derived from remote sensing and thus drives the model, but previous quantitative experiments did not show significant enhancement of the outputs; thus, we calculated it here from the GAI and two parameters (see Table 5).
The water module is impacted by the vegetation module through the calculation of FCOVER (which considers the GAI), and in turn, the vegetation module is impacted by the water module through the water stress function (WtSV).

3.2. Model Parameters

In addition to requiring meteorological data, the SAFYE-CO2 model requires crop-specific parameters. These parameters came from several sources. Most of them could be obtained from the literature, some could be estimated via in situ measurements, and the rest were obtained from calibration (see Table 5). The model could run using literature and in situ parameters only, but when the last category of parameters is not constrained by remote sensing, the model cannot cope with changing conditions in the fields, which renders the outputs inaccurate. The calibration procedure based on the satellite-derived GAI allows optimization of the phenological parameters and the A parameter of the light-use efficiency function. Once these parameters are calibrated, the outputs are representative of the GAI extraction area (e.g., pixels, fields, farms). To calculate the plot-scale C budget, data supplied by farmers, concerning the amount of C imported into fields as organic amendments, are eventually needed.

3.3. Model Parametrization and Calibration

The calibration procedure incorporating remotely sensed GAI was detailed in the report by [18]. Six phenological parameters (PlA, PlB, SENA, SENB, SLA, D0) as well as the A parameter of the light-use efficiency (see Table 5) were calibrated during this process for each field and each year independently so that the simulated GAI matched the remotely sensed GAI. This procedure allows the model to accurately reproduce crop growth. This inversion process is adapted from the Nelder–Mead simplex method [61], in which a priori values and minimum/maximum boundaries are set for each retrieved parameter. A quadratic difference cost function—root mean square error (RMSE)—is minimized using the adapted simplex between the remotely sensed GAI and the simulated GAI. For minimization, a global approach based on 30 different a priori samples is used to avoid local minima while preventing unnecessary runs. Once completed, the inversion process results in a set of parameters, allowing accurate reproduction of the remotely sensed GAI while avoiding biomass estimates as well as water and CO2 flux estimates (see Section 3.1).

3.4. Model Validation and Sensitivity Analysis

The validation strategy adopted in this study is summarized in Figure 4. It relies on several types of datasets. First, we evaluated the ability of SAFYE-CO2 to simulate crop growth and development (leaf area index (LAI) and DAM), YLD, and CO2 (GPP, RECO and NEE) and water (Ev, TR) fluxes. To accomplish this, the model was run for the two crop seasons of sunflower at the FR-Aur site (CY-07 and CY-16) using in situ (IS) meteorological data. Additionally, in situ measurements were used to estimate θFC and θWP and the soil depth. These results are presented in Section 4.1.
For upscaling of the model, the impact of the soil products (GSM and SG were used to estimate θFC, θWP, and soil depth) on the model’s performance must be evaluated. For this purpose, the SAFYE-CO2 model was run on the same two crop years at FR-Aur using the GSM or the SG products to estimate soil water retention parameters. The outputs were compared with those obtained with in situ measurements instead of soil products. Additionally, the results obtained with SAFYE-CO2 were compared to those from the original SAFY-CO2 version to assess whether the coupling of SAFY-CO2 with the water module improves the model’s outputs or not; Section 4.2.1 describes those results. Finally, the SAFY-CO2 and SAFYE-CO2 models were used to simulate DAM and YLD over a large area. In addition to the two crop years at FR-Aur, the fields from the 2013, 2014, and 2015 campaigns were simulated, and in situ measurements of DAM and YLD were used for validation. In situ meteorological variables were used to simulate the crop years at FR-Aur, while data from SAFRAN were used to simulate the other fields/crop years. The findings related to this exercise are described in Section 4.2.2. All simulations were performed from 15 November to 14 November of the following year, usually representing one cropping season in the area of study.

4. Results

4.1. Local Assessment of the Model’s Performances

The SAFYE-CO2 model was first tested against the in situ data of the FR-Aur site. Figure 5 presents the correlation between the simulated and observed GAI, DAM, and CO2 (GPP, RECO, and NEE) and water (ETR) fluxes for the two crop years of sunflower. Only data corresponding to the crop cycle (i.e., from emergence to harvest) were included to compute the GPP statistics, while the entire simulation period was considered to compute the model performances when estimating the other flux variables.
Since the model is calibrated by minimizing the quadratic difference between simulated and remotely sensed GAI, the estimation of GAI is used to verify that the inversion procedure is correctly performed rather than being used to validate the modeled GAI. Consequently, the statistical results related to the GAI estimation are very good (R2 = 0.99, RMSE = 0.099 m2·m−2, and rRMSE = 6.9%), indicating that the sunflower’s GAI dynamics are accurately reproduced by the model. The DAM is also robustly estimated in terms of the correlation (R2 = 0.96) and error (RMSE = 66 g·m−2 and rRMSE = 16.2%) data, and the yield (only 2 values at FR-Aur) is also estimated with a low error (RMSE = 53.8 g·m−2). Accurate estimation of biomass, and thus yield, is a prerequisite for assessing NECB with high precision. CO2 fluxes are correctly simulated, although compared with CY-07 CO2 fluxes, CY-16 CO2 fluxes are estimated with lower accuracy due to the lack of satellite images throughout the agricultural season. Indeed, for CY-16, only four images were available during sunflower development (it was the beginning of Sentinel’s mission). In particular, no images were available from 30 July to 28 September (see Figure 3), leading to a senescence period that was poorly documented by the remotely sensed GAI and a model-simulated senescence stage that was too early. The poor representation of the senescence phase in 2016 led to an underestimation of GPP and RECO. This period is easily identifiable on the scatter plots (Figure 5C,D) and corresponds to the green points underestimated by the model. Indeed, the mean bias of the simulated GPP and RECO is higher for CY-16 (respectively −2.21 and −0.69 gC·m−2·d−1) than for CY-07 (respectively −0.11 and −0.16 gC·m−2·d−1) as well as the slopes associated to these data (respectively 0.75 and 0.65 for CY-07 and 1.22 and 1.00 for CY-16). Nevertheless, when enough images are available to describe the crop dynamics (notably key dates, such as those corresponding to emergence, maximum vegetation, start of senescence, harvest), the crop production as well as the CO2 and water fluxes are estimated with acceptable accuracy (see [18]).
GPP and RECO are estimated with relatively low errors (RMSE of 2.36 gC·m−2·d−1 and 0.96 gC·m−2·d−1, respectively, and rRMSE of 39.5% and 40.9%, respectively) and good correlations (R2 of 0.81 and 0.83, respectively).
In the proposed approach, the errors associated with GPP and RECO compensate for each other when estimating NEE (which is calculated from the difference between the two variables). This resulted in an accurate estimation of the NEE (RMSE of 0.97 gC·m−2·d−1 and R² of 0.83) despite the lack of images during the 2016 senescence phase.
The model is also able to reproduce the ETR (Ev+TR) dynamics of the two crop years with good accuracy (RMSE of 0.6 mm·d−1, rRMSE = 40%, and R2 of 0.79). For the CO2 fluxes, the ETR is underestimated by the model for CY-16 during the senescence period, which led to a decrease in TR before the end of the observed crop cycle (Figure 5F).
In the proposed approach, RH depends only on the soil temperature, while it is known to be impacted by many other factors (e.g., AWC and soil organic C content). In their study, [18] discussed taking into account AWC in the estimation of RH. We thus estimated RH according to the [62] Equation III.5 that accounts for AWC. During the two simulated crop years at FR-Aur (data not shown), the temperature-dependent-only formulation performed better (RMSE = 0.5 gC·m−2·d−1 and R2 = 0.36) than the combined temperature- and AWC-dependent equation (RMSE = 0.73 gC·m−2·d−1 and R2 = 0.06), even when considering in situ AWC to estimate RH (RMSE = 0.61 gC·m−2·d−1 and R2 = 0.15). We therefore kept the first-order exponential equation that depends only on soil temperature for the rest of the study.

4.2. Large-Scale Application

4.2.1. Impact of Soil Data and of Coupling with the Water Module on Model Performance

The two sunflower crop years at FR-Aur were simulated with the SAFYE-CO2 model using field capacity, the permanent wilting point, and soil depth estimated from (i) the in situ measurements, (ii) the GSM products, and (iii) the SG products. These two crop years were also simulated with (iv) the SAFY-CO2 model to evaluate the benefits of coupling with the water module compared to the previous version. The different evaluated configurations of the model are referred to as the following: (i) SAFYE-CO2IS, (ii) SAFYE-CO2GSM, (iii) SAFYE-CO2SG, and (iv) SAFY-CO2.
Figure 6 summarizes the statistics (RMSE and R2) obtained when testing the different configurations of the model in simulating the target variables (GAI, DAM, GPP, RECO, NEE, and ETR) for CY-07 and CY-16. There are no statistical data on ETR for SAFY-CO2, as this version does not simulate water flux. The differences observed between the different tested configurations are mainly due to water stress. In addition to simulating water fluxes (Ev, TR), the water module impacts vegetation photosynthesis and growth through a water stress function (WtSV). Therefore, correctly estimating AWC is essential for accurate simulation of crop growth.
Regarding CY-07, the GAI and RECO are estimated with nearly the same accuracy for each model configuration. Concerning the GAI, the good performance of each model configuration indicates that all the configurations were able to accurately fit the remotely sensed GAI regardless of the soil data used. The errors in RECO are also similar between the different configurations for two reasons: (i) The heterotrophic component of RECO is estimated in the same way for all versions; and (ii) the autotrophic component, depending on the GPP, differs between versions only at the end of the crop cycle when water stress occurs. However, greater differences occur for the DAM, GPP, NEE, and ETR estimates, where SAFYE-CO2GSM and SAFYE-CO2SG show low accuracy. This is due to the overestimation of the water stress at the end of the season, which in turn underestimates the water fluxes and GPP and thus the DAM. Indeed, TWC is underestimated for both SAFYE-CO2GSM and SAFYE-CO2SG during this period (see Section 2.2.4), leading to an underestimation of AWC, which induces stress at the end of the cropping season.
For CY-16, the same observation can be made concerning the estimations of the GAI. Each version is able to accurately reproduce GAI dynamics. However, SAFYE-CO2IS largely outperformed SAFYE-CO2GSM and SAFYE-CO2SG for all the other variables. SAFYE-CO2SG and SAFYE-CO2GSM show high errors when estimating GPP and ETR because, like CY-07, those versions largely overestimated water stress, inducing early senescence and thus an underestimation of the CO2 and water fluxes.
Despite the observed water stress during both simulated years (at the end of the crop cycle), the best performances in estimating CO2 fluxes were achieved by the SAFY-CO2 model that did not simulate water stress, indicating that the water stress is already taken into account through the remotely sensed GAI. Additionally, the water module adds an additional level of complexity, while the performance of SAFY-CO2 is already very high.
Figure 7 shows the components of the NECB estimated with the different configurations of the model compared with observed values: The net annually cumulated CO2 fluxes (net annual ecosystem productivity (NEP)), the C input as seeds (CINP), and the C exported at harvest (CEXP). Because of the poor estimates of TWC deduced from GSM and SG products, SAFYE-CO2GSM and SAFYE-CO2SG overestimate NEP and underestimate CEXP for both simulated years. Both models’ configurations overestimate water stress, especially at the end of the season when AWC is low, which leads to premature senescence and thus an early harvest. In the absence of vegetation, the only CO2 flux is RH (flux towards the atmosphere), so shortening of the crop cycle by the two model configurations leads to an overestimation of NEP. The same reason explains the underestimation of CEXP. The more stressed the vegetation is, the lower the total biomass production and the yield are. On the other end, both SAFYE-CO2IS and SAFY-CO2 were able to accurately reproduce the NECB for CY-07 (−4.3 and −10.8% for SAFYE-CO2IS and SAFY-CO2, respectively) and for CY-16 (+15 and −6.2% for SAFYE-CO2IS and SAFYE-CO2, respectively), as both configurations were able to accurately simulate the components of the C budget.

4.2.2. Biomass and Yield Estimates

Apart from the two sunflower crop years available at FR-Aur, the model was evaluated for DAM and YLD over numerous fields in 2013, 2014, and 2015 (see Section 2.2.2) to assess its ability to reproduce those variables under a variety of climatic, soil, and management conditions. Since no soil in situ measurements were available at the plot scale in 2013, 2014, and 2015, the SAFYE-CO2IS configuration could not be tested.
Figure 8 presents DAM and YLD simulated by SAFYE-CO2GSM (A,B), SAFYE-CO2SG (C,D), and SAFY-CO2 (E,F). Although SAFYE-CO2GSM and SAFYE-CO2SG showed good overall performance at estimating DAM over the whole dataset (R² of 0.78 and 0.74, respectively; RMSE of 110 and 119 g·m−2, respectively), these two models’ configurations had difficulty accurately estimating the final DAM for the 2013 field campaign. This could be explained by the water stress simulated at the end of the season, which slows down crop growth and thus underestimates final biomass. This inaccurate estimation of water stress was the result of poor estimation of the TWC.
In the previous sections, TWC derived from GSM and SG were compared to the value derived from in situ measurements performed at FR-Aur (Table 3). The values associated with the global soil databases were lower than those observed in situ, and only a few variations were observed in terms of texture (particularly for SG, Figure 2). The TWC distributions derived from GSM and SG for all the studied fields are presented in Figure 9 and were compared to values derived from in situ measurements collected nearby (i.e., 45 fields located in the northwestern part of the study area—experiment described by [18]). The values of TWC derived from the global soil databases are low compared to those obtained from the previous studied fields, and very little spatial variation is present. The TWC derived from GSM ranges between 0.11 and 0.14 m3·m−3, while the TWC derived from SG is even less variable—between 0.12 and 0.13 m3·m−3. Considering the in situ measurements, soils in the study area present more variation, with the TWC between 0.10 and 0.16 m3·m−3. These low and quite spatially constant TWC values derived from the global soil databases undoubtedly lead to an overestimation of water stress and thus to an underestimation of the final DAM. Nevertheless, SAFY-CO2, which does not simulate water flux or water stress, better estimated DAM (R2 of 0.91 and RMSE of 73 g·m−2) and especially better reproduced the final highest DAM values (red dots in Figure 8E), confirming that SAFYE-CO2GSM and SAFYE-CO2SG simulated too much water stress in 2013.
Regarding YLD estimates (Figure 8B,D,F), all versions have difficulty reproducing the observed values accurately. This is possibly due to several factors. First, YLD is calculated from the maximum value of DAM using a constant HI. The HI of sunflower is known to vary with the length of the grain filling period and with temperature [63] but also with water stress [64]; moreover, there is some additional variations due to crop variety [65]. These factors are not yet accounted for in the proposed approach. Additionally, the observed yields were provided by the farmers, and may have high uncertainty. Indeed, compared to the “true” yield, those data are subject to significant grain losses at harvest.

4.3. C Budget Assessment at the Landscape Scale

Since the SAFYE-CO2 model has difficulty reproducing sunflower growth if no accurate information on soil TWC is available (i.e., SAFYE-CO2GSM and SAFYE-CO2SG, which is the case at the landscape scale), the SAFY-CO2 model was used to simulate the components of the C budget of CY-15 at the landscape scale (Figure 9) for more than 200 fields (black dashed square in Figure 1). Figure 10 presents the corresponding NEP (A), CEXP (B), and resulting NECB (C).
Considering all 200 fields, the simulated NEP is between −91 (sink) and 232 (source) gC·m−2·yr−1, while the amount of C exported at harvest is between 68 and 155 gC·m−2·yr−1. This results in an NECB between 64 and 300 gC·m−2·yr−1, which means that each simulated field behaves as a C source whose flux is towards the atmosphere. The topography seems to slightly affect sunflower development. This tendency can be seen in Figure 10A; plots where the C assimilation is higher (lowest NEP) are located mainly in plains and valleys. The same observation can be made for the CEXP term, meaning that the highest yields are observed in plains or in the lowest parts of hillsides. However, high yields are also observed in hilly areas, and poor yields are also observed in plains. Concerning the NECB, the tendency is less clear as the C budget is conditioned by the NEP and the CEXP term, whose signs may be opposite, with each favorably affected by good soil conditions (not short of water). Additionally, NEP can be affected by spontaneous regrowth events or weeds, increasing the CO2 net assimilation, albeit with no effect on yield. Differences in the availability of water resources for crop management can also explain the spatial variability in the NECB and its components. All plots simulated here are sources of C because the length of the sunflower cycle is short, which leaves a long fallow period that releases C from RH to the atmosphere. This observation is consistent with the results of a meta-analysis conducted for more than 15 flux crop sites in Europe by [66].

5. Discussion

5.1. SAFYE-CO2 Overall Performances

The SAFYE-CO2 model for sunflower was validated in this study. Section 4.1 presents the results achieved with this approach using in situ soil measurements to calculate SWC parameters and soil depth. The performances for simulating the DAM and yield are within the range of those of other studies. As an example, [67] compared the performances of three crop models (AquaCrop, CropSyst, and WOFOST) for simulating the aboveground biomass and yield of sunflower. Those authors showed that under rainfed conditions (similar to the conditions in the present study), the accuracy of the models ranged from 63 (WOFOST) to 204 (AquaCrop) g·m−2, while SAFYE-CO2 estimated DAM at FR-Aur with an error of 66 g·m−2. Concerning yield, SAFYE-CO2 (using in situ measurements to estimate the TWC) provided estimates with an accuracy (error of 54 g·m−2) similar to that shown by [67], who reported errors ranging from 70 (AquaCrop) to 94 g·m−2 (CropSyst). However, the performances for estimating yield with SAFYE-CO2 should be taken with caution because of the small number of comparisons with field observations (only 2 corresponding to the Auradé site in 2007 and 2016), meaning that additional validation exercises are needed.
CO2 fluxes are also robustly estimated by the SAFYE-CO2 model and are within the range of other models. The authors of [68] evaluated the CERES-EGC model against CO2 fluxes over several European flux sites, including FR-Aur. The model reproduced the net CO2 fluxes over the growing period of CY-07 with an error of 30 gC·m−2, while SAFYE-CO2IS (using in situ measurements for soil description) estimated NEP (throughout the whole cropping year) with an error of −12.1 gC·m−2. Additionally, in 2016, despite the lack of satellite images covering the period of senescence at FR-Aur, the model provided accurate estimates of NEE (mainly because the errors in GPP and RECO compensate for each other). This accuracy in estimating CO2 fluxes makes possible the realistic estimation of the NECB at the plot scale with the proposed approach. SAFYE-CO2IS was able to reproduce NECB with very good accuracy (relative errors of −4.3% and +15% for CY-07 and CY-16, respectively). SAFYE-CO2IS also performed well at estimating ETR. During the two simulated crop years at FR-Aur, the model reproduced ETR with a low error (RSME = 0.68 mm·d−1), meaning that our approach can also be used to calculate annual water budgets for sunflower plots.
At the landscape scale, coupling with a soil water module tended to lower biomass and yield estimates. The poor estimates of TWC derived from soil products (GSM and SG) led to an underestimation of the final biomass and, thus, yield. In contrast, the biomass was accurately estimated by the SAFY-CO2 model (R2 = 0.91 and RMSE = 73 g·m−2), and the estimations were within the range of those of models specific to sunflower: Using the OILCROP-SUN model, [21] estimated DAM with an error of 170 g·m−2, while using the SUNFLO model, [22] simulated DAM with an error of 114 g·m−2. In contrast, the yield was poorly estimated by SAFY-CO2 in terms of its correlation (R2 = 0.31), but the error was still within the range or better than those other studies. The error achieved by SAFY-CO2 in estimating YLD was 51 g·m−2, while errors with OILCROP-SUN and SUNFLO were 80 and 50 g·m−2, respectively (in [21,22]). The poor results achieved for yield with SAFY-CO2 can be explained in part by differences in the HI between sunflower plots. Indeed, sunflower HI is known to change with genotype, management, and environment [63,64] and depends on many factors that are not yet accounted for in the proposed approach. Accounting for varietal potentials [65] and for abiotic stress after flowering as in CropSyst [19] is an avenue for improvement with our approach.

5.2. Prerequisite to Compute Water Fluxes Over a Wide Area

Estimating water fluxes and thus water stress over croplands requires an accurate description of soil properties (θFC, θWP, and soil depth). This need has been highlighted in the present study. Soil parameters were derived from in situ measurements when available or from soil texture products for upscaling. The two existing soil products tested in this study (GSM and SG) did not meet the needs for the expected accuracy. Indeed, we demonstrated (see Section 4.2) that the use of soil products led to an underestimation of the soil TWC, which in turn led to an overestimation of the water stress and an underestimation of the biomass and yield compared to the in situ data. However, this conclusion concerns only our study area, which is not representative of the spatial scale of GSM and SG products. Indeed, those products provide, among other things, soil texture and depth at a fine resolution (SG: 250 m; GSM: 90 m) all across the globe. As such, it would be difficult to expect those products to accurately describe differences in soil properties between plots over a small area of study. Additionally, taking water fluxes into account in our approach does not improve the estimation of biomass or CO2 fluxes, even under water-stressed conditions (see Section 4.2.2). Both sunflower crop years at FR-Aur experienced water stress (especially at the end of the season), and SAFY-CO2 was still able to accurately reproduce the crop dynamics and production as well as the CO2 flux, indicating that the plant stresses (i.e., water and nitrogen) were already correctly taken into account through the remotely sensed GAI. Vegetation stress slows down (or even stops) growth, which is observable through GAI. Therefore, we believe that coupling SAFY-CO2 with a soil water module is useful for estimating ETR and water budgets or for assessing crop water-use efficiency [1] when safeguarding soil and water resources, but this coupling does not seem to improve the simulations of the annual C budget components. During the two sunflower crop years at FR-Aur, simulations were performed using SAFYE-CO2, with no feedback of the SWC on GPP computation (WtSV = 1, Equation (2)). The accuracy in estimating biomass and CO2 flux was similar to that of SAFY-CO2; thus, the accuracy was slightly better than that of simulations considering the water stress function. Concerning the estimation of water fluxes, the same accuracy was obtained with or without the consideration of effects of water stress on vegetation development. We thus argue that SAFYE-CO2 should be used in a decoupled way between vegetation growth and SWC since information about water stress is already provided by the remotely sensed GAI. However, further analysis that includes additional water flux data and the simulation of other crop species should reveal whether the water module should be used systematically in a decoupled way, allowing us to estimate water fluxes in the absence of vegetation stress.

5.3. Limitations and Potential Improvements

In addition to simulating crop dynamics and production, our approach fills the gap in assessing the C and water budget simultaneously at the plot scale over large areas. Notably, to achieve this, information on both organic fertilization and eventually irrigation at the plot scale are needed. This information is not yet available at a landscape spatial scale, which can lead, for example, to great uncertainties in the C budget calculations in areas with livestock farming. However, this drawback could be compensated by the combination of our approach with the use of data in farm management information systems (FMISs). Additionally, the representation of RH is quite simplistic in our approach and would require improvement. As shown in the present study and as discussed by [18], this simple method allows for the accurate assessment of the annual C budget, but when considering relatively long time scales (e.g., crop rotation), the effects of other variables, such as soil organic matter content, or processes, such as priming effects, will have to be taken into account. For those reasons, the SAFYE-CO2 model should be coupled with a soil organic matter model (e.g., AMG [66,67,68,69]; RothC [70]) in the future.
Another limitation of this approach concerns the availability of remote sensing images. To date, we used only optical data in this approach, which is thus limited in the case of cloudy sky conditions. Our method could even be inoperable if long periods of cloudy sky conditions occur during crop development. This constraint is partially mitigated because of the ongoing Sentinel missions that will provide high spatial- (10 m) and temporal-resolution (5 days) images all across the globe. Additionally, the assimilation of the C-band synthetic-aperture radar (SAR) from Sentinel-1 should be considered to lift this constraint, as shown in many studies [71,72,73,74]. Finally, thermal-infrared remote sensing, which allows the estimation of vegetation water stress [75], could be considered to better constrain the model.

6. Conclusions

The newly established SAFYE-CO2 model was validated at local and landscape scales for sunflower crops against a wide range of in situ measurements. The model was able to accurately reproduce biomass, yield, CO2 flux, and ETR at the FR-Aur flux site when using in situ measurements of texture to deduce soil parameters (field capacity and permanent wilting point) and soil depth. For upscaling the model to the landscape scale, the sensitivity of the model to those soil parameters estimated from existing state-of-the-art soil databases (GlobalSoilMap and SoilGrids) was tested. Biomass, over a population of fields, was correctly estimated by SAFY-CO2 but underestimated by SAFYE-CO2. This is possibly due to the underestimation of the total soil water capacity derived from the soil products, which leads to overestimations of water stress. Moreover, the low heterogeneity and spatial variability of the derived soil parameters over the study area, compared to observed values, indicates that the current accuracy of soil products (GSM and SG) is not sufficient to provide an accurate description of TWC, especially for approaches developed at the field scale. Since SAFYE-CO2 requires information about soil water capacity that is not yet available at the plot scale over wide areas with sufficient accuracy, in the case where only a C budget is needed, we suggest using SAFY-CO2 for large-scale simulation, as the assimilation of frequent HSTR GAI data enables the indirect incorporation of the integrated impact of water stress. For the same reason, we argue that crop growth in SAFYE-CO2 can be decoupled from soil water dynamics by omitting the water stress function to obtain water flux estimates as well as accurate CO2 flux estimates. However, this approach needs further validation using other flux sites to ensure the robustness of the model for simulating various soil and climatic conditions and management practices.

Author Contributions

Conceptualization, G.P., R.F. and E.C.; Data curation, P.D. and T.T.; Formal analysis, G.P., R.F., P.D., A.A.B. and E.C.; Funding acquisition, E.C.; Investigation, G.P., R.F. and A.A.B.; Methodology, G.P.; Software, A.A.B.; Supervision, G.P., R.F. and E.C.; Validation, G.P., R.F. and E.C.; Visualization, G.P. and A.A.B.; Writing—original draft, G.P; Writing—review & editing, G.P., R.F., P.D., A.A.B., T.T. and E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Climate-KIC 2020 soil carbon farming project (75000142), the Agence De l’Environnement et de la Maîtrise de l’Energie (ADEME), which financed the project “Couverts Intermédiaires pour l’atténuation du Changement Climatique” (CICC) and half of Gaétan’s dissertation research, the Centre National d’Etudes Spatiales (CNES), which financed half of Gaétan’s dissertation research, and the Agence de l’Eau Adour Garonne (AEAG), which financed the Bag’ages project.

Acknowledgments

This work was made possible through the support of the Agence De l’Environnement et de la Maîtrise de l’Energie (ADEME), which financed the project “Couverts Intermédiaires pour l’atténuation du Changement Climatique” (CICC) and half of Gaétan’s dissertation research, the Centre National d’Etudes Spatiales (CNES), which financed half of Gaétan’s dissertation research, and the Agence de l’Eau Adour Garonne (AEAG), which financed the Bag’ages project. The data acquisition at FR-Aur was funded mainly by the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique (CNRS-INSU) through the ICOS and OSR SW observatories. Facility and staff were also funded and supported by the University Toulouse III—Paul Sabatier, the CNES and IRD (Institut de Recherche pour le Développement). We are grateful to Franck Granouillac, Nicole Claverie, Baptiste Lemaire and Bartosz Zawilski for their technical support, advice, and valuable assistance in the field and with site management. We extend our special thanks to Andréoni (farmer) for accommodating our measurement devices in his field, the FR-Aur site. We are also very grateful to Jean-François Dejoux, Hervé Gibrin, Luc Champolivier, etc., for granting and facilitating our access to the 2013 campaign dataset. Data from 2014 and 2015 campaigns were collected during the CasDAR 1311 project, which was supported by the French Ministry of Agriculture.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tallec, T.; Béziat, P.; Jarosz, N.; Rivalland, V.; Ceschia, E. Crops’ water use efficiencies in temperate climate: Comparison of stand, ecosystem and agronomical approaches. Agric. For. Meteorol. 2013, 168, 69–81. [Google Scholar] [CrossRef]
  2. De Witt, T.; Brouwer, R. The simulation of photosynthetic systems. In Proceedings of the IBP/PP Technical Meeting, Trebon, Czech Republic, 14–21 September 1969; p. 24. [Google Scholar]
  3. Arkin, G.F.; Vanderlip, R.L.; Ritchie, J.T. A Dynamic Grain Sorghum Growth Model. Trans. ASAE 1976, 19, 622–626. [Google Scholar] [CrossRef]
  4. Wilkerson, G.G.; Jones, J.W.; Boote, K.J.; Ingram, K.T.; Mishoe, J.W. Modeling Soybean Growth for Crop Management. Trans. ASAE 1983, 26, 63–73. [Google Scholar] [CrossRef]
  5. Swaney, D.P.; Jones, J.W.; Boggess, W.G.; Wilkerson, G.G.; Mishoe, J.W. Real-Time Irrigation Decision Analysis Using Simulation. Trans. ASAE 1983, 26, 562–568. [Google Scholar] [CrossRef]
  6. McKinion, J.M.; Baker, D.N.; Whisler, F.D.; Lambert, J.R. Application of the GOSSYM/COMAX system to cotton crop management. Agric. Syst. 1989, 31, 55–65. [Google Scholar] [CrossRef]
  7. Bouman, B.A.M.; van Keulen, H.; van Laar, H.H.; Rabbinge, R. The “School of de Wit” crop growth simulation models: A pedigree and historical overview. Agric. Syst. 1996, 52, 171–198. [Google Scholar] [CrossRef]
  8. Van Laar, H.H. Simulation of Crop Growth for Potential and Water-Limited Production Situations: As Applied to Spring Wheat. 1992. Available online: https://library.wur.nl/WebQuery/wurpubs/fulltext/359573 (accessed on 10 September 2020).
  9. Yin, X.; van Laar, H.H. Crop Systems Dynamics: An Ecophysiological Model of Genotype-by-Environment Interactions (GECROS); Wageningen Academic Publishers: Wageningen, The Netherlands, 2005. [Google Scholar]
  10. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop—The FAO Crop Model to Simulate Yield Response to Water: I. Concepts and Underlying Principles. Agron. J. 2009, 101, 426. [Google Scholar] [CrossRef] [Green Version]
  11. Constantin, J.; Willaume, M.; Murgue, C.; Lacroix, B.; Therond, O. The soil-crop models STICS and AqYield predict yield and soil water content for irrigated crops equally well with limited data. Agric. For. Meteorol. 2015, 206, 55–68. [Google Scholar] [CrossRef]
  12. Monteith, J.L.; Moss, C.J. Climate and the Efficiency of Crop Production in Britain and Discussion. Philos. Trans. R. Soc. B Biol. Sci. 1977, 281, 277–294. [Google Scholar] [CrossRef]
  13. Jones, C.A.; Kiniry, J.R.; Dyke, P.T. CERES-Maize: A Simulation Model of Maize Growth and Development; Texas A&M University Press: College Station, TX, USA, 1986. [Google Scholar]
  14. Williams, J.R.; Jones, C.A.; Kiniry, J.R.; Spanel, D.A. The EPIC Crop Growth Model. Trans. ASAE 1989, 32, 497–511. [Google Scholar] [CrossRef]
  15. Brisson, N.; Gary, C.; Justes, E.; Roche, R.; Mary, B.; Ripoche, D.; Zimmer, D.; Sierra, J.; Bertuzzi, P.; Bussière, F.; et al. An overview of the crop model stics. Eur. J. Agron. 2003, 18, 309–332. [Google Scholar] [CrossRef]
  16. Duchemin, B.; Maisongrande, P.; Boulet, G.; Benhadj, I. A simple algorithm for yield estimates: Evaluation for semi-arid irrigated winter wheat monitored with green leaf area index. Environ. Model. Softw. 2008, 23, 876–892. [Google Scholar] [CrossRef] [Green Version]
  17. Duchemin, B.; Fieuzal, R.; Rivera, M.A.; Ezzahar, J.; Jarlan, L.; Rodriguez, J.C.; Hagolle, O.; Watts, C. Impact of Sowing Date on Yield and Water Use Efficiency of Wheat Analyzed through Spatial Modeling and FORMOSAT-2 Images. Remote Sens. 2015, 7, 5951–5979. [Google Scholar] [CrossRef] [Green Version]
  18. Pique, G.; Fieuzal, R.; Al Bitar, A.; Veloso, A.; Tallec, T.; Brut, A.; Ferlicoq, M.; Zawilski, B.; Gibrin, H.; Ceschia, E.; et al. Estimation of daily CO2 fluxes and of the components of the carbon budget for winter wheat by the assimilation of Sentinel 2-like remote sensing data into a crop model. Geoderma 2020, 376, 114428. [Google Scholar] [CrossRef]
  19. Stöckle, C.O.; Donatelli, M.; Nelson, R. CropSyst, a cropping systems simulation model. Eur. J. Agron. 2003, 18, 289–307. [Google Scholar] [CrossRef]
  20. Diepen, C.A.; Wolf, J.; Keulen, H.; Rappoldt, C. WOFOST: A simulation model of crop production. Soil Use Manag. 1989, 5, 16–24. [Google Scholar] [CrossRef]
  21. Villalobos, F.J.; Hall, A.J.; Ritchie, J.T.; Orgaz, F. OILCROP-SUN: A Development, Growth, and Yield Model of the Sunflower Crop. Agron. J. 1996, 88, 403. [Google Scholar] [CrossRef]
  22. Casadebaig, P.; Guilioni, L.; Lecoeur, J.; Christophe, A.; Champolivier, L.; Debaeke, P. SUNFLO, a model to simulate genotype-specific performance of the sunflower crop in contrasting environments. Agric. For. Meteorol. 2011, 151, 163–178. [Google Scholar] [CrossRef]
  23. Battude, M.; Al Bitar, A.; Morin, D.; Cros, J.; Huc, M.; Sicre, C.M.; Le Dantec, V.; Demarez, V. Estimating maize biomass and yield over large areas using high spatial and temporal resolution Sentinel-2 like remote sensing data. Remote Sens. Environ. 2016, 184, 668–681. [Google Scholar] [CrossRef]
  24. Sibley, A.M.; Grassini, P.; Thomas, N.E.; Cassman, K.G.; Lobell, D.B. Testing Remote Sensing Approaches for Assessing Yield Variability among Maize Fields. Agron. J. 2014, 106, 24. [Google Scholar] [CrossRef]
  25. Sus, O.; Heuer, M.W.; Meyers, T.P.; Williams, M. A data assimilation framework for constraining upscaled cropland carbon flux seasonality and biometry with MODIS. Biogeosciences 2013, 10, 2451–2466. [Google Scholar] [CrossRef] [Green Version]
  26. Wu, S.; Huang, J.; Liu, X.; Fan, J.; Ma, G.; Zou, J. Assimilating MODIS-LAI into Crop Growth Model with EnKF to Predict Regional Crop Yield. In Computer and Computing Technologies in Agriculture V; Li, D., Chen, Y., Eds.; Springer: Berlin/Heidelberg, Germany, 2012; Volume 370, pp. 410–418. [Google Scholar]
  27. Running, S.W.; Thornton, P.E.; Nemani, R.; Glassy, J.M. Global Terrestrial Gross and Net Primary Productivity from the Earth Observing System. In Methods in Ecosystem Science; Sala, O.E., Jackson, R.B., Mooney, H.A., Howarth, R.W., Eds.; Springer: New York, NY, USA, 2000; pp. 44–57. [Google Scholar]
  28. Cicek, H.; Sunohara, M.; Wilkes, G.; McNairn, H.; Pick, F.; Topp, E.; Lapen, D.R. Using vegetation indices from satellite remote sensing to assess corn and soybean response to controlled tile drainage. Agric. Water Manag. 2010, 98, 261–270. [Google Scholar] [CrossRef]
  29. Panda, S.S.; Ames, D.P.; Panigrahi, S. Application of Vegetation Indices for Agricultural Crop Yield Prediction Using Neural Network Techniques. Remote Sens. 2010, 2, 673–696. [Google Scholar] [CrossRef] [Green Version]
  30. Shanahan, J.F.; Schepers, J.S.; Francis, D.D.; Varvel, G.E.; Wilhelm, W.W.; Tringe, J.M.; Schlemmer, M.R.; Major, D.J. Use of Remote-Sensing Imagery to Estimate Corn Grain Yield. Agron. J. 2001, 93, 583–589. [Google Scholar] [CrossRef] [Green Version]
  31. Claverie, M.; Demarez, V.; Duchemin, B.; Hagolle, O.; Ducrot, D.; Marais-Sicre, C.; Dejoux, J.-F.; Huc, M.; Keravec, P.; Fieuzal, R.; et al. Maize and sunflower biomass estimation in southwest France using high spatial and temporal resolution remote sensing data. Remote Sens. Environ. 2012, 124, 844–857. [Google Scholar] [CrossRef]
  32. Ferrant, S.; Gascoin, S.; Veloso, A.; Salmon-Monviola, J.; Claverie, M.; Rivalland, V.; Dedieu, G.; Demarez, V.; Ceschia, E.; Durand, P.; et al. Agro-hydrology and multi-temporal high-resolution remote sensing: Toward an explicit spatial processes calibration. Hydrol. Earth Syst. Sci. 2014, 18, 5219–5237. [Google Scholar] [CrossRef] [Green Version]
  33. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; Volume 300, p. D05109. [Google Scholar]
  34. Battude, M.; Al Bitar, A.; Brut, A.; Tallec, T.; Huc, M.; Cros, J.; Weber, J.-J.; Lhuissier, L.; Simonneaux, V.; Demarez, V. Modeling water needs and total irrigation depths of maize crop in the south west of France using high spatial and temporal resolution satellite imagery. Agric. Water Manag. 2017, 189, 123–136. [Google Scholar] [CrossRef]
  35. Champolivier, L.; Debaeke, P.; Dejoux, J.F.; Dizien, C.; Micheneau, A.; Colombet, C.; Gibrin, H.; Pontet, C.; Al Bitar, A.; Ansart, A.; et al. Outils, références et méthodes pour la construction d’un simulateur pour la prévision du rendement et de la qualité du tournesol à l’échelle territoriale mobilisant la télédétection satellitaire. Innov. Agron. 2019, 71, 15–34. [Google Scholar] [CrossRef]
  36. Champolivier, L.; Debaeke, P.; Merrien, A. Pourquoi irriguer le tournesol, une culture réputée tolérante à la sécheresse? Innov. Agron. 2011, 14, 151–164. [Google Scholar]
  37. Béziat, P.; Ceschia, E.; Dedieu, G. Carbon balance of a three crop succession over two cropland sites in South West France. Agric. For. Meteorol. 2009, 149, 1628–1645. [Google Scholar] [CrossRef] [Green Version]
  38. Durand, Y.; Brun, E.; Merindol, L.; Guyomarc’h, G.; Lesaffre, B.; Martin, E. A meteorological estimation of relevant parameters for snow models. Ann. Glaciol. 1993, 18, 65–71. [Google Scholar] [CrossRef] [Green Version]
  39. Aubinet, M.; Vesala, T.; Papale, D. (Eds.) Eddy Covariance: A Practical Guide to Measurement and Data Analysis; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  40. Baldocchi, D.D. Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: Past, present and future. Glob. Chang. Biol. 2003, 9, 479–492. [Google Scholar] [CrossRef] [Green Version]
  41. Moncrieff, J.B.; Massheder, J.M.; De Bruin, H.; Elbers, J.; Friborg, T.; Heusinkveld, B.; Kabat, P.; Scott, S.; Soegaard, H.; Verhoef, A. A system to measure surface fluxes of momentum, sensible heat, water vapour and carbon dioxide. J. Hydrol. 1997, 188–189, 589–611. [Google Scholar] [CrossRef]
  42. Reichstein, M.; Falge, E.; Baldocchi, D.; Papale, D.; Aubinet, M.; Berbigier, P.; Bernhofer, C.; Buchmann, N.; Gilmanov, T.; Grünwald, T.; et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: Review and improved algorithm. Glob. Chang. Biol. 2005, 11, 1424–1439. [Google Scholar] [CrossRef]
  43. Dobarco, M.R.; Bourennane, H.; Arrouays, D.; Saby, N.P.A.; Cousin, I.; Martin, M.P. Uncertainty assessment of GlobalSoilMap soil available water capacity products: A French case study. Geoderma 2019, 344, 14–30. [Google Scholar] [CrossRef]
  44. Hengl, T.; Mendes de Jesus, J.; Heuvelink, G.B.; Ruiperez Gonzalez, M.; Kilibarda, M.; Blagotić, A.; Shangguan, W.; Wright, M.N.; Geng, X.; Guevara, M.A.; et al. SoilGrids250m: Global gridded soil information based on machine learning. PLoS ONE 2017, 12, e0169748. [Google Scholar] [CrossRef] [Green Version]
  45. Minasny, B.; McBratney, A.B. Methodologies for Global Soil Mapping. In Digital Soil Mapping; Boettinger, J.L., Howell, D.W., Moore, A.C., Hartemink, A.E., Kienast-Brown, S., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 429–436. [Google Scholar]
  46. Dobarco, M.R.; Cousin, I.; le Bas, C.; Martin, M.P. Pedotransfer functions for predicting available water capacity in French soils, their applicability domain and associated uncertainty. Geoderma 2019, 336, 81–95. [Google Scholar] [CrossRef]
  47. Liu, J.; Pattey, E.; Miller, J.R.; McNairn, H.; Smith, A.; Hu, B. Estimating crop stresses, aboveground dry biomass and yield of corn using multi-temporal optical data combined with a radiation use efficiency model. Remote Sens. Environ. 2010, 114, 1167–1177. [Google Scholar] [CrossRef]
  48. Lobell, D.B.; Asner, G.P.; Ortiz-Monasterio, J.I.; Benning, T.L. Remote sensing of regional crop production in the Yaqui Valley, Mexico: Estimates and uncertainties. Agric. Ecosyst. Environ. 2003, 94, 205–220. [Google Scholar] [CrossRef] [Green Version]
  49. Hagolle, O.; Sylvander, S.; Huc, M.; Claverie, M.; Clesse, D.; Dechoz, C.; Lonjou, V.; Poulain, V. SPOT-4 (Take 5): Simulation of Sentinel-2 Time Series on 45 Large Sites. Remote Sens. 2015, 7, 12242–12264. [Google Scholar] [CrossRef] [Green Version]
  50. Lafrance, B.; Lenot, X.; Ruffel, C.; Cao, P.; Rabaute, T. Outils de prétraitements des images optiques KALIDEOS. Rev. Fr. Photogramm. Télédétect. 2012, 197, 10–16. [Google Scholar]
  51. Baret, F.; Hagolle, O.; Geiger, B.; Bicheron, P.; Miras, B.; Huc, M.; Berthelot, B.; Weiss, M.; Samain, O.; Roujean, J.L.; et al. LAI, fAPAR and fCover CYCLOPES global products derived from VEGETATION. Remote Sens. Environ. 2007, 110, 275–286. [Google Scholar] [CrossRef] [Green Version]
  52. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarco-Tejada, P.J.; Asner, G.P.; François, C.; Ustin, S.L. PROSPECT+SAIL models: A review of use for vegetation characterization. Remote Sens. Environ. 2009, 113, S56–S66. [Google Scholar] [CrossRef]
  53. Demarez, V.; Duthoit, S.; Baret, F.; Weiss, M.; Dedieu, G. Estimation of leaf area and clumping indexes of crops with hemispherical photographs. Agric. For. Meteorol. 2008, 148, 644–655. [Google Scholar] [CrossRef] [Green Version]
  54. McCree, K.J. Equations for the Rate of Dark Respiration of White Clover and Grain Sorghum, as Functions of Dry Weight, Photosynthetic Rate, and Temperature1. Crop Sci. 1974, 14, 509. [Google Scholar] [CrossRef]
  55. Van’t Hoff, J.H. Lectures on Theoretical and Physical Chemistry. Part I. Chemical Dynamics; Edward Arnold: London, UK, 1898; pp. 224–229. [Google Scholar]
  56. Amthor, J.S. Respiration and Crop Productivity; Springer: New York, NY, USA, 1989. [Google Scholar]
  57. Choudhury, B.J. A sensitivity analysis of the radiation use efficiency for gross photosynthesis and net carbon accumulation by wheat. Agric. For. Meteorol. 2000, 101, 217–234. [Google Scholar] [CrossRef]
  58. Baret, F.; Olioso, A.; Luciani, J.L. Root biomass fraction as a function of growth degree days in wheat. Plant Soil 1992, 140, 137–144. [Google Scholar] [CrossRef]
  59. Ma, T.; Zeng, W.; Li, Q.; Yang, X.; Wu, J.; Huang, J. Shoot and Root Biomass Allocation of Sunflower Varying with Soil Salinity and Nitrogen Applications. Agron. J. 2017, 109, 2545. [Google Scholar] [CrossRef]
  60. Soares, J.V.; Bernard, R.; Taconet, O.; Vidal-Madjar, D.; Weill, A. Estimation of bare soil evaporation from airborne measurements. J. Hydrol. 1988, 99, 281–296. [Google Scholar] [CrossRef]
  61. Lagarias, J.C.; Reeds, J.A.; Wright, M.H.; Wright, P.E. Convergence Properties of the Nelder–Mead Simplex Method in Low Dimensions. SIAM J. Optim. 1998, 9, 112–147. [Google Scholar] [CrossRef] [Green Version]
  62. Delogu, E. Modélisation de la Respiration du Sol dans les Agrosystèmes; Université Paul Sabatier: Toulouse, France, 2017. [Google Scholar]
  63. Bange, M.P.; Hammer, G.L.; Rickert, K.G. Temperature and Sowing Date Affect the Linear Increase of Sunflower Harvest Index. Agron. J. 1998, 90, 324. [Google Scholar] [CrossRef]
  64. Gholinezhad, E.; Aynaband, A.; Hassanzade, A.; Noormohamadi, G.; Bernousi, I. Study of the Effect of Drought Stress on Yield, Yield Components and Harvest Index of Sunflower Hybrid Iroflor at Different Levels of Nitrogen and Plant Population. Not. Bot. Horti Agrobot. 2009, 37, 85–94. [Google Scholar]
  65. Debaeke, P.; Casadebaig, P.; Haquin, B.; Mestries, E.; Palleau, J.-P.; Salvi, F. Simulation de la réponse variétale du tournesol à l’environnement à l’aide du modèle SUNFLO. Ol. Corps Gras Lipides 2010, 17, 143–151. [Google Scholar] [CrossRef]
  66. Ceschia, E.; Béziat, P.; Dejoux, J.F.; Aubinet, M.; Bernhofer, C.; Bodson, B.; Buchmann, N.; Carrara, A.; Cellier, P.; Elbers, J.A.; et al. Management effects on net ecosystem carbon and GHG budgets at European crop sites. Agric. Ecosyst. Environ. 2010, 139, 363–383. [Google Scholar] [CrossRef] [Green Version]
  67. Todorovic, M.; Albrizio, R.; Zivotic, L.; Saab, M.-T.A.; Stöckle, C.; Steduto, P. Assessment of AquaCrop, CropSyst, and WOFOST Models in the Simulation of Sunflower Growth under Different Water Regimes. Agron. J. 2009, 101, 509. [Google Scholar] [CrossRef]
  68. Lehuger, S.; Gabrielle, B.; Cellier, P.; Loubet, B.; Roche, R.; Béziat, P.; Ceschia, E.; Wattenbach, M. Predicting the net carbon exchanges of crop rotations in Europe with an agro-ecosystem model. Agric. Ecosyst. Environ. 2010, 139, 384–395. [Google Scholar] [CrossRef] [Green Version]
  69. Saffih-Hdadi, K.; Mary, B. Modeling consequences of straw residues export on soil organic carbon. Soil Biol. Biochem. 2008, 40, 594–607. [Google Scholar] [CrossRef]
  70. Coleman, K.; Jenkinson, D.S. RothC-26.3—A Model for the turnover of carbon in soil. In Evaluation of Soil Organic Matter Models; Powlson, D.S., Smith, P., Smith, J.U., Eds.; Springer: Berlin/Heidelberg, Germany, 1996; pp. 237–246. [Google Scholar]
  71. Ameline, M.; Fieuzal, R.; Betbeder, J.; Berthoumieu, J.-F.; Baup, F. Estimation of Corn Yield by Assimilating SAR and Optical Time Series into a Simplified Agro-Meteorological Model: From Diagnostic to Forecast. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 4747–4760. [Google Scholar] [CrossRef]
  72. Baup, F.; Ameline, M.; Fieuzal, R.; Frappart, F.; Corgne, S.; Berthoumieu, J.-F. Temporal Evolution of Corn Mass Production Based on Agro-Meteorological Modelling Controlled by Satellite Optical and SAR Images. Remote Sens. 2019, 11, 1978. [Google Scholar] [CrossRef] [Green Version]
  73. Fieuzal, R.; Sicre, C.M.; Baup, F. Estimation of corn yield using multi-temporal optical and radar satellite data and artificial neural networks. Int. J. Appl. Earth Obs. Geoinf. 2017, 57, 14–23. [Google Scholar] [CrossRef]
  74. Revill, A.; Sus, O.; Barrett, B.; Williams, M. Carbon cycling of European croplands: A framework for the assimilation of optical and microwave Earth observation data. Remote Sens. Environ. 2013, 137, 84–93. [Google Scholar] [CrossRef] [Green Version]
  75. Labbé, S.; Lebourgeois, V.; Jolivot, A.; Marti, R. Thermal infra-red remote sensing for water stress estimation in agriculture. Options Méditerr. Sér. B Etudes Rech. 2012, 67, 175–184. [Google Scholar]
Figure 1. Top-left: Fr-Aur field and the location of the measurement station; Top-right: Location of the study area in France; Bottom: Location of the fields in which biomass and yield data were collected in 2013, 2014, and 2015, as well as the selected area for the model’s upscaling (see Section 4.3).
Figure 1. Top-left: Fr-Aur field and the location of the measurement station; Top-right: Location of the study area in France; Bottom: Location of the fields in which biomass and yield data were collected in 2013, 2014, and 2015, as well as the selected area for the model’s upscaling (see Section 4.3).
Remotesensing 12 02967 g001
Figure 2. Soil type distribution derived from GlobalSoilMap (left) and SoilGrids (right) products based on the United States Department of Agriculture (USDA) classification system. The red, black, and blue circles respectively correspond to fields in the 2013, 2014, and 2015 seasons.
Figure 2. Soil type distribution derived from GlobalSoilMap (left) and SoilGrids (right) products based on the United States Department of Agriculture (USDA) classification system. The red, black, and blue circles respectively correspond to fields in the 2013, 2014, and 2015 seasons.
Remotesensing 12 02967 g002
Figure 3. Timeline of images for each year of study and each sensor.
Figure 3. Timeline of images for each year of study and each sensor.
Remotesensing 12 02967 g003
Figure 4. Summary of the validation strategy in terms of upscaling the SAFY-CO2 and SAFYE-CO2 models. IS is in situ, GSM is GlobalSoilMap, and SG is SoilGrids.
Figure 4. Summary of the validation strategy in terms of upscaling the SAFY-CO2 and SAFYE-CO2 models. IS is in situ, GSM is GlobalSoilMap, and SG is SoilGrids.
Remotesensing 12 02967 g004
Figure 5. Scatter plots of the observed and simulated GAI (A), DAM (B), GPP (C), RECO (D), NEE (E), and ETR (F) for the two crop years of sunflower at FR-Aur and mean statistic over the 2 crop years.
Figure 5. Scatter plots of the observed and simulated GAI (A), DAM (B), GPP (C), RECO (D), NEE (E), and ETR (F) for the two crop years of sunflower at FR-Aur and mean statistic over the 2 crop years.
Remotesensing 12 02967 g005
Figure 6. RMSE and R2 values for the main model’s output for CY-07 (a) and CY-16 (b) at FR-Aur via SAFY-CO2 (orange bars) or SAFYE-CO2 run with in situ soil data (SAFYE-CO2IS, green bars), GlobalSoilMap data (SAFYE-CO2GSM, red bars), and SoilGrids data (SAFYE-CO2SG, blue bars).
Figure 6. RMSE and R2 values for the main model’s output for CY-07 (a) and CY-16 (b) at FR-Aur via SAFY-CO2 (orange bars) or SAFYE-CO2 run with in situ soil data (SAFYE-CO2IS, green bars), GlobalSoilMap data (SAFYE-CO2GSM, red bars), and SoilGrids data (SAFYE-CO2SG, blue bars).
Remotesensing 12 02967 g006
Figure 7. Components of the net ecosystem carbon budget for CY-07 (a) and CY-16 (b) at FR-Aur for the SAFYE-CO2IS (green bars), SAFYE-CO2GSM (red bars), SAFYE-CO2SG (blue bars), and SAFY-CO2 (orange bars) models compared to the observed value (black bars).
Figure 7. Components of the net ecosystem carbon budget for CY-07 (a) and CY-16 (b) at FR-Aur for the SAFYE-CO2IS (green bars), SAFYE-CO2GSM (red bars), SAFYE-CO2SG (blue bars), and SAFY-CO2 (orange bars) models compared to the observed value (black bars).
Remotesensing 12 02967 g007
Figure 8. Scatter plots of simulated and observed DAM (left) and YLD (right) obtained with the SAFYE-CO2GSM (A,B), SAFYE-CO2SG (C,D), and SAFY-CO2 (E,F) model.
Figure 8. Scatter plots of simulated and observed DAM (left) and YLD (right) obtained with the SAFYE-CO2GSM (A,B), SAFYE-CO2SG (C,D), and SAFY-CO2 (E,F) model.
Remotesensing 12 02967 g008
Figure 9. TWC of the fields simulated in this study derived from GSM and SG and of 45 fields within the study area estimated with in situ (IS) measurements.
Figure 9. TWC of the fields simulated in this study derived from GSM and SG and of 45 fields within the study area estimated with in situ (IS) measurements.
Remotesensing 12 02967 g009
Figure 10. (AC) The annual carbon budget (NECB) and its component, NEP and Cexp, simulated for 200 fields with SAFY-CO2 inside the study area for CY-15. Digital elevation model from Shuttle Radar Topography Mission is displayed in the background.
Figure 10. (AC) The annual carbon budget (NECB) and its component, NEP and Cexp, simulated for 200 fields with SAFY-CO2 inside the study area for CY-15. Digital elevation model from Shuttle Radar Topography Mission is displayed in the background.
Remotesensing 12 02967 g010
Table 1. Sowing and harvest dates of sunflower cultivated at FR-Aur and climatic variables averaged from the 15th of April to the 15th of September for all crop years. Bold values are averages over the 2007–2016 period.
Table 1. Sowing and harvest dates of sunflower cultivated at FR-Aur and climatic variables averaged from the 15th of April to the 15th of September for all crop years. Bold values are averages over the 2007–2016 period.
Cropping
Year
SowingHarvestMean
Temperature [°C]
Mean Radiation
[MJ·m−2]
Cumulative
Precipitation [mm]
200710th of April20th of September18.019.3339
2013--17.320.1364
2014--17.619.8178
2015--19.020.2275
201620th of April23th of September18.120.8294
Average 2007–2016--18.220.3267
Table 2. Dataset of dry aboveground mass (DAM) and yield as well as their measurement scale.
Table 2. Dataset of dry aboveground mass (DAM) and yield as well as their measurement scale.
YearNumber of Fields MonitoredDry Aboveground MassFinal YieldWorkable Final YieldCO2 and H2O FluxesScale
2007FR-Aur711YesField/ESU *
201393699NoESU
201470-3621NoField
2015116-7949NoField
2016FR-Aur311YesESU
* Elementary Sampling Units.
Table 3. Soil characteristics derived from in situ measurements and GlobalSoilMap and SoilGrids data at FR-Aur corresponding to 3 layers simulated by the model.
Table 3. Soil characteristics derived from in situ measurements and GlobalSoilMap and SoilGrids data at FR-Aur corresponding to 3 layers simulated by the model.
Field Capacity
[m3·m−3]
Wilting Point
[m3·m−3]
Available Water Capacity [m3·m3]Soil Depth [cm]
Source of DataSurface LayerRoot LayerDeep LayerSurface LayerRoot LayerDeep LayerSurface LayerRoot LayerDeep Layer-
In situ0.310.390.390.130.190.190.180.20.2(100–200)
GSM0.290.30.290.180.170.160.110.130.13159
SoilGrids0.340.330.330.250.210.210.090.120.12188
Table 4. Main characteristics of the satellite sensors used in this study and the bands used to derive the green area index (GAI).
Table 4. Main characteristics of the satellite sensors used in this study and the bands used to derive the green area index (GAI).
SPOTFormosat-2DEIMOS-1Landsat-8Sentinel-2A
Spatial resolution20 m (SPOT2-4) 10 m (SPOT5)8 m22 m30 m10 m
Revisit frequencyBetween 1 and 3 days1 dayBetween 3 and 4 days16 days5 days
Footprint60 km × 60 km24 km × 24 kmUp to 625 km170 km × 185 km290 km
Wavelength (μm)Green: [0.5–0.59]; Red: [0.61–0.68]; NIR: [0.79–0.89]; SWIR: [1.53–1.75] Green: [0.52–0.6]; Red: [0.63–0.69]; NIR: [0.79–0.9]Green: [0.52–0.6]; Red: [0.63–0.69]; NIR: [0.77–0.9]Green: [0.53–0.59]; Red: [0.64–0.67]; NIR: [0.85–0.88]; SWIR1: [1.57–1.65]; SWIR2: [2.11–2.29] Green: [0.54–0.58];
Red: [0.65–0.68];
VRE1: [0.7–0.71];
VRE2: [0.73–0.75];
VRE3: [0.77–0.79];
NIR: [0.72–0.95];
Narrow NIR [0.85–0.87];
SWIR1: [1.57–1.66];
SWIR2: [2.11–2.29]
Viewing angleConstantConstantVariableConstantConstant
Table 5. List of the SAFYE-CO2 model parameters for sunflower, values, or a priori intervals of the parameters, method for setting them, and relevant sources.
Table 5. List of the SAFYE-CO2 model parameters for sunflower, values, or a priori intervals of the parameters, method for setting them, and relevant sources.
ModelNameNotationUnitValueMethodSource
SAFY modelClimatic efficiencyεC-0.48Literature Varlet grancher 1982
Light interception coefficientKEXT-0.8Literature Flénet 96/Steer 96
Minimal temperature for growthTMIN°C8Literature STICS
Maximal temperature for growthTMAX°C42Literature STICS
Optimal temperature for growthTOPT°C25Literature STICS
Polynomial degree of thermal stress functionβ-2Literature Duchemin 2008
Harvest indexHI-0.29Literature Andrade 1995
Specific leaf areaSLAm2·g−1(10−3–4.10−2)Calibration-
Partition-to-leaf function parameter APlA-(0.01–0.5)Calibration-
Partition-to-leaf function parameter BPlB-(10−4–10−2)Calibration-
Sum of temperature for senescenceSENA°C(200–1800)Calibration-
Rate of senescenceSENB°C.d−1(102–3.104)Calibration-
Day of plant emergenceEmergday(100–200)Calibration-
Water moduleDepth of evaporative layerDEVAPm0.1Literature FAO-56 Allen 1998
Maximum depth of the root systemDROOTSm1.4Literature Stone 2001
Soil depthDm-Soil Data-
Humidity at field capacityθFCm3·m−3-Soil Data-
Humidity at the permanent wilting pointθWPm3·m−3-Soil Data-
Root mean growth rateVRm·°C−10.0025Literature STICS
Evaporation coefficient AEvAmm−10.09MeasurementsIn situ ETR measurements
Evaporation coefficient BEvB-1.09MeasurementsIn situ ETR measurements
Evaporation coefficient CKBS-1.2MeasurementsIn situ ETR measurements
Maximal transpiration coefficientKcbMAX-1Literature FAO-56 Allen 98
Critical humidity parameterDfe-0.41Literature Casadebaig 2008
Maximal vegetation fraction coverKCOV-0.92MeasurementsICOS site data
Fraction cover exponential coefficientECOV-0.62MeasurementsICOS site data
Transpiration reduction coefficientKTRP1.01MeasureMeasurementsIn situ ETR measurements
Carbon moduleEffective light-use efficiency parameter AELUEAgC·MJ−1[0.75–0.95]Calibration-
Effective light-use efficiency parameter BELUEB-1.14MeasurementsIn-situ measurments
Q10 maintenance respiration parameterQ10-2Literature Amthor 2000
Reference maintenance respiration at 10 °CR10gC·gDM−10.0036Literature Szaniawski 1982
Growth respiration conversion efficiency parameterYG-0.76Literature Szaniawski 1982
Root fraction parametersf0−f−c-0.3/0.1/1.5Literature Baret 1992/Ma 2017
Carbon content coefficientCVEGgC·gVEG−10.46Literature Béziat 2009
Reference heterotrophic respiration at 0 °CRhREFgC·m−2·d−10.34Literature Suleau 2011
Q10 heterotrophic respiration parameterQ10_H-2.3Literature Suleau 2011
Conversion factor of Tair in Tsoilt-1.07MeasurementsICOS site data

Share and Cite

MDPI and ACS Style

Pique, G.; Fieuzal, R.; Debaeke, P.; Al Bitar, A.; Tallec, T.; Ceschia, E. Combining High-Resolution Remote Sensing Products with a Crop Model to Estimate Carbon and Water Budget Components: Application to Sunflower. Remote Sens. 2020, 12, 2967. https://doi.org/10.3390/rs12182967

AMA Style

Pique G, Fieuzal R, Debaeke P, Al Bitar A, Tallec T, Ceschia E. Combining High-Resolution Remote Sensing Products with a Crop Model to Estimate Carbon and Water Budget Components: Application to Sunflower. Remote Sensing. 2020; 12(18):2967. https://doi.org/10.3390/rs12182967

Chicago/Turabian Style

Pique, Gaétan, Rémy Fieuzal, Philippe Debaeke, Ahmad Al Bitar, Tiphaine Tallec, and Eric Ceschia. 2020. "Combining High-Resolution Remote Sensing Products with a Crop Model to Estimate Carbon and Water Budget Components: Application to Sunflower" Remote Sensing 12, no. 18: 2967. https://doi.org/10.3390/rs12182967

APA Style

Pique, G., Fieuzal, R., Debaeke, P., Al Bitar, A., Tallec, T., & Ceschia, E. (2020). Combining High-Resolution Remote Sensing Products with a Crop Model to Estimate Carbon and Water Budget Components: Application to Sunflower. Remote Sensing, 12(18), 2967. https://doi.org/10.3390/rs12182967

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop