Loading [MathJax]/jax/output/SVG/jax.js
Research article

Modeling the temporal dynamics of chlordecone in the profile of tropical polluted soils as affected by land use change

  • Received: 26 April 2021 Accepted: 29 June 2021 Published: 08 July 2021
  • The insecticide chlordecone (CLD) was applied from 1972 to 1993 to banana fields in the French Antilles, which resulted in the long-term pollution of soils and the contamination of crops and water resources. We coupled two biophysical models describing CLD and soil organic carbon (SOC) dynamics to determine the impact of a change of cropping system from banana to vegetable crops on the temporal pattern of CLD content in different soil types, and to assess how this might impact crop contamination and environmental pollution. The results indicated that a change of the cropping system when the CLD content in the topsoil (0–0.3 m) drops below the threshold (1 mg CLD kg-1) established by local authorities to allow the cultivation of vegetable crops (e.g., cucumber, melon, watermelon, pumpkin), might cause crop contamination due to the presence of relatively high CLD levels in the subsurface layer (i.e., 0.3–0.6 m) of nitisols and ferralsols. The impact of changing the cropping system on the risk of environmental pollution depends on the time of that change, and it is much greater for vegetable crop systems established in the early 1990s following a financial crisis affecting the banana sector. This is linked to the progressive decline of SOC stocks caused by vegetable crop systems, which reduces CLD retention in the soil and increases CLD leaching. Overall, this study highlights the urgent need to include the monitoring of CLD in the subsurface layer and the dynamics of SOC stocks in current soil testing campaigns carried out on polluted soils in French Antilles.

    Citation: Jorge Sierra, Antoine Richard. Modeling the temporal dynamics of chlordecone in the profile of tropical polluted soils as affected by land use change[J]. AIMS Environmental Science, 2021, 8(4): 304-320. doi: 10.3934/environsci.2021020

    Related Papers:

    [1] Thierry Woignier, Florence Clostre, Philippe Cattan, Magalie Lesueur-Jannoyer . Pollution of soils and ecosystems by a permanent toxic organochlorine pesticide: chlordecone—numerical simulation of allophane nanoclay microstructure and calculation of its transport properties. AIMS Environmental Science, 2015, 2(3): 494-510. doi: 10.3934/environsci.2015.3.494
    [2] M.A. Rahim, M.G. Mostafa . Impact of sugar mills effluent on environment around mills area. AIMS Environmental Science, 2021, 8(1): 86-99. doi: 10.3934/environsci.2021006
    [3] Dongyan Mu, John Hawks, Andrew Diaz . Impacts on vegetable yields, nutrient contents and soil fertility in a community garden with different compost amendments. AIMS Environmental Science, 2020, 7(4): 350-365. doi: 10.3934/environsci.2020023
    [4] Carlos Garcia, Teresa Hernandez, Maria D Coll, Sara Ondoño . Organic amendments for soil restoration in arid and semiarid areas: a review. AIMS Environmental Science, 2017, 4(5): 640-676. doi: 10.3934/environsci.2017.5.640
    [5] Jerry R. Miller, John P. Gannon, Kyle Corcoran . Concentrations, mobility, and potential ecological risks of selected metals within compost amended, reclaimed coal mine soils, tropical South Sumatra, Indonesia. AIMS Environmental Science, 2019, 6(4): 298-325. doi: 10.3934/environsci.2019.4.298
    [6] Parthiban Rajasekaran, Harikishan Kannan, Smruti Das, Mikaeel Young, Swadeshmukul Santra . Comparative analysis of copper and zinc based agrichemical biocide products: materials characteristics, phytotoxicity and in vitro antimicrobial efficacy. AIMS Environmental Science, 2016, 3(3): 439-455. doi: 10.3934/environsci.2016.3.439
    [7] Ioannis Panagopoulos, Athanassios Karayannis, Georgios Gouvalias, Nikolaos Karayannis, Pavlos Kassomenos . Chromium and nickel in the soils of industrial areas at Asopos river basin. AIMS Environmental Science, 2016, 3(3): 420-438. doi: 10.3934/environsci.2016.3.420
    [8] Lawrence H. Tanner, Megan M. Vandewarker . Significance of vegetation cover differences on albedo and soil carbon on a basaltic sandplain in southern Iceland. AIMS Environmental Science, 2019, 6(6): 435-444. doi: 10.3934/environsci.2019.6.435
    [9] Mahidin, Asri Gani, Saiful, Muhammad Irham, Wulan Windari, Erdiwansyah . An overview of the potential risks, sources, and analytical methods for microplastics in soil. AIMS Environmental Science, 2022, 9(2): 185-216. doi: 10.3934/environsci.2022013
    [10] Maja Radziemska, Agnieszka Bęś, Zygmunt M. Gusiatin, Jerzy Jeznach, Zbigniew Mazur, Martin Brtnický . Novel combined amendments for sustainable remediation of the Pb-contaminated soil. AIMS Environmental Science, 2020, 7(1): 1-12. doi: 10.3934/environsci.2020001
  • The insecticide chlordecone (CLD) was applied from 1972 to 1993 to banana fields in the French Antilles, which resulted in the long-term pollution of soils and the contamination of crops and water resources. We coupled two biophysical models describing CLD and soil organic carbon (SOC) dynamics to determine the impact of a change of cropping system from banana to vegetable crops on the temporal pattern of CLD content in different soil types, and to assess how this might impact crop contamination and environmental pollution. The results indicated that a change of the cropping system when the CLD content in the topsoil (0–0.3 m) drops below the threshold (1 mg CLD kg-1) established by local authorities to allow the cultivation of vegetable crops (e.g., cucumber, melon, watermelon, pumpkin), might cause crop contamination due to the presence of relatively high CLD levels in the subsurface layer (i.e., 0.3–0.6 m) of nitisols and ferralsols. The impact of changing the cropping system on the risk of environmental pollution depends on the time of that change, and it is much greater for vegetable crop systems established in the early 1990s following a financial crisis affecting the banana sector. This is linked to the progressive decline of SOC stocks caused by vegetable crop systems, which reduces CLD retention in the soil and increases CLD leaching. Overall, this study highlights the urgent need to include the monitoring of CLD in the subsurface layer and the dynamics of SOC stocks in current soil testing campaigns carried out on polluted soils in French Antilles.



    The organochlorine insecticide chlordecone (CLD) was used in French Antilles (Guadeloupe and Martinique) from 1972 to 1993 to control banana weevil (Cosmopolites sordidus), which resulted in diffuse soil pollution that induced environmental damage, crop contamination and severe problems for human health [1]. Chlordecone is a very stable and recalcitrant molecule, with a strong affinity for soil organic matter (SOM) and very low volatility [2]. Cabidoche et al. (2009) were the first to report that the persistence of CLD differs between soil types, ranging from decades for nitisols, centuries for ferralsols to about a millennium for andosols; this is linked to their different mineralogy and SOM contents that markedly affect their sorption properties [3]. Some 30 years after these treatments ended, one third of the agricultural land in French Antilles is still polluted by CLD, and both animal breeding and the cultivation of many crops for human consumption (mainly tuber and root crops) have been banned on these soils [4]. In this context, the Guadeloupe Department of Food, Agriculture and Forestry established three thresholds for soil CLD content in the topsoil of polluted land for the cultivation of vegetables (<1 mg CLD kg-1) and tuber crops (<0.1 mg CLD kg-1), and for animal breeding (<0.003 mg CLD kg-1) [5] that would enable the level of contamination of harvested organs and meat not to exceed the maximum residue limit (MRL) defined by European Union health regulations [6]. At present, only banana, orchards and some vegetable crops such as tomato and pepper can be cultivated in soils containing > 1 mg CLD kg-1, because CLD is no longer detected in the harvested fruits [5].

    Despite the high stability of CLD, some authors have reported evidence of limited CLD degradation in laboratory experiments under aerobic [7] and anaerobic conditions [8], and its relatively high degradation in the field under anaerobic conditions in a nitisol amended with zero valent iron [1]. Although the CLD content in topsoil had fallen by up to 30% three months after the application of the treatment in the latter study, maintaining anaerobic conditions for several months under real-life soil management by farmers is extremely difficult because of the high water drainage that characterizes the volcanic soils of French Antilles [9]. In a recent study, Cattan et al. (2019) used a model of CLD and metabolite sorption, leaching and degradation to estimate the annual rate of CLD degradation in Martinique soils from data on the CLD and metabolite contents in surface water and groundwater [10]. They found that the estimated annual rate of CLD degradation ranged from 0.02% to 0.15% of the soil CLD stock, which was much less than the rate of CLD leaching. Based on the information available concerning CLD dynamics it can therefore be assumed that leaching is the principal process affecting CLD stocks in the volcanic soils of French Antilles [3,10].

    Most studies of soil CLD have only accounted for its dynamics or status in the topsoil layer [3,11,12]. Clostre et al. (2014) analyzed the vertical heterogeneity of CLD levels in six agricultural fields in French Antilles, and observed that although the topsoil layer (e.g., 0–0.3 m) generally contained higher levels than the subsoil layer (e.g., 0.3–0.6 m), the differences between layers were very small for some soils (e.g., 3%), and in one case more CLD was found in the subsoil layer [13]. These authors also noted that within each soil type analyzed (andosols and ferralsols) the CLD content in both layers displayed considerable variability (e.g., the coefficients of variation ranged from 30% to 45%). Such variability may be associated with differences in the total rate of CLD spreading, management of the cropping system and the tillage practices applied by farmers [11]. Cabidoche et al. (2009) also found that the capacity of CLD sorption varies strongly within each soil type, which might differently affect the ratio of CLD content between soil layers since the time of the last CLD application [3]. Taken together, these results suggest that CLD levels in the subsurface layer cannot be predicted directly from observed data on the topsoil layer. Indeed, knowledge of CLD dynamics in the subsurface layer might be crucial when the CLD content in the topsoil approaches the established threshold that permits the cultivation of vegetable crops. This would contribute to reducing the risk of crop contamination for vegetables with a relatively deep root system such as cucumber, pumpkin, melon and watermelon (e.g., 0.6 m depth) [14], which are currently banned on soils with CLD levels > 1 mg kg-1 in the 0–0.3 m layer [5].

    Another factor involved in the cultivation of vegetable crops is the impact of these cropping systems on the level of soil organic carbon (SOC) stocks, which may affect the sorption properties of the soil. Sierra et al. (2015, 2017) found that although SOC under the current banana and orchards systems are near or at the steady-state, most vegetable cropping systems in French Antilles induce a reduction in SOC stocks at rates that range from 0.5% to 1% per year [15,16]. Sierra et al. (2015) proposed that such differences between systems are mainly due to the frequency and intensity of soil tillage, which is much higher for vegetable crops [15]. It could be hypothesized that the cultivation of vegetable crops after the topsoil reaches the critical threshold may enhance CLD desorption and leaching, and thus increase the risk of environmental pollution. The same could have occurred since 1993 when a financial crisis affecting the banana sector generated a 32% reduction in the area cultivated with this crop that was partially replaced by vegetable and food crops [17]. It is interesting to note that this crisis was not linked to soil pollution by CLD, which was revealed for the first time in the early 2000s [18]. In a recent study, Sabatier et al. (2021) observed that changes in farming practices on polluted soils in French Antilles, including use of the herbicide glyphosate from the late 1990s, induced a rise in CLD levels in marine sediment through the cumulated erosion of contaminated bare soils [19]. To our knowledge, no study has been performed to assess the overall effect of changes in land use and farming practices on CLD dynamics in soil in terms of SOC and CLD interactions. This is necessary to better manage the risk of crop contamination and to prevent any changes to the management of the cropping systems that might cause a further release of soil CLD towards other environmental compartments.

    The aim of this work was therefore to determine the impact of changes in land use on the temporal pattern of CLD content in the topsoil and subsurface layers of different soil types in French Antilles (nitisols, ferralsols and andosols). To achieve this, we applied a modeling approach coupling the WISORCH model that describes CLD dissipation [3] and the MorGwanik model that describes SOC dynamics in different soil types and under different cropping systems [15]. Simulations were carried out over the 1972–2045 period, thus including the time of CLD spreading (1972–1993) and the subsequent phase characterized by CLD losses. We analyzed two cropping systems involving banana in monoculture and banana followed by a monoculture of vegetable crops.

    The soils analyzed during this study are located in Basse-Terre, the main island of the Guadeloupe archipelago (French Antilles) located in the eastern Caribbean (16°05′ N, 61°40′ W). Basse-Terre Island covers 848 km2 and is dominated by a mountain chain oriented northwest to southeast where the crest stands at 1467 m.a.s.l. (La Soufrière volcano). Although the land west of the crest slopes steeply toward the Caribbean Sea (20–30% slopes), the eastern side of the chain slopes more gently toward the Atlantic Ocean (10% slopes). The soils affected by CLD pollution are located in the southern part of the island where banana cultivated as a monoculture is the main cropping system. The soils are nitisols and ferralsols in the lowlands and andosols in the uplands (FAO classification) [15]. The mean air temperature and mean annual rainfall are respectively 25.0 ℃ and 2200 mm yr-1 in the lowlands, and 23.9 ℃ and 3800 mm yr-1 in the uplands. Nitisols are rich in halloysite clay (around 50%) and have developed on young ash deposits; the soil pH ranges from 5.0 to 6.5, and the cation exchange capacity from 15 cmol kg-1 to 35 cmol kg-1. Ferralsols are rich in kaolinite (around 35%) and aluminum and iron hydrous oxides (around 40%) and have developed on old ash deposits. The soil pH ranges from 4.5 to 5.5, while their cation exchange capacity ranges from 10 cmol kg-1 to 20 cmol kg-1. Andosols are characterized by high amorphous clay content (around 75%) and have developed on young ash deposits. They are very porous and have high aluminum content; the soil pH ranges from 5.0 to 6.5 and the cation exchange capacity from 30 cmol kg-1 to 50 cmol kg-1. SOC stocks are relatively high in all the soils and vary markedly between soil types (Table 1).

    Table 1.  Characteristics of the 0–0.3 m layer of the soils analyzed and parameter values for the WISORCH model. Data were taken from the study of Cabidoche et al. (2009) [3].
    Soils Chlordecone spreading Rainfall**
    (m yr-1)
    Bulk density
    (Mg m-3)
    SOC stock#
    (Mg C ha-1)
    KOC##
    (m3 kg-1)
    Period Total rate*
    (kg ha-1)
    Andosol 1 1972–78 & 1982–93 57 6.3 0.40 163 18.03
    Andosol 2 1982–93 36 3.3 0.79 83 17.50
    Ferralsol 1 1972–93 66 3.6 0.92 58 8.28
    Ferralsol 2 1982–87 18 2.8 0.91 60 9.14
    Nitisol 1 1972–79 & 1981–93 63 2.5 0.90 51 3.53
    Nitisol 2 1972–79 & 1982–87 42 2.5 0.95 40 3.78
    *The annual rate was 3 kg ha-1 yr-1 for all analyzed soils.
    **Mean value for the 1970–2015 period.
    #Initial value for the topsoil layer in 1972.
    ## Partitioning coefficient of chlordecone between SOC and water.

     | Show Table
    DownLoad: CSV

    For this study, we selected two soils of each type included in the database generated during the study by Cabidoche et al. (2009) [3]. The soils selected within each type corresponded to plots that had received extreme rates of CLD applications between 1972 and 1993. Some characteristics of the selected soils are presented in Table 1. Note that the level of rainfall varies between soil types and within each type as a function of the microclimate surrounding the plots.

    We used the WISORCH model proposed by Cabidoche et al. (2009) to assess the dynamics of CLD in the topsoil layer [3]. Taking account of the fact that CLD volatilization is nil [2], and losses of CLD by crop uptake [6] and runoff [10] are negligible when compared to CLD stocks in the soil, this model only considers desorption and leaching as the main soil processes affecting CLD stocks. We extended the original model to account for CLD dynamics in the subsurface layer. We also included a term dealing with CLD degradation by soil microorganisms to assess the potential impact of such process on changes in CLD dynamics, using the approach proposed by Cattan et al (2019) [10].

    The basic equations of the model are as follows: for the topsoil layer (i.e., 0–0.3 m in this study)

    CLDn(T)=[(CLDn-1(T)+Sn)×(1Deg(T))]×f(Drain,KOC(T),SOCn(T)) (1)

    for the subsurface layer (i.e., 0.3–0.6 m in this study)

    CLDn(S)=[(CLDn-1(S)+CLDleach, n(T))×(1Deg(S))]×f(Drain,KOC(S),SOCn(S)) (2)

    with

    f(Drain,KOC(L),SOC(L))={exp[(10×Drain)÷(KOC(L)×SOCn(L))]} (3)

    where the suffixes (T) and (S) indicate respectively the topsoil and subsurface layers, the suffix (L) in Eq 3 indicates (T) or (S), CLDn(L) and CLDn-1(L) (kg CLD ha-1) are the CLD stocks in both layers in years n and n-1, respectively, Sn (kg CLD ha-1) is the CLD spread on the soil surface in year n, Deg(L) (yr-1) is the annual rate of CLD degradation for each layer, Drain (m) is the water drainage in year n, KOC(L) (m3 kg-1) is the partitioning coefficient of CLD between SOC and water for each layer, SOCn(L) (Mg SOC ha-1) is the stock of SOC for each layer in year n, and CLDleach, n(T) (kg CLD ha-1) is the CLD leached from the topsoil and entering the subsurface layer in year n, which is calculated using Eq 1. The factor 10 in Eq 3 arises from the units chosen and allowed us to express CLD stocks in kg ha-1. The exponential term in Eq 3 represents the CLD desorption process which is affected by KOC(L) and SOC(L), and the leaching of the CLD released through water drainage. Note that KOC, as expressed in Eq 3, also depends on soil mineralogy [3,20]. This is further discussed in section 3.1. Eq 2 implies that all the CLD leached from the topsoil layer was a CLD input for the subsurface layer. As mentioned above, this assumption is based on the fact that other CLD losses are negligible in the soils being analyzed. The model considers two types of water drainage: i- forced drainage linked to the stemflow induced by the banana plant, which affects CLD leaching during the five years after the last CLD application; it is calculated as 1.2 × annual rainfall, and ii- after that period, homogeneous drainage estimated according to the water balance [3]. The forced drainage was included in the model to take into account that rainfall redistribution by plant canopy induces strong water drainage at the foot of the banana stem, downstream from the stemflow, which crosses the foot spreading of CLD during the application period (i.e., 1972–1993) [3].

    The corresponding CLD gravimetric content (CLDcont, in mg CLD kg-1) of each soil layer was calculated as:

    CLDcontn(L)=CLDn(L)÷(BD(L)×Z(L)×10) (4)

    where BD(L) (Mg m-3) and Z(L) (m) are respectively the bulk density and thickness of each soil layer. The factor 10 arises from the units chosen and allowed us to express CLDcontn(L) in mg CLD kg-1.

    The original WISORCH model accounts for the mechanical dilution of CLD between layers when the depth of soil tillage (e.g., every 5 yr for the banana cropping systems) exceeds the bottom of the topsoil layer [3]. However, it is not clear from the formalisms of the model how the subsurface layer contributes to CLD dilution in the topsoil, and whether SOC dilution between layers is also accounted for. The latter could be important because changes in SOC may affect the sorption properties of the soil (Eq 3). In the present study we did not include any effect of soil tillage on CLD dilution by assuming that tillage only concerned the topsoil layer (i.e., tillage depth equal to 0.3 m). This is further discussed in section 3.1.

    We used the MorGwanik model to assess the impact of changing the cropping system from banana to vegetable crops on the SOC stocks of the topsoil layer, and to set the SOCn(T) value in Eq 1. This model was calibrated and tested using data of the time course of SOC stocks obtained from 253 plots covering most soil types, climates and cropping systems present in the Caribbean, and has been applied successfully to assessing SOC dynamics in many Guadeloupe cropping systems [15,16]. MorGwanik simulates SOC dynamics as a function of annual C inputs and outputs at the plot scale. The basic equation of the model is:

    SOCn(T)=SOCn-1(T)+(Cres, n×hres)+(Came, n×hame)[SOCn-1(T)×(ksoil×ktill)] (5)

    where Cres, n (Mg C ha-1) is the C input from crop residues (aboveground and roots) in year n, Came, n (Mg C ha-1) is the C input from organic amendments in year n, hres and hame (unitless) are respectively the humification coefficients of crop residues and organic amendments, ksoil (yr-1) is the SOC mineralization rate constant for each soil type, and ktill (unitless) is the coefficient accounting for the effect of soil tillage on SOC mineralization. Although ksoil is specific to each soil type, ktill is specific to each cropping system (e.g., higher for annual crops than perennial crops). In this study we did not consider the application of organic amendments and the only C input corresponded to crop residues.

    As the SOC stock values in the subsurface layer were lacking for the soils analyzed in this study, we used experimental data reported by other authors to set the initial value of SOCn(S) in Eq 2. In this way, SOC stocks in the subsurface layer represented 56% of the SOC stock in the topsoil layer of andosols [20], 81% of ferralsols [21], and 85% of nitisols [22]. As MorGwanik only accounts for the topsoil layer, we assumed that annual changes in SOC stocks in the subsurface layer corresponded to a fraction of those occurring in the topsoil layer:

    SOCn(S)=SOCn-1(S)+0.1×(SOCn(T)SOCn-1(T)) (6)

    where 0.1 is the fraction of the change in the topsoil that was applied to the subsurface layer (i.e., 10%). That value is an average of the changes observed by Poeplau et al. (2011) and Balesdent et al. (2018) in studies dealing with the temporal dynamics of SOC as a function of soil depth [23,24]. Note that SOCn(S) may increase, decrease or remain constant as a function of changes in the topsoil.

    A first series of simulations was performed from 1972 to 2045 to assess the effect of soil type on CLD dynamics in the topsoil and subsurface layers. Simulations were made for soils cultivated with a banana monoculture, which corresponds to the more common cropping system observed on polluted soils in French Antilles [10]. This was performed with and without CLD degradation in order to test the impact of this process on changes in CLD content. Data of KOC for the subsurface layer of the analyzed soils was lacking, but the information available indicated that KOC would remain quite stable up to 0.6 m for soils in Guadeloupe [3,20]. We therefore used the same KOC value for both layers in each soil. The same could reasonably be assumed for BD(L) in Eq 4 [21,25]. For the simulations carried out to test the impact of CLD degradation, Deg(T) in Eq 1 was set at 0.0015 yr-1 (i.e., 0.15% of the CLD stock), which is the median value reported by Cattan et al. (2019) for natural CLD degradation in the soils of Martinique (French Antilles) [10]. These authors estimated CLD degradation in soil from data on the pollutant contents in surface water and groundwater at the watershed scale, which involved several soil types. Because data of CLD degradation for each soil type is lacking, in the present study we considered the same Deg value for all the soils and layers analyzed. This simplification is further discussed in section 3.2. The parameter values for the WISORCH model are presented in Table 1. The values of soil and banana parameters for the MorGwanik model are presented in Tables S1, S2 and S3 in Appendix A. Both models were initialized with the SOC stocks presented in Table 1.

    A second series of simulations was carried out to assess the impact on CLD and SOC dynamics of changing the cropping system from a banana monoculture to vegetable crops. This was performed in two steps. First, we tested the impact of such a change as from the next year when the CLD content was < 1mg CLD kg-1 in the topsoil layer, which is the maximum CLD level in topsoil established by the Guadeloupe Department of Food, Agriculture and Forestry to permit the cultivation of vegetable crops [5]. In the second step we set the change of cropping system at 1993 for all soils. As mentioned above, the crisis experienced by the banana sector that year induced an increase in the area cultivated with vegetables at the expense banana crops. All these simulations were carried out considering a banana monoculture from 1972 to the year of the change, and then annual vegetable crops up to 2045. The values of soil and vegetable crop parameters for the MorGwanik model are presented in Tables S1 and S2, respectively. Changes in CLD stocks were estimated as for the first simulation series using the parameter values presented in Table 1.

    In this section we present the results obtained excluding CLD degradation.The topsoil layer presented a first phase corresponding to the period of CLD spreading, which was characterized by a rapid increase in CLD levels, followed by a second phase when the CLD content diminished at different rates as a function of the soil type (Figure 1). The changes of CLD content in the subsurface layer depended on the soil type: a continuous increase for andosols (Figure 1a, b) and an increase followed by a decrease for ferralsols and nitisols (Figure 1c–f). The slight temporal fall in CLD levels in the topsoil of andosol 1 (Figure 1a), nitisol 1 (Figure 1e) and nitisol 2 (Figure 1f) during the late 1970s reflected the non-application of CLD during two or three years within this period (Table 1). At the time of the maximum CLD content (e.g., 1987 for ferralsol 2 and nitisol 2, and 1993 for the other soils), the correlation between soil CLD stocks up to 0.6 m and the total CLD rate was small and not significant (i.e., R2 = 0.29, P < 0.27). Indeed, at that time soil CLD stocks represented 93% of the CLD applied in andosols, 77% in ferralsols and only 41% in nitisols, which indicates that a large quantity of CLD was lost by leaching below 0.6 m during the period of CLD spreading, and that CLD losses by leaching were markedly dependent on soil type.

    Figure 1.  Temporal pattern of chlordecone (CLD) content in the topsoil and subsurface layers of the six soils analyzed under banana monoculture. Horizontal lines at 1 mg kg−1 and 0.1 mg kg−1 indicate the maximum thresholds established by local authorities to allow the cultivation of vegetable and tuber crops, respectively.

    During the 1973–1993 period, the rate of CLD leaching varied strongly between soil types: nitisols (average 1.6 kg ha-1 yr-1), ferralsols (0.7 kg ha-1 yr-1), and andosols (0.2 kg ha-1 yr-1) (Figure 2). Differences in CLD content and leaching between soil types can be explained by differences in KOC and SOC stocks, which were the highest in andosols (Table 1) inducing greater CLD sorption and longer residence time in the solid phase of the soil [3]. Cabidoche et al. (2009) had observed that KOC may differ markedly between soils of different types but with similar SOC stocks. This is due to the influence of the different mineralogy of volcanic soils on chlordecone sorption [3]. This was subsequently confirmed by Fernández-Bayo et al. (2013) in a laboratory study [20]. This is particularly true in andosols where the allophane microstructure favors CLD and SOC sequestration [26], thus reducing CLD desorption [27] and SOC mineralization [16]. It is interesting to note that rainfall levels were not the principal factor affecting CLD leaching during the first phase. For example, although CLD leaching was 3.6 times higher for nitisol 1 than for andosol 1 (Figure 2a, e), rainfall was 2.5 times lower for the former (Table 1). Indeed, the high rates of CLD leaching during the 1972–1993 period contributed to the high pollution of rivers, groundwater and crops, as revealed by surveys carried out by the French Department of Health of Martinique in the early 2000s [18].

    Figure 2.  Temporal pattern of the rate of chlordecone leaching below 0.6 m for the six soils analyzed under banana monoculture.

    Figure 1 shows that the CLD content in topsoil declined rapidly over the five years after the last application (i.e., until the late 1990s), which was due to both a lack of CLD input and the effects of forced drainage induced by stemflow. During that period, the rate of CLD leaching was relatively high, after which it fell sharply (Figure 2). Because the negative balance between CLD inputs and outputs was less for the subsurface layer of ferralsols and nitisols, the CLD content in these soils was gradually higher in that layer than in the topsoil (Figure 1c–f). For andosols, the balance for the subsurface layer was positive but CLD contents were always higher in the topsoil (Figure 1a, b). Once again, this was linked to the greater CLD retention in the topsoil of andosols that limited downwards CLD transfer. This implies that the proposal made by Cabidoche et al. (2009) [3] concerning long-term CLD pollution of the topsoil of andosols is also applicable to their subsurface layers. The levels of CLD contents and stocks estimated during the present study for the different soil types were close to the experimental data obtained by several authors between 2001 and 2010 [3,11,28] and the data in a GIS database including more recent CLD measurements [29]. In particular, the GIS database shows a regular decrease in CLD stocks in soil from the early 2000s, which is more pronounced for nitisols and ferralsols than for andosols. Moreover, CLD content in some nitisols in the late 2010s is only slightly higher than the threshold of 1 mg CLD kg-1 established by local authorities.

    The rates of CLD leaching over the 2020–2040 period differed markedly from those mentioned above for 1973–1993 and decreased in the order: ferralsols (average 0.24 kg ha-1 yr-1) > andosols (0.16 kg ha-1 yr-1) > nitisols (0.08 kg ha-1 yr-1). In other words, although nitisols were the major source of environmental CLD pollution during the period of spreading, at present the main sources would be ferralsols followed by andosols. Indeed, this is linked to CLD stocks remaining in soils after the period of the highest CLD losses. In this sense, Table 2 shows that in 2045, at the end of the period simulated, CLD stocks will account for < 1% of the total CLD applied in nitisols, 22% in ferralsols and 70% in andosols. At some sites affected by polluted soils, the shallow water table fluctuates between 0.6 and 3.4 m depending on the soil type and rainfall levels, and appears to be the principal contributor to stream contamination [30]. For this reason our calculations of CLD stocks focused only on the first 0.6 m of soils.

    Table 2.  Total leaching and soil stock (0–0.6 m) of chlordecone at the end of the simulated period in 2045 for the six soils analyzed under banana monoculture.
    Soils Total leaching Soil stock
    kg ha-1 % rate* kg ha-1 % rate*
    Andosol 1 16 28 41 72
    Andosol 2 12 33 24 67
    Ferralsol 1 57 86 9 14
    Ferralsol 2 12 69 6 31
    Nitisol 1 63 99 0.4 1
    Nitisol 2 42 99 0.1 1
    * % of the total rate of chlordecone applied in the 1972–1993 period (see Table 1).

     | Show Table
    DownLoad: CSV

    The inclusion of an effect of CLD degradation in the model had little effect on the general pattern of soil CLD dynamics described above. Throughout the simulated period, CLD degradation averaged 3 kg CLD ha-1 for andosols, 2 kg CLD ha-1 for ferralsols and 1 kg CLD ha-1 for nitisols. These values represented 7% of total applied amount of CLD for andosols, 4% for ferralsols and 1% for nitisols (Figure 3). It is clear that in both absolute and relative terms, CLD degradation was dependent on the residence time of CLD in soils, which is directly associated to the sorption capacity of each soil type (e.g., andosols > ferralsols > nitisols). Total CLD leaching at the end of the simulated period only decreased slightly when compared to the results obtained without CLD degradation: 4% for andosols, 3% for ferralsols and 1% for nitisols (data not shown). CLD degradation therefore appears to be sufficient to induce the presence of transformation products in surface water and groundwater [10], but it is insufficient to markedly modify soil CLD stocks. As a consequence, when estimated considering CLD degradation, the temporal patterns of CLD content and leaching were similar to those presented in Figures 1 and 2 for the simulations excluding degradation (data not shown). For andosols in French Antilles, Woignier et al. (2019) reported that allophane microporosity can strongly limit the ability of microorganisms to decompose sequestrated CLD because of mechanical constraints and poor transport properties [27]. It is thus possible that the CLD degradation in andosols calculated during our study overestimated the true degradation in these soils. In the remainder of this paper we will focus on the results obtained excluding CLD degradation.

    Figure 3.  Relative impact of chlordecone degradation on chlordecone losses (0–0.6 m) at the end of the simulated period in 2045 for the six soils analyzed under banana monoculture. Results correspond to the simulations performed considering the effect of chlordecone degradation.

    The CLD content in topsoil was equal to the critical threshold of 1 mg CLD kg-1 in 1993 for nitisol 2, in 2005 for nitisol 1 and in 2036 for ferralsol 2 (Figure 1). At these dates, CLD content in the subsurface layer were 1.8 mg CLD kg-1 for nitisol 2, 2.1 mg CLD kg-1 for nitisol 1, and 1.5 mg CLD kg-1 for ferralsol 2. The earlier date observed for nitisol 2 reflects its low rate of CLD spreading and the fact that it was last applied in 1987 (Table 1). The CLD content in the topsoil of both andosols and ferralsol 1 did not reach that threshold within the simulated period (Figure 1). Only nitisols reached the second threshold at 0.1 mg CLD kg-1 in topsoil (i.e., in 2016 for nitisol 2 and in 2030 for nitisol 1; Figure 1). Some studies carried out in Guadeloupe have shown that the roots of annual crops such as tropical maize can reach a depth of 0.5 m [21]. In addition, it is well known that vegetables such as cucumber, watermelon and pumpkin have a deep rooting system that can reach 0.6 m [14]. These crops are currently banned in soils with > 1 mg CLD kg-1 in the 0–0.3 m layer [5]. On the other hand, studies performed in French Antilles to assess crop CLD contamination in polluted soils were carried out in the field or under greenhouse conditions using homogenized soils [6,12]. During these studies, only the topsoil was analyzed for CLD content so that any potential contribution of the subsurface layer to crop contamination remained unknown. Our results suggest that the subsurface layer could contribute to attaining crop CLD contents that exceed the MRL when vegetables with deep root systems are cultivated in nitisols and ferralsols just after the topsoil reaches the threshold of 1 mg CLD kg-1. Indeed, further research is necessary to determine the actual CLD content in the subsurface layer of these soils and establish the effects of this layer on crop contamination in the field.

    Clostre et al (2014) observed that deeper tillage (i.e., > 0.4 m in their study) causes homogenization and dilution of the CLD content in the soil profile, which reduces the differences between the topsoil and subsurface layers [13]. Soil tillage in the more intensive banana systems in French Antilles includes subsoiling at ≥0.4 m [31]. Unlike deep plowing and deep mixing, subsoiling aims to reduce soil compaction without turning or mixing soil horizons [32]. Although some parts of the bottom of the topsoil and the top of the subsurface layer may be partially mixed, the precise mixed fraction in each layer is hard to estimate. Therefore, although it is clear that deeper tillage may induce mechanical dilution of CLD and high CLD content in the subsurface layer, as observed by Clostre et al. (2014) [13], such an effect cannot be estimated accurately from the data currently available on the effects of subsoiling. This information could be obtained during future research by estimating the fraction of the mixed layers from data on the CLD and SOC contents before and after subsoiling in soils that differ in terms of their mineralogy and structure.

    As described above, during the first step of this analysis we tested the impact of changes of the cropping systems as from the year after the CLD content was < 1mg CLD kg-1 in the topsoil layer; this concerned three soils: nitisol 1 from 2006, nitisol 2 from 1994, and ferralsol 2 from 2037. As reported in previous studies performed in French Antilles [15,16], a change in the cropping system may strongly affect SOC stocks, which could modify the sorption properties of the soil and then CLD leaching. The results showed that SOC stocks remained constant during the period under banana monoculture (Figure 4a), thus reflecting that this system is at the steady-state in terms of SOC stocks in French Antilles [15]. This is associated to the fact that C outputs (mineralization) and C inputs (crop residues) converge in the long-term under the same system management (e.g., more than five decades for banana systems in French Antilles) [3,10,15]. Following the change of cropping system, SOC stocks in the topsoil under vegetable crops fell by 1% yr-1 on average. This was in line with the data observed in Guadeloupe and reported by Sierra et al. (2017) for monocultures of vegetables when organic amendments were not applied by farmers [16]. Because C inputs and outputs were the same in both nitisols (Table S1 and S2), their SOC stocks converged towards the same new equilibrium at around 36 Mg C ha-1 (Figure 4a). High SOC losses under vegetables crops are mainly caused by the high intensity of soil tillage (e.g., up to six tillage operations per year) that enhances SOC mineralization [16]. The ktill parameter in Eq 5 reflects this effect, which is about double with vegetables than with banana systems where tillage is only applied every 5 years [15] (Table S2).

    Figure 4.  Impact of changing the cropping system from banana monoculture to vegetable crops on (a) and (b) SOC stocks in the 0–0.3 m layer, and (c) the increase in total chlordecone leaching below 0.6 m at the end of the simulated period in 2045. In (a), the change of cropping system was set from the year after the chlordecone content was < 1mg kg-1 in the topsoil layer, which only concerned ferralsol 2 (from 2037), nitisol 1 (from 2006) and nitisol 2 (from 1994). In (b) and (c), the change of cropping system was fixed in 1993 for the six soils analyzed. In (c), the increase in chlordecone leaching is related to the values obtained for the banana monoculture system (see Table 2). In (a) and (b), vertical lines indicate the time of the change of cropping system. In (c), the increase in chlordecone leaching is related to the values obtained for the banana monoculture system (see Table 2).

    Despite the relatively large fall in SOC stocks under vegetable crops, affecting the sorption properties of the soils, its impact on CLD content and leaching was negligible (data not shown). For example, total CLD leaching at the end of the simulated period rose less than 1% in the three soils when compared to the values obtained with banana monoculture and presented in Table 2. Thus the data presented in Figures 1 and 2 for these soils were not substantially modified by the change of cropping system when CLD content in the topsoil reached the threshold of 1 mg kg-1. Indeed, at that time, CLD stocks were already too small to further affect CLD dynamics and leaching. From these results, it therefore appears that the effect of changing the cropping system when only the topsoil reached the threshold might be more important in terms of crop contamination than environmental pollution.

    In the second step of this analysis we fixed the change of cropping system in 1993 following the financial crisis experienced by the banana sector. Figure 4b shows that the fall in SOC stocks under vegetable crops was relatively marked in all soils. At the end of the simulated period, SOC losses averaged 32% of the initial SOC stocks. As an index of the effect of SOC changes on CLD dynamics Figure 4c presents the increase in total CLD leaching compared to those estimated under banana monoculture (i.e., Table 2). On average, total leaching increased by 6% in andosols, 5% in ferralsols and 0.1% in nitisols. This increase was negligible in nitisols because their sorption capacity is too small to be further affected by a gradual decline in SOC stocks under vegetable crops. By contrast, the increase in CLD leaching was more notable in andosols and ferralsols with a relatively high sorption capacity and CLD stocks at the time of the change in the cropping system. In these cases, a reduction of CLD sorption through the decrease in SOC enhanced CLD release and leaching. For this reason, the increase in CLD leaching was the highest for andosol 1 (Figure 4c) which had the highest SOC and KOC values (Table 1). These results suggest that the early substitution of banana systems by vegetable crops following the financial crisis might have contributed to increasing environmental pollution and the risk of crop contamination during the past two decades.

    Although the partial replacement of banana by vegetables operated in 1993 mainly involved smallholder farmers [17], the land area cultivated with vegetables and food crops on polluted soil has not yet been fully identified. This is a strategic priority for local authorities [4]. Undoubtedly, to perform CLD analyses on each plot in the territory would be a very long and expensive operation. In a recent study, Rochette et al. (2020) proposed a method for the demarcation of CLD-polluted areas based on surface water analyses and the delimitation of hydrological units at the watershed and sub-watershed scales [29]. Such an approach might be helpful to identify plots that should be targeted as a priority during soil testing campaigns, which should also include an analysis of the subsurface layer.

    Coupling the WISORCH model of CLD dynamics and the MorGwanik model accounting for SOC evolution enabled us to highlight the fact that a change of cropping system when the topsoil reaches the threshold for CLD content established by local authorities might increase the risk of crop contamination because of the contribution of the subsurface layer. Although the change in cropping system might not markedly impact CLD fluxes towards other environmental compartments at that time, the impact would much higher for vegetable cropping systems installed in the early 1990s following the financial crisis that affected the banana sector. Under these systems, SOC stocks decline progressively due to intensive soil tillage that reduces the sorption capacity of soils and increases CLD leaching. This would be particularly important for nitisols and ferralsols with a mineralogy inducing a low CLD sorption, which favors CLD leaching.

    Our results indicate that in order to reduce the risk of crop contamination and environmental pollution, current soil testing campaigns should include CLD monitoring of not only the topsoil but also of the subsurface layer, paying special attention to temporal changes of SOC stocks in soils cultivated with vegetables. This implies that continuous long-term observations, rather that one-off sampling, are essential to test model predictions for different soil type × cropping system scenarios and improve the accuracy of decision tools at the regional scale. Meanwhile, until decontamination techniques become available, polluted soils should be managed carefully to prevent further CLD and transformation products leaching linked to changes in land use and CLD degradation. To achieve this, farming practices such as reduced soil tillage and the use of organic amendments should be promoted in order to maintain or increase SOC stocks and the sorption properties of polluted soils.

    This study formed part of the TropEmis Project funded by the Reacctif Program of the French Agency for Ecological Transition (ADEME) (Grant 410-00159), the European Regional Development Fund (FEDER) (Grant 410-00160) and the Regional Council of Guadeloupe (Grant 410-00161). We would like to thank Victoria Hawken for reviewing the English manuscript.

    The authors declare no conflicts of interest in this paper.



    [1] Mouvet C, Collet B, Gaude JM, et al. (2020) Physico-chemical and agronomic results of soil remediation by in situ chemical reduction applied to a chlordecone-contaminated nitisol at plot scale in a French Caribbean banana plantation. Environ Sci Pollut Res 27: 41063-41092. doi: 10.1007/s11356-020-07603-z
    [2] Dromard CR, Devault DA, Bouchon-Navaro Y, et al. (2019) Environmental fate of chlordecone in coastal habitats: recent studies conducted in Guadeloupe and Martinique (Lesser Antilles). Environ Sci Pollut Res
    [3] Cabidoche YM, Achard R, Cattan P, et al. (2009) Long-term pollution by chlordecone of tropical volcanic soils in the French West Indies: a simple leaching model accounts for current residue. Environ Pollut 157: 1697-1705. doi: 10.1016/j.envpol.2008.12.015
    [4] Letchimy S, Benin J (2019) Impact économique, sanitaire et environnemental de l'utilisation du chlordécone et du paraquat comme insecticides agricoles dans les territoires de Guadeloupe et de Martinique, sur les responsabilités publiques et privées dans la prolongation de leur autorisation et évaluant la nécessité et les modalités d'une indemnisation des préjudices des victimes et de ces territoires. Report to the National Assembly of France, No 2440. Available from: https://www.assemblee-nationale.fr/dyn/15/rapports/cechlordec/l15b2440-ti_rapport-enquete.pdf (in French).
    [5] DAAF (2019) Actualisation de la carte de contamination des sols par la chlordécone. Available from: https://daaf.guadeloupe.agriculture.gouv.fr/Actualisation-de-la-carte-de (in French).
    [6] Cabidoche YM, Lesueur-Jannoyer M (2012) Contamination of harvested organs in root crops grown on chlordecone-polluted soils. Pedosphere 22: 562-571. doi: 10.1016/S1002-0160(12)60041-1
    [7] Fernández-Bayo JD, Saison C, Voltz M, et al. (2013) Chlordecone fate and mineralisation in a tropical soil (andosol) microcosm under aerobic conditions. Sci Total Environ 463-464: 395-403. doi: 10.1016/j.scitotenv.2013.06.044
    [8] Mouvet C, Dictor MC, Bristeau S, et al. (2017) Remediation by chemical reduction in laboratory mesocosms of three chlordecone-contaminated tropical soils. Environ Sci Pollut Res 24: 25500-25512. doi: 10.1007/s11356-016-7582-4
    [9] Sansoulet J, Cabidoche YM, Cattan P (2007) Adsorption and transport of nitrate and potassium in an Andosol under banana (Guadeloupe, French West Indies). Eur J Soil Sci 58: 478-489. doi: 10.1111/j.1365-2389.2007.00904.x
    [10] Cattan P, Charlier JB, Clostre F, et al. (2019) A conceptual model of organochlorine fate from a combined analysis of spatial and mid- to long-term trends of surface and ground water contamination in tropical areas (FWI). Hydrol Earth Syst Sci 23: 691-709. doi: 10.5194/hess-23-691-2019
    [11] Levillain J, Cattan P, Colin F, et al. (2012) Analysis of environmental and farming factors of soil contamination by a persistent organic pollutant, chlordecone, in a banana production area of French West Indies. Agric Ecosyst Environ 159: 123-132. doi: 10.1016/j.agee.2012.07.005
    [12] Clostre F, Letourmy P, Lesueur-Jannoyer M (2017) Soil thresholds and a decision tool to manage food safety of crops grown in chlordecone polluted soil in the French West Indies. Environ Pollut 223: 357-366. doi: 10.1016/j.envpol.2017.01.032
    [13] Clostre F, Lesueur-Jannoyer M, Achard R, et al. (2014) Decision support tool for soil sampling of heterogeneous pesticide (chlordecone) pollution. Environ Sci Pollut Res 21: 1980-1992. doi: 10.1007/s11356-013-2095-x
    [14] Lott DE, Hammond VE (2013) Vegetable and fruit production. NebGuide No G2189. University of Nebraska. Available from: https://extensionpublications.unl.edu/assets/pdf/g2189.pdf
    [15] Sierra J, Causeret F, Diman JL, et al. (2015) Observed and predicted changes in soil carbon stocks under export and diversified agriculture in the Caribbean. The case study of Guadeloupe. Agric Ecosyst Environ 213: 252-264.
    [16] Sierra J, Causeret F, Chopin P (2017) A framework coupling farm typology and biophysical modelling to assess the impact of vegetable crop-based systems on soil carbon stocks. Application in the Caribbean. Agric Syst 153: 172-180.
    [17] Tillieut O (2006) Cartographie de la pollution des sols de Guadeloupe par la Chlordecone. Rapport technique 2005-2006. DAF & INRA Guadeloupe. Available from: https://daaf.guadeloupe.agriculture.gouv.fr/IMG/pdf/071005_DAF-SPV-OT_cartographie-pollution-sols-Guadeloupe-chlordecone_cle8da565.pdf (in French).
    [18] DSDS (2001) Pesticides et alimentation en eau potable en Martinique. Etat des lieux et position sanitaire. Bilan actualisé en octobre 2001. Direction de la Santé et du Développement Social de la Martinique, France. (in French).
    [19] Sabatier P, Mottes C, Cottin N, et al. (2021) Evidence of chlordecone resurrection by glyphosate in French West Indies. Environ Sci Technol 55: 2296-2306. doi: 10.1021/acs.est.0c05207
    [20] Fernández-Bayo JD, Saison C, Geniez C, et al. (2013) Sorption characteristics of chlordecone and cadusafos in tropical agricultural soils. Curr Org Chem 2013: 2976-2984. doi: 10.2174/13852728113179990121
    [21] Sierra J, Ozier-Lafontaine H, Dufour L, et al. (2006) Nutrient and assimilate partitioning in two tropical maize cultivars in relation to their tolerance to soil acidity. Field Crop Res 95: 234-249. doi: 10.1016/j.fcr.2005.03.002
    [22] Raphaël L, Sierra J, Recous S, et al. (2012) Soil turnover of crop residues from the banana (Musa AAA cv. Petite-Naine) mother plant and simultaneous uptake by the daughter plant of released nitrogen. Eur J Agron 38: 117-123.
    [23] Poeplau C, Don A, Vesterdal L, et al. (2011) Temporal dynamics of soil organic carbon after land-use change in the temperate zone - carbon response functions as a model approach. Global Change Biol 17: 2415-2427. doi: 10.1111/j.1365-2486.2011.02408.x
    [24] Balesdent J, Basile-Doelsch I, Chadoeuf J, et al. (2018) Atmosphere-soil carbon transfer as a function of soil depth. Nature 559: 599-602. doi: 10.1038/s41586-018-0328-3
    [25] Dorel M, Lakhia S, Pététin C, et al. (2010) No-till banana planting on crop residue mulch: effect on soil quality and crop functioning. Fruits 65: 55-68. doi: 10.1051/fruits/20010001
    [26] Woignier T, Fernandes P, Soler A, et al. (2013) Soil microstructure and organic matter: keys for chlordecone sequestration. J Hazard Mater 262: 357-364. doi: 10.1016/j.jhazmat.2013.08.070
    [27] Woignier T, Rangon L, Clostre F, et al. (2019) Physical limitation of pesticides (chlordecone) decontamination in volcanic soils: fractal approach and numerical simulation. Environ Sci Pollut Res 27: 40980-40991. doi: 10.1007/s11356-019-05899-0
    [28] Crabit A, Cattan P, Colin F, et al. (2016) Soil and river contamination patterns of chlordecone in a tropical volcanic catchment in the French West Indies (Guadeloupe). Environ Pollut 212: 615-626. doi: 10.1016/j.envpol.2016.02.055
    [29] Rochette R, Bonnal V, Andrieux P, et al. (2020) Analysis of surface water reveals land pesticide contamination: an application for the determination of chlordecone-polluted areas in Guadeloupe, French West Indies. Environ Sci Pollut Res 27: 41132-41142. doi: 10.1007/s11356-020-10718-y
    [30] Charlier JB, Cattan P, Voltz M, et al. (2009) Transport of a nematicide in surface and groundwaters in a tropical volcanic catchment. J Environ Qual 38: 1031-1041. doi: 10.2134/jeq2008.0355
    [31] Chambre d'Agriculture de Martinique (2014) Fiche d'itinéraire technique Banane (export). Available from: https://martinique.chambre-agriculture.fr/fileadmin/user_upload/National/FAL_commun/publications/Outre-Mer/FIT_Banane_Export_2014.pdf (in French).
    [32] Schneider F, Dona A, Hennings I, et al. (2017) The effect of deep tillage on crop yield - What do we really know? Soil Till Res 174: 193-204.
  • Environ-08-04-20-s1.pdf
  • This article has been cited by:

    1. Annelise Tran, Victor Dufleit, Vincent Bonnal, Pascal Degenne, Jeremy Lavarenne, Lucie Lecat, Philippe Cattan, Agri-STAMP: A spatial model representing the use and fate of pesticides over the long term, 2025, 959, 00489697, 178225, 10.1016/j.scitotenv.2024.178225
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2480) PDF downloads(131) Cited by(1)

Figures and Tables

Figures(4)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog