# Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe

Céline Schoving, Claudio Osvaldo Stockle, Céline Colombet, Luc

Champolivier, Philippe Debaeke, Pierre Maury

► To cite this version:

Céline Schoving, Claudio Osvaldo Stockle, Céline Colombet, Luc Champolivier, Philippe Debaeke, et al.. Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe. Frontiers in Plant Science, 2020, 10, pp.1755. 10.3389/fpls.2019.01755. [hal-02619374](https://hal.inrae.fr/hal-02619374)

HAL Id: hal-02619374

https://hal.inrae.fr/hal-02619374

Submitted on 5 Feb 2024

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

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

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

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

Check for
updates

# Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe

Céline Schoving¹*, Claudio Osvaldo Stöckle², Céline Colombet¹, Luc Champolivier³, Philippe Debaeke¹ and Pierre Maury⁴

**OPEN ACCESS**

**Edited by:**
Luis A. N. Aguirrezabal,
National University of Mar del Plata,
Argentina

**Reviewed by:**
Ralf Uptmoor,
University of Rostock, Germany
Alan Severini,
Instituto Nacional de Tecnología
Agropecuaria (Argentina), Argentina

*Correspondence:
Céline Schoving
celine.schoving@inrae.fr

**Specialty section:**
This article was submitted to
*Crop and Product Physiology,*
a section of the journal
*Frontiers in Plant Science*

**Received:** 24 May 2019
**Accepted:** 16 December 2019
**Published:** 29 January 2020

**Citation:**
Schoving C, Stöckle CO, Colombet C, Champolivier L, Debaeke P and Maury P (2020) Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe.
Front. Plant Sci. 10:1755.
doi: 10.3389/fpls.2019.01755

¹ Université de Toulouse, INRAE, UMR AGIR, Castanet-Tolosan, France, ² Biological Systems Engineering, Washington State University, Pullman, WA, United States, ³ Terres Inova–Institut technique des oléagineux, des protéagineux et du chanvre, Paris, France, ⁴ Université de Toulouse, INRAE, INP-ENSAT Toulouse, UMR AGIR, Castanet-Tolosan, France

Developing new cropping strategies (very early sowing, crop expansion at higher latitudes, double cropping) to improve soybean production in Europe under climate change needs a good prediction of phenology under different temperature and photoperiod conditions. For that purpose, a simple phenology algorithm (SPA) was developed and parameterized for 10 contrasting soybean cultivars (maturity group 000 to II). Two experiments were carried out at INRA Toulouse (France) for parameterization: 1) Phenological monitoring of plants in pots on an outdoor platform with 6 planting dates. 2) Response of seed germination to temperature in controlled conditions. Multi-location field trials including 5 sites, 4 years, 2 sowing dates, and 10 cultivars were used to evaluate the SPA phenology predictions. Mean cardinal temperatures (minimum, optimum, and maximum) for germination were ca. 2, 30, and 40°C, respectively with significant differences among cultivars. The photoperiod sensitivity coefficient varied among cultivars when fixing Popt and Pct, optimal and critical photoperiods respectively, by maturity group. The parameterized algorithm showed an RMSE of less than 6 days for the prediction of crop cycle duration (i.e. cotyledons stage to physiological maturity) in the field trials including 75 data points. Flowering (R1 stage), and beginning of grain filling (R5 stage) dates were satisfactorily predicted with RMSEs of 8.2 and 9.4 days respectively. Because SPA can be also parameterized using data from field experiments, it can be useful as a plant selection tool across environments. The algorithm can be readily applied to species other than soybean, and its incorporation into cropping systems models would enhance the assessment of the performance of crop cultivars under climate change scenarios.

**Keywords:** planting date, glycine max, varietal phenotyping, photoperiod, temperature
<figure><img src="image_2.png" /></figure>

Check for
updates

# Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe

Céline Schoving¹*, Claudio Osvaldo Stöckle², Céline Colombet¹, Luc Champolivier³, Philippe Debaeke¹ and Pierre Maury⁴

**OPEN ACCESS**

**Edited by:**
Luis A. N. Aguirrezabal,
National University of Mar del Plata,
Argentina

**Reviewed by:**
Ralf Uptmoor,
University of Rostock, Germany
Alan Severini,
Instituto Nacional de Tecnología
Agropecuaria (Argentina), Argentina

*Correspondence:
Céline Schoving
celine.schoving@inrae.fr

**Specialty section:**
This article was submitted to
*Crop and Product Physiology,*
a section of the journal
*Frontiers in Plant Science*

**Received:** 24 May 2019
**Accepted:** 16 December 2019
**Published:** 29 January 2020

**Citation:**
Schoving C, Stöckle CO, Colombet C, Champolivier L, Debaeke P and Maury P (2020) Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe.
Front. Plant Sci. 10:1755.
doi: 10.3389/fpls.2019.01755

¹ Université de Toulouse, INRAE, UMR AGIR, Castanet-Tolosan, France, ² Biological Systems Engineering, Washington State University, Pullman, WA, United States, ³ Terres Inova–Institut technique des oléagineux, des protéagineux et du chanvre, Paris, France, ⁴ Université de Toulouse, INRAE, INP-ENSAT Toulouse, UMR AGIR, Castanet-Tolosan, France

Developing new cropping strategies (very early sowing, crop expansion at higher latitudes, double cropping) to improve soybean production in Europe under climate change needs a good prediction of phenology under different temperature and photoperiod conditions. For that purpose, a simple phenology algorithm (SPA) was developed and parameterized for 10 contrasting soybean cultivars (maturity group 000 to II). Two experiments were carried out at INRA Toulouse (France) for parameterization: 1) Phenological monitoring of plants in pots on an outdoor platform with 6 planting dates. 2) Response of seed germination to temperature in controlled conditions. Multi-location field trials including 5 sites, 4 years, 2 sowing dates, and 10 cultivars were used to evaluate the SPA phenology predictions. Mean cardinal temperatures (minimum, optimum, and maximum) for germination were ca. 2, 30, and 40°C, respectively with significant differences among cultivars. The photoperiod sensitivity coefficient varied among cultivars when fixing Popt and Pct, optimal and critical photoperiods respectively, by maturity group. The parameterized algorithm showed an RMSE of less than 6 days for the prediction of crop cycle duration (i.e. cotyledons stage to physiological maturity) in the field trials including 75 data points. Flowering (R1 stage), and beginning of grain filling (R5 stage) dates were satisfactorily predicted with RMSEs of 8.2 and 9.4 days respectively. Because SPA can be also parameterized using data from field experiments, it can be useful as a plant selection tool across environments. The algorithm can be readily applied to species other than soybean, and its incorporation into cropping systems models would enhance the assessment of the performance of crop cultivars under climate change scenarios.

**Keywords:** planting date, glycine max, varietal phenotyping, photoperiod, temperature
# INTRODUCTION

EU-28 needs to import about 40 million tons of soybean (seed equivalent) for feeding livestock, and humans to a lesser extent. In recent years, soybean production in Europe has largely increased reaching 2.7 million tons in 2018 (USDA, 2019), but remains insufficient for the increasing demand of conventional and organic grains and soyfoods. In the context of climate change, southern areas in Europe are expected to experience more frequent and severe drought and heat waves while northern areas may benefit from higher temperatures in late season (IPCC, 2014). Three cropping strategies could be proposed to grow soybean under these conditions: (i) early sowing to use less irrigation water in summer as escaping strategy (Maury et al., 2015); (ii) northward cultivation of the crop to extend its cropping area toward regions becoming more suitable with increasing temperatures (Olesen et al., 2011); or (iii) double cropping after cereals under late sowing to fully use the thermal time window (Seifert and Lobell, 2015). These new crop management systems using either early or late sowing and new cultivation areas will expose the crop to a wider range of thermal and photoperiodic conditions under higher latitudes (from 43° to 52° N).

Soybean is a thermophilic plant whose development rate is governed by three cardinal temperatures: minimum temperature (T0, also referred to as base temperature), below which soybean does not develop; optimum temperature (Topt), at which soybean develops at the highest rate; and maximum temperature (Tmax), above which development is stopped. It is difficult to determine a definitive T0 temperature from the literature. Reported T0 values differ greatly, ranging from 2.5 to 13.2°C (Grimm et al., 1993; Sinclair et al., 2005; Salem et al., 2007; Setiyono et al., 2007) with a potential variation of T0 during the cycle (Setiyono et al., 2007). These variations are explained by the diversity of studied cultivars and methods used. Optimum temperature has been reported between 22 and 32°C depending on the studies (Grimm et al., 1993; Boote et al., 1998a; Setiyono et al., 2007; Boote, 2011; Parent and Tardieu, 2012) while maximum temperature has been less explored, with reported values between 40 and 47°C (Salem et al., 2007; Setiyono et al., 2007).

The change in sowing dates, mostly tested on late maturity groups from III to VII (Bowers, 1995; Hu and Wiatrak, 2012; Salmeron et al., 2014), can place the early maturing soybean cultivars of interest in this study under suboptimal temperature conditions. Indeed, early planting increases the likelihood of exposing soybean to temperatures lower than T0 during early developmental stages (Meyer and Badaruddin, 2001; Maury et al., 2015), whereas late sowings for double cropping will lead to exposure to suboptimal temperatures for development at the end of the cycle and a risk of frost before harvest (Mahieu and Brinkman, 1990; Shapiro et al., 1992). The accumulation of thermal time for development can be calculated in different ways: (i) a linear equation, where daily thermal time is calculated by subtracting T0 from daily average temperature and then accumulated over the period of interest; (ii) a bilinear equation, which takes into account both T0 and Topt to obtain an upper bound temperature above which thermal time is not

accumulated (Brisson et al., 2003); another bilinear triangular equation adds Tmax to impact the temperature accumulation between Topt and Tmax (Jones et al., 2003; Keating et al., 2003), and finally (iii) a nonlinear equation that uses the three cardinal temperatures and more closely represents temperature effect on biological functioning (Wang and Engel, 1998).

The characterization of cardinal temperatures is a basic step to parameterize thermal time equations, and the use of fast and reproducible methods for their determination is desirable. Parent et al. (2010) have shown that all the development processes respond in the same way to temperature. Therefore, it is convenient to choose a simple process for subsequent application to the whole plant. Salem et al. (2007) have chosen to characterize cardinal temperatures using pollen germination percentage and tube length, but it is also possible to use the germination process (Andreucci et al., 2016), which has the advantage of requiring only an incubator.

Soybean development rate is also strongly impacted by photoperiod due to its short-day plant behavior. Photoperiod regulates the onset of flowering (Garner and Allard, 1920; Hadley et al., 1984; Caffaro et al., 1988) and the duration of the phenological phases (Raper and Thomas, 1978; Kantolic and Slafer, 2001; Nico et al., 2015). An increase in day length above a given threshold, mainly related to the maturity group (Mcblain and Bernard, 1987; Xia et al., 2012; Lu et al., 2017), will slow down the development process. This characteristic has a particular impact on the post-flowering phase, which is more sensitive to photoperiod in late cultivars than in early ones (Major et al., 1975; Kantolic and Slafer, 2001). The type of growth (determinate vs. indeterminate) also results in different responses to photoperiod, especially between the onset of flowering (R1) and beginning seed fill (R5) stages where the final number of nodes is fixed for indeterminate cultivars. For these cultivars, long photoperiods will increase the length of R1 to R5 phase (Piper et al., 1996; Han et al., 2006; Kantolic et al., 2013).

Due to this response to photoperiod, a change in sowing dates will lead to development differences in soybean. Earlier planting will tend to initiate an early flowering, at a less advanced vegetative stage, since day length will be short enough to trigger it. The R1-R5 phase, which takes place when the day lengths are the longest, will be lengthened, allowing soybean to produce more pods and seeds (Kantolic and Slafer, 2001; Kantolic and Slafer, 2007; Kantolic et al., 2013; Nico et al., 2015). On the other hand, a late planting date will have the opposite effect, often with limitations of yield potential compared to conventional planting dates (Zhang et al., 2010; Hu, 2013; Clovis et al., 2015). Growth type, determinate or indeterminate, also plays a significant role. An indeterminate cultivar will have to be sown sooner, since the final number of nodes is reached at R5 (Piper et al., 1996).

In this study, we used a simple phenology algorithm accounting for temperature and photoperiod along with a protocol to determine the necessary parameters from controlled or field experiments. Models with different requirements of parameterization have been proposed that consider: (i) the total cycle (Piper et al., 1996), (ii) two phases
# INTRODUCTION

EU-28 needs to import about 40 million tons of soybean (seed equivalent) for feeding livestock, and humans to a lesser extent. In recent years, soybean production in Europe has largely increased reaching 2.7 million tons in 2018 (USDA, 2019), but remains insufficient for the increasing demand of conventional and organic grains and soyfoods. In the context of climate change, southern areas in Europe are expected to experience more frequent and severe drought and heat waves while northern areas may benefit from higher temperatures in late season (IPCC, 2014). Three cropping strategies could be proposed to grow soybean under these conditions: (i) early sowing to use less irrigation water in summer as escaping strategy (Maury et al., 2015); (ii) northward cultivation of the crop to extend its cropping area toward regions becoming more suitable with increasing temperatures (Olesen et al., 2011); or (iii) double cropping after cereals under late sowing to fully use the thermal time window (Seifert and Lobell, 2015). These new crop management systems using either early or late sowing and new cultivation areas will expose the crop to a wider range of thermal and photoperiodic conditions under higher latitudes (from 43° to 52° N).

Soybean is a thermophilic plant whose development rate is governed by three cardinal temperatures: minimum temperature (T0, also referred to as base temperature), below which soybean does not develop; optimum temperature (Topt), at which soybean develops at the highest rate; and maximum temperature (Tmax), above which development is stopped. It is difficult to determine a definitive T0 temperature from the literature. Reported T0 values differ greatly, ranging from 2.5 to 13.2°C (Grimm et al., 1993; Sinclair et al., 2005; Salem et al., 2007; Setiyono et al., 2007) with a potential variation of T0 during the cycle (Setiyono et al., 2007). These variations are explained by the diversity of studied cultivars and methods used. Optimum temperature has been reported between 22 and 32°C depending on the studies (Grimm et al., 1993; Boote et al., 1998a; Setiyono et al., 2007; Boote, 2011; Parent and Tardieu, 2012) while maximum temperature has been less explored, with reported values between 40 and 47°C (Salem et al., 2007; Setiyono et al., 2007).

The change in sowing dates, mostly tested on late maturity groups from III to VII (Bowers, 1995; Hu and Wiatrak, 2012; Salmeron et al., 2014), can place the early maturing soybean cultivars of interest in this study under suboptimal temperature conditions. Indeed, early planting increases the likelihood of exposing soybean to temperatures lower than T0 during early developmental stages (Meyer and Badaruddin, 2001; Maury et al., 2015), whereas late sowings for double cropping will lead to exposure to suboptimal temperatures for development at the end of the cycle and a risk of frost before harvest (Mahieu and Brinkman, 1990; Shapiro et al., 1992). The accumulation of thermal time for development can be calculated in different ways: (i) a linear equation, where daily thermal time is calculated by subtracting T0 from daily average temperature and then accumulated over the period of interest; (ii) a bilinear equation, which takes into account both T0 and Topt to obtain an upper bound temperature above which thermal time is not

accumulated (Brisson et al., 2003); another bilinear triangular equation adds Tmax to impact the temperature accumulation between Topt and Tmax (Jones et al., 2003; Keating et al., 2003), and finally (iii) a nonlinear equation that uses the three cardinal temperatures and more closely represents temperature effect on biological functioning (Wang and Engel, 1998).

The characterization of cardinal temperatures is a basic step to parameterize thermal time equations, and the use of fast and reproducible methods for their determination is desirable. Parent et al. (2010) have shown that all the development processes respond in the same way to temperature. Therefore, it is convenient to choose a simple process for subsequent application to the whole plant. Salem et al. (2007) have chosen to characterize cardinal temperatures using pollen germination percentage and tube length, but it is also possible to use the germination process (Andreucci et al., 2016), which has the advantage of requiring only an incubator.

Soybean development rate is also strongly impacted by photoperiod due to its short-day plant behavior. Photoperiod regulates the onset of flowering (Garner and Allard, 1920; Hadley et al., 1984; Caffaro et al., 1988) and the duration of the phenological phases (Raper and Thomas, 1978; Kantolic and Slafer, 2001; Nico et al., 2015). An increase in day length above a given threshold, mainly related to the maturity group (Mcblain and Bernard, 1987; Xia et al., 2012; Lu et al., 2017), will slow down the development process. This characteristic has a particular impact on the post-flowering phase, which is more sensitive to photoperiod in late cultivars than in early ones (Major et al., 1975; Kantolic and Slafer, 2001). The type of growth (determinate vs. indeterminate) also results in different responses to photoperiod, especially between the onset of flowering (R1) and beginning seed fill (R5) stages where the final number of nodes is fixed for indeterminate cultivars. For these cultivars, long photoperiods will increase the length of R1 to R5 phase (Piper et al., 1996; Han et al., 2006; Kantolic et al., 2013).

Due to this response to photoperiod, a change in sowing dates will lead to development differences in soybean. Earlier planting will tend to initiate an early flowering, at a less advanced vegetative stage, since day length will be short enough to trigger it. The R1-R5 phase, which takes place when the day lengths are the longest, will be lengthened, allowing soybean to produce more pods and seeds (Kantolic and Slafer, 2001; Kantolic and Slafer, 2007; Kantolic et al., 2013; Nico et al., 2015). On the other hand, a late planting date will have the opposite effect, often with limitations of yield potential compared to conventional planting dates (Zhang et al., 2010; Hu, 2013; Clovis et al., 2015). Growth type, determinate or indeterminate, also plays a significant role. An indeterminate cultivar will have to be sown sooner, since the final number of nodes is reached at R5 (Piper et al., 1996).

In this study, we used a simple phenology algorithm accounting for temperature and photoperiod along with a protocol to determine the necessary parameters from controlled or field experiments. Models with different requirements of parameterization have been proposed that consider: (i) the total cycle (Piper et al., 1996), (ii) two phases
—before and after a given stage, either R1 (Keating et al., 2003; Sinclair et al., 2005) or R5 (Brisson et al., 2003), (iii) a finer differentiation within phases (Setiyono et al., 2007; Salmerón and Purcell, 2016). Phenotyping these parameters is tedious, and mainly involves experimentation in pots inside culture chambers or greenhouses (Cregan and Hartwig, 1984; Ellis et al., 1992; Wang et al., 1998; Cober et al., 2014) or by extension or shortening of the photoperiod in the field through lights or shading systems (Kantolic and Slafer, 2001; Kantolic et al., 2013; Xu et al., 2013; Nico et al., 2015). It is also possible to use the natural variation of day length by testing several planting dates.

Depending on their responses to temperature and photoperiod, soybean cultivars are classified into different maturity groups (MGs). Today, thirteen maturity groups are included in the international classification system for soybean (Poehlman, 1987), classified from “000” for the very early to “X” (Roman numeral) for very late ones. Commonly cropped cultivars in Europe belong to maturity groups from 000 to II (Kurasch et al., 2017) while cultivars in the USA range from 00 to IX (Scott and Aldrich, 1970; Zhang et al., 2007), in Brazil from V to IX (Alliprandini et al., 2009), and in Argentina from II to IX (Dardanelli et al., 2006). The distribution of MGs for these regions depends mainly on the latitude, i.e. early cultivars are usually sown in high latitudes in northern hemisphere. MGs range from 000 to IX in China but are more related to eco-regions based on climatic and geographical conditions, cropping systems and season sowing types than to latitude (Yuesheng et al., 2006). The main difficulty encountered for extending soybean cropping to northern Europe and/or for modifying the sowing period is relative to the short-day behavior of the plant. As long photoperiods are usual in Europe, it is important to understand the temperature-photoperiod interactions of current and future cultivars.

The aim of this study was three-fold: (i) Develop a simplified phenotyping method for phenology in (semi-)controlled conditions (EXP1, EXP2), (ii) parameterize a simple phenology algorithm (SPA) using EXP1 and EXP2, and (iii) evaluate SPA using field experiments (EXP3).

# MATERIAL AND METHODS

## Experimental Design and Measurements
### Cultivars' Characteristics

Ten soybean cultivars representing diverse maturity groups and growth “types” (more or less indeterminate) used in Western Europe cropping systems were selected for this study. CVs ‘Isidor,’ ‘Santana,’ ‘Blancas,’ and ‘Ecudor,’ respectively from maturity groups I, I/II, and II, are usually grown under irrigation in South-West of France, which is the main French soybean production area. CVs ‘Klaxon,’ ‘Sultana,’ ‘RGT_Shouna,’ ‘ES_Mentor,’ and ‘Sigalia’ are considered as “very early” cultivars belonging to maturity groups from 000 to 00. They are usually grown as rainfed crops in the North-East of France, and in the North of Europe (Netherlands, Germany, Poland). Despite their small number, these cultivars have been

chosen in order to represent the widest range of phenological characteristics (MG, growth type) among elite cultivars cultivated in Europe (Table 1). For example, ‘ES_Pallador,’ which belongs to the same MG as ‘Isidor,’ has a totally different type of growth, being very indeterminate. CV ‘Blancas’ has variable classification depending on location, being classified as MG I/II in Italy and II in South-West of France.

## Germination Rate Experiment to Determine Cardinal Temperatures—EXP1

This experimentation was conducted in 2018. Seeds of the 10 cultivars were incubated under eleven temperatures (3, 6.5, 10, 15, 20, 25, 30, 35, 37.5, 40, and 43°C). To do this, hundred seeds of each cultivar were placed in four 90 mm diameter Petri dishes, resulting in a total of 25 seeds by plate. To maintain sufficient moisture inside the plate, the seeds were placed between two filter papers and kept moist with 8 ml of purified (osmosis) water. Petri dishes were arranged on four trays (one for each replicate) in the incubator or the cold chamber depending on the temperature tested. Temperature was monitored using a sensor (THERMOCHRON DS1922T, Embedded Data Systems). Germinated seeds were counted and removed two or three times a day. One seed was considered as germinated when it presented a radicle, outside of the tegument, of at least 3 mm long (McDonald et al., 1988). Data produced during EXP1 are presented in the first tab of **Data Sheet 1 (Supplementary Material)**. To determine the cardinal temperatures for each cultivar from EXP1, we followed the method presented by

<table><caption>TABLE 1 | Characteristics of the 10 cultivars used in this study.</caption><thead><tr><th>Cultivar</th><th>Breeder</th><th>Maturity group</th><th>Growth type</th><th>Leaf type</th></tr></thead><tbody><tr><td>KLAXON</td><td>RAGT 2N</td><td>000(0)</td><td>indeterminate</td><td>pointed-oval</td></tr><tr><td>RGT_SHOUNA</td><td>RAGT 2N</td><td>000</td><td>semi-indeterminate</td><td>pointed-oval</td></tr><tr><td>SULTANA</td><td>RAGT 2N</td><td>000</td><td>indeterminate</td><td>pointed-oval</td></tr><tr><td>ES_MENTOR</td><td>EURALIS SEMENCES</td><td>00</td><td>semi-indeterminate</td><td>rounded-oval</td></tr><tr><td>SIGALIA</td><td>RAGT 2N</td><td>00</td><td>semi-indeterminate to indeterminate</td><td>rounded-oval</td></tr><tr><td>ES_PALLADOR</td><td>EURALIS SEMENCES</td><td>I</td><td>indeterminate</td><td>lanceolate</td></tr><tr><td>ISIDOR</td><td>EURALIS SEMENCES</td><td>I</td><td>semi-determinate</td><td>big round leaves</td></tr><tr><td>SANTANA</td><td>RAGT 2N</td><td>I/II</td><td>indeterminate</td><td>rounded-oval</td></tr><tr><td>BLANCAS</td><td>CAUSSADE SEMENCES</td><td>II</td><td>indeterminate</td><td>rounded-oval</td></tr><tr><td>ECUDOR</td><td>EURALIS SEMENCES</td><td>II</td><td>indeterminate</td><td>big round leaves</td></tr></tbody></table>

They are ordered by maturity group (MG), from the earliest ('Klaxon'—MG 000) to the latest maturing one ('Ecudor'—MG II). 'Klaxon' was registered as MG 000 but according to the breeder's expertise it was classified in this study as MG 0000. The studied germplasm comes only from French breeders. At least two contrasted cultivars by MG in terms of growth type and leaf type were kept in order to maximize the difference regarding photoperiod-temperature interactions.
—before and after a given stage, either R1 (Keating et al., 2003; Sinclair et al., 2005) or R5 (Brisson et al., 2003), (iii) a finer differentiation within phases (Setiyono et al., 2007; Salmerón and Purcell, 2016). Phenotyping these parameters is tedious, and mainly involves experimentation in pots inside culture chambers or greenhouses (Cregan and Hartwig, 1984; Ellis et al., 1992; Wang et al., 1998; Cober et al., 2014) or by extension or shortening of the photoperiod in the field through lights or shading systems (Kantolic and Slafer, 2001; Kantolic et al., 2013; Xu et al., 2013; Nico et al., 2015). It is also possible to use the natural variation of day length by testing several planting dates.

Depending on their responses to temperature and photoperiod, soybean cultivars are classified into different maturity groups (MGs). Today, thirteen maturity groups are included in the international classification system for soybean (Poehlman, 1987), classified from “000” for the very early to “X” (Roman numeral) for very late ones. Commonly cropped cultivars in Europe belong to maturity groups from 000 to II (Kurasch et al., 2017) while cultivars in the USA range from 00 to IX (Scott and Aldrich, 1970; Zhang et al., 2007), in Brazil from V to IX (Alliprandini et al., 2009), and in Argentina from II to IX (Dardanelli et al., 2006). The distribution of MGs for these regions depends mainly on the latitude, i.e. early cultivars are usually sown in high latitudes in northern hemisphere. MGs range from 000 to IX in China but are more related to eco-regions based on climatic and geographical conditions, cropping systems and season sowing types than to latitude (Yuesheng et al., 2006). The main difficulty encountered for extending soybean cropping to northern Europe and/or for modifying the sowing period is relative to the short-day behavior of the plant. As long photoperiods are usual in Europe, it is important to understand the temperature-photoperiod interactions of current and future cultivars.

The aim of this study was three-fold: (i) Develop a simplified phenotyping method for phenology in (semi-)controlled conditions (EXP1, EXP2), (ii) parameterize a simple phenology algorithm (SPA) using EXP1 and EXP2, and (iii) evaluate SPA using field experiments (EXP3).

# MATERIAL AND METHODS

## Experimental Design and Measurements
### Cultivars' Characteristics

Ten soybean cultivars representing diverse maturity groups and growth “types” (more or less indeterminate) used in Western Europe cropping systems were selected for this study. CVs ‘Isidor,’ ‘Santana,’ ‘Blancas,’ and ‘Ecudor,’ respectively from maturity groups I, I/II, and II, are usually grown under irrigation in South-West of France, which is the main French soybean production area. CVs ‘Klaxon,’ ‘Sultana,’ ‘RGT_Shouna,’ ‘ES_Mentor,’ and ‘Sigalia’ are considered as “very early” cultivars belonging to maturity groups from 000 to 00. They are usually grown as rainfed crops in the North-East of France, and in the North of Europe (Netherlands, Germany, Poland). Despite their small number, these cultivars have been

chosen in order to represent the widest range of phenological characteristics (MG, growth type) among elite cultivars cultivated in Europe (Table 1). For example, ‘ES_Pallador,’ which belongs to the same MG as ‘Isidor,’ has a totally different type of growth, being very indeterminate. CV ‘Blancas’ has variable classification depending on location, being classified as MG I/II in Italy and II in South-West of France.

## Germination Rate Experiment to Determine Cardinal Temperatures—EXP1

This experimentation was conducted in 2018. Seeds of the 10 cultivars were incubated under eleven temperatures (3, 6.5, 10, 15, 20, 25, 30, 35, 37.5, 40, and 43°C). To do this, hundred seeds of each cultivar were placed in four 90 mm diameter Petri dishes, resulting in a total of 25 seeds by plate. To maintain sufficient moisture inside the plate, the seeds were placed between two filter papers and kept moist with 8 ml of purified (osmosis) water. Petri dishes were arranged on four trays (one for each replicate) in the incubator or the cold chamber depending on the temperature tested. Temperature was monitored using a sensor (THERMOCHRON DS1922T, Embedded Data Systems). Germinated seeds were counted and removed two or three times a day. One seed was considered as germinated when it presented a radicle, outside of the tegument, of at least 3 mm long (McDonald et al., 1988). Data produced during EXP1 are presented in the first tab of **Data Sheet 1 (Supplementary Material)**. To determine the cardinal temperatures for each cultivar from EXP1, we followed the method presented by

<table><caption>TABLE 1 | Characteristics of the 10 cultivars used in this study.</caption><thead><tr><th>Cultivar</th><th>Breeder</th><th>Maturity group</th><th>Growth type</th><th>Leaf type</th></tr></thead><tbody><tr><td>KLAXON</td><td>RAGT 2N</td><td>000(0)</td><td>indeterminate</td><td>pointed-oval</td></tr><tr><td>RGT_SHOUNA</td><td>RAGT 2N</td><td>000</td><td>semi-indeterminate</td><td>pointed-oval</td></tr><tr><td>SULTANA</td><td>RAGT 2N</td><td>000</td><td>indeterminate</td><td>pointed-oval</td></tr><tr><td>ES_MENTOR</td><td>EURALIS SEMENCES</td><td>00</td><td>semi-indeterminate</td><td>rounded-oval</td></tr><tr><td>SIGALIA</td><td>RAGT 2N</td><td>00</td><td>semi-indeterminate to indeterminate</td><td>rounded-oval</td></tr><tr><td>ES_PALLADOR</td><td>EURALIS SEMENCES</td><td>I</td><td>indeterminate</td><td>lanceolate</td></tr><tr><td>ISIDOR</td><td>EURALIS SEMENCES</td><td>I</td><td>semi-determinate</td><td>big round leaves</td></tr><tr><td>SANTANA</td><td>RAGT 2N</td><td>I/II</td><td>indeterminate</td><td>rounded-oval</td></tr><tr><td>BLANCAS</td><td>CAUSSADE SEMENCES</td><td>II</td><td>indeterminate</td><td>rounded-oval</td></tr><tr><td>ECUDOR</td><td>EURALIS SEMENCES</td><td>II</td><td>indeterminate</td><td>big round leaves</td></tr></tbody></table>

They are ordered by maturity group (MG), from the earliest ('Klaxon'—MG 000) to the latest maturing one ('Ecudor'—MG II). 'Klaxon' was registered as MG 000 but according to the breeder's expertise it was classified in this study as MG 0000. The studied germplasm comes only from French breeders. At least two contrasted cultivars by MG in terms of growth type and leaf type were kept in order to maximize the difference regarding photoperiod-temperature interactions.
Tribouillois et al. (2016). This method consists of using the Beta function described by Yin et al. (1995), which is fitted to experimental data to obtain the three cardinal temperatures (Equation 1). The Beta function adjustment was performed by simultaneously optimizing the three function parameters (mu, alpha, beta) with the Excel solver tool. A constraint was imposed on the T0 parameter so that it must be greater than or equal to 0 and Tmax were fixed according observations. The method used was the non-linear GRG with the goal of minimizing the RMSE between the simulated and the observed germination rate values.

Equation 1:

$$ \frac{1}{t(G)} = \exp(\mu) \times (T - T0)^{\alpha} \times (Tmax - T)^{\beta} $$

Where μ, α and β are the model parameters; 1/t(G), the germination rate; T0, T, and Tmax, the minimum, actual and maximum temperatures respectively.

## Response of Phenology to Temperature and Photoperiod in an Outdoor Pot Platform—EXP2

The 10 cultivars were grown on an automated outdoor pot platform (Heliaphen) designed for high-throughput plant phenotyping as previously described in Gosseau et al. (2019). As this experiment was set up to understand and model the temperature by photoperiod interactions, six planting dates were used: 17th March (D1), 6th April (D2), 16th May (D3), 26th June (D4), 24th July (D5), and 1st September 2017 (D6). These dates were chosen to cover the whole soybean growing season and to exhibit contrasting photoperiods and temperatures. Earlier (D1) and late (D6) planting dates experienced damaging cold stress and were not included for further analysis. In southwestern France, field planting dates generally occur between April 10th and May 20th, corresponding to the planting dates D2 and D3 of this experiment. D1 and D5 represent new management practices. D1 begins to be tested in areas in order to escape water deficit around flowering and seed filling time (Maury et al., 2015). D5 is typical of double cropping systems, where soybean is planted after an early harvested crop (winter barley, pea) (Quinsac et al., 2014).

Seeds were sown in a greenhouse at a mean temperature of 20°C with natural day length to standardize the initial growth of seedlings. The most vigorous seedlings were transplanted in pots and transferred to the platform. This was also a way to prevent emergence losses on the platform due to possible poor seed quality. Before transferring the seedlings, an acclimation regime at a temperature around 10°C during 7 days was applied for the two first planting dates to prevent a risk of thermal shock. Seedlings which reached VC stage (first unifoliate leaf) were transplanted in individual pots filled with 15 L of P.A.M.2 potting substrate (Proveen, distributed by Soprimex, Chateaurenard, Bouches-du-Rhône, France) and 125 g of extended released fertilizer (Osmocote Exact High K 5-6M, ICL Specialty Fertilizers, distributed by Agri Garonne, Castelginest, Haute-Garonne, France). For each cultivar and planting date, three pots of five plants were prepared. The experiment consisted of a split-plot design with three blocks. The main plot consisted of planting date treatments and subplots

of cultivars. The pots were organized by planting date to avoid shading effects on the youngest plants (cf. experimental design in **Figure A—Supplementary Material**). The 10 cultivars were randomized within each treatment-block combination. Pots were well-watered using drips. Hourly temperatures were recorded by an automatic weather station located within 150 m from the platform.

Plant phenology was recorded two times a week using Fehr and Caviness (1977) scale for vegetative and reproductive stages (description in **Table A—Supplementary Material**). Records began around VC stage and ended at R7 stage, which is the beginning of maturity (Setiyono et al., 2007). In total, the phenology of 900 plants was recorded and the day of appearance of each stage was calculated using the five plants in a pot (mean, standard deviation). Climatic and phenological data produced in EXP2 are summarized in the second and third tabs of **Data Sheet 1 (Supplementary Material)**.

## Field Experiments—EXP3

Eight field experiments were conducted in 2013, 2014, 2017, and 2018 in different experimental sites in the Toulouse region (43°N 1°E), namely Mondonville (2013, 2014), Rivières (2013, 2014), En Crambade (2013, 2014), and Auzeville (2017, 2018). The cultivars (CVs 'Ecudor,' 'Isidor,' 'Santana') were tested during the four seasons, but three more were added in 2017 (CVs 'RGT_Shouna,' 'Sultana,' 'Blancas'), and three more in 2018 (CVs 'Sigalia,' 'ES_Mentor,' 'ES_Pallador'). Two sowing dates (early and early + 1 to 1,5 month) were tested, except in 2018 where normal and late sowings were tested because of the wet spring. All the experiments were optimally irrigated using a decision tool (Terres Inovia, 2019), resulting in water applications from 65 to 161 mm. Daily temperatures were recorded by an automatic weather station located on each experimental site. At least four phenological stages were recorded when 50 percent of the plants in one plot had reached the stage: VE (emergence), R1 (first flower), R5 (beginning seed filling), R7 (beginning maturity). This information is summarized in **Table 2**. Climatic and agronomic data produced in EXP3 are summarized in the last two tabs of **Data Sheet 1 (Supplementary Material)**.

## Simple Phenology Algorithm Description

A simple phenology algorithm (SPA) was developed to predict the occurrence of soybean development stages, accounting for diverse photoperiod and temperature conditions based on concepts introduced by Robertson (1968), and expanded by Wang and Engel (1998). Under non-limiting conditions of photoperiod and temperature, the optimal physiological development days (PDDoptc,p) to complete a given phenological phase and the number of calendar days are the same (Wang and Engel, 1998). When the conditions are limiting, the actual physiological development day (PDD) accumulated each day is less than 1, and the number of calendar days to complete the phase is larger. The accumulation of PDD required to complete a phenological phase for each cultivar (PDDc,p) is
Tribouillois et al. (2016). This method consists of using the Beta function described by Yin et al. (1995), which is fitted to experimental data to obtain the three cardinal temperatures (Equation 1). The Beta function adjustment was performed by simultaneously optimizing the three function parameters (mu, alpha, beta) with the Excel solver tool. A constraint was imposed on the T0 parameter so that it must be greater than or equal to 0 and Tmax were fixed according observations. The method used was the non-linear GRG with the goal of minimizing the RMSE between the simulated and the observed germination rate values.

Equation 1:

$$ \frac{1}{t(G)} = \exp(\mu) \times (T - T0)^{\alpha} \times (Tmax - T)^{\beta} $$

Where μ, α and β are the model parameters; 1/t(G), the germination rate; T0, T, and Tmax, the minimum, actual and maximum temperatures respectively.

## Response of Phenology to Temperature and Photoperiod in an Outdoor Pot Platform—EXP2

The 10 cultivars were grown on an automated outdoor pot platform (Heliaphen) designed for high-throughput plant phenotyping as previously described in Gosseau et al. (2019). As this experiment was set up to understand and model the temperature by photoperiod interactions, six planting dates were used: 17th March (D1), 6th April (D2), 16th May (D3), 26th June (D4), 24th July (D5), and 1st September 2017 (D6). These dates were chosen to cover the whole soybean growing season and to exhibit contrasting photoperiods and temperatures. Earlier (D1) and late (D6) planting dates experienced damaging cold stress and were not included for further analysis. In southwestern France, field planting dates generally occur between April 10th and May 20th, corresponding to the planting dates D2 and D3 of this experiment. D1 and D5 represent new management practices. D1 begins to be tested in areas in order to escape water deficit around flowering and seed filling time (Maury et al., 2015). D5 is typical of double cropping systems, where soybean is planted after an early harvested crop (winter barley, pea) (Quinsac et al., 2014).

Seeds were sown in a greenhouse at a mean temperature of 20°C with natural day length to standardize the initial growth of seedlings. The most vigorous seedlings were transplanted in pots and transferred to the platform. This was also a way to prevent emergence losses on the platform due to possible poor seed quality. Before transferring the seedlings, an acclimation regime at a temperature around 10°C during 7 days was applied for the two first planting dates to prevent a risk of thermal shock. Seedlings which reached VC stage (first unifoliate leaf) were transplanted in individual pots filled with 15 L of P.A.M.2 potting substrate (Proveen, distributed by Soprimex, Chateaurenard, Bouches-du-Rhône, France) and 125 g of extended released fertilizer (Osmocote Exact High K 5-6M, ICL Specialty Fertilizers, distributed by Agri Garonne, Castelginest, Haute-Garonne, France). For each cultivar and planting date, three pots of five plants were prepared. The experiment consisted of a split-plot design with three blocks. The main plot consisted of planting date treatments and subplots

of cultivars. The pots were organized by planting date to avoid shading effects on the youngest plants (cf. experimental design in **Figure A—Supplementary Material**). The 10 cultivars were randomized within each treatment-block combination. Pots were well-watered using drips. Hourly temperatures were recorded by an automatic weather station located within 150 m from the platform.

Plant phenology was recorded two times a week using Fehr and Caviness (1977) scale for vegetative and reproductive stages (description in **Table A—Supplementary Material**). Records began around VC stage and ended at R7 stage, which is the beginning of maturity (Setiyono et al., 2007). In total, the phenology of 900 plants was recorded and the day of appearance of each stage was calculated using the five plants in a pot (mean, standard deviation). Climatic and phenological data produced in EXP2 are summarized in the second and third tabs of **Data Sheet 1 (Supplementary Material)**.

## Field Experiments—EXP3

Eight field experiments were conducted in 2013, 2014, 2017, and 2018 in different experimental sites in the Toulouse region (43°N 1°E), namely Mondonville (2013, 2014), Rivières (2013, 2014), En Crambade (2013, 2014), and Auzeville (2017, 2018). The cultivars (CVs 'Ecudor,' 'Isidor,' 'Santana') were tested during the four seasons, but three more were added in 2017 (CVs 'RGT_Shouna,' 'Sultana,' 'Blancas'), and three more in 2018 (CVs 'Sigalia,' 'ES_Mentor,' 'ES_Pallador'). Two sowing dates (early and early + 1 to 1,5 month) were tested, except in 2018 where normal and late sowings were tested because of the wet spring. All the experiments were optimally irrigated using a decision tool (Terres Inovia, 2019), resulting in water applications from 65 to 161 mm. Daily temperatures were recorded by an automatic weather station located on each experimental site. At least four phenological stages were recorded when 50 percent of the plants in one plot had reached the stage: VE (emergence), R1 (first flower), R5 (beginning seed filling), R7 (beginning maturity). This information is summarized in **Table 2**. Climatic and agronomic data produced in EXP3 are summarized in the last two tabs of **Data Sheet 1 (Supplementary Material)**.

## Simple Phenology Algorithm
### Description

A simple phenology algorithm (SPA) was developed to predict the occurrence of soybean development stages, accounting for diverse photoperiod and temperature conditions based on concepts introduced by Robertson (1968), and expanded by Wang and Engel (1998). Under non-limiting conditions of photoperiod and temperature, the optimal physiological development days (PDDopt,c,p) to complete a given phenological phase and the number of calendar days are the same (Wang and Engel, 1998). When the conditions are limiting, the actual physiological development day (PDD) accumulated each day is less than 1, and the number of calendar days to complete the phase is larger. The accumulation of PDD required to complete a phenological phase for each cultivar (PDDc,p) is
<table><caption>**TABLE 2** | Summary of field experiments used for the validation of the simple phenology algorithm.</caption><thead><tr><th>Year</th><th>Site</th><th>Planting type</th><th>Irrigation (mm)</th><th>Cultivars</th><th>Planting date</th><th>VE date</th><th>R1 date</th><th>R5 date</th><th>R7 date</th><th>Harvest date</th></tr></thead><tbody><tr><td>2013</td><td>En Crambade</td><td>Conventional</td><td>145</td><td>Ecudor, Isidor, Santana</td><td>25-Apr</td><td>10-May</td><td>9-Jul</td><td>6-Aug</td><td>21-Sep</td><td>7-Oct</td></tr><tr><td>2013</td><td>En Crambade</td><td>Very Early</td><td>143</td><td>Ecudor, Isidor, Santana</td><td>15-Mar</td><td>13-Apr</td><td>29-Jun</td><td>23-Jul</td><td>12-Sep</td><td>7-Oct</td></tr><tr><td>2013</td><td>Mondonville</td><td>Conventional</td><td>86</td><td>Ecudor, Isidor, Santana</td><td>27-May</td><td>13-Jun</td><td>22-Jul</td><td>21-Aug</td><td>27-Sep</td><td>17-Oct</td></tr><tr><td>2013</td><td>Mondonville</td><td>Very Early</td><td>65</td><td>Ecudor, Isidor, Santana</td><td>22-Mar</td><td>16-Apr</td><td>17-Jun</td><td>12-Aug</td><td>15-Sep</td><td>14-Oct</td></tr><tr><td>2013</td><td>Rivieres</td><td>Conventional</td><td>125</td><td>Ecudor, Isidor, Santana</td><td>6-May</td><td>29-May</td><td>9-Jul</td><td>15-Aug</td><td>27-Sep</td><td>19-Oct</td></tr><tr><td>2013</td><td>Rivieres</td><td>Very Early</td><td>161</td><td>Ecudor, Isidor, Santana</td><td>22-Mar</td><td>17-Apr</td><td>15-Jun</td><td>5-Aug</td><td>8-Sep</td><td>2-Oct</td></tr><tr><td>2014</td><td>En Crambade</td><td>Conventional</td><td>120</td><td>Ecudor, Isidor, Santana</td><td>30-Apr</td><td>11-May</td><td>27-Jun</td><td>29-Jul</td><td>10-Sep</td><td>22-Sep</td></tr><tr><td>2014</td><td>En Crambade</td><td>Very Early</td><td>129</td><td>Ecudor, Isidor, Santana</td><td>14-Mar</td><td>7-Apr</td><td>11-Jun</td><td>14-Jul</td><td>30-Aug</td><td>13-Sep</td></tr><tr><td>2014</td><td>Mondonville</td><td>Conventional</td><td>0</td><td>Ecudor, Isidor, Santana</td><td>6-May</td><td>23-May</td><td>29-Jun</td><td>3-Aug</td><td>6-Sep</td><td>6-Oct</td></tr><tr><td>2014</td><td>Rivieres</td><td>Conventional</td><td>140</td><td>Ecudor, Isidor, Santana</td><td>6-May</td><td>16-May</td><td>1-Jul</td><td>1-Aug</td><td>16-Sep</td><td>9-Oct</td></tr><tr><td>2014</td><td>Rivieres</td><td>Very Early</td><td>100</td><td>Ecudor, Isidor, Santana</td><td>18-Mar</td><td>9-Apr</td><td>13-Jun</td><td>9-Jul</td><td>30-Aug</td><td>26-Sep</td></tr><tr><td>2017</td><td>Auzeville</td><td>Conventional</td><td>112</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana</td><td>10-May</td><td>17-May</td><td>21-Jun</td><td>29-Jul</td><td>1-Sep</td><td>15-Sep</td></tr><tr><td>2017</td><td>Auzeville</td><td>Very Early</td><td>76</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana</td><td>21-Mar</td><td>7-Apr</td><td>2-Jun</td><td>4-Jul</td><td>13-Aug</td><td>5-Sep</td></tr><tr><td>2018</td><td>Auzeville</td><td>Conventional</td><td>76</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana, ES_Pallador, Sigalia, ES_Mentor</td><td>24-Apr</td><td>10-May</td><td>22-Jun</td><td>24-Jul</td><td>26-Aug</td><td>10-Sep</td></tr><tr><td>2018</td><td>Auzeville</td><td>Late</td><td>111</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana, ES_Pallador, Sigalia, ES_Mentor</td><td>4-Jun</td><td>11-Jun</td><td>19-Jul</td><td>13-Aug</td><td>10-Sep</td><td>21-Sep</td></tr></tbody></table>

The four sites (En Crambade, Mondonville, Rivières, Auzeville) were located in the South-West of France, around Toulouse. All experiments were irrigated and had two planting dates, except Mondonville in 2014. The same three cultivars ('Ecudor,' 'Santana,' 'Isidor') were tested in 2013, 2014, 2017, and 2018. Respectively three and six cultivars were added in 2017 and 2018. Mean phenological stages needed for the study are reported for each experiment.

calculated using equation 2 (detailed in **Supplementary Material**). PDD<sub>c,p</sub> was calculated for four development phases: VC-R7, VC-R1, R1-R5, and R5-R7.

Equation 2:

$$PDD_{c,p} = PDDopt_{c,p} / \sum_{i=1}^{PDD_{c,p}} (f(Ti) \times f(Pi) / PDD_{c,p})$$

Where:

PDD<sub>c,p</sub>: Physiological Development Days of cultivar c and phase p (in calendar day)

f(T) and f(P) are the temperature and photoperiod functions, respectively.

c: cultivar

p: phase

i: day number

## Temperature Function

The temperature function (Equation 3) was described by Wang and Engel (1998).

Equation 3:

$$f(T) = \frac{2(T-T_0)^\alpha (T_{opt}-T_0)^\alpha - (T-T_0)^{2\alpha}}{(T_{opt}-T_0)^{2\alpha}} \quad \text{if } T_0 \le T \le T_{\text{max}}$$

$$f(T) = 0 \quad \text{if } T < T_0 \text{ or } T > T_{\text{max}}$$

Where:

T: average temperature of the day (°C).

T0, Topt, Tmax: cultivar-specific cardinal temperatures (°C).

The parameter α is calculated using Equation 4:
Equation 4:

$$\alpha = \log(2)/\log[(T_{\text{max}} - T_0)/(T_{\text{opt}} - T_0)]$$

## Photoperiod Function

The photoperiod function f(P) was calculated as shown in Equation 5.

Equation 5:

$$f(P) = 1 - \left[ \frac{(P - Popt)}{(Pcrt - Popt)} \right]^S \quad \text{if } Popt \le P \le Pcrt$$

$$f(P) = 1 \quad \text{if } Popt < P$$

$$f(P) = 0 \quad \text{if } Popt > P$$

Where Pcrt is the day length above which development rate is zero, Popt is the day length below which development rate is not limited, P is the day length including civil twilight, all in hours, and S is the sensitivity coefficient of the cultivar to photoperiod (S=0 for highest sensitivity; >100 for no sensitivity).
<table><caption>**TABLE 2** | Summary of field experiments used for the validation of the simple phenology algorithm.</caption><thead><tr><th>Year</th><th>Site</th><th>Planting type</th><th>Irrigation (mm)</th><th>Cultivars</th><th>Planting date</th><th>VE date</th><th>R1 date</th><th>R5 date</th><th>R7 date</th><th>Harvest date</th></tr></thead><tbody><tr><td>2013</td><td>En Crambade</td><td>Conventional</td><td>145</td><td>Ecudor, Isidor, Santana</td><td>25-Apr</td><td>10-May</td><td>9-Jul</td><td>6-Aug</td><td>21-Sep</td><td>7-Oct</td></tr><tr><td>2013</td><td>En Crambade</td><td>Very Early</td><td>143</td><td>Ecudor, Isidor, Santana</td><td>15-Mar</td><td>13-Apr</td><td>29-Jun</td><td>23-Jul</td><td>12-Sep</td><td>7-Oct</td></tr><tr><td>2013</td><td>Mondonville</td><td>Conventional</td><td>86</td><td>Ecudor, Isidor, Santana</td><td>27-May</td><td>13-Jun</td><td>22-Jul</td><td>21-Aug</td><td>27-Sep</td><td>17-Oct</td></tr><tr><td>2013</td><td>Mondonville</td><td>Very Early</td><td>65</td><td>Ecudor, Isidor, Santana</td><td>22-Mar</td><td>16-Apr</td><td>17-Jun</td><td>12-Aug</td><td>15-Sep</td><td>14-Oct</td></tr><tr><td>2013</td><td>Rivieres</td><td>Conventional</td><td>125</td><td>Ecudor, Isidor, Santana</td><td>6-May</td><td>29-May</td><td>9-Jul</td><td>15-Aug</td><td>27-Sep</td><td>19-Oct</td></tr><tr><td>2013</td><td>Rivieres</td><td>Very Early</td><td>161</td><td>Ecudor, Isidor, Santana</td><td>22-Mar</td><td>17-Apr</td><td>15-Jun</td><td>5-Aug</td><td>8-Sep</td><td>2-Oct</td></tr><tr><td>2014</td><td>En Crambade</td><td>Conventional</td><td>120</td><td>Ecudor, Isidor, Santana</td><td>30-Apr</td><td>11-May</td><td>27-Jun</td><td>29-Jul</td><td>10-Sep</td><td>22-Sep</td></tr><tr><td>2014</td><td>En Crambade</td><td>Very Early</td><td>129</td><td>Ecudor, Isidor, Santana</td><td>14-Mar</td><td>7-Apr</td><td>11-Jun</td><td>14-Jul</td><td>30-Aug</td><td>13-Sep</td></tr><tr><td>2014</td><td>Mondonville</td><td>Conventional</td><td>0</td><td>Ecudor, Isidor, Santana</td><td>6-May</td><td>23-May</td><td>29-Jun</td><td>3-Aug</td><td>6-Sep</td><td>6-Oct</td></tr><tr><td>2014</td><td>Rivieres</td><td>Conventional</td><td>140</td><td>Ecudor, Isidor, Santana</td><td>6-May</td><td>16-May</td><td>1-Jul</td><td>1-Aug</td><td>16-Sep</td><td>9-Oct</td></tr><tr><td>2014</td><td>Rivieres</td><td>Very Early</td><td>100</td><td>Ecudor, Isidor, Santana</td><td>18-Mar</td><td>9-Apr</td><td>13-Jun</td><td>9-Jul</td><td>30-Aug</td><td>26-Sep</td></tr><tr><td>2017</td><td>Auzeville</td><td>Conventional</td><td>112</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana</td><td>10-May</td><td>17-May</td><td>21-Jun</td><td>29-Jul</td><td>1-Sep</td><td>15-Sep</td></tr><tr><td>2017</td><td>Auzeville</td><td>Very Early</td><td>76</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana</td><td>21-Mar</td><td>7-Apr</td><td>2-Jun</td><td>4-Jul</td><td>13-Aug</td><td>5-Sep</td></tr><tr><td>2018</td><td>Auzeville</td><td>Conventional</td><td>76</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana, ES_Pallador, Sigalia, ES_Mentor</td><td>24-Apr</td><td>10-May</td><td>22-Jun</td><td>24-Jul</td><td>26-Aug</td><td>10-Sep</td></tr><tr><td>2018</td><td>Auzeville</td><td>Late</td><td>111</td><td>Ecudor, Isidor, Santana, Blancas, RGT_Shouna, Sultana, ES_Pallador, Sigalia, ES_Mentor</td><td>4-Jun</td><td>11-Jun</td><td>19-Jul</td><td>13-Aug</td><td>10-Sep</td><td>21-Sep</td></tr></tbody></table>

The four sites (En Crambade, Mondonville, Rivières, Auzeville) were located in the South-West of France, around Toulouse. All experiments were irrigated and had two planting dates, except Mondonville in 2014. The same three cultivars ('Ecudor,' 'Santana,' 'Isidor') were tested in 2013, 2014, 2017, and 2018. Respectively three and six cultivars were added in 2017 and 2018. Mean phenological stages needed for the study are reported for each experiment.

calculated using equation 2 (detailed in **Supplementary Material**). PDD<sub>c,p</sub> was calculated for four development phases: VC-R7, VC-R1, R1-R5, and R5-R7.

Equation 2:

$$ PDD_{c,p} = PDDopt_{c,p} / \sum_{i=1}^{PDD_{c,p}} (f(Ti) \times f(Pi) / PDD_{c,p}) $$

Where:

PDD<sub>c,p</sub>: Physiological Development Days of cultivar c and phase p (in calendar day)

f(T) and f(P) are the temperature and photoperiod functions, respectively.

c: cultivar

p: phase

i: day number

## Temperature Function

The temperature function (Equation 3) was described by Wang and Engel (1998).

Equation 3:

$$ f(T) = \frac{2(T-T_0)^\alpha (T_{opt}-T_0)^\alpha - (T-T_0)^{2\alpha}}{(T_{opt}-T_0)^{2\alpha}} \quad \text{if } T_0 \le T \le T_{max} $$

$$ f(T) = 0 \quad \text{if } T < T_0 \text{ or } T > T_{max} $$

Where:

T: average temperature of the day (°C).

T0, Topt, Tmax: cultivar-specific cardinal temperatures (°C).

The parameter α is calculated using Equation 4:
Equation 4:

$$ \alpha = \log(2)/\log[(T_{max} - T_0)/(T_{opt} - T_0)] $$

## Photoperiod Function

The photoperiod function f(P) was calculated as shown in Equation 5.

Equation 5:

$$ f(P) = 1 - \left[ \frac{(P - Popt)}{(Pcrt - Popt)} \right]^S \quad \text{if } Popt \le P \le Pcrt $$

$$ f(P) = 1 \quad \text{if } Popt < P $$

$$ f(P) = 0 \quad \text{if } Popt > P $$

Where Pcrt is the day length above which development rate is zero, Popt is the day length below which development rate is not limited, P is the day length including civil twilight, all in hours, and S is the sensitivity coefficient of the cultivar to photoperiod (S=0 for highest sensitivity; >100 for no sensitivity).
<table><caption>**TABLE 3** | Parameters used in Simple Phenology Algorithm (SPA).</caption><thead><tr><th>Name</th><th>Unit</th><th>Abbreviation</th><th>Origin of data</th></tr></thead><tbody><tr><td>Minimum temperature of development</td><td>(°C)</td><td>T0</td><td>EXP1</td></tr><tr><td>Optimal temperature of development</td><td>(°C)</td><td>Topt</td><td>EXP1</td></tr><tr><td>Maximum temperature of development</td><td>(°C)</td><td>Tmax</td><td>EXP1</td></tr><tr><td>Optimal daylength for maximum development</td><td>(h)</td><td>Popt</td><td>Literature</td></tr><tr><td>Critical daylength for zero development</td><td>(h)</td><td>Pcrt</td><td>Literature</td></tr><tr><td>Physiological development days in optimum conditions of temperature and photoperiod</td><td>d</td><td>PDDopt<sub>c,p</sub></td><td>Optimization</td></tr><tr><td>Sensitivity coefficient to the photoperiod</td><td>unitless</td><td>S</td><td>Optimization</td></tr></tbody></table>

Plant parameters are divided into three types depending on their method of acquisition. Cardinal temperatures (T0, Topt, Tmax) were directly produced in EXP1 at a cultivar level. Popt and Pcrt have been adopted from (Setiyono et al., 2007), based on the cultivar's MG. PDDopt<sub>c,p</sub> and sensitivity coefficient to the photoperiod were optimized during the parameterization phase.

## Parameterization

Plant parameters needed in SPA are presented in **Table 3**, obtained from three methods of acquisition: experimental data, literature, and optimization.

Once all other parameters were determined (i.e. cardinal temperatures, Popt and Pcrt—**Table 3**), PDDopt<sub>c,p</sub> and S were set by simultaneous optimization over the entire life cycle VC-R7. The optimization was performed with the Excel solver using the non-linear GRG method (Newton Global). It was carried out in two steps: 1) simultaneous optimization of the minimum number of days for the VC-R7 phase (PDDopt<sub>c,p</sub>) and of the sensitivity to the photoperiod (S); 2) optimization of the minimum number of days for the R1-R5 and R5-R7 phases while maintaining the value of S as determined in (1). The objective of the optimization was to minimize the RMSE between predicted and actual day of appearance of R7 for all planting dates of a given cultivar.

## Evaluation of SPA

SPA was evaluated based on a set of field data corresponding to 75 cases including 5 sites, 4 years, 2 sowing dates, and 10 cultivars (cf **Table 3**). The parameterized SPA was used to predict the phenology in these independent conditions. As we did not have a VC rating in 2013 and 2014, we estimated it by considering that it was reached at 50°C.day⁻¹ after VE (as observed in 2017 and 2018). Indeed, only thermal effects play a role at this stage of development (Hodges and French, 1985).

## Statistical Analysis

The effects of cultivars, planting dates, and their interactions were identified by analysis of variance using a two-way ANOVA performed with R software (version 3.2.2; [R Core Team 2015]) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna]). Comparison of means was done with post-hoc Student-Newman-Keuls (SNK) tests with a significance level of p < 0.05.

<table><caption>**TABLE 4** | Performance indicators used for the evaluation of SPA. Equations and units are reported.</caption><thead><tr><th>Performance Indicator</th><th>Equation</th><th>Unit</th></tr></thead><tbody><tr><td>Root Mean Square Error (RMSE)</td><td>RMSE = √(1/n ∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - y<sub>i</sub>)<sup>2</sup>)</td><td>days</td></tr><tr><td>Mean Bias Error (MBE)</td><td>MBE = ∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - y<sub>i</sub>) / n</td><td>days</td></tr><tr><td>Relative RMSE (RRMSE)</td><td>RRMSE = RMSE / ŷ</td><td>unitless</td></tr><tr><td>Efficiency (EF)</td><td>[EF = 1 - (∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - ŷ̂<sub>i</sub>)<sup>2</sup>) / (∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - ŷ̂<sub>i</sub>)<sup>2</sup>)] × 100</td><td>%</td></tr></tbody></table>

Where ŷ<sub>i</sub> is the observed value for the i observation, ŷ̂<sub>i</sub> is the predicted value, and ŷ the average of observations.

## Indicators of Algorithm Performance

The performance of SPA was evaluated with several indicators including RMSE (Root Mean Square Error), MBE (Mean Bias Error), RRMSE (Relative Root Mean Square Error), and EF (Efficiency). Equations of performance indicators are presented in **Table 4**.

# RESULTS

## Cardinal Temperatures of Germination

The Yin model fitted well the observation data; these are available in **Data Sheet 1 (Supplementary Material)**. Determination coefficients R² ranked from 0.95 to 1. The germination rate progressively increased from the minimum temperature to the optimum temperature and then decreased. The maximum germination rates varied according to the cultivars (**Figure 1**).

<figure><img src="image_2.png" /><figcaption>**FIGURE 1** | Graphical display of the Yin adjustments for all cultivars at 50% of germination. Dots represent mean observed data and lines the adjusted Yin function. Position of cardinal temperatures (T0, Topt, Tmax) are illustrated on the graph.</figcaption></figure>
<table><caption>**TABLE 3** | Parameters used in Simple Phenology Algorithm (SPA).</caption><thead><tr><th>Name</th><th>Unit</th><th>Abbreviation</th><th>Origin of data</th></tr></thead><tbody><tr><td>Minimum temperature of development</td><td>(°C)</td><td>T0</td><td>EXP1</td></tr><tr><td>Optimal temperature of development</td><td>(°C)</td><td>Topt</td><td>EXP1</td></tr><tr><td>Maximum temperature of development</td><td>(°C)</td><td>Tmax</td><td>EXP1</td></tr><tr><td>Optimal daylength for maximum development</td><td>(h)</td><td>Popt</td><td>Literature</td></tr><tr><td>Critical daylength for zero development</td><td>(h)</td><td>Pcrt</td><td>Literature</td></tr><tr><td>Physiological development days in optimum conditions of temperature and photoperiod</td><td>d</td><td>PDDopt<sub>c,p</sub></td><td>Optimization</td></tr><tr><td>Sensitivity coefficient to the photoperiod</td><td>unitless</td><td>S</td><td>Optimization</td></tr></tbody></table>

Plant parameters are divided into three types depending on their method of acquisition. Cardinal temperatures (T0, Topt, Tmax) were directly produced in EXP1 at a cultivar level. Popt and Pcrt have been adopted from (Setiyono et al., 2007), based on the cultivar's MG. PDDopt<sub>c,p</sub> and sensitivity coefficient to the photoperiod were optimized during the parameterization phase.

## Parameterization

Plant parameters needed in SPA are presented in **Table 3**, obtained from three methods of acquisition: experimental data, literature, and optimization.

Once all other parameters were determined (i.e. cardinal temperatures, Popt and Pcrt—**Table 3**), PDDopt<sub>c,p</sub> and S were set by simultaneous optimization over the entire life cycle VC-R7. The optimization was performed with the Excel solver using the non-linear GRG method (Newton Global). It was carried out in two steps: 1) simultaneous optimization of the minimum number of days for the VC-R7 phase (PDDopt<sub>c,p</sub>) and of the sensitivity to the photoperiod (S); 2) optimization of the minimum number of days for the R1-R5 and R5-R7 phases while maintaining the value of S as determined in (1). The objective of the optimization was to minimize the RMSE between predicted and actual day of appearance of R7 for all planting dates of a given cultivar.

## Evaluation of SPA

SPA was evaluated based on a set of field data corresponding to 75 cases including 5 sites, 4 years, 2 sowing dates, and 10 cultivars (cf **Table 3**). The parameterized SPA was used to predict the phenology in these independent conditions. As we did not have a VC rating in 2013 and 2014, we estimated it by considering that it was reached at 50°C.day⁻¹ after VE (as observed in 2017 and 2018). Indeed, only thermal effects play a role at this stage of development (Hodges and French, 1985).

## Statistical Analysis

The effects of cultivars, planting dates, and their interactions were identified by analysis of variance using a two-way ANOVA performed with R software (version 3.2.2; [R Core Team {2015} R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna]). Comparison of means was done with post-hoc Student-Newman-Keuls (SNK) tests with a significance level of p < 0.05.

<table><caption>**TABLE 4** | Performance indicators used for the evaluation of SPA. Equations and units are reported.</caption><thead><tr><th>Performance Indicator</th><th>Equation</th><th>Unit</th></tr></thead><tbody><tr><td>Root Mean Square Error (RMSE)</td><td>RMSE = √(1/n ∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - y<sub>i</sub>)<sup>2</sup>)</td><td>days</td></tr><tr><td>Mean Bias Error (MBE)</td><td>MBE = ∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - y<sub>i</sub>) / n</td><td>days</td></tr><tr><td>Relative RMSE (RRMSE)</td><td>RRMSE = RMSE / ŷ</td><td>unitless</td></tr><tr><td>Efficiency (EF)</td><td>[EF = 1 - (∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - ŷ̂<sub>i</sub>)<sup>2</sup>) / (∑<sub>i=1</sub><sup>n</sup>(ŷ<sub>i</sub> - ŷ̂<sub>i</sub>)<sup>2</sup>)] × 100</td><td>%</td></tr></tbody></table>

Where ŷ<sub>i</sub> is the observed value for the i observation, ŷ̂<sub>i</sub> is the predicted value, and ŷ the average of observations.

## Indicators of Algorithm Performance

The performance of SPA was evaluated with several indicators including RMSE (Root Mean Square Error), MBE (Mean Bias Error), RRMSE (Relative Root Mean Square Error), and EF (Efficiency). Equations of performance indicators are presented in **Table 4**.

# RESULTS

## Cardinal Temperatures of Germination

The Yin model fitted well the observation data; these are available in **Data Sheet 1 (Supplementary Material)**. Determination coefficients R² ranked from 0.95 to 1. The germination rate progressively increased from the minimum temperature to the optimum temperature and then decreased. The maximum germination rates varied according to the cultivars (**Figure 1**).

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Graphical display of the Yin adjustments for all cultivars at 50% of germination. Dots represent mean observed data and lines the adjusted Yin function. Position of cardinal temperatures (T0, Topt, Tmax) are illustrated on the graph.</figcaption></figure>
<table><caption>**TABLE 5** | Estimates (and standard errors (±)) of optimum (Topt) and maximum (Tmax) temperatures obtained by the Yin function for the germination of the 10 soybean cultivars studied.</caption><thead><tr><th>MG</th><th>Cultivars</th><th>Topt (°C)</th><th></th><th>Tmax (°C)</th><th></th></tr></thead><tbody><tr><td>000(0)</td><td> Klaxon</td><td>29.8 ± 0.11</td><td>c</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>000</td><td> RGT_Shouna</td><td>29.8 ± 0.23</td><td>c</td><td>40.7 ± 0.38</td><td>bc</td></tr><tr><td>000</td><td> Sultana</td><td>29.6 ± 0.48</td><td>c</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>00</td><td> ES_Mentor</td><td>32.1 ± 0.76</td><td>a</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>00</td><td> Sigalia</td><td>29.2 ± 0.50</td><td>c</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>I</td><td> ES_Pallador</td><td>30.7 ± 0.21</td><td>b</td><td>40.3 ± 0.00</td><td>c</td></tr><tr><td>I</td><td> Isidor</td><td>29.8 ± 0.46</td><td>c</td><td>40.3 ± 0.00</td><td>c</td></tr><tr><td>I/II</td><td> Santana</td><td>30.5 ± 0.40</td><td>b</td><td>40.3 ± 0.00</td><td>c</td></tr><tr><td>II</td><td> Blancas</td><td>29.7 ± 0.54</td><td>c</td><td>41.0 ± 0.00</td><td>b</td></tr><tr><td>II</td><td> Ecudor</td><td>28.1 ± 0.64</td><td>d</td><td>42.0 ± 0.00</td><td>a</td></tr><tr><td></td><td> Average</td><td>29.9 ± 0.43</td><td></td><td>40.6 ± 0.16</td><td></td></tr></tbody></table>

T0 was set to 2°C. In a column, means followed by common letter(s) are not significantly different at *p* < 0.05 level by Student-Newman-Keuls test.

The average T0, Topt, Tmax of all the 10 cultivars were 2, 28.1, 40.6°C respectively. The lowest temperature tested was 3°C and having no point of reference we chose to set T0 to 2 °C. The final germination percentage was not significantly affected in the 6.5 to 40 °C temperature range under these experimental conditions (Table B—Supplementary Material).

Cultivar differences were significant for optimum and maximum temperatures (Table 5). The optimum temperature was the most discriminating, and allowed to separate four significantly different groups at the *p* < 0.05 threshold. 'ES_Mentor' is the cultivar which takes the best advantage of high temperatures (Topt = 32.1°C) while 'Ecudor' requires a lower temperature (28.1°C). In contrast to Topt, Tmax for 'Ecudor' was the highest of the cultivar group tested (42°C), with a majority of cultivars being between 40.3 and 41°C.

## Temperature and Photoperiod Effects on Plant Phenology in the Outdoor Platform Temperature and Photoperiod Conditions

The different planting dates did not experience the same photoperiod and temperature conditions. For example, D1, D2, and D3 met a varying photoperiod, increasing at the beginning and decreasing at the end of the cycle. D4, D5, and D6 only met a decreasing photoperiod. The temperature and photoperiod conditions encountered during EXP2 are presented in Figure B (Supplementary Material).

### Impact of Planting Date and Cultivars on the Duration of Phenological Phases

The variance analysis of cycle duration (from VC to R7 stage) in calendar days (CD VC-R7) showed a highly significant effect (*p* < 0.001) of planting date, cultivar, and the interaction of these two factors (Table 6). Block effect and Block X Cultivar interaction were not significant. Total cycle duration showed an amplitude from 100 to 155 days across planting dates, all cultivars combined. Results are available in Data Sheet 1 (Supplementary Material).

<table><caption>**TABLE 6** | Analysis of variance of cycle duration from VC to R7 stage in calendar days (CD VC-R7) of 10 cultivars of soybean for five planting dates in outdoor pot experiment (EXP2).</caption><thead><tr><th>Source of variation</th><th>Df</th><th>CD VC-R7</th><th>MS</th></tr></thead><tbody><tr><td>Planting date</td><td>4</td><td></td><td>5118***</td></tr><tr><td>Block</td><td>2</td><td></td><td>48.</td></tr><tr><td>Cultivar</td><td>9</td><td></td><td>900***</td></tr><tr><td>Planting date × Cultivar</td><td>36</td><td></td><td>59***</td></tr><tr><td>Planting date × Block</td><td>8</td><td></td><td>58**</td></tr><tr><td>Block × Cultivar</td><td>18</td><td></td><td>20</td></tr><tr><td>Residuals</td><td>66</td><td></td><td>17</td></tr></tbody></table>

Planting dates analyzed were D1 to D5, as D6 did not reach R7 stage.

(.), (**), and (***), significant at 0.1, 0.01, 0.001 levels, respectively.

Df: Degree of freedom; MS: Mean Square.

## Simple Phenology Algorithm Parameterization of SPA

The parameterization led to satisfactory prediction of the dates of stages. Figure 2 shows predicted and observed occurrence (day of year) of R7 stage for planting dates D2 to D5. Despite a slight overestimation for D2, predicted values were close to those observed, and the quality of adjustment was satisfactory with an efficiency of 81% and an RMSE of 6.42 days for the appearance of R7 stage since the day of planting. Parameter values are presented in Tables 5, 6, 7.

### When Does Photoperiod and Temperature Affect Soybean Phenology?

The parameterization of SPA made it possible to highlight the phases impacted by limiting temperatures or photoperiods for each cultivar. As an example, Figure 3 shows individually the effects of f(T) and f(P) and their multiplicative effects on the

<figure><img src="image_3.png" /><figcaption>**FIGURE 2** | Comparison of fitted (FIT) and observed (OBS) day of appearance of R7 stage, for algorithm parameterization on outdoor pot experiment data (EXP2). Colors represent the four planting dates D2 to D5.</figcaption></figure>
<table><caption>**TABLE 5** | Estimates (and standard errors (±)) of optimum (Topt) and maximum (Tmax) temperatures obtained by the Yin function for the germination of the 10 soybean cultivars studied.</caption><thead><tr><th>MG</th><th>Cultivars</th><th>Topt (°C)</th><th></th><th>Tmax (°C)</th><th></th></tr></thead><tbody><tr><td>000(0)</td><td> Klaxon</td><td>29.8 ± 0.11</td><td>c</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>000</td><td> RGT_Shouna</td><td>29.8 ± 0.23</td><td>c</td><td>40.7 ± 0.38</td><td>bc</td></tr><tr><td>000</td><td> Sultana</td><td>29.6 ± 0.48</td><td>c</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>00</td><td> ES_Mentor</td><td>32.1 ± 0.76</td><td>a</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>00</td><td> Sigalia</td><td>29.2 ± 0.50</td><td>c</td><td>40.4 ± 0.31</td><td>c</td></tr><tr><td>I</td><td> ES_Pallador</td><td>30.7 ± 0.21</td><td>b</td><td>40.3 ± 0.00</td><td>c</td></tr><tr><td>I</td><td> Isidor</td><td>29.8 ± 0.46</td><td>c</td><td>40.3 ± 0.00</td><td>c</td></tr><tr><td>I/II</td><td> Santana</td><td>30.5 ± 0.40</td><td>b</td><td>40.3 ± 0.00</td><td>c</td></tr><tr><td>II</td><td> Blancas</td><td>29.7 ± 0.54</td><td>c</td><td>41.0 ± 0.00</td><td>b</td></tr><tr><td>II</td><td> Ecudor</td><td>28.1 ± 0.64</td><td>d</td><td>42.0 ± 0.00</td><td>a</td></tr><tr><td></td><td> Average</td><td>29.9 ± 0.43</td><td></td><td>40.6 ± 0.16</td><td></td></tr></tbody></table>

T0 was set to 2°C. In a column, means followed by common letter(s) are not significantly different at *p* < 0.05 level by Student-Newman-Keuls test.

The average T0, Topt, Tmax of all the 10 cultivars were 2, 28.1, 40.6°C respectively. The lowest temperature tested was 3°C and having no point of reference we chose to set T0 to 2 °C. The final germination percentage was not significantly affected in the 6.5 to 40 °C temperature range under these experimental conditions (Table B—Supplementary Material).

Cultivar differences were significant for optimum and maximum temperatures (Table 5). The optimum temperature was the most discriminating, and allowed to separate four significantly different groups at the *p* < 0.05 threshold. ‘ES_Mentor’ is the cultivar which takes the best advantage of high temperatures (Topt = 32.1°C) while ‘Ecudor’ requires a lower temperature (28.1°C). In contrast to Topt, Tmax for ‘Ecudor’ was the highest of the cultivar group tested (42°C), with a majority of cultivars being between 40.3 and 41°C.

## Temperature and Photoperiod Effects on Plant Phenology in the Outdoor Platform Temperature and Photoperiod Conditions

The different planting dates did not experience the same photoperiod and temperature conditions. For example, D1, D2, and D3 met a varying photoperiod, increasing at the beginning and decreasing at the end of the cycle. D4, D5, and D6 only met a decreasing photoperiod. The temperature and photoperiod conditions encountered during EXP2 are presented in **Figure B (Supplementary Material).**

### Impact of Planting Date and Cultivars on the Duration of Phenological Phases

The variance analysis of cycle duration (from VC to R7 stage) in calendar days (CD VC-R7) showed a highly significant effect (*p* < 0.001) of planting date, cultivar, and the interaction of these two factors (Table 6). Block effect and Block X Cultivar interaction were not significant. Total cycle duration showed an amplitude from 100 to 155 days across planting dates, all cultivars combined. Results are available in **Data Sheet 1 (Supplementary Material).**

