SPATIAL DISTRIBUTION OF SPONTANEOUS PRODUCTION OF Attalea maripa (Aubl) Mart. IN THE LOW TOCANTINS

DISTRIBUIÇÃO ESPACIAL DA PRODUÇÃO ESPONTÂNEA DE Attalea maripa (Aubl) Mart. DO BAIXO TOCANTINS

ABSTRACT

The production of biofuels from native Amazonian palm trees, belonging to the family of Arecaceae, like Attalea maripa (Aubl) Mart. (Maripa palm), are objects of studies of alternative energy sources for the biofuels production. Thus, the aim of this research was to evaluate the spatial dynamics of agronomic variables of Maripa palm crop production associated with the characteristics of the sampling sites in a spontaneous growth area in the Low Tocantins region. 100 palms of A. maripa were selected, of spontaneous growth and in production phase, on a 10 ha area in the municipality of Cametá / PA, in which the morpho-agronomical variables were collected (MIN, FIN, ABN, LN, DBH, AWAB and AFW), from January 2014 to June 2014. The semivariograms were produced for each variable, which presented a spatial dependence range of 75 meters to MIN, FIN, AWAB and AFW, 80 meters to LN and DBH and 98 meters to ABN. The area of reach found was 17,671.5 m² for MIN, FIN, AWAB and AFW variables, 20,106.2 m² for LN and DBH and 30,171.9 m² for ABN. The data of the variables were adjusted to the spherical model presenting determination coefficient (R²) that ranged from 0.57 to 0.99. The spatial dependence index (SDI) for MIN and FIN variables was 0.49 and 0.31, respectively, indicating moderate spatial dependence. For ABN, LN, DBH and AWAB variables ranged from 0.17 to 0.22, resulting in poor spatial dependence. The AFW variable presented SDI of 0.77. For the visual analysis of the spatial distribution of the studied variables, kriging maps were created.

Index terms
Spatial dependence; Maripa palm; Kriging; Semivariogram; Variables morphoagronomics

RESUMO

A produção de biocombustível a partir de palmeiras nativas da Amazônia, pertencentes à família das Arecaceae, como a Attalea maripa (Aubl) Mart. (inajá), são objetos de estudos de fontes energéticas alternativas para a produção de biocombustíveis. Dessa forma, o objetivo desta pesquisa foi avaliar a dinâmica espacial das variáveis morfoagronômicas de produção da cultura do inajazeiro associado às características dos locais de amostragem, em uma área de crescimento espontâneo, na região do Baixo Tocantins. Foram selecionadas 100 palmeiras de A. maripa, de crescimento espontâneo e em fase de produção, em uma área de 10 ha, no município de Cametá-PA, nas quais foram realizados os levantamentos das variáveis morfoagronômicas (NIM, NIF, NCA, NF, DAP, PMCA e PMF), no período de janeiro de 2014 a junho de 2014. Foram elaborados os semivariogramas para cada variável estudada, que apresentaram um alcance de dependência espacial de 75 metros para NIM, NIF, PMCA e PMF; 80 metros para NF e DAP, e 98 metros para NCA. Já a área de alcance encontrada foi de 17.671,5 m² para as variáveis NIM, NIF, PMCA e PMF; 20.106,2 m² para NF e DAP, e 30.171,9 m² para NCA. Os dados das variáveis se ajustaram-se ao modelo esférico, apresentando coeficiente de determinação (R²) que variou de 0,57 a 0,99. O índice de dependência espacial (IDE) para as variáveis NIM e NIF foi de 0,49 e 0,31, respectivamente, que indica dependência espacial moderada. Para as variáveis NCA, NF, DAP e PMCA, variaram de 0,17 a 0,22, resultando em dependência espacial fraca. Já a variável PMF apresentou IDE com valor 0,77. Para análise visual da distribuição espacial das variáveis estudas, foram elaborados mapas de krigagem.

Termos para indexação
Dependência espacial; Inajá; Krigagem; Semivariograma; Variáveis morfoagronômicas

INTRODUCTION

Negative environmental impacts in Amazonia, caused mainly by disordered deforestation, inadequate agricultural practices, livestock and illegal occupation of public land, have led to the disappearance of large areas of primary forests and strong soil degradation, resulting in the emergence of large almost homogeneous populations of palm trees, mainly from Attalea maripa (Aubl) Mart. (MIRANDA; RABELO, 2008).

Many species of the Arecaceae family have fruits and seeds with high oil content, with economic potential and a perspective of important appreciation for the region (PALLET, 2002).According to Miranda and Rabelo (2008), this oleaginous species, mainly destined to the feeding of traditional populations, can also form the basis of a self-sustaining technological development model.

Thus, among the palm trees known as “invasive”, deserves attention the Maripa palm, from the Arecaceae family, it is distributed in the tropical rainforest of the north of South America, mainly in Brazil, occurring in the states of Acre, Maranhão, Amapá, Amazonas, Pará, Rondônia and Roraima (Anderson et al., 2005). In the adult phase, the Maripa palm is an erect, simple and cylindrical palm tree with no tillers, solitary, reaching in the adult phase 13 m of height, almost 100 cm in diameter at breast height (DBH), presents composite leaf, inflorescence, interfoliar bunch protected by persistent spatula that are modified leaves, fruit of the drupe type, containing from one to three seeds (MATOS, 2010). From the fruit is extracted a clear and transparent oil, with excellent physicochemical properties for cooking, as well as raw material for biodiesel production (MOTA; FRANÇA, 2008).

The state of Pará has many areas of native populations of Maripa palm with great representativeness located mainly in the northeast region of Pará (Municipality of Cametá/PA).

However, in spite of the oleaginous potential of the fruits (Maripa palm), their natural populations are being gradually eliminated during the cleaning of pastures and by the practice of slash and burning of sponteneous secondary vegetation, used in the preparation of areas for agricultural crops planting.

Therefore, in order to provide the economic use of the Maripa palm and its incorporation in the list of commercial non-timber products, it is necessary to expand the research on the characteristics of bunches and fruits, as well as the production and productivity of this promising palm tree. Thus, the aim of this research was to evaluate the spatial dynamics of morpho-agronomical variables of the Maripa palm crop production associated to the characteristics of the sampling sites in an area of spontaneous growth in the Low Tocantins region.

MATERIALS AND METHODS

The studied area was located in the municipality of Cametá/PA, belonging to the Mesoregion of the Northeast of the state of Para and the Low Tocantins Region, with the geographic coordinates 02º44’54” south latitude and 49º03’42” west longitude of the Greenwich Meridian. The climate is classified as AMI, according to Köppen (1948), with an annual average temperature of approximately 25ºC and rainfall around 2,250 mm annually, with a higher concentration from January to June, representing around 80% of the total. The average air relative humidity is around 85%.

The study of the Maripa palm populations was carried out at the Continental Farm, which has an approximate area of 30 ha of spontaneous Maripa palm trees, with an estimated of 40 palms/ha, being selected an area of 10 ha and within that area were evaluated 100 palm trees of Maripa palms, free of pollination, with approximately 20 years old.

The palms were selected at random (Figure 1), obeying a minimum distance between plants of approximately 6 m within the selected area. Then, they were spatially located through the collection of their geographic coordinates (UTM).

The Morpho-agronomical parameters were evaluated in Maripa palms in areas with different uses, pasture, sponteneous secondary vegetation, flood and forest. The morpho-agronomical parameters studied were diameter at breast height (DBH), male inflorescence number (MIN) and female (FIN), leaves number (LN), androgynous bunch number (ABN), average weight of the androgynous bunch (AWAB) and average fruit weight (AFW).

To evaluate the morpho-agronomical variables spatially, a map of the sampling sites of the plants was elaborated (Figure 3A), and the types of environment observed in these areas were determined. The observed environments refer to flooded areas or areas subject to floods, pasture, sponteneous secondary vegetation and forest.

The DBH of the stem was measured at 1.30 m from the soil, using a tape measure, and then the circumference was measured. For the determination of the MIN and FIN, the inflorescences were collected with the aid of a sickle of the Malaysian type, sexed and checked.

The leaf number (LN) evaluation was carried out with the aid of a reinforced fiberglass ladder of two throws, measuring six meters, with all the leaves of each plant counted.

For the determination of the ABN and AWAB, the androgynous bunches were collected manually with the aid of a adjustable rod for six meters with a Malayan type sickle at the top and then transported by tractors to the premises of the factory of the green oils of Amazon plant, in which all bunches were counted, individually weighed on an English scale, then the length and circumference of the apex, middle and base of the bunch were measured with the aid of a millimeter scale.

After collecting the bunches, all of them were placed together in an industrial harvester designed for Maripa palm and the separation of the fruits from the ráqui and ráquilas (components of the empty bunch), the loose fruits were weighed in an English scale to estimate the average weight of the fruits.

In addition, the chemical characteristics of the soil were evaluated, and four simple soil samples were collected at 0-0.2 m depth, which subsequently formed a composite sample for each studied area. The phosphorus (P) and the potassium (K) were extracted by Mehlich-1 and determined by spectrophotometry and flame photometry respectively. The calcium (Ca), magnesium (Mg) and aluminum (Al) were extracted with KCl 1 mol L-1 and determined by titration.

The pH, in potentiometer in the soil: water relation of 1:2.5. The total organic carbon (TOC), by the volumetric method of oxidation with K2Cr2O7 and titration with ammonium ferrous sulfate, the organic matter was estimated by multiplying the TOC by 1.72 (EMBRAPA, 2011).

For the spatial distribution analysis of the production, the geostatistics was used, from the semivariogram modeling and the creation of kriging maps. For this procedure, the position information of each sample and the values observed at each point were used. Each plant had its known location, with its coordinates (UTM) in units of measures that enabled the calculation of the distance between plants (VIEIRA et al., 1983).

Then, semivariograms were elaborated, which defined the parameters referring to the semivariogram: the threshold (C or C0 + C1), the spatial variance (C1), the nugget effect (C0) and the range (a), through the estimation of the semivariogram, according to the equation below.

e01 $${?}^{*}\left(h\right)=\frac{1}{2N\left(h\right)}\underset{i=1}{\overset{N\left(h\right)}{?}}{\left[Z\left(x_i\right)-Z\left(x_i+h\right)\right]}^2>$$

Where N(h) is the number of pairs of measured values Z (xi) Z (xi + h), separated by a vector h. ?*(h) is the semivariance given in function of the distance, so it depends on the direction of h.

The data were tested to the experimental models with threshold and adjusted according to their behaviors. The models were determined according to the determination coefficient (R2), selecting those that presented higher values. The experimental models used were spherical, exponential, Gaussian and random (pure nugget effect).

From the adjustments, the parameters of the experimental semivariograms were calculated using the Spatial Dependency Indexes (SDI), k parameter, which is the relation between C1/(C0+C1) to determine the spatial dependence of each evaluation, and classified according to the intervals proposed by Zimback (2001) where the SDI values <0.25 is considered weak dependence, values of 0.25=SDI= 0.75 are considered moderate and, values> 0.75 are considered with strong spatial dependence.

For the creation of the kriging maps and evaluation of the study, it was necessary to interpolate the points sampled to obtain an estimate, given by the following equation:

e02 $$Z^{*}\left(x_0\right)=\underset{i=1}{\overset{N}{?}}{?}_{i}Z\left(x_i\right)$$

Where: z* = estimate; x0= linear combination of values of neighboring measures; N = number of measured values involved in the estimate z (xi); ?i = weight associated with each measured value.

For the elaboration, adjustments of semivariograms and the creation of kriging maps the Surfer 8.0 software was used.

RESULTS AND DISCUSSION

The different soil use systems showed low levels of Ca, Mg, P and K, (Table 1).The low concentration of exchangeable bases and high acidity, pH ranging from 4.3 to 5.3 are related to the leaching of significant quantities of the soil exchangeable bases caused by the high precipitation of the region.

In addition, the soils of the region are poor in primary minerals responsible for the replacement of exchangeable bases in the soil (CRAVO et al., 2007).

The spatial dependence of the variables was studied through semivariograms, presenting better adjust of the data to the spherical model, which provided the highest determination coefficient for all evaluations (Figure 2).The parameters of the semivariograms adjusted for the variables under study can be observed in Table 2.

Semivariograms were then elaborated for each studied variable (Figure 2A to G), which presented a spatial dependence range of 75 meters for MIN, FIN, AWAB and AFW, 80 meters for LN and DBH and 98 meters for ABN (Table 2).

The aggregation radius of the variables in the field is provided by the spatial dependence range, which allows making an analysis of the production within the studied area. From the radius of variables aggregation, it was estimated the area of reach for each variable which results in the influence of the same for the plants. The range area was 17,671.5 m² for MIN, FIN, AWAB and AFW, 20,106.2 m² for LN and DBH and 30,171.9 m² for ABN (Table 2).

From the range area, it was possible to identify, through the maps, which morpho-agronomical variables are related to each other and those with the soil chemical attributes content in each type in each type of soil use studied, being able to indicate in what type of soil use there is higher plant yield.

Regarding the spatial dependence index (SDI), for the MIN and FIN variables, moderate spatial dependence was observed, with indexes of 0.49 and 0.31, respectively. For the ABN, LN, DBH and AWAB variables ranged from 0.17 to 0.22, resulting in weak spatial dependence. The AFW variable presented SDI with a value of 0.77, with strong spatial dependence (Table 2).

The semivariograms adjusted to the models are shown in Figure 2, which presented adjustments to the spherical model, presenting determination coefficient (R²) ranging from 0.57 to 0.99, which justifies the aggregate distribution of the variables in the studied area.

Through the adjusted models to the semivariograms, it was possible to estimate the sampled values by the ordinary kriging method to construct the maps of the studied areas and the spatial distribution of morpho-agronomical characteristics (Figure 3), which allows the identification of homogeneous subareas, thus enabling an adequate management of the Maripa palm in the different areas.

According to the kriging maps, it was verified from the interpolation, that all variables formed coppices inside the area (Figure 3).

From the visual analysis of the sampling area maps and the kriging maps of the morphoagronomical variables, we observed that in areas subject to flooding, the MIN was higher, while in the pasture environment these values were reduced by half (Figure 3B), and may be related to the availability of water and/or light of the area.

According to Salm et al. (2015) the A. maripa show low tolerance to floods, obtaining a greater affinity for dry land, fact that may have influenced the greater emission of male inflorescences in this research.

For the FIN, the forest area provided higher values, but in the flooded and sponteneous secondary vegetation areas there were isolated sites that also presented high values. However, in the flooded area, isolated sites with reduced FIN were also observed (Figure 3C). These isolated sites may be related to the genotypic characteristics of certain palms and/ or such sites are close to other forest areas and/or nutrient availability.

For the ABN, in flooded areas, there is a significant decrease in the production of androgynous bunches, and in other places, these values were satisfactory (Figure 3D). However, the occurrence of isolated sites in the sponteneous secondary vegetation and pasture areas with low production of androgynous bunches was observed, which may be correlated with plant nutrition in the area. This fact is proven when the values obtained from AWAB and AFW were also observed, which were also low for flooded areas and at these isolated sites in the areas of sponteneous secondary vegetation and pasture, ranging from 17 to 42 kg for AWAB and 12 to 52 kg for AFW in all the areas.

In addition, the lowest average yields occurred in areas of floods (51.52 kg) and pasture (56.36 kg), when compared to areas of forest (71.60 kg) and sponteneous secondary vegetation (72.81 kg). According to Duarte (2008), the satisfactory index for the production of androgynous bunches varies from 3 to 5 bunches per palm tree per year with an average weight of 35 kg, corroborating with the results obtained in this research.

For the LN, the most significant values are observed in pasture and forest areas, and less significant in flooded and sponteneous secondary vegetation areas (Figure 3E). The lowest values for DBH also occurred in flooded and sponteneous secondary vegetation areas, but it was not generalized for all areas but for isolated sites of plants (Figure 3F). Beside that, the DBH AWAB and AFW variables presented higher values in the forest area (Figure 3F, Figure 3 G and Figure 3H) and pasture, only for the DBH variable (Figure 3F). According to Clark et al. (1995), the palm production is strongly influenced by environmental conditions and may lead to different responses due to environmental changes, thus, it can be stated that the low emission of leaves, DBH, AWAB and AFW in the areas of sponteneous secondary vegetation and flooded was influenced by soil and climate conditions.

Among the main management and maintenance activities of the plantation of some palm trees are: the crowning of the plants, skimming of the interlines, elimination of invasive plants, plant health, pruning, harvest and collection of the bunches and, the supervised pollination (BERTHAUD et al., 2000).Thus, this lack of management and maintenance of the A. maripa spontaneous palms in the areas of sponteneous secondary vegetation and flooded were possibly determinants facts that influenced the less expressive values in these areas. Thus, the adoption of cultural practices in a correct way and at the appropriate time would be of fundamental importance for the good development and production of these palm trees.

From the soils chemical characteristics at the four sampling sites and associating them with the low yields of the spontaneous populations found at some sampling sites in this region (Table 1 and Figure 3), we can assume that the mineral nutrition may also have exerted an influence on the production components of A. maripa. It is observed that other palm trees such as the palm oil crop have a significant demand for nutrients (SOH, 1997).Beside that, other factors such as nutritional requirement vary widely, depending on the yield required, the type of genetic material used, spacing, plant age, soil type and environmental factors. Environmental factors can significantly influence soil properties (CHAPLOT et al., 2001; OLLINGER et al., 2002) and reflect this on palm production. Kusnu et al. (1996) found an increase in the bunch weight when evaluating the growing effect of potassium doses on the palm oil crop.

The results presented in the spontaneous Maripa palm trees show that this material is in good reproductive state in pasture and forest areas (Figure 3). This population, for characteristics like AWAB, ABN, and LN, meets the selection standards used (DUARTE, 2008).

FIGURE 1
Location of the georeferenced palm trees. The points represent the spatial location of each plant within the studied area.

FIGURE 2
Semivariograms of an area of 10 ha of Maripa palm plantation (Attalea maripa (Aubl) Mart.) in the municipality of Cametá, PA, from January to March, 2014.

FIGURE 3
(A) Map of the sampling sites ( and maps of the spatial distribution of (B) male inflorescence number (MIN), (C) female inflorescence (FIN), (D) androgynous bunches (ABN), (E) leaves number (LN), (F) diameter at breast height (DBH), (G) average weight of androgynous bunches (AWAB) and (H) average fruit weight (AFW) of an area of 10 ha of maripa palm tree plantation (Attalea maripa (Aubl) Mart.) in the municipality of Cametá, PA, from January 2014 to March 2014.

TABLE 1
Result of the analysis of soil chemical attributes in the different environments within the study area.

• ¹
  Organic Matter, 2Cation Exchange Capacity, 3Base Sum, 4Base Saturation

TABLE 2
Semivariogram parameters adjusted to the models, determination coefficient (R²) and k parameter for the geostatistical analysis in an area of 10 ha of Maripa palm plantation (Attalea maripa (Aubl) Mart.) in the municipality of Cametá, PA, in the period from January 2014 to March 2014.

• ^a
  Calculated by p.r2, where p = 3.14 and r = a.
• ^b
  Relation between C1/ (C?+ C1).

CONCLUSION

The spherical model best represents the spatial distribution of the Maripa palm tree production variables, which indicates an aggregate distribution pattern of the morpho-agronomical variables in the field.

Spontaneous Maripa palm tree populations, although they are older in the field, have a good reproductive status, in relation to the Average Weight of Androgynous Bunches (AWAB), Average Fruit Weight (AFW) and Leaves Number (LN).

The pasture area provides higher values for diameter under the bifurcation point, leaves number and number of androgynous bunches, which indicates that the area is more suitable for a better production of Maripa palm tree.

The association of low yields of some morpho-agronomical characteristics with the sampling sites indicates the need to prioritize the selection of genotypes to increase the production of these palms in this region.

• ANDERSON, L. O.; SHIMABUKURO, Y. E.; LATORRE, M. L.; ESPÍRITO-SANTO, F. DEL B.; LIMA, A.; DUARTE, V.; ARAI, E.; DEFRIES, R. S.; MORTON, D.; HANSEN, M.; JASINSKY, E. Coleta de dados para estudos de cobertura arbórea, mapeamento da cobertura da terra, desmatamento e avaliação de unidades de conservação/territórios indígenas: metodologias de trabalhos de campo. Brasília: Instituto Nacional de Pesquisa Espacial, 2005. 37p.
• BERTHAUD, A.; NUNES, C. D. M.; BARCELOS, E.; CUNHA, R. N. V. da. Implantação e exploração da cultura do dendezeiro. In: VIÉGAS, I. de J. M.; MÜLLER, A. A. (Ed.). A cultura do dendezeiro na Amazônia Brasileira. Belém: Embrapa Amazônia Oriental, Manaus: Embrapa Amazônia Ocidental, 2000. p.193-227.
• CHAPLOT, V.; BERNOUX, M.; WATLER, C.; CURMI, P.; HERPIN, U. Soil carbon storage prediction in temperate hydromorphic soils using a morphologic index and digital elevation model. Soil Science, Baltimore, v.166, p.48–60, 2001.
• CLARK, D.A.; CLARK, D.B.; SANDOVAL, M.R.;CASTRO, C.M.V. Edaphic and human effects on landscape-scale distributions of tropical rain forest palms. Ecology, Brooklyn, v.76, p.2581-2594, 1995.
• CRAVO, M.S.; VIÉGAS, I.J.M.; BRASIL, E.C. (Ed.). Recomendações de adubação e calagem para o estado do Pará. Belém: Embrapa Amazônia Oriental, 2007. 262p.
• DUARTE, O. R. Avaliação quantitativa e análise dos parâmetros biológicos, químicos e físico-químicos de frutos de Maximiliana maripa (Aubl.) Drude (Inajá) como subsídio do estudo do potencial oleífero de populações promissoras para o Estado de Roraima. 2008. 145 f. Tese (Doutorado em Química) - Universidade Federal do Amazonas, Instituto Nacional de Pesquisas da Amazônia, Manaus, 2008.
• EMBRAPA - Empresa Brasileira de Pesquisa Agropecuária. Manual de análises químicas de solos, plantas e fertilizantes. 2. ed. Brasília, DF, 2011. 627p.
• KÖPPEN, W. Climatologia: con um estúdio de los climas de la Tierra. México: Fondo de Cultura Economica, 1948. 478p.
• KUSNU, M.; SIAHAAN, M. M.; POELOENGAN, Z. Effects of N, P, K, and Mg fertilizer on the growth and yield of oil palm on typic paleudult. In: SUSTAINABILITY OF OIL PALM PLANTATIONS. AGRONOMIC AND ENVIRONMENTAL PERSPECTIVES, 1996, Kuala Lumpu... 14 p.
• MATOS, A.K.M.G. Biometria e morfologia de Attalea maripa (Aubl) Mart. (INAJÁ) em sistema silvipastoril no Nordeste Paraense. 2010. 80 f. Dissertação (Mestrado em Ciências Florestal) - Universidade Federal Rural da Amazônia, Manaus, 2010.
• MIRANDA, I. P. A.; RABELO, A. Guia de identificação das palmeiras de um fragmento florestal urbano. Manaus: INPA, EDUA, 2008. 228p.
• MOTA, R.V.; FRANÇA, L.F. de. Estudo das características da Ucuuba (Virola surinamensis) e do Inajá (Maximiliana regia) com vistas à produção de biodiesel. Revista Científica da UFPA, Belém, v.6, n.1, 2007. Disponível em: http://www.cultura.ufpa.br/rcientifica/artigos_cientificos/ed_08/pdf/rafael_vitti.pdf Acesso em: 29 set. 2015.
  » http://www.cultura.ufpa.br/rcientifica/artigos_cientificos/ed_08/pdf/rafael_vitti.pdf Acesso em: 29 set. 2015.
• OLLINGER, S. V.; SMITH, M. L.; MARTIN, M. E.; HALLETT, R. A.; GOODALE, C. L.; ABER, J. D. Regional variation in foliar chemistry and N cycling among forests of diverse history and composition. Ecology, Brooklyn, v.83, n. 2, p.339–355, 2002.
• PALLET, D. Perspectivas de valorização dos frutos amazônicos obtidos por extrativismo. Montpellier: Cirad. Outubro de 2002. São Paulo: Folha, 2002.
• SALM, R.; PRATES, A.; SIMÕES, N. R.; FEDER, L. Palm community transitions along a topographic gradient from floodplain to terra firme in the eastern Amazon. Acta amazônica, v. 45, n. 1, p. 65-74, 2015.
• SOH, K. G. Fertilizer use by crop: a global assessment. Beijing: Fertilizer Demand Working Group Meeting, 1997. 6p.
• VIEIRA, S. R.; HATFIELD, T. L.; NIELSEN, D. R.; BIGGAR, J. W. Geostatistical theory and application to variability of some agronomical properties. Hilgardia, Berkeley, v.51, n.3, p.1-75, 1983.
• ZIMBACK, C.R.L. Análise espacial de atributos químicos de solos para fins de mapeamento da fertilidade. 2001. 114 f. Tese (Livre-Docência) – Universidade Estadual de São Paulo, São Paulo, 2001.

Publication Dates

• Publication in this collection
  Aug 2017

History

• Received
  03 Nov 2015
• Accepted
  15 Aug 2016