The importance of climate scenarios in setting the productivity indexes in the Pampa Arenosa sub-region of the province of Buenos Aires, Argentina

Importância dos cenários climáticos na determinação de índices de produtividade na sub-região Pampa Arenosa da província de Buenos Aires, Argentina

Abstract

Starting in the 1970's, the Pampa Arenosa sub-region experienced an increasing water regime that generated an increased area for dryland farming. With increasing agricultural activity, examining the land is an essential strategic tool for its planning. The objective of this work was to highlight the importance of considering several climate scenarios when setting the productivity index in the sector of longitudinal dunes in the Pampa Arenosa sub-region in the province of Buenos Aires, Argentina. The climate scenarios were set in relation to the shifts in average rain values, according to the Pettitt Test. The land were classified according to their productivity index. It was found that the productivity index of the land increased with the rain, reaching its highest value in the period right after the abrupt shift. Land of moderate productive capacity whose productivity index values are between 65 and 51 comprised the majority of the area of our study.

Keywords:
agro climatology; land evaluation; Pettitt test; rainfall.

Resumo

A partir da década de 70, a sub-região Pampa Arenosa experimentou um aumento no seu regime de água, o que permitiu o crescimento da área dedicada à agricultura de secas. Em vista da intensificação da atividade agrícola, uma avaliação da terra é uma ferramenta estratégica essencial para o planejamento. O objetivo do trabalho foi destacar a importância de considerar diferentes cenários climáticos, na determinação do Índice de Produtividade, no setor de dunas longitudinais da sub-região Pampa Arenosa na Província de Buenos Aires, Argentina. Os cenários climáticos foram estabelecidos de acordo com as mudanças nos valores médios de precipitação, pelo teste Pettitt. As terras foram classificadas pelo índice de produtividade. Verificou-se que o índice de produtividade do solo aumentou com o aumento da precipitação, atingindo sua máxima expressão climática no período após a mudança abrupta. As terras de capacidade produtiva moderada com valores de Índice de Produtividade entre 65-51 ocuparam a maior área da área de estudo.

Palavras-chave:
agroclimatologia; avaliação da terra; precipitação; teste de Pettitt.

1. INTRODUCTION

The Pampa region, located in the east-central part of the country (30°S - 40°S y 56°O - 65°O), is Argentina’s main agricultural area. The Pampa Region has a humid temperate climate (Cf in the Köppen-Geiger classification, as revised by Kottek et al., 2006). East winds predominate, driven by a semi-permanent anticyclone from the coast of Brazil. After being drawn across the Brazilian coastline, maritime subtropical air moves southeast, reaching up to 40° latitude in summer and about 30° latitude in winter. In this way, the Pampa Region receives sea winds throughout the year, with a moisture gradient decreasing from east to west (Pérez et al., 2015).

Several studies indicate that the westward advance of the agricultural frontier in the Pampas during the last quarter of the twentieth century (Viglizzo et al., 1995) was partly a consequence of increased rainfall (Castañeda and Barros 1994; Pérez et al., 2011; 2015; Pérez and Sierra, 2012). This increase in precipitation acted synergistically with an increasing demand from international markets (Trigo, 2005) and technological innovations (Satorre, 2005).

Some authors believe that the above-mentioned increase in rainfall is permanent. They attribute it to increased energy in the climate system caused by global warming. In their view, this has led to an increased thermal regime throughout the country, affecting the whole of its climate (Carril et al., 1997; Minetti et al., 2003; Barros, 2004). In contrast, others (Suriano and Ferpozzi, 1993; Roberto et al., 1994; Pérez et al., 2003; 2011; Sierra and Pérez, 2006, Pérez et al., 2015) suggest that these changes are reversible. In their view, the Pampas have a long-term water cycle with wet and dry phases separated by transition periods during which the agricultural frontier either advances or retreats.

Starting in the 1970’s, the Pampa Arenosa sub-region experienced an increasing water regime that generated an increasing area for dryland farming.

With increasing agricultural activity, examining the land becomes an essential strategic tool for its planning.

Several methods are used when creating land evaluation systems in terms of adaptability and/or vulnerability (De la Rosa et al., 2004). In Argentina, the most frequent quantitative indirect classification is the productivity index (PI). This system was developed initially by FAO (Riquier et al., 1970) and adopted by the National Institute of Agricultural Technology (in Spanish, Instituto Nacional de Tecnología Agropecuaria) according to the local agro-ecological conditions. The productivity index relates the property values or levels that have some influence on land productivity (Irigoin, 2011).

The objective of this work was to highlight the importance of considering several climate scenarios when setting the productivity index in the sector of longitudinal dunes in the Pampa Arenosa sub-region in the province of Buenos Aires, Argentina.

2. MATERIAL AND METHODS

The area of study comprises the second-level administrative subdivisions of Nueve de Julio, Carlos Casares, Pehuajó and Trenque Lauquen, which are located in the province of Buenos Aires, Argentina, in the sector of longitudinal dunes of the Pampa Arenosa sub-region. Annual rainfall data from four towns in the Pampa Arenosa sub-region for the period 1918-2011 was used (Table 1). Data was provided by the National Weather Service and the National Institute of Agricultural Technology.

The homogeneity of the precipitation series was tested using Alexandersson and Moberg’s (1997) Standard Normal Homogeneity Test (SNHT) on AnClim software (Štĕpánek, 2006). The test was applied to series of annual values, using the average annual rainfall of each sub-region as a reference series.

Table 1.
Locations of towns in the Pampa Arenosa sub-region.

Town/City	Latitude (S)	Longitude (W)	Altitude (msl)
Nueve de Julio	35° 27´	60° 53´	76
Carlos Casares	35° 55´	61° 20´	88
Pehuajó	35° 48´	61° 54´	85
Trenque Lauquen	35° 58´	62° 43´	80

With the homogeneous precipitation series of each town, shifts in the average values were found using the Pettitt test. The climate scenarios were defined using the temporal series of rainfall for the periods before and after the shifts in the average values, identified as the dry period (DP) and the humid period (HP), respectively.

The productivity index (PI) was obtained by using the changes introduced by the National Institute of Agricultural Technology (INTA) for the Pampa region (Sobral and Nakama 1988; Maccarini, 1990; Sobral et al., 2010). This parametric index comprises the following factors: Macro-climate condition (H); drain (D); effective depth (Pe); surface texture (Ta); sub-surface texture (Tb); salinity (Sa); sodicity (Na); organic matter (Mo); cation-exchange capacity (T), current and potential water and wind erosion (E).

The quantitative expression of the productivity index is calculated by multiplying the factors listed above. This multiplication intends to highlight the influence of each factor on the final value of the index (Equation 1):

$\mathbf{P}\mathbf{I}=\mathbf{}\mathbf{H}\mathbf{}\mathbf{*}\mathbf{}\mathbf{D}\mathbf{}\mathbf{*}\mathbf{}\mathbf{P}\mathbf{e}\mathbf{}\mathbf{*}\mathbf{}\mathbf{T}\mathbf{a}\mathbf{}\mathbf{*}\mathbf{}\mathbf{T}\mathbf{b}\mathbf{}\mathbf{*}\mathbf{}\mathbf{S}\mathbf{a}\mathbf{}\mathbf{*}\mathbf{}\mathbf{N}\mathbf{a}\mathbf{}\mathbf{*}\mathbf{M}\mathbf{o}\mathbf{}\mathbf{*}\mathbf{}\mathbf{T}\mathbf{}\mathbf{*}\mathbf{}\mathbf{E},$ (1)

The results of the PI calculation are positive values, the highest of which is equal to 100, where the highest value represents the greater productive capacity of the land.

2.1. Climate factor of the PI

The climate factor H of the PI was calculated following the INTA methodology (Sobral and Nakama 1988; Maccarini, 1990).

The climate factor was calculated using Thornthwaite’s global humidity index (Im), the soil moisture regimes (Moscatelli, 1991; Van Wambeke and Scoppa, 1980) and the 16°C isotherm (Table 2).

El Im (Thornthwaite, 1948; 1955) was obtained from the difference between the annual water excess (humidity index, Ih) and the deficiencies (aridity index, Ia) as the result of the water balance for maximum water storage capacity in soil of 100 mm and its relation with the potential evapotranspiration (ETP) estimated by Thornthwaite, where Im= (Ih-Ia)/ETP (Castellví Sentí and Castillo, 2001). The value for Im for the area as a whole was obtained from the average of the value for Im for each area.

Table 2.
Conversion of climate characteristics into coefficient H values (Macroclimate condition) (Maccarini 1990).

Global humidity index (Im)	Average annual temperature (°C)	Macro-climate condition (H)
> 10	≥16	100
10 to 0	≥16	95
0 to -10	≥16	90
-10 at ustic/aridic	≥16	75

2.2. Edaphics factors of the PI

The edaphic information for the different factors that make up the PI is the following: drain, alkalinity, salinity, texture, cation-exchange capacity, content of organic matter, current and potential water and wind erosion. All this data was obtained from soil charts (INTA, 1990).

The PI value was obtained for each taxonomic unit defined at a series level, from the use of conversion tables (Irigoin, 2011), where each edaphic factor participating in the parametric equation is classified into categories or ranges with their respective coefficients with values between 0 and 1 (Sobral et al., 2010, Sobral and Nakama 1988).

The PI value of the cartographic unit was calculated by weighing the PI values of each taxonomic unit and in relation to the percentage of area covered by each of them.

The PI values of the cartographic units were divided into six productive categories: very high (>80), high (80-66), moderate (65-51), moderately low (50-36), low (35-20) and very low (<20) (Irigoin, 2011).

3. RESULTS AND DISCUSSION

3.1. Homogeneity Test

According to the Homogeneity Test (SNHT), two of the annual precipitation series available (Table 3) showed a T value lower than the critical value (Khaliq and Quarda, 2007), and can be considered homogeneous at the level of significance α = 0.05. The other two series showed T values higher than the critical value and were thus considered non-homogeneous. For the analysis, the non-homogeneous series were adjusted according to the method proposed by Alexandersson (1986).

Table 3.
Test results of the SNHT applied to annual precipitation series for the Pampa Arenosa sub-region.

Town/City	Shift year	T Value	Adjusted shift year	Adjusted T value
Nueve de Julio	1984	7.107
Carlos Casares	2010	4.157
Pehuajo	1935	10.747 *	1954	2.191
Trenque Lauquen	1972	14.080 *	2010	4.647

• (* indicates that the T value exceeds 95%).

3.2. Pettitt Test

The shifts in the values of the annual precipitation average according to the Pettitt Test (Pettitt, 1979) (Table 4) were produced in 1962 and 1965, thus defining the dry scenarios, before the shift (1918-1962/65) and wet, after the shift (1962/65-2005).

Table 4.
Results after the Pettitt Test for the annual precipitation of the Pampa Arenosa sub-region.

Town/City	Shift year	Average annual rainfall (mm) Before the shift	Rainfall increase (%) After the shift
Nueve de Julio	1962	893	16.8
Carlos Casares	1965	787	22.5
Pehuajo	1965	820	19.3
Trenque Lauquen	1962	708	38.3

The average of Nueve de Julio (Figure 1) showed a positive abrupt change, with the average annual precipitation rising from 893 mm during the 1918-1962 sub-period, to 1043 mm during the 1963-2011 sub-period.

Figure 1.
Annual precipitation and means for sub-periods in the Nueve de Julio by Pettitt’s method.

The average of Carlos Casares (Figure 2) showed a positive abrupt change, with the average annual precipitation rising from 787 mm during the 1918-1965 sub-period, to 964 mm during the 1966-2011 sub-period.

Figure 2.
Annual precipitation and means for sub-periods in the Carlos Casares by Pettitt’s method.

The average of Pehuajo (Figure 3) showed a positive abrupt change, with the average annual precipitation rising from 820 mm during the 1918-1965 sub-period, to 978.2 mm during the 1966-2011 sub-period.

Figure 3.
Annual precipitation and means for sub-periods in the Pehuajo by Pettitt’s method

The average of Trenque Lauquen (Figure 4) showed a positive abrupt change, with the average annual precipitation rising from 708 mm during the 1918-1962 sub-period, to 979.2 mm during the 1963-2011 sub-period.

Figure 4.
Annual precipitation and means for sub-periods in the Trenque Lauquen by Pettitt’s method.

3.3. Productivity Index

Table 5 shows the productivity indices in the sector of longitudinal dunes in the Pampa Arenosa sub-region in the province of Buenos Aires, Argentina. In the study area, approximately 458,700 hectares have values of productivity index between 65 and 51, indicating that these lands have some type of permanent limitations for the production of common crops. On the other hand, these IP values would indicate that the feasible yield to obtain will correspond between 65 and 50% of the optimum yield for the region (Sobral and Nakama, 1988). A similar area (460,720 hectares) corresponds to lands of moderately low productive capacity with IP values of 50-36. These two IP categories represent more than half of the study area, approximately 54.7%. At the same time, lands with severe and very severe limitations for agriculture occupy 668,320 hectares with IP values of less than 35.

Table 5.
Classification of the lands according to the Productivity Index (PI) for the study area and for each Town/City, expressed in % of occupation and km².

		Productivity Index (PI)
Town/City	Occupation	100-81	80-66	65-51	50-36	35-21	20-0
Nueve de Julio	km2	88.5	623.1	1717.8	215.2	1026.2	614.9
	%	2.1	14.5	40.1	5.0	23.9	14.4
Carlos Casares	km2	0.0	200.0	830.4	28.3	943.5	516.5
	%	0.0	7.9	33.0	1.1	37.5	20.5
Pehuajó	km2	0.0	3.4	1312.6	1381.7	740.6	1081.8
	%	0.0	0.1	29.0	30.6	16.4	23.9
Trenque Lauquen	km2	0.0	0.0	726.2	2982.0	1340.1	419.6
	%	0.0	0.0	13.3	54.5	24.5	7.7
Study area	km2	88.5	826.5	4587.0	4607.2	4050.4	2632.8
	%	0.5	4.8	27.3	27.4	24.1	15.7

According to the PI values for both climate scenarios and their location, we can observe a shift in the productivity of the land as the consequence of shifts in the climate factor.

The PI weighted average obtained for the period right after the abrupt shift for the Nueve de Julio area was 48.5, for Carlos Casares was 43.7, for Pehuajó was 39.9 and for Trenque Lauquen was 46.1 (Figure 5). On the other hand, the PI values for the period before were 39.4, 36.1, 36.9 and 34.4, respectively, according to the transect of these areas, from East to West (Figure 6).

Figure 5.
Distribution pattern of the classes of productivity index of the land (PI) for the period after the abrupt shift in the area of Nueve de Julio (A), Carlos Casares (B), Pehuajó (C) and Trenque Lauquen (D).

Figure 6.
Distribution pattern of the classes of productivity index of the land (PI) for the period before the abrupt shift in the area of Nueve de Julio (A), Carlos Casares (B), Pehuajó (C) and Trenque Lauquen (D).

In this way, the land resulted of moderately low productive capacity (50-36), both for the period after the abrupt shift and for the period before, thus showing that such land has permanent edaphic limitations for the production of conventional crops. The productivity index is interpreted as a proportion of the potential maximum yield of the most common crops in the region, eco-typically adapted, under a determined level of management.

4. CONCLUSION

The Pettitt Test (Pettitt, 1979), applied to the annual precipitation series of the Pampa Arenosa sub-region, showed that the sub-region had abrupt positive shifts during the second half of the 20^th century. This sub-region was considerably affected by such shifts. Indeed, they changed from a semi-arid steppe climate (Bs in the Köppen Climate Classification System) (Köppen, 1948) during dry phases to a dry-summer subtropical climate (Cw in the Köppen classification) during wet phases, which significantly affects agricultural capacity.

It was found that the PI of the lands increased with the rain, reaching its highest value in the period right after the abrupt shift for all the studied areas. The Trenque Lauquen area was the area most positively affected by increasing precipitation.

5. REFERENCES

• ALEXANDERSSON, H. A homogeneity test applied to precipitation data. International Journal of Climatology, v. 6, n. 6, p. 661-675, 1986. https://doi.org/10.1002/joc.3370060607
  » https://doi.org/10.1002/joc.3370060607
• ALEXANDERSSON, H.; MOBERG, A. Homogenization of Swedish temperature data. Part I: homogeneity test for linear trends. International Journal of Climatology, v. 17, p. 25-34, 1997. https://doi.org/10.1002/(SICI)1097-0088(199701)17:1%3C25::AID-JOC103%3E3.0.CO;2-J
  » https://doi.org/10.1002/(SICI)1097-0088(199701)17:1%3C25::AID-JOC103%3E3.0.CO;2-J
• BARROS, V. El cambio climático global. [S.l.]: Del Zorzal, 2004.
• CARRIL, A. F.; MENÉNDEZ, C. G.; NUÑEZ, M. M. Climate change scenarios over the South American region: an intercomparison of coupled general atmosphere-ocean circulation models. International Journal of Climatology, v. 17, n. 15, p. 1613-1633, 1997. https://doi.org/10.1002/(SICI)1097-0088(199712)17:15%3C1613::AID-JOC209%3E3.0.CO;2-8
  » https://doi.org/10.1002/(SICI)1097-0088(199712)17:15%3C1613::AID-JOC209%3E3.0.CO;2-8
• CASTAÑEDA, M. E.; BARROS, V. Las tendencias de la precipitación en el Cono Sur de América al este de los Andes. Meteorológica, v. 19, n. 1-2, p. 23-32, 1994.
• CASTELLVÍ SENTÍS, F.; CASTILLO, F. Agrometeorología. [S.l.]: Mundi-Prensa, 2001, 517 p.
• DE LA ROSA, D.; MAYOL, F.; DÍAZ-PEREIRA, E.; FERNÁNDEZ, M.; DE LA ROSA, D. J. R. A land evaluation decision support system (MicroLEIS DSS) for agricultural soil protection with special reference to Mediterranen region. Environmental Modelling & Sofware, v. 19, p. 929-942, 2004. https://doi.org/10.1016/j.envsoft.2003.10.006
  » https://doi.org/10.1016/j.envsoft.2003.10.006
• INTA. Atlas de Suelos de la República Argentina. Escala 1: 500.000 y 1: 1.000.000. Tomo I y II. Proyecto PNUD ARG. 85/019. Buenos Aires, 1990, 667 p.
• IRIGOIN, J. Sistemas de evaluación de tierras y elaboración de modelos de aptitud de uso agrícola, para distintos escenarios climáticos, en un sector de la subregión Pampa Arenosa (Provincia de Buenos Aires, Argentina). 2011. 160 p. Tesis (Maestría de Agronomía) - Universidad de Buenos Aires, Buenos Aires, 2011.
• KHALIQ, M. N.; QUARDA, T. B. M. J. On the critical values on the standard normal homogeneity test (SNHT). International Journal of Climatology, v. 27, p. 681-687, 2007. https://doi.org/10.1002/joc.1438
  » https://doi.org/10.1002/joc.1438
• KÖPPEN, W. Climatología. Mexico: Fondo de Cultura Económica, 1948. 478 p.
• KOTTEK, M.; GRIESER, J.; BECK, C.; RUDOLF, B.; RUBEL, F. World Map of the Köppen-Geiger climate classification updated. Meteorologische Zeitschrift, v. 15, n. 3, p. 259-263, 2006. https://doi.org/10.1127/0941-2948/2006/0130
  » https://doi.org/10.1127/0941-2948/2006/0130
• MACCARINI, G. D. Utilización del Método de evaluación de tierras paramétrico de Bramao y Riquier y su adaptación por INTA en un sector de la Región Pampeana Húmeda. 1990. Tesis (Maestría en Ciencias del Suelo) - Escuela de Posgrado, Facultad de Agronomía, Universidad de Buenos Aires, Buenos Aires,1990.
• MINETTI, J. L.; VARGAS, W. M.; POBLETE, A. G.; ACUNA, L. G.; CASAGRANDE, G. Non-linear trends and low frequency oscillations in annual precipitation over Argentina and Chile, 1931- 1999. Atmosfera, v. 16, p. 119-135, 2003.
• MOSCATELLI, G. Los suelos de la Región Pampeana. In: BARSKY, O. (Ed.). El desarrollo agropecuario pampeano. Buenos Aires: INDEC-INTA-IICA, 1991. p. 1-76.
• PÉREZ, S.; SIERRA, E. M.; CASAGRANDE, G.; VERGARA, G.; BERNAL, F. Comportamiento de las precipitaciones (1918/2000) en el centro oeste de la provincia de Buenos Aires (Argentina). Revista De la Facultad de Agronomía de la Universidad Nacional de La Pampa, v. 14, n. 1-2, p. 39-46, 2003.
• PÉREZ, S.; SIERRA, E.; LÓPEZ, E.; NIZZERO, G.; MOMO, F.; MASSOBRIO, M. Abrupt changes in rainfall in the Eastern area of La Pampa Province, Argentina. Theoretical and Applied Climatology, v. 103, p. 159-165, 2011. https://doi.org/10.1007/s00704-010-0290-y
  » https://doi.org/10.1007/s00704-010-0290-y
• PÉREZ, S.; SIERRA, E. Changes in rainfall patterns in the eastern area of La Pampa province, Argentina. Revista Ambiente & Agua, v. 7, n. 1, p. 24-35, 2012. http://dx.doi.org/10.4136/ambi-agua.692
  » http://dx.doi.org/10.4136/ambi-agua.692
• PÉREZ, S.; SIERRA, E.; MOMO, F.; MASSOBRIO, M. Changes in Average Annual Precipitation in Argentina’s Pampa Region and Their Possible Causes. Climate, v. 3, p. 150-167, 2015. http://d.doi.org/10.3390/cli3010150
  » http://d.doi.org/10.3390/cli3010150
• PETTITT, A. N. A non-parametric approach to the change point problem. Applied Statistics, v. 28, n. 2, p. 126-135, 1979. http://dx.doi.org/10.2307/2346729
  » http://dx.doi.org/10.2307/2346729
• RIQUIER, J.; BRAMAO, L.; CORNET, S. P. A new system of soil appraisal in terms of actual and potential productivity. Roma: FAO, 1970.
• ROBERTO, Z. E.; CASAGRANDE, G.; VIGLIZZO, E. Lluvias en la Pampa Central: tendencia y variaciones del siglo. Cambio Climático y Agricultura Sustentable en la Región Pampeana. Boletine INTA Centro Regional La Pampa-San Luis, n. 2, 1994.
• SATORRE, E. H. Cambios tecnológicos en la agricultura argentina actual. Ciencia Hoy, v. 15, n. 87, p. 24-31, 2005.
• SIERRA, E. M.; PÉREZ, S. P. Tendencia del régimen de precipitación y el manejo sustentable de los agroecosistemas: estudio de un caso en el noroeste de la provincia de Buenos Aires, Argentina. Revista de Climatología, v. 6, p. 1-12, 2006.
• SOBRAL, R. E.; NAKAMA, V. Índices de productividad, método paramétrico para evaluación de tierras. In: CONGRESO ARGENTINO DE LA CIENCIA DEL SUELO, 12., 12 al 16 septiembre 1988, Corrientes. Actas […] [S.l.]: Comisión mineralogía, génesis, clasificación y cartografía de suelos, 1988. 259 p.
• SOBRAL, R. E.; NAKAMA , V.; DE ANTUENO, L. Actualización de los índices de productividad de los suelos de la provincia de Buenos Aires. In: CONGRESO ARGENTINO DE LA CIENCIA DEL SUELO, 22., 31 de mayo al 4 de Junio 2010, Rosario, Argentina. Actas […] [S.l.]: Comisión mineralogía, génesis, clasificación y cartografía de suelos, 2010.
• SURIANO, J. M.; FERPOZZI, L. H. Los cambios climáticos en la Pampa también son historia. Todo es Historia, n. 306, p. 8-25, 1993.
• ŠTĔPÁNEK, P. AnClim Software for time series analysis. Masaryk: Masaryk University, 2006.
• THORNTHWAITE, C. W. An approach toward a rational classification of climate. The Geografical Review, v. 38, p. 55-94, 1948. http://dx.doi.org/10.2307/210739
  » http://dx.doi.org/10.2307/210739
• THORNTHWAITE, C. W; HARE, F. K. Climate classification in Forestry. Unasylva, v. 9, n. 2, 1955.
• TRIGO, E. Consecuencias económicas de la transformación agrícola. Ciencia Hoy, v. 15, n. 87, p. 46-51, 2005.
• VAN WAMBEKE, A.; SCOPPA, C. Las taxas climáticas de los suelos argentinos. Castelar: Intacirn, 1980.
• VIGLIZZO, E. F.; ROBERTO, Z. E.; FILIPPIN, M. C.; PORDOMINGO, A. J. Climate variability and agroecological change in the Central Pampas of Argentina. Agriculture Ecosystems and Environment, v. 55, p. 7-16, 1995. https://doi.org/10.1016/0167-8809(95)00608-U
  » https://doi.org/10.1016/0167-8809(95)00608-U

Publication Dates

• Publication in this collection
  11 Mar 2019
• Date of issue
  2019

History

• Received
  04 June 2018
• Accepted
  18 Dec 2018