<table><caption>**TABLE 6** | Analysis of variance of cycle duration from VC to R7 stage in calendar days (CD VC-R7) of 10 cultivars of soybean for five planting dates in outdoor pot experiment (EXP2).</caption><thead><tr><th>Source of variation</th><th>Df</th><th>CD VC-R7</th><th>MS</th></tr></thead><tbody><tr><td>Planting date</td><td>4</td><td></td><td>5118***</td></tr><tr><td>Block</td><td>2</td><td></td><td>48.</td></tr><tr><td>Cultivar</td><td>9</td><td></td><td>900***</td></tr><tr><td>Planting date × Cultivar</td><td>36</td><td></td><td>59***</td></tr><tr><td>Planting date × Block</td><td>8</td><td></td><td>58**</td></tr><tr><td>Block × Cultivar</td><td>18</td><td></td><td>20</td></tr><tr><td>Residuals</td><td>66</td><td></td><td>17</td></tr></tbody></table>

Planting dates analyzed were D1 to D5, as D6 did not reach R7 stage.

(.), (**), and (***), significant at 0.1, 0.01, 0.001 levels, respectively.

Df: Degree of freedom; MS: Mean Square.

## Simple Phenology Algorithm Parameterization of SPA

The parameterization led to satisfactory prediction of the dates of stages. **Figure 2** shows predicted and observed occurrence (day of year) of R7 stage for planting dates D2 to D5. Despite a slight overestimation for D2, predicted values were close to those observed, and the quality of adjustment was satisfactory with an efficiency of 81% and an RMSE of 6.42 days for the appearance of R7 stage since the day of planting. Parameter values are presented in **Tables 5, 6, 7**.

### When Does Photoperiod and Temperature Affect Soybean Phenology?

The parameterization of SPA made it possible to highlight the phases impacted by limiting temperatures or photoperiods for each cultivar. As an example, **Figure 3** shows individually the effects of f(T) and f(P) and their multiplicative effects on the

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Comparison of fitted (FIT) and observed (OBS) day of appearance of R7 stage, for algorithm parameterization on outdoor pot experiment data (EXP2). Colors represent the four planting dates D2 to D5.</figcaption></figure>
<table><caption>**TABLE 7** | Parameters used for each phase. The temperatures T0, Topt, and Tmax were those determined during EXP2 (Table 5).</caption><thead><tr><th rowspan="2">Cultivar</th><th rowspan="2">GM</th><th rowspan="2">Pcrt (h)</th><th rowspan="2">Popt (h)</th><th rowspan="2">S</th><th colspan="4">PDDopt<sub>cp</sub> (d)</th></tr><tr><th>VC-R7</th><th>VC-R1</th><th>R1-R5</th><th>R5-R7</th></tr></thead><tbody><tr><td>Klaxon</td><td>000<br>(0)</td><td>20,5</td><td>13,0</td><td>1,29</td><td>41,9</td><td>11,4</td><td>12,7</td><td>19,2</td></tr><tr><td>RGT_Shouna</td><td>000</td><td>20,0</td><td>13,0</td><td>1,16</td><td>42,8</td><td>10,4</td><td>13,2</td><td>19,9</td></tr><tr><td>Sultana</td><td>000</td><td>20,0</td><td>13,0</td><td>1,18</td><td>41,9</td><td>10,5</td><td>11,9</td><td>20,4</td></tr><tr><td>ES_Mentor</td><td>00</td><td>19,50</td><td>12,50</td><td>1,47</td><td>31,3</td><td>7,1</td><td>9,8</td><td>15,3</td></tr><tr><td>Sigalia</td><td>00</td><td>19,50</td><td>12,50</td><td>1,06</td><td>42,8</td><td>10,0</td><td>11,6</td><td>21,9</td></tr><tr><td>ES_Pallador</td><td>I</td><td>19,00</td><td>12,50</td><td>1,00</td><td>32,7</td><td>8,3</td><td>11,3</td><td>13,6</td></tr><tr><td>Isidor</td><td>I</td><td>19,00</td><td>12,50</td><td>1,23</td><td>41,9</td><td>10,6</td><td>12,8</td><td>19,0</td></tr><tr><td>Santana</td><td>I/II</td><td>19,0</td><td>12,5</td><td>1,32</td><td>42,9</td><td>10,2</td><td>16,4</td><td>17,2</td></tr><tr><td>Blancas</td><td>II</td><td>18,3</td><td>12,0</td><td>1,50</td><td>42,1</td><td>10,3</td><td>15,7</td><td>16,7</td></tr><tr><td>Ecudor</td><td>II</td><td>18,3</td><td>12,0</td><td>1,35</td><td>46,3</td><td>10,1</td><td>16,6</td><td>18,9</td></tr></tbody></table>

Pcrt and Popt were calculated from the bibliography depending on the maturity group of the cultivars. S and PDDopt<sub>c,p</sub> have been optimized for the VC-R7 phase. For the three intermediate phases (VC-R1, R1-R5, R5-R7), only PDDopt<sub>c,p</sub> changed.

development of CV 'Ecudor' for planting date D3 on the outdoor platform. The photoperiod acts on development at a time when the temperature would not be limiting (f[T] ~ 1) around DOY 170, the value of f(P) being at its minimum at ~0.4. For earlier planting around DOY 80, soybean development would be more impacted by temperature with a value of f(T) close to 0.4 and a f (P) of 0.8.

Sensitivity coefficient (S—Table 7) ranged from 1 for 'ES_Pallador' to 1.50 for 'Blancas.' However, this sensitivity was also impacted by Popt and Pcrt, which were determined according to the maturity group. Figure 4 shows the photoperiod function for all cultivars with this calibration. Overall, there is a

trend for the response to photoperiod to decrease from late to early cultivars.

## Evaluation of SPA With Independent Field Data

Figure 5A shows the evaluation results for VE-R7 phase and all years combined. The algorithm predicted the occurrence of the R7 stage with 5.6 days of error on average and a prediction efficiency of 94%. The prediction was satisfactory, even for early planting resulting in a date of appearance of stage R7 outside the conventional range, here before calendar day 210. Data were a little more scattered during 2014 testing year. The same analysis was performed by phases, and the results are shown in Figure 5B. Predictions by phases gave variable levels of performance, but always lower than for the entire cycle. Performance indicators for each phase are reported in the table below the two graphs. VE-R1 phase tended to be underestimated for short phase length (before calendar day 160) while R1-R5 phase was the only stage overestimated by SPA, with a MBE of 4.02 days.

## DISCUSSION

The aim of this study was to develop and apply a phenology phenotyping method coupled with a simple phenology algorithm for soybean cultivars. Phenotyping was divided into two complementary protocols to determine a) cardinal temperatures (EXP1) and b) development response to temperature and photoperiod (EXP2). Data produced in these experiments were used to parameterize SPA (Simple Phenology Algorithm), which was evaluated with independent data from field experiments.

<figure><img src="image_4.png" /><figcaption>**FIGURE 3** | Graphical representation of the temperature function f(T), dotted blue line; photoperiod function f(P), dashed yellow line; and their multiplicative effect f(T) × f(P), solid green line, with parameterization parameters of 'Ecudor' cultivar for planting date D3 (10th may). The planting date as well as the stages R1, R5, and R7 are reported on the abscissa axis.</figcaption></figure>
<table><caption>**TABLE 7** | Parameters used for each phase. The temperatures T0, Topt, and Tmax were those determined during EXP2 (Table 5).</caption><thead><tr><th rowspan="2">Cultivar</th><th rowspan="2">GM</th><th rowspan="2">Pcrt (h)</th><th rowspan="2">Popt (h)</th><th rowspan="2">S</th><th colspan="4">PDDopt<sub>cp</sub> (d)</th></tr><tr><th>VC-R7</th><th>VC-R1</th><th>R1-R5</th><th>R5-R7</th></tr></thead><tbody><tr><td>Klaxon</td><td>000<br>(0)</td><td>20,5</td><td>13,0</td><td>1,29</td><td>41,9</td><td>11,4</td><td>12,7</td><td>19,2</td></tr><tr><td>RGT_Shouna</td><td>000</td><td>20,0</td><td>13,0</td><td>1,16</td><td>42,8</td><td>10,4</td><td>13,2</td><td>19,9</td></tr><tr><td>Sultana</td><td>000</td><td>20,0</td><td>13,0</td><td>1,18</td><td>41,9</td><td>10,5</td><td>11,9</td><td>20,4</td></tr><tr><td>ES_Mentor</td><td>00</td><td>19,50</td><td>12,50</td><td>1,47</td><td>31,3</td><td>7,1</td><td>9,8</td><td>15,3</td></tr><tr><td>Sigalia</td><td>00</td><td>19,50</td><td>12,50</td><td>1,06</td><td>42,8</td><td>10,0</td><td>11,6</td><td>21,9</td></tr><tr><td>ES_Pallador</td><td>I</td><td>19,00</td><td>12,50</td><td>1,00</td><td>32,7</td><td>8,3</td><td>11,3</td><td>13,6</td></tr><tr><td>Isidor</td><td>I</td><td>19,00</td><td>12,50</td><td>1,23</td><td>41,9</td><td>10,6</td><td>12,8</td><td>19,0</td></tr><tr><td>Santana</td><td>I/II</td><td>19,0</td><td>12,5</td><td>1,32</td><td>42,9</td><td>10,2</td><td>16,4</td><td>17,2</td></tr><tr><td>Blancas</td><td>II</td><td>18,3</td><td>12,0</td><td>1,50</td><td>42,1</td><td>10,3</td><td>15,7</td><td>16,7</td></tr><tr><td>Ecudor</td><td>II</td><td>18,3</td><td>12,0</td><td>1,35</td><td>46,3</td><td>10,1</td><td>16,6</td><td>18,9</td></tr></tbody></table>

Pcrt and Popt were calculated from the bibliography depending on the maturity group of the cultivars. S and PDDopt<sub>c,p</sub> have been optimized for the VC-R7 phase. For the three intermediate phases (VC-R1, R1-R5, R5-R7), only PDDopt<sub>c,p</sub> changed.

development of CV 'Ecudor' for planting date D3 on the outdoor platform. The photoperiod acts on development at a time when the temperature would not be limiting (f[T] ~ 1) around DOY 170, the value of f(P) being at its minimum at ~0.4. For earlier planting around DOY 80, soybean development would be more impacted by temperature with a value of f(T) close to 0.4 and a f (P) of 0.8.

Sensitivity coefficient (S—Table 7) ranged from 1 for 'ES_Pallador' to 1.50 for 'Blancas.' However, this sensitivity was also impacted by Popt and Pcrt, which were determined according to the maturity group. Figure 4 shows the photoperiod function for all cultivars with this calibration. Overall, there is a

trend for the response to photoperiod to decrease from late to early cultivars.

## Evaluation of SPA With Independent Field Data

**Figure 5A** shows the evaluation results for VE-R7 phase and all years combined. The algorithm predicted the occurrence of the R7 stage with 5.6 days of error on average and a prediction efficiency of 94%. The prediction was satisfactory, even for early planting resulting in a date of appearance of stage R7 outside the conventional range, here before calendar day 210. Data were a little more scattered during 2014 testing year. The same analysis was performed by phases, and the results are shown in **Figure 5B**. Predictions by phases gave variable levels of performance, but always lower than for the entire cycle. Performance indicators for each phase are reported in the table below the two graphs. VE-R1 phase tended to be underestimated for short phase length (before calendar day 160) while R1-R5 phase was the only stage overestimated by SPA, with a MBE of 4.02 days.

## DISCUSSION

The aim of this study was to develop and apply a phenology phenotyping method coupled with a simple phenology algorithm for soybean cultivars. Phenotyping was divided into two complementary protocols to determine a) cardinal temperatures (EXP1) and b) development response to temperature and photoperiod (EXP2). Data produced in these experiments were used to parameterize SPA (Simple Phenology Algorithm), which was evaluated with independent data from field experiments.

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Graphical representation of the temperature function f(T), dotted blue line; photoperiod function f(P), dashed yellow line; and their multiplicative effect f(T) × f(P), solid green line, with parameterization parameters of 'Ecudor' cultivar for planting date D3 (10th may). The planting date as well as the stages R1, R5, and R7 are reported on the abscissa axis.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**FIGURE 4** | Photoperiod function representation for all cultivars using the calibration parameters. Maturity groups of the cultivars are the following: 000: 'Klaxon,' 'Sultana,' 'RGT_Shouna'; 00: 'Sigalia,' 'ES_Mentor'; I: 'ES_Pallador,' 'Isidor'; I/II: 'Sultana'; II: 'Blancas,' 'Ecudor'.</figcaption></figure>

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

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

<table><caption>**FIGURE 5** | Observed (OBS) and predicted (SIM) days after planting for maturity (R7), all cultivars combined (A) and for three successive phenological phases (B). In graph A, colors indicate years from 2013 to 2018 on data produced at INRA Auzeville, En Crambade, Rivières, and Mondonville. In graph B, colors indicate dates appearance of R1, R5, and R7 with simulations beginning respectively at VE, R1, and R5. Line 1:1 is reported on both graphs A and B. Table below shows the indicators of performance of SPA during its field evaluation.</caption><thead><tr><td>Phase</td><td>MBE (d)</td><td>RMSE (d)</td><td>RRMSE</td><td>EF (%)</td></tr></thead><tbody><tr><td>VE-R7</td><td>-0,70</td><td>5,61</td><td>0,05</td><td>94</td></tr><tr><td>VE-R1</td><td>-0,71</td><td>8,22</td><td>0,17</td><td>75</td></tr><tr><td>R1-R5</td><td>4,02</td><td>9,40</td><td>0,13</td><td>80</td></tr><tr><td>R5-R7</td><td>-3,54</td><td>7,37</td><td>0,06</td><td>90</td></tr></tbody></table>

## Is the Maturity Group Sufficient to Characterize the Phenology of a Cultivar?

The 10 cultivars tested showed differences in phenology in response to the planting dates tested in EXP2. To compare cultivar cycle lengths in a way familiar to breeders, we also calculated cumulative thermal times based on cardinal temperatures and f(T). The analysis of variance showed a

significant effect of planting date, cultivars, and interaction of both on thermal time accumulation for $p < 0.001$. **Table 8** lists the cumulative thermal time to reach R7 for each planting date and cultivar. Differences in thermal time accumulation across cultivars (cultivar amplitude) decreased when delaying planting date; it reached its maximum at date D2 and its lowest at D5. The MG classification of cultivars in Europe is usually performed
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Photoperiod function representation for all cultivars using the calibration parameters. Maturity groups of the cultivars are the following: 000: 'Klaxon,' 'Sultana,' 'RGT_Shouna'; 00: 'Sigalia,' 'ES_Mentor'; I: 'ES_Pallador,' 'Isidor'; I/II: 'Sultana'; II: 'Blancas,' 'Ecudor'.</figcaption></figure>

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

<table><caption>**FIGURE 5** | Observed (OBS) and predicted (SIM) days after planting for maturity (R7), all cultivars combined (A) and for three successive phenological phases (B). In graph A, colors indicate years from 2013 to 2018 on data produced at INRA Auzeville, En Crambade, Rivières, and Mondonville. In graph B, colors indicate dates appearance of R1, R5, and R7 with simulations beginning respectively at VE, R1, and R5. Line 1:1 is reported on both graphs A and B. Table below shows the indicators of performance of SPA during its field evaluation.</caption><thead><tr><td>Phase</td><td>MBE (d)</td><td>RMSE (d)</td><td>RRMSE</td><td>EF (%)</td></tr></thead><tbody><tr><td>VE-R7</td><td>-0,70</td><td>5,61</td><td>0,05</td><td>94</td></tr><tr><td>VE-R1</td><td>-0,71</td><td>8,22</td><td>0,17</td><td>75</td></tr><tr><td>R1-R5</td><td>4,02</td><td>9,40</td><td>0,13</td><td>80</td></tr><tr><td>R5-R7</td><td>-3,54</td><td>7,37</td><td>0,06</td><td>90</td></tr></tbody></table>

## Is the Maturity Group Sufficient to Characterize the Phenology of a Cultivar?

The 10 cultivars tested showed differences in phenology in response to the planting dates tested in EXP2. To compare cultivar cycle lengths in a way familiar to breeders, we also calculated cumulative thermal times based on cardinal temperatures and f(T). The analysis of variance showed a

significant effect of planting date, cultivars, and interaction of both on thermal time accumulation for $p < 0.001$. **Table 8** lists the cumulative thermal time to reach R7 for each planting date and cultivar. Differences in thermal time accumulation across cultivars (cultivar amplitude) decreased when delaying planting date; it reached its maximum at date D2 and its lowest at D5. The MG classification of cultivars in Europe is usually performed
<table><caption>TABLE 8 | Thermal time accumulation from VC (unifoliate leaves) to R7 (beginning maturity stage), by cultivar and planting date.</caption><thead><tr><th>Maturity group</th><th>Cultivar</th><th>D1</th><th>D2</th><th>D3</th><th>D4</th><th>D5</th><th>Mean</th><th>Planting date amplitude</th></tr></thead><tbody><tr><td>000(0)</td><td>Klaxon</td><td>1691</td><td>1717</td><td>1881</td><td>1794</td><td>1513</td><td>1719 d</td><td>367</td></tr><tr><td>000</td><td>RGT_Shouna</td><td>1787</td><td>1891</td><td>2174</td><td>1845</td><td>1518</td><td>1843 cd</td><td>656</td></tr><tr><td>000</td><td>Sultana</td><td>1769</td><td>1818</td><td>2143</td><td>1735</td><td>1514</td><td>1796 d</td><td>630</td></tr><tr><td>00</td><td>ES_Mentor</td><td>1311</td><td>1412</td><td>1705</td><td>1280</td><td>1063</td><td>1354 e</td><td>642</td></tr><tr><td>00</td><td>Sigalia</td><td>2360</td><td>2544</td><td>2609</td><td>2088</td><td>1767</td><td>2274 a</td><td>842</td></tr><tr><td>I</td><td>ES_Pallador</td><td>1993</td><td>2007</td><td>2294</td><td>1729</td><td>1407</td><td>1886 c</td><td>886</td></tr><tr><td>I</td><td>Isidor</td><td>2265</td><td>2210</td><td>2414</td><td>1831</td><td>1537</td><td>2051 b</td><td>877</td></tr><tr><td>I/II</td><td>Santana</td><td>2134</td><td>2346</td><td>2296</td><td>1829</td><td>1598</td><td>2041 b</td><td>748</td></tr><tr><td>II</td><td>Blancas</td><td>2415</td><td>2670</td><td>2487</td><td>2193</td><td>1801</td><td>2313 a</td><td>869</td></tr><tr><td>II</td><td>Ecudor</td><td>2503</td><td>2591</td><td>2525</td><td>2153</td><td>1798</td><td>2314 a</td><td>793</td></tr><tr><td></td><td>Mean</td><td>2023 c</td><td>2121 b</td><td>2253 a</td><td>1848 d</td><td>1552 e</td><td></td><td></td></tr><tr><td></td><td>Cultivar amplitude</td><td>1192</td><td>1258</td><td>904</td><td>913</td><td>737</td><td></td><td></td></tr></tbody></table>

Values followed by common letter(s) are not significantly different at *p* < 0.001 level by SNK test. Planting date and genotypic amplitude are the calculation of maximum minus minimum values from each line and each column, respectively.

using D3 planting date (10th may). As the cultivars are normally not tested for late planting (D5), which would be useful in double cropping decisions, breeders may choose to recommend an early and less productive cultivar while a later one could take advantage of the effect of the photoperiod to reach maturity. We showed in this study that with a planting date around 20th July (D5), a cultivar from MG I could reach beginning maturity (R7) at the same cumulative thermal time than a cultivar from MG 000 (Table 8). This behavior was neither expected nor yet studied in our area. It would be interesting to characterize the cultivars not by their cycle length in days for a given planting date, but by their PDDoptc,p, which would allow breeders to have another criteria for selection. Subsequently, cultivar suitability maps could be developed and characterized, as described by Kurasch et al. (2017). In addition, photoperiod sensitivity of cultivars appears variable among maturity groups. Indeed, two cultivars of the same maturity group can be classified into two different groups according to the SNK test when considering the cumulative thermal time needed to complete the cycle.

## Simple Phenotyping Method

EXP1 was used to determine the cardinal temperatures of cultivars. As presented earlier, fitting the non-linear function described by Yin et al. (1995) led to setting T0 = 0°C for all cultivars tested. However, by fitting a linear function on our data, we found T0 values between 4.5 and 6.1°C (data not shown), which agree with Covell et al. (1986), who also used a linear model. During EXP1 all cultivars germinated at 3°C, indicating that soybean can develop at lower temperature than reported T0 values based on linear fitting. It also illustrates that T0 is not truly a physiological parameter but rather influenced by the fitting method used. Taking into account the observed values, we therefore decided to set T0 at 2°C to get closer to biological reality. Topt, which exhibited a significant difference among cultivars, did not follow any MG-related order. A cultivar-dependent value of Topt has not been reported to our knowledge for soybean. According to Parent and Tardieu (2012), Topt would depend only on the cultivated species,

these authors finding no differences within maize, wheat, and rice cultivars and lines. In order to confirm this finding, it would be interesting to study the response to temperature for a different process, for example the development phase VC-V3, as used by Parent and Tardieu (2012) in their study. The determined Tmax are in the range of what has already been published (Salem et al., 2007; Setiyono et al., 2007). We carried out a supplementary analysis with an nls model performed with R software, with or without genotypic effect in order to test the genotypic effect of cardinal temperatures. The values obtained according to Aikake's information criterion were -3417 and -3323 respectively, indicating the interest in keeping the information at the cultivar level. In addition, genotypic differences for the different parameters were highly significant (*p* < 0.001, Table 5). The calculation of Topt and Tmax by germination percentage and their comparison by a post-hoc test gave the same rankings as the method initially used. The method we used (fitting a model per cultivar) has the advantage of determining the cardinal temperatures of a cultivar independently of the other cultivars, which can be useful for phenotyping a new cultivar. Using the cardinal temperatures determined during the germination experiment for the entire cycle was preferred for the sake of simplification (experimentation, modeling). However, it was theoretically justified by the work of Parent and Tardieu (2012); these authors showed that the cardinal temperatures remained essentially the same regardless of the development process (at least until flowering). By extension, we have applied it to the different phases of the cycle, including the post-flowering period. This assumption was not refuted when applying the model to the different phases then comparing the resulting simulations to field observations. However, other models (e.g. Setiyono et al., 2007) use different cardinal temperatures depending on the phases of the cycle. It may be possible to improve phase prediction by doing this with SPA but this would need specific experiments in controlled conditions or very wide experimental field network to determine them.

Using EXP2, we have phenotyped the response to photoperiod of contrasting cultivars grown in Europe for the
<table><caption>TABLE 8 | Thermal time accumulation from VC (unifoliate leaves) to R7 (beginning maturity stage), by cultivar and planting date.</caption><thead><tr><th>Maturity group</th><th>Cultivar</th><th>D1</th><th>D2</th><th>D3</th><th>D4</th><th>D5</th><th>Mean</th><th>Planting date amplitude</th></tr></thead><tbody><tr><td>000(0)</td><td>Klaxon</td><td>1691</td><td>1717</td><td>1881</td><td>1794</td><td>1513</td><td>1719 d</td><td>367</td></tr><tr><td>000</td><td>RGT_Shouna</td><td>1787</td><td>1891</td><td>2174</td><td>1845</td><td>1518</td><td>1843 cd</td><td>656</td></tr><tr><td>000</td><td>Sultana</td><td>1769</td><td>1818</td><td>2143</td><td>1735</td><td>1514</td><td>1796 d</td><td>630</td></tr><tr><td>00</td><td>ES_Mentor</td><td>1311</td><td>1412</td><td>1705</td><td>1280</td><td>1063</td><td>1354 e</td><td>642</td></tr><tr><td>00</td><td>Sigalia</td><td>2360</td><td>2544</td><td>2609</td><td>2088</td><td>1767</td><td>2274 a</td><td>842</td></tr><tr><td>I</td><td>ES_Pallador</td><td>1993</td><td>2007</td><td>2294</td><td>1729</td><td>1407</td><td>1886 c</td><td>886</td></tr><tr><td>I</td><td>Isidor</td><td>2265</td><td>2210</td><td>2414</td><td>1831</td><td>1537</td><td>2051 b</td><td>877</td></tr><tr><td>I/II</td><td>Santana</td><td>2134</td><td>2346</td><td>2296</td><td>1829</td><td>1598</td><td>2041 b</td><td>748</td></tr><tr><td>II</td><td>Blancas</td><td>2415</td><td>2670</td><td>2487</td><td>2193</td><td>1801</td><td>2313 a</td><td>869</td></tr><tr><td>II</td><td>Ecudor</td><td>2503</td><td>2591</td><td>2525</td><td>2153</td><td>1798</td><td>2314 a</td><td>793</td></tr><tr><td></td><td>Mean</td><td>2023 c</td><td>2121 b</td><td>2253 a</td><td>1848 d</td><td>1552 e</td><td></td><td></td></tr><tr><td></td><td>Cultivar amplitude</td><td>1192</td><td>1258</td><td>904</td><td>913</td><td>737</td><td></td><td></td></tr></tbody></table>

Values followed by common letter(s) are not significantly different at *p* < 0.001 level by SNK test. Planting date and genotypic amplitude are the calculation of maximum minus minimum values from each line and each column, respectively.

using D3 planting date (10th may). As the cultivars are normally not tested for late planting (D5), which would be useful in double cropping decisions, breeders may choose to recommend an early and less productive cultivar while a later one could take advantage of the effect of the photoperiod to reach maturity. We showed in this study that with a planting date around 20th July (D5), a cultivar from MG I could reach beginning maturity (R7) at the same cumulative thermal time than a cultivar from MG 000 (Table 8). This behavior was neither expected nor yet studied in our area. It would be interesting to characterize the cultivars not by their cycle length in days for a given planting date, but by their PDDopt<sub>c,p</sub>, which would allow breeders to have another criteria for selection. Subsequently, cultivar suitability maps could be developed and characterized, as described by Kurasch et al. (2017). In addition, photoperiod sensitivity of cultivars appears variable among maturity groups. Indeed, two cultivars of the same maturity group can be classified into two different groups according to the SNK test when considering the cumulative thermal time needed to complete the cycle.

## Simple Phenotyping Method

EXP1 was used to determine the cardinal temperatures of cultivars. As presented earlier, fitting the non-linear function described by Yin et al. (1995) led to setting T0 = 0°C for all cultivars tested. However, by fitting a linear function on our data, we found T0 values between 4.5 and 6.1°C (data not shown), which agree with Covell et al. (1986), who also used a linear model. During EXP1 all cultivars germinated at 3°C, indicating that soybean can develop at lower temperature than reported T0 values based on linear fitting. It also illustrates that T0 is not truly a physiological parameter but rather influenced by the fitting method used. Taking into account the observed values, we therefore decided to set T0 at 2°C to get closer to biological reality. Topt, which exhibited a significant difference among cultivars, did not follow any MG-related order. A cultivar-dependent value of Topt has not been reported to our knowledge for soybean. According to Parent and Tardieu (2012), Topt would depend only on the cultivated species,

these authors finding no differences within maize, wheat, and rice cultivars and lines. In order to confirm this finding, it would be interesting to study the response to temperature for a different process, for example the development phase VC-V3, as used by Parent and Tardieu (2012) in their study. The determined Tmax are in the range of what has already been published (Salem et al., 2007; Setiyono et al., 2007). We carried out a supplementary analysis with an nls model performed with R software, with or without genotypic effect in order to test the genotypic effect of cardinal temperatures. The values obtained according to Aikake's information criterion were -3417 and -3323 respectively, indicating the interest in keeping the information at the cultivar level. In addition, genotypic differences for the different parameters were highly significant (*p* < 0.001, Table 5). The calculation of Topt and Tmax by germination percentage and their comparison by a post-hoc test gave the same rankings as the method initially used. The method we used (fitting a model per cultivar) has the advantage of determining the cardinal temperatures of a cultivar independently of the other cultivars, which can be useful for phenotyping a new cultivar. Using the cardinal temperatures determined during the germination experiment for the entire cycle was preferred for the sake of simplification (experimentation, modeling). However, it was theoretically justified by the work of Parent and Tardieu (2012); these authors showed that the cardinal temperatures remained essentially the same regardless of the development process (at least until flowering). By extension, we have applied it to the different phases of the cycle, including the post-flowering period. This assumption was not refuted when applying the model to the different phases then comparing the resulting simulations to field observations. However, other models (e.g. Setiyono et al., 2007) use different cardinal temperatures depending on the phases of the cycle. It may be possible to improve phase prediction by doing this with SPA but this would need specific experiments in controlled conditions or very wide experimental field network to determine them.

Using EXP2, we have phenotyped the response to photoperiod of contrasting cultivars grown in Europe for the
developmental phase VC-R7. Phenotyping phenology in a pot experiment over the complete life cycle and including several planting dates is not common. In most cases, under these conditions, experiments are limited to a maximum of R1 stage (Wang et al., 1998; Sinclair et al., 2005). This experiment made it possible to position the soybean cycle over day length ranging from 9 h (D5) to 16.5 h (D2 and D3). Pot cultivation simplified stage ratings, which could have been laborious in the field. The tested density of 5 plants per pot of 28 cm diameter corresponds to a density of 83 plants.m⁻² uniformly distributed on the platform. For comparison, in the field, the sowing rate is close to 50 plants.m⁻² spaced 50 cm apart and therefore arranged tighter on the same row. We increased the density of plants in the pots mainly to increase the number of repetitions per cultivar while maintaining a density that does not have a known negative impact on the agronomic performance of the crop in the field (Edwards and Purcell, 2005). Certainly, EXP2 did not allow soybean to be exposed to continuous optimal photoperiod and temperature conditions throughout the cycle, which would have made it possible to directly determine PDDoptc,p of each phase, therefore not requiring optimization. Experiments in a growth chamber would allow collecting this information, but this requires the necessary equipment and such a long cycle (VC-R7) has never been tested in this type of experiments to our knowledge.

## Simple Algorithm of Phenology

The parameterization of the algorithm gave satisfactory results, even if the early planting date D2 cycle length was overestimated by 5.5 days on average (data not shown). The performance during evaluation was satisfactory for VE-R7 phase, with an efficiency of 94% and a RMSE of 5.61 days. This compares well with the performance of the more detailed and harder to parameterize approach of Setiyono et al. (2007) (RMSE of 3 days) and soil-crop models like STICS and CropGro (RMSE ~10 days) (Boote et al., 1998b; Sojamip research project, http://www6.inra.fr/sojamip). These results show that it is possible to use a semi-controlled phenotyping method such as an automated outdoor platform to parameterize an algorithm that can be used to predict phenological stages in the field. This could be useful as a plant breeding tool. Because the algorithm has been developed based on a limited range of MGs, it is not possible to check its suitability for later MGs with our data. Nevertheless, maturity groups tested correspond to most MGs grown in Europe.

The duration of R1-R5 phase was overestimated for early plantings by 4 days on average. Longer photoperiods were met during this stage, the summer solstice taking place on DOY 172 with 16.5 h a day at a latitude of 43°N. This means that the photoperiod function f(P) applied had a strong impact on the development rate (f(P) values around 0.4) on this phase. In SPA, we chose to use a non-linear photoperiod function, unlike other models such as STICS (Brisson et al., 2003). The effect of f(P) in SPA is therefore greater for photoperiods close to Pcrt for long photoperiods compared to a linear formalism. Popt and Pcrt have been adopted from Setiyono et al. (2007), but it would

be desirable to re-evaluate these thresholds with experiments such as the one presented by Daba et al. (2016). The optimization of the photoperiod sensitivity coefficient (S) made it possible to highlight differences in sensitivity within the same maturity group, which have the same Popt and Pcrt. This element must be taken into consideration by breeders when adapting cultivars to new management practices or new latitudes.

To test if the approach could be simplified by replacing cultivar-specific Topt by soybean-specific Topt, we calculated an average value of Topt across the 10 cultivars then re-estimated the values of PDDoptc,p, and S parameters. The prediction quality of VC-R7 duration was compared at field level (75 data points) for the two sets of parameters, either cultivar- or soybean-specific (Table C—Supplementary Material). RMSE (6.20 vs. 5.61 days) and MBE (-2.42 vs. -0.70) indicators performances were lessened when using soybean- instead of cultivar-specific Topt values respectively. Therefore using a cultivar-specific calibration in this study was justified. However, due to the slight decrease of the performance of SPA with average Topt, end-users could prefer adopt a soybean-specific Topt value for a direct and simpler application of the phenological algorithm.

One important goal of this study was to develop and evaluate a simple algorithm of phenology. Its simplicity comes from the low number (5) of parameters (photoperiodic thresholds and cardinal temperatures), observed or adopted from literature, plus two additional parameters obtained by optimization (PDDoptc,p and S) based on observed phenology data. Other models can predict phenology well, but they are more difficult to parameterize (Setiyono et al., 2007) or they reach up to R1 only although with more detailed analysis (Sinclair et al., 2005).

## CONCLUSION

In this study we provide a path for soybean cultivar phenology × environment × crop management (sowing date) advice by: (1) developing a simple algorithm of phenology; and (2) proposing a simple method in semi-controlled conditions for phenotyping cultivar phenology and determining parameters for the algorithm. SPA allows users to quantify the effects of temperature and photoperiod on soybean phenology that affect the possible extension of the crop to northern Europe. It will also help to define new cultivation areas and management practices (i.e. early or late sowing) according to cultivar characteristics in the context of climate change. Because SPA is based on observations, it is possible to parameterize SPA directly from phenology data produced in field experiments, enhancing the value of these experiments and assisting plant selection efforts. For the same reason, SPA can be parameterized for other crop species without modification. Finally, integration of SPA into cropping system models can enhance the performance of these models for assessing the adaptation of cultivars to climate change scenarios and changing management options.
developmental phase VC-R7. Phenotyping phenology in a pot experiment over the complete life cycle and including several planting dates is not common. In most cases, under these conditions, experiments are limited to a maximum of R1 stage (Wang et al., 1998; Sinclair et al., 2005). This experiment made it possible to position the soybean cycle over day length ranging from 9 h (D5) to 16.5 h (D2 and D3). Pot cultivation simplified stage ratings, which could have been laborious in the field. The tested density of 5 plants per pot of 28 cm diameter corresponds to a density of 83 plants.m⁻² uniformly distributed on the platform. For comparison, in the field, the sowing rate is close to 50 plants.m⁻² spaced 50 cm apart and therefore arranged tighter on the same row. We increased the density of plants in the pots mainly to increase the number of repetitions per cultivar while maintaining a density that does not have a known negative impact on the agronomic performance of the crop in the field (Edwards and Purcell, 2005). Certainly, EXP2 did not allow soybean to be exposed to continuous optimal photoperiod and temperature conditions throughout the cycle, which would have made it possible to directly determine PDDoptc,p of each phase, therefore not requiring optimization. Experiments in a growth chamber would allow collecting this information, but this requires the necessary equipment and such a long cycle (VC-R7) has never been tested in this type of experiments to our knowledge.

## Simple Algorithm of Phenology

The parameterization of the algorithm gave satisfactory results, even if the early planting date D2 cycle length was overestimated by 5.5 days on average (data not shown). The performance during evaluation was satisfactory for VE-R7 phase, with an efficiency of 94% and a RMSE of 5.61 days. This compares well with the performance of the more detailed and harder to parameterize approach of Setiyono et al. (2007) (RMSE of 3 days) and soil-crop models like STICS and CropGro (RMSE ~10 days) (Boote et al., 1998b; Sojamip research project, http://www6.inra.fr/sojamip). These results show that it is possible to use a semi-controlled phenotyping method such as an automated outdoor platform to parameterize an algorithm that can be used to predict phenological stages in the field. This could be useful as a plant breeding tool. Because the algorithm has been developed based on a limited range of MGs, it is not possible to check its suitability for later MGs with our data. Nevertheless, maturity groups tested correspond to most MGs grown in Europe.

The duration of R1-R5 phase was overestimated for early plantings by 4 days on average. Longer photoperiods were met during this stage, the summer solstice taking place on DOY 172 with 16.5 h a day at a latitude of 43°N. This means that the photoperiod function f(P) applied had a strong impact on the development rate (f(P) values around 0.4) on this phase. In SPA, we chose to use a non-linear photoperiod function, unlike other models such as STICS (Brisson et al., 2003). The effect of f(P) in SPA is therefore greater for photoperiods close to Pcrt for long photoperiods compared to a linear formalism. Popt and Pcrt have been adopted from Setiyono et al. (2007), but it would

be desirable to re-evaluate these thresholds with experiments such as the one presented by Daba et al. (2016). The optimization of the photoperiod sensitivity coefficient (S) made it possible to highlight differences in sensitivity within the same maturity group, which have the same Popt and Pcrt. This element must be taken into consideration by breeders when adapting cultivars to new management practices or new latitudes.

To test if the approach could be simplified by replacing cultivar-specific Topt by soybean-specific Topt, we calculated an average value of Topt across the 10 cultivars then re-estimated the values of PDDoptc,p, and S parameters. The prediction quality of VC-R7 duration was compared at field level (75 data points) for the two sets of parameters, either cultivar- or soybean-specific (Table C—Supplementary Material). RMSE (6.20 vs. 5.61 days) and MBE (-2.42 vs. -0.70) indicators performances were lessened when using soybean- instead of cultivar-specific Topt values respectively. Therefore using a cultivar-specific calibration in this study was justified. However, due to the slight decrease of the performance of SPA with average Topt, end-users could prefer adopt a soybean-specific Topt value for a direct and simpler application of the phenological algorithm.

One important goal of this study was to develop and evaluate a simple algorithm of phenology. Its simplicity comes from the low number (5) of parameters (photoperiodic thresholds and cardinal temperatures), observed or adopted from literature, plus two additional parameters obtained by optimization (PDDoptc,p and S) based on observed phenology data. Other models can predict phenology well, but they are more difficult to parameterize (Setiyono et al., 2007) or they reach up to R1 only although with more detailed analysis (Sinclair et al., 2005).

## CONCLUSION

In this study we provide a path for soybean cultivar phenology × environment × crop management (sowing date) advice by: (1) developing a simple algorithm of phenology; and (2) proposing a simple method in semi-controlled conditions for phenotyping cultivar phenology and determining parameters for the algorithm. SPA allows users to quantify the effects of temperature and photoperiod on soybean phenology that affect the possible extension of the crop to northern Europe. It will also help to define new cultivation areas and management practices (i.e. early or late sowing) according to cultivar characteristics in the context of climate change. Because SPA is based on observations, it is possible to parameterize SPA directly from phenology data produced in field experiments, enhancing the value of these experiments and assisting plant selection efforts. For the same reason, SPA can be parameterized for other crop species without modification. Finally, integration of SPA into cropping system models can enhance the performance of these models for assessing the adaptation of cultivars to climate change scenarios and changing management options.
# DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the article/
**Supplementary Material**.

# AUTHOR CONTRIBUTIONS

CS, COS, LC, PD, and PM, designed the research. PD, LC, and PM acquired the funding for the research project. CS and CC performed the experiments. COS designed the algorithm. CS, COS, PD, and PM analyzed the data. CS, COS, PD, and PM wrote the manuscript.

# FUNDING

The authors thank the Occitanie / Pyrénées-Méditerranée region and Terres Inovia for financial support of this research.

# REFERENCES

Alliprandini, L. F., Abatti, C., Bertagnolli, P. F., Cavassim, J. E., Gabe, H. L., Kurek, A., et al. (2009). Understanding soybean maturity groups in brazil: environment, cultivar classification, and stability. *Crop Sci.* 49, 801–808. doi: 10.2135/cropsci2008.07.0390

Andreucci, M. P. P., Moot, D. J. J., Black, A. D. D., and Sedcole, R. (2016). A comparison of cardinal temperatures estimated by linear and nonlinear models for germination and bulb growth of forage brassicas. *Eur. J. Agron.* 81, 52–63. doi: 10.1016/j.eja.2016.08.010

Boote, K. J., Jones, J. W., and Hoogenboom, G. (1998a). "Simulation of crop growth: CROPGRO Model," in *Agricultural Systems Modeling and Simulation* (New York: Marcel Dekker Inc.), 651–692.

Boote, K. J., Jones, J. W., Hoogenboom, G., and Pickering, N. B. (1998b). The CROPGRO model for grain legumes. *Underst. options Agric. Prod.* 1983, 99–128. doi: 10.1007/978-94-017-3624-4_6

Boote, K. J. (2011). "Improving Soybean Cultivars for Adaptation to Climate Change and Climate Variability," in *Crop Adaptation to Climate Change* (Oxford, UK: Wiley-Blackwell), 370–395. doi: 10.1002/9780470960929.ch26

Bowers, G. R. (1995). An early soybean production system for drought avoidance. *J. Prod. Agric.* 8, 112. doi: 10.2134/jpa1995.fpage

Brisson, N., Gary, C., Justes, E., Roche, R., Mary, B., Ripoche, D., et al. (2003). An overview of the crop model STICS. *Eur. J. Agron.* 18, 309–332. doi: 10.1016/S1161-0301(02)00110-7

Caffaro, S. V., Martigmone, R. A., Torres, R., and Nakayama, F. (1988). Photoperiod regulation of vegetative growth and meristem behavior toward flower initiation of an indeterminate Soybean. *Bot. Gaz.* 149, 311–316. doi: 10.1086/337720

Clovis, P. J., Jackson, K., Marcelo, B., Murilo, V. D. C., Leandro, M., Pierozan Junior, C., et al. (2015). Phenological and quantitative plant development changes in soybean cultivars caused by sowing date and their relation to yield. *Afr. J. Agric. Res.* 10, 515–523. doi: 10.5897/AJAR2014.9325

Cober, E., Curtis, D., Stewart, D., and Morrison, M. (2014). Quantifying the effects of photoperiod, temperature and daily irradiance on flowering time of soybean isolines. *Plants* 3, 476–497. doi: 10.3390/plants3040476

Covell, S., Ellis, R. H., Roberts, E. H., and Summerfield, R. J. (1986). The influence of temperature on seed germination rate in grain legumes: I. a comparison of chickpea, lentil, soyabean and cowpea at constant temperatures. *J. Exp. Bot.* 37, 705–715. doi: 10.1093/jxb/37.5.705

# ACKNOWLEDGMENTS

The works presented are part of the “Sojamip” research project (http://www6.inra.fr/sojamip) involving various partners (Euralis Semences, RAGT2n, Caussade Semences, Terres Inovia, Terres Univia, EIP Purpan, Toulouse INP, ENSAT, INRA). We thank Toulouse INP for the invitation of Prof. COS as “Visiting Professor” at UMR AGIR. We would like to thank Nicolas Capoul, Michel Labarrère, Damien Marchand, and Hugo Quintana for their technical assistance in carrying out the agronomic tests and Elie Maza for his expert advice in statistical analysis.

# SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01755/full#supplementary-material

Cregan, P. B., and Hartwig, E. E. (1984). Characterization of flowering response to photoperiod in diverse soybean genotypes. *Crop Sci.* 24, 659–662. doi: 10.2135/cropsci1984.0011183X002400040008x

Daba, K., Warkentin, T. D., Bueckert, R., Todd, C. D., and Tar'an, B. (2016). Determination of photoperiod-sensitive phase in chickpea (Cicer arietinum L.). *Front. Plant Sci.* 7, 478. doi: 10.3389/fpls.2016.00478

Dardanelli, J. L., Balzarini, M., Martínez, M. J., Cuniberti, M., Resnik, S., Ramunda, S. F., et al. (2006). Soybean maturity groups, environments, and their interaction define mega-environments for seed composition in Argentina. *Crop Sci.* 46, 1939–1947. doi: 10.2135/cropsci2005.12.0480

Edwards, J. T., and Purcell, L. C. (2005). Soybean yield and biomass responses to increasing plant population among diverse maturity groups. *Crop Sci.* 45, 1770. doi: 10.2135/cropsci2004.0564

Ellis, R. H., Collinson, S. T., Hudson, D., and Patefield, W. M. (1992). The analysis of reciprocal transfer experiments to estimate the durations of the photoperiod-sensitive and photoperiod-insensitive phases of plant development: An example in soya bean. *Ann. Bot.* 70, 87–92. doi: 10.1093/oxfordjournals.aob.a088443

Fehr, W. R., and Caviness, C. E. (1977). Stages of soybean development. *Iowa Coop. Ext. Serv. Iowa Agric. Home Econ. Exp. Stn. Spec. Rep.* 80, 11.

Garner, W. W., and Allard, H. A. (1920). Effect of the relative length of day and night and other factors of the environment on growth and reproduction in plants. *Mon. Weather Rev.* 48, 415–415. doi: 10.1175/1520-0493(1920)48<415b:EOTRLO>2.0.CO;2

Gosseau, F., Blanchet, N., Varès, D., Burger, P., Campergue, D., Colombet, C., et al. (2019). Heliaphen, an outdoor high-throughput phenotyping platform designed to integrate genetics and crop modeling. *Front. Plant Sci.* 9, 1–26. doi: 10.3389/fpls.2018.01908

Grimm, S. S., Jones, J. W., Boote, K. J., and Hesketh, J. D. (1993). Parameter estimation for predicting flowering date of soybean cultivars. *Crop Sci.* 33, 137–144. doi: 10.2135/cropsci1993.0011183X003300010025x

Hadley, P., Roberts, E. H., Summerfield, R. J., and Minchin, F. R. (1984). Effects of temperature and photoperiod on flowering in soybean [Glycine max (L.) Merrill]: a quantitative model. *Ann. Bot.* 53, 669–681. doi: 10.1093/oxfordjournals.aob.a086732

Han, T., Wu, C., Tong, Z., Mentreddy, R. S., Tan, K., and Gai, J. (2006). Postflowering photoperiod regulates vegetative growth and reproductive development of soybean. *Environ. Exp. Bot.* 55, 120–129. doi: 10.1016/j.envexbot.2004.10.006

Hodges, T., and French, V. (1985). Soyphen: soybean growth stages modeled from temperature, daylength, and water availability. *Agron. J.* 77, 500–505. doi: 10.2134/agronj1985.00021962007700030031x
# DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the article/
**Supplementary Material**.

# AUTHOR CONTRIBUTIONS

CS, COS, LC, PD, and PM, designed the research. PD, LC, and PM acquired the funding for the research project. CS and CC performed the experiments. COS designed the algorithm. CS, COS, PD, and PM analyzed the data. CS, COS, PD, and PM wrote the manuscript.

# FUNDING

The authors thank the Occitanie / Pyrénées-Méditerranée region and Terres Inovia for financial support of this research.

# REFERENCES

Alliprandini, L. F., Abatti, C., Bertagnolli, P. F., Cavassim, J. E., Gabe, H. L., Kurek, A., et al. (2009). Understanding soybean maturity groups in brazil: environment, cultivar classification, and stability. *Crop Sci.* 49, 801–808. doi: 10.2135/cropsci2008.07.0390

Andreucci, M. P. P., Moot, D. J. J., Black, A. D. D., and Sedcole, R. (2016). A comparison of cardinal temperatures estimated by linear and nonlinear models for germination and bulb growth of forage brassicas. *Eur. J. Agron.* 81, 52–63. doi: 10.1016/j.eja.2016.08.010

Boote, K. J., Jones, J. W., and Hoogenboom, G. (1998a). "Simulation of crop growth: CROPGRO Model," in *Agricultural Systems Modeling and Simulation* (New York: Marcel Dekker Inc.), 651–692.

Boote, K. J., Jones, J. W., Hoogenboom, G., and Pickering, N. B. (1998b). The CROPGRO model for grain legumes. *Underst. options Agric. Prod.* 1983, 99–128. doi: 10.1007/978-94-017-3624-4_6

Boote, K. J. (2011). "Improving Soybean Cultivars for Adaptation to Climate Change and Climate Variability," in *Crop Adaptation to Climate Change* (Oxford, UK: Wiley-Blackwell), 370–395. doi: 10.1002/978047 0960929.ch26

Bowers, G. R. (1995). An early soybean production system for drought avoidance. *J. Prod. Agric.* 8, 112. doi: 10.2134/jpa1995.fpage

Brisson, N., Gary, C., Justes, E., Roche, R., Mary, B., Ripoche, D., et al. (2003). An overview of the crop model STICS. *Eur. J. Agron.* 18, 309–332. doi: 10.1016/ S1161-0301(02)00110-7

Caffaro, S. V., Martigmone, R. A., Torres, R., and Nakayama, F. (1988). Photoperiod regulation of vegetative growth and meristem behavior toward flower initiation of an indeterminate Soybean. *Bot. Gaz.* 149, 311–316. doi: 10.1086/337720

Clovis, P. J., Jackson, K., Marcelo, B., Murilo, V. D. C., Leandro, M., Pierozan Junior, C., et al. (2015). Phenological and quantitative plant development changes in soybean cultivars caused by sowing date and their relation to yield. *Afr. J. Agric. Res.* 10, 515–523. doi: 10.5897/AJAR2014. 9325

Cober, E., Curtis, D., Stewart, D., and Morrison, M. (2014). Quantifying the effects of photoperiod, temperature and daily irradiance on flowering time of soybean isolines. *Plants* 3, 476–497. doi: 10.3390/plants3040476

Covell, S., Ellis, R. H., Roberts, E. H., and Summerfield, R. J. (1986). The influence of temperature on seed germination rate in grain legumes: I. a comparison of chickpea, lentil, soyabean and cowpea at constant temperatures. *J. Exp. Bot.* 37, 705–715. doi: 10.1093/jxb/37.5.705

# ACKNOWLEDGMENTS

The works presented are part of the “Sojamip” research project (http://www6.inra.fr/sojamip) involving various partners (Euralis Semences, RAGT2n, Caussade Semences, Terres Inovia, Terres Univia, EIP Purpan, Toulouse INP, ENSAT, INRA). We thank Toulouse INP for the invitation of Prof. COS as “Visiting Professor” at UMR AGIR. We would like to thank Nicolas Capoul, Michel Labarrère, Damien Marchand, and Hugo Quintana for their technical assistance in carrying out the agronomic tests and Elie Maza for his expert advice in statistical analysis.

# SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01755/full#supplementary-material

Cregan, P. B., and Hartwig, E. E. (1984). Characterization of flowering response to photoperiod in diverse soybean genotypes. *Crop Sci.* 24, 659–662. doi: 10.2135/ cropsci1984.0011183X002400040008x

Daba, K., Warkentin, T. D., Bueckert, R., Todd, C. D., and Tar'an, B. (2016). Determination of photoperiod-sensitive phase in chickpea (Cicer arietinum L.). *Front. Plant Sci.* 7, 478. doi: 10.3389/fpls.2016.00478

Dardanelli, J. L., Balzarini, M., Martínez, M. J., Cuniberti, M., Resnik, S., Ramunda, S. F., et al. (2006). Soybean maturity groups, environments, and their interaction define mega-environments for seed composition in Argentina. *Crop Sci.* 46, 1939–1947. doi: 10.2135/cropsci2005.12.0480

Edwards, J. T., and Purcell, L. C. (2005). Soybean yield and biomass responses to increasing plant population among diverse maturity groups. *Crop Sci.* 45, 1770. doi: 10.2135/cropsci2004.0564

Ellis, R. H., Collinson, S. T., Hudson, D., and Patefield, W. M. (1992). The analysis of reciprocal transfer experiments to estimate the durations of the photoperiod-sensitive and photoperiod-insensitive phases of plant development: An example in soya bean. *Ann. Bot.* 70, 87–92. doi: 10.1093/ oxfordjournals.aob.a088443

Fehr, W. R., and Caviness, C. E. (1977). Stages of soybean development. *Iowa Coop. Ext. Serv. Iowa Agric. Home Econ. Exp. Stn. Spec. Rep.* 80, 11.

Garner, W. W., and Allard, H. A. (1920). Effect of the relative length of day and night and other factors of the environment on growth and reproduction in plants. *Mon. Weather Rev.* 48, 415–415. doi: 10.1175/1520-0493(1920) 48<415b:EOTRLO>2.0.CO;2

Gosseau, F., Blanchet, N., Varès, D., Burger, P., Campergue, D., Colombet, C., et al. (2019). Heliaphen, an outdoor high-throughput phenotyping platform designed to integrate genetics and crop modeling. *Front. Plant Sci.* 9, 1–26. doi: 10.3389/fpls.2018.01908

Grimm, S. S., Jones, J. W., Boote, K. J., and Hesketh, J. D. (1993). Parameter estimation for predicting flowering date of soybean cultivars. *Crop Sci.* 33, 137–144. doi: 10.2135/cropsci1993.0011183X003300010025x

Hadley, P., Roberts, E. H., Summerfield, R. J., and Minchin, F. R. (1984). Effects of temperature and photoperiod on flowering in soybean [Glycine max (L.) Merrill]: a quantitative model. *Ann. Bot.* 53, 669–681. doi: 10.1093/ oxfordjournals.aob.a086732

Han, T., Wu, C., Tong, Z., Mentreddy, R. S., Tan, K., and Gai, J. (2006). Postflowering photoperiod regulates vegetative growth and reproductive development of soybean. *Environ. Exp. Bot.* 55, 120–129. doi: 10.1016/j.envexbot.2004.10.006

Hodges, T., and French, V. (1985). Soyphen: soybean growth stages modeled from temperature, daylength, and water availability. *Agron. J.* 77, 500–505. doi: 10.2134/agronj1985.00021962007700030031x
Hu, M., and Wiatrak, P. (2012). Effect of planting date on soybean growth, yield, and grain quality: review. *Agron. J.* 104, 785-790. doi: 10.2134/agronj2011.0382

Hu, M. (2013). Effects of late planting dates, maturity groups and management systems on growth, development and yield of soybean in South. *M.Sc. thesis*, 109, Clemson University, Clemson, South Carolina USA.

IPCC. (2014). Climate Change 2014 synthesis report summary chapter for policymakers. IPCC, 31. doi: 10.1017/CBO9781107415324

Jones, J. W., Hoogenboom, G., Porter, C. H., Boote, K. J., Batchelor, W. D., Hunt, L. A., et al. (2003). "The DSSAT cropping system model," *Eur. J. Agron.* 18, 235-265. doi: 10.1016/S1161-0301(02)00107-7

Kantolic, A. G., and Slafer, G. A. (2001). Photoperiod sensitivity after flowering and seed number determination in indeterminate soybean cultivars. *F. Crop Res.* 72, 109-118. doi: 10.1016/S0378-4290(01)00168-X

Kantolic, A. G., and Slafer, G. A. (2007). Development and seed number in indeterminate soybean as affected by timing and duration of exposure to long photoperiods after flowering. *Ann. Bot.* 99, 925-933. doi: 10.1093/aob/mcm033

Kantolic, A. G., Peralta, G. E., and Slafer, G. A. (2013). Seed number responses to extended photoperiod and shading during reproductive stages in indeterminate soybean. *Eur. J. Agron.* 51, 91-100. doi: 10.1016/j.eja.2013.07.006

Keating, B. A., Carberry, P. S., Hammer, G. L., Probert, M. E., Robertson, M. J., Holzworth, D., et al. (2003). "An overview of APSIM, a model designed for farming systems simulation," *Eur. J. Agro.* 18, 267-288. doi: 10.1016/S1161-0301(02)00108-9

Kurasch, A. K., Hahn, V., Leiser, W. L., Vollmann, J., Schori, A., Bétrix, C. A., et al. (2017). Identification of mega-environments in Europe and effect of allelic variation at maturity E loci on adaptation of European soybean. *Plant Cell Environ.* 40, 765-778. doi: 10.1111/pce.12896

Lu, S., Zhao, X., Hu, Y., Liu, S., Nan, H., Li, X., et al. (2017). Natural variation at the soybean J locus improves adaptation to the tropics and enhances yield. *Nat. Genet.* 49, 773-779. doi: 10.1038/ng.3819

Mahieu, P. J., and Brinkman, M. A. (1990). Double-cropping soybean after harvesting small grains as forage in the north central USA. *J. Prod. Agric.* 3, 385-389. doi: 10.2134/jpa1990.0385

Major, D. J., Johnson, D. R., Tanner, J. W., and Anderson, I. C. (1975). Effects of daylength and temperature on soybean development. *Crop Sci.* 15, 174. doi: 10.2135/cropsci1975.0011183X001500020009x

Maury, P., Andrianasolo, F. N., Alric, F., Berger, M., Beugniet, G., Chambert, C., et al. (2015). Le semis très précoce : une stratégie agronomique pour améliorer les performances du soja en France ? *OCL* 22, D503. doi: 10.1051/ocl/2015028

Mcblain, B. A., and Bernard, R. L. (1987). A new gene affecting the time of flowering and maturity in soybeans. *J. Hered.* 78, 160-162. doi: 10.1093/oxfordjournals.jhered.a110349

McDonald, M. B. B., Vertucci, C. W. W., and Roos, E. E. E. (1988). Soybean seed imbibition: water absorption by seed parts. *Crop Sci.* 28, 993-997. doi: 10.2135/cropsci1988.0011183X002800060026x

Meyer, D. W., and Badaruddin, M. (2001). Frost tolerance of ten seedling legume species at four growth stages. *Crop Sci.* 41, 1838-1842. doi: 10.2135/cropsci2001.1838

Nico, M., Miralles, D. J., and Kantolic, A. G. (2015). Post-flowering photoperiod and radiation interaction in soybean yield determination: direct and indirect photoperiodic effects. *F. Crop Res.* 176, 45-55. doi: 10.1016/j.fcr.2015.02.018

Olesen, J. E., Trnka, M., Kersebaum, K. C., Skjelvåg, A. O., Seguin, B., Peltonen-Sainio, P., et al. (2011). Impacts and adaptation of European crop production systems to climate change. *Eur. J. Agron.* 34, 96-112. doi: 10.1016/j.eja.2010.11.003

Parent, B., and Tardieu, F. (2012). Temperature responses of developmental processes have not been affected by breeding in different ecological areas for 17 crop species. *New Phytol.* 194, 760-774. doi: 10.1111/j.1469-8137.2012.04086.x

Parent, B., Turc, O., Gibon, Y., Stitt, M., and Tardieu, F. (2010). Modelling temperature-compensated physiological rates, based on the co-ordination of responses to temperature of developmental processes. *J. Exp. Bot.* 61, 2057-2069. doi: 10.1093/jxb/erq003

Piper, E. L., Boote, K. J., Jones, J. W., and Grimm, S. S. (1996). Comparison of two phenology models for predicting flowering and maturity date of soybean. *Crop Sci.* 36, 1606-1614. doi: 10.2135/cropsci1996.0011183X003600060033x

Poehlman, J. M. (1987). "Breeding Soybeans," in *Breeding Field Crops* (Dordrecht: Springer Netherlands), 421-450. doi: 10.1007/978-94-015-7271-2_17

Quinsac, A., Lecomte, V., Bourrel, C., Labalette, F., Ledoux, S., and Jouffret, P. (2014). Panorama et futur de la filière du soja français. *Oléagineux Corps. gras Lipides.* 17, 345-355. doi: 10.1051/ocl.2010.0345

Raper, C. D., and Thomas, J. F. (1978). Photoperiodic alteration of dry matter partitioning and seed yield in soybeans. *Crop Sci.* 18, 654-656. doi: 10.2135/cropsci1978.0011183X001800040032x

Robertson, G. W. (1968). A biometeorological time scale for a cereal crop involving day and night temperatures and photoperiod. *Int. J. Biometeorol.* 12, 191-223. doi: 10.1007/BF01553422

Salem, M. A., Kakani, V. G., Koti, S., and Reddy, K. R. (2007). Pollen-based screening of soybean genotypes for high temperatures. *Crop Sci.* 47, 219-231. doi: 10.2135/cropsci2006.07.0443

Salmerón, M., and Purcell, L. C. (2016). Simplifying the prediction of phenology with the DSSAT-CROPGRO-soybean model based on relative maturity group and determinacy. *Agric. Syst.* 148, 178-187. doi: 10.1016/j.agsy.2016.07.016

Salmeron, M., Gbur, E. E., Bourland, F. M., Buehring, N. W., Earnest, L., Fritschi, F. B., et al. (2014). Soybean maturity group choices for early and late plantings in the midsouth. *Agron. J.* 106, 1893-1901. doi: 10.2134/agronj14.0222

Scott, W. O., and Aldrich, S. R. (1970). *Modern soybean production* (Champaign, Illinois: S & A Publications.). Available at: https://www.cabdirect.org/cabdirect/abstract/19790781104 [Accessed March 19, 2018].

Seifert, C. A., and Lobell, D. B. (2015). Response of double cropping suitability to climate change in the United States. *Environ. Res. Lett.* 10 (2), 024002. doi: 10.1088/1748-9326/10/2/024002

Setiyono, T. D., Weiss, A., Specht, J., Bastidas, A. M., Cassman, K. G., and Dobermann, A. (2007). Understanding and modeling the effect of temperature and daylength on soybean phenology under high-yield conditions. *F. Crop Res.* 100, 257-271. doi: 10.1016/j.fcr.2006.07.011

Shapiro, B. I., Brorsen, B. W., and Doster, D. H. (1992). Adoption of double-cropping soybeans and wheat. *South. J. Agric. Econ.* 24, 33-40. doi: 10.1017/S0081305200018355

Sinclair, T. R., Neumaier, N., Farias, J. R. B., and Nepomuceno, A. L. (2005). Comparison of vegetative development in soybean cultivars for low-latitude environments. *F. Crop Res.* 92, 53-59. doi: 10.1016/j.fcr.2004.08.008

Terres Inovia. (2019). Outil IRRIsoja. Available at: https://www.terresinovia.fr/p/irrisoja

Tribouillois, H., Dürr, C., Demilly, D., Wagner, M.-H., and Justes, E. (2016). Determination of germination response to temperature and water potential for a wide range of cover crop species and related functional groups. *PloS One* 11 (8), e0161185. doi: 10.1371/journal.pone.0161185

USDA. (2019). *PSD_Oilseeds. PSD_Oilseeds*. Available at: https://apps.fas.usda.gov/psdonline/app/index.html#/app/downloads.

Wang, E., and Engel, T. (1998). Simulation of phenological development of wheat crops. *Agric. Syst.* 58, 1-24. doi: 10.1016/S0308-521X(98)00028-6

Wang, Z., Reddy, V. R., and Acock, M. C. (1998). Testing for early photoperiod insensitivity in soybean. *Agron. J.* 90, 389-392. doi: 10.2134/agronj1998.00021962009000030012x

Xia, Z., Zhai, H., Liu, B., Kong, F., Yuan, X., Wu, H., et al. (2012). Molecular identification of genes controlling flowering time, maturity, and photoperiod response in soybean. *Plant Syst. Evol.* 298, 1217-1227. doi: 10.1007/s00606-012-0628-2

Xu, M., Xu, Z., Liu, B., Kong, F., Tsubokura, Y., Watanabe, S., et al. (2013). Genetic variation in four maturity genes affects photoperiod insensitivity and PHYA-regulated post-flowering responses of soybean. *BMC Plant Biol.* 13 (1), 91. doi: 10.1186/1471-2229-13-91

Yin, X., Kropff, M. J., McLaren, G., and Visperas, R. M. (1995). A nonlinear model for crop development as a function of temperature. *Agric. For. Meteorol.* 77, 1-16. doi: 10.1016/0168-1923(95)02236-Q

Yuesheng, W., Jianbing, Q., Junyi, G., and Guangyuang, H. (2006). Classification and characteristic of maturity groups of Chinese landraces of soybean [Glycine max (L.) Merr.]. *Genet. Resour. Crop Evol.* 53, 803-809. doi: 10.1007/s10722-004-5731-y

Zhang, L. X., Kyei-Boahen, S., Zhang, J., Zhang, M. H., Freeland, T. B., Watson, C. E., et al. (2007). Modifications of optimum adaptation zones for soybean maturity groups in the USA. *Crop Manage.* 6, 1-10. doi: 10.1094/CM-2007-0927-01RS Abstract
Hu, M., and Wiatrak, P. (2012). Effect of planting date on soybean growth, yield, and grain quality: review. *Agron. J.* 104, 785–790. doi: 10.2134/agronj2011.0382

Hu, M. (2013). Effects of late planting dates, maturity groups and management systems on growth, development and yield of soybean in South. *M.Sc. thesis*, 109, Clemson University, Clemson, South Carolina USA.

IPCC. (2014). Climate Change 2014 synthesis report summary chapter for policymakers. IPCC, 31. doi: 10.1017/CBO9781107415324

Jones, J. W., Hoogenboom, G., Porter, C. H., Boote, K. J., Batchelor, W. D., Hunt, L. A., et al. (2003). "The DSSAT cropping system model," *Eur. J. Agron.* 18, 235–265. doi: 10.1016/S1161-0301(02)00107-7

Kantolic, A. G., and Slafer, G. A. (2001). Photoperiod sensitivity after flowering and seed number determination in indeterminate soybean cultivars. *F. Crop Res.* 72, 109–118. doi: 10.1016/S0378-4290(01)00168-X

Kantolic, A. G., and Slafer, G. A. (2007). Development and seed number in indeterminate soybean as affected by timing and duration of exposure to long photoperiods after flowering. *Ann. Bot.* 99, 925–933. doi: 10.1093/aob/mcm033

Kantolic, A. G., Peralta, G. E., and Slafer, G. A. (2013). Seed number responses to extended photoperiod and shading during reproductive stages in indeterminate soybean. *Eur. J. Agron.* 51, 91–100. doi: 10.1016/j.eja.2013.07.006

Keating, B. A., Carberry, P. S., Hammer, G. L., Probert, M. E., Robertson, M. J., Holzworth, D., et al. (2003). "An overview of APSIM, a model designed for farming systems simulation," *Eur. J. Agro.* 18, 267–288. doi: 10.1016/S1161-0301(02)00108-9

Kurasch, A. K., Hahn, V., Leiser, W. L., Vollmann, J., Schori, A., Bétrix, C. A., et al. (2017). Identification of mega-environments in Europe and effect of allelic variation at maturity E loci on adaptation of European soybean. *Plant Cell Environ.* 40, 765–778. doi: 10.1111/pce.12896

Lu, S., Zhao, X., Hu, Y., Liu, S., Nan, H., Li, X., et al. (2017). Natural variation at the soybean J locus improves adaptation to the tropics and enhances yield. *Nat. Genet.* 49, 773–779. doi: 10.1038/ng.3819

Mahieu, P. J., and Brinkman, M. A. (1990). Double-cropping soybean after harvesting small grains as forage in the north central USA. *J. Prod. Agric.* 3, 385–389. doi: 10.2134/jpa1990.0385

Major, D. J., Johnson, D. R., Tanner, J. W., and Anderson, I. C. (1975). Effects of daylength and temperature on soybean development. *Crop Sci.* 15, 174. doi: 10.2135/cropsci1975.0011183X001500020009x

Maury, P., Andrianasolo, F. N., Alric, F., Berger, M., Beugniet, G., Chambert, C., et al. (2015). Le semis très précoce : une stratégie agronomique pour améliorer les performances du soja en France ? *OCL* 22, D503. doi: 10.1051/ocl/2015028

Mcblain, B. A., and Bernard, R. L. (1987). A new gene affecting the time of flowering and maturity in soybeans. *J. Hered.* 78, 160–162. doi: 10.1093/oxfordjournals.jhered.a110349

McDonald, M. B. B., Vertucci, C. W. W., and Roos, E. E. E. (1988). Soybean seed imbibition: water absorption by seed parts. *Crop Sci.* 28, 993–997. doi: 10.2135/cropsci1988.0011183X002800060026x

Meyer, D. W., and Badaruddin, M. (2001). Frost tolerance of ten seedling legume species at four growth stages. *Crop Sci.* 41, 1838–1842. doi: 10.2135/cropsci2001.1838

Nico, M., Miralles, D. J., and Kantolic, A. G. (2015). Post-flowering photoperiod and radiation interaction in soybean yield determination: direct and indirect photoperiodic effects. *F. Crop Res.* 176, 45–55. doi: 10.1016/j.fcr.2015.02.018

Olesen, J. E., Trnka, M., Kersebaum, K. C., Skjelvåg, A. O., Seguin, B., Peltonen-Sainio, P., et al. (2011). Impacts and adaptation of European crop production systems to climate change. *Eur. J. Agron.* 34, 96–112. doi: 10.1016/j.eja.2010.11.003

Parent, B., and Tardieu, F. (2012). Temperature responses of developmental processes have not been affected by breeding in different ecological areas for 17 crop species. *New Phytol.* 194, 760–774. doi: 10.1111/j.1469-8137.2012.04086.x

Parent, B., Turc, O., Gibon, Y., Stitt, M., and Tardieu, F. (2010). Modelling temperature-compensated physiological rates, based on the co-ordination of responses to temperature of developmental processes. *J. Exp. Bot.* 61, 2057–2069. doi: 10.1093/jxb/erq003

Piper, E. L., Boote, K. J., Jones, J. W., and Grimm, S. S. (1996). Comparison of two phenology models for predicting flowering and maturity date of soybean. *Crop Sci.* 36, 1606–1614. doi: 10.2135/cropsci1996.0011183X003600060033x

Poehlman, J. M. (1987). "Breeding Soybeans," in *Breeding Field Crops* (Dordrecht: Springer Netherlands), 421–450. doi: 10.1007/978-94-015-7271-2_17

Quinsac, A., Lecomte, V., Bourrel, C., Labalette, F., Ledoux, S., and Jouffret, P. (2014). Panorama et futur de la filière du soja français. *Oléagineux Corps. gras Lipides.* 17, 345–355. doi: 10.1051/ocl.2010.0345

Raper, C. D., and Thomas, J. F. (1978). Photoperiodic alteration of dry matter partitioning and seed yield in soybeans. *Crop Sci.* 18, 654–656. doi: 10.2135/cropsci1978.0011183X001800040032x

Robertson, G. W. (1968). A biometeorological time scale for a cereal crop involving day and night temperatures and photoperiod. *Int. J. Biometeorol.* 12, 191–223. doi: 10.1007/BF01553422

Salem, M. A., Kakani, V. G., Koti, S., and Reddy, K. R. (2007). Pollen-based screening of soybean genotypes for high temperatures. *Crop Sci.* 47, 219–231. doi: 10.2135/cropsci2006.07.0443

Salmerón, M., and Purcell, L. C. (2016). Simplifying the prediction of phenology with the DSSAT-CROPGRO-soybean model based on relative maturity group and determinacy. *Agric. Syst.* 148, 178–187. doi: 10.1016/j.agsy.2016.07.016

Salmeron, M., Gbur, E. E., Bourland, F. M., Buehring, N. W., Earnest, L., Fritschi, F. B., et al. (2014). Soybean maturity group choices for early and late plantings in the midsouth. *Agron. J.* 106, 1893–1901. doi: 10.2134/agronj14.0222

Scott, W. O., and Aldrich, S. R. (1970). *Modern soybean production* (Champaign, Illinois: S & A Publications.). Available at: https://www.cabdirect.org/cabdirect/abstract/19790781104 [Accessed March 19, 2018].

Seifert, C. A., and Lobell, D. B. (2015). Response of double cropping suitability to climate change in the United States. *Environ. Res. Lett.* 10 (2), 024002. doi: 10.1088/1748-9326/10/2/024002

Setiyono, T. D., Weiss, A., Specht, J., Bastidas, A. M., Cassman, K. G., and Dobermann, A. (2007). Understanding and modeling the effect of temperature and daylength on soybean phenology under high-yield conditions. *F. Crop Res.* 100, 257–271. doi: 10.1016/j.fcr.2006.07.011

Shapiro, B. I., Brorsen, B. W., and Doster, D. H. (1992). Adoption of double-cropping soybeans and wheat. *South. J. Agric. Econ.* 24, 33–40. doi: 10.1017/S0081305200018355

Sinclair, T. R., Neumaier, N., Farias, J. R. B., and Nepomuceno, A. L. (2005). Comparison of vegetative development in soybean cultivars for low-latitude environments. *F. Crop Res.* 92, 53–59. doi: 10.1016/j.fcr.2004.08.008

Terres Inovia. (2019). Outil IRRIsoja. Available at: https://www.terresinovia.fr/p/irrisoja

Tribouillois, H., Dürr, C., Demilly, D., Wagner, M.-H., and Justes, E. (2016). Determination of germination response to temperature and water potential for a wide range of cover crop species and related functional groups. *PloS One* 11 (8), e0161185. doi: 10.1371/journal.pone.0161185

USDA. (2019). *PSD_Oilseeds. PSD_Oilseeds*. Available at: https://apps.fas.usda.gov/psdonline/app/index.html#/app/downloads.

Wang, E., and Engel, T. (1998). Simulation of phenological development of wheat crops. *Agric. Syst.* 58, 1–24. doi: 10.1016/S0308-521X(98)00028-6

Wang, Z., Reddy, V. R., and Acock, M. C. (1998). Testing for early photoperiod insensitivity in soybean. *Agron. J.* 90, 389–392. doi: 10.2134/agronj1998.00021962009000030012x

Xia, Z., Zhai, H., Liu, B., Kong, F., Yuan, X., Wu, H., et al. (2012). Molecular identification of genes controlling flowering time, maturity, and photoperiod response in soybean. *Plant Syst. Evol.* 298, 1217–1227. doi: 10.1007/s00606-012-0628-2

Xu, M., Xu, Z., Liu, B., Kong, F., Tsubokura, Y., Watanabe, S., et al. (2013). Genetic variation in four maturity genes affects photoperiod insensitivity and PHYA-regulated post-flowering responses of soybean. *BMC Plant Biol.* 13 (1), 91. doi: 10.1186/1471-2229-13-91

Yin, X., Kropff, M. J., McLaren, G., and Visperas, R. M. (1995). A nonlinear model for crop development as a function of temperature. *Agric. For. Meteorol.* 77, 1–16. doi: 10.1016/0168-1923(95)02236-Q

Yuesheng, W., Jianbing, Q., Junyi, G., and Guangyuang, H. (2006). Classification and characteristic of maturity groups of Chinese landraces of soybean [Glycine max (L.) Merr.]. *Genet. Resour. Crop Evol.* 53, 803–809. doi: 10.1007/s10722-004-5731-y

Zhang, L. X., Kyei-Boahen, S., Zhang, J., Zhang, M. H., Freeland, T. B., Watson, C. E., et al. (2007). Modifications of optimum adaptation zones for soybean maturity groups in the USA. *Crop Manage.* 6, 1–10. doi: 10.1094/CM-2007-0927-01- RS.Abstract
Zhang, Q.-Y., Gao, Q.-L., Herbert, S. J., Li, Y.-S., and Hashemi, A. M. (2010). Influence of sowing date on phenological stages, seed growth and marketable yield of four vegetable soybean cultivars in North-eastern USA. *Afr. J. Agric. Res.* 5, 2556–2562.

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

Copyright © 2020 Schoving, Stöckle, Colombet, Champolivier, Debaeke and Maury. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Zhang, Q.-Y., Gao, Q.-L., Herbert, S. J., Li, Y.-S., and Hashemi, A. M. (2010). Influence of sowing date on phenological stages, seed growth and marketable yield of four vegetable soybean cultivars in North-eastern USA. *Afr. J. Agric. Res.* 5, 2556–2562.

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

Copyright © 2020 Schoving, Stöckle, Colombet, Champolivier, Debaeke and Maury. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.