Hydraulic characterization from porous aquifers of the Brazilian Federal District Caracterização hidráulica dos aquíferos porosos do Distrito Federal

Hydraulic conductivity (K) in unsaturated soil is a key input parameter for modeling subsurface water and solute movements. K-values are also important to better define the potential of aquifers and to optimize water resources management activities. Since K-values are usually not readily available, different techniques are applied to estimate them. This study aimed to estimate unsaturated K-values from porous aquifers found in the Federal District of Brazil. Infiltration tests were conducted in different soil types using the open-end-hole approach and the permeability test using shallow boreholes with specific depths, as reported by Heitfeld in 1979. Soil structure was taken into consideration in such estimations. In order to consider important soil properties such as soil texture and bulk density, K-values were also estimated by means of pedotransfer functions (PTFs). Soil texture was determined in the laboratory and used as input parameter for PTFs. Results from open-end-hole method and permeability test compared to those obtained from pedotransfer functions. K-values from four different shallow porous aquifers systems encountered in the Federal District varied from 10-8 ms-1 to 10-6ms-1. Highest K-values were found in Oxisols while the lowest rates were found in Inceptisols. Decreasing conductivity trend was found with increasing depth due to the increase of loamy soils. Variations in the rate of hydraulic conductivity indicated heterogeneity of porous aquifers due to differences in textural and structural characteristics of the soils.


INTRODUCTION
Groundwater demand has grown significantly in the Federal District of Brazil in the past decades.According to the 2010 demographic census carried out by IBGE (Brazilian Institute of Geography and Statistics), 2.5 million inhabitants were living in the central Federal District by the year 2012, far beyond the originally planned 600,000 people for the year 2000.The study developed by Donald Belcher & Associates in 1950, known as "Belcher Report" (CODEPLAN 1995), pointed out that the water resources in this region were environmentally fragile because of the dominance of high plateaus and the presence of important headwaters, small streams and rivers.
Nowadays, the basic services of water and energy supply, and wastewater treatment in the Federal District are insufficient to meet the demand of a growing population.Water supply in Brasilia (capital of Brazil) and its surrounding satellite cities is mainly dependent on two reservoirs (Santa Maria and Descoberto), some small streams, and groundwater.Major groundwater demand comes from more than 500 urbanizationsbuilt after 1990 (SEDUH 2006).Most of these residential communities are located in high plateaus and thick soils, which play an important role in recharging the aquifers.Urban occupation often seals surface soils, making rainfall infiltration more restricted.Consequently, there is a decrease in the groundwater and surface water availability.Residential communitiesoften establish their own water supply projects, based solely on their needs.
Such situation has motivated researchers to improve the knowledge of hydraulic conditions of local aquifers, which is necessary in order to better control and manage the groundwater utilization by local inhabitants or water suppliers.Campos and Freitas-Silva (1998), Souza (2001), Lousada and Campos (2005), Almeida et al. (2006) and Fiori et al.(2010), among others, classified aquifers based mainly on unsaturated hydraulic conductivity (K), although the results are still incipient because of limited data concerning the potential and exploitation of the analyzed aquifers.
Estimating K-values from the Federal District shallow aquifers is important, because they are key parameters for the differentiation of aquifer bodies and their yield.Hydraulic tests such as infiltration tests support the delineation of different groundwater systems.These shallow aquifers play a crucial rolein recharging deeper aquifers, regulating stream flows, and contaminant flows, among other controls.Until now, there is no long-term monitoring system to create a consistent database necessary for the correct management of these aquifers.K-values from soils and saprolites can be measured by direct field methods or by indirect methods that estimate infiltration potential through other parameters.Most of the studies developed in the region used direct methods such as the double ring and open-end-hole.Field direct methods are time-consuming and expensive.As an alternative, pedotransfer functions (PTFs) are often used to estimate K-values indirectly, based on soil attributes such as texture, morphological structure, organic matter content and density (Bouma & Van Lanen 1987).These methods were developed for soils from temperate regions.Therefore, adequation for tropical soils, mainly Oxisols (the most extensive soil type from tropics), is necessary.
The K-values showed in this study refer to vertical hydraulic conductivity in a specific depth.As the soil can be considered homogeneous in each horizon, it can be treated as an isotropic hydrogeologic medium, so that the vertical and horizontal hydraulic conductivity may be seen as the same order of magnitude.The estimated hydraulic conductivities are related to the unsaturated zone of the aquifers, since the in situ infiltration tests were developed up to 200 cm and the soil samples were collected in the shallow horizons.
Specific PTFs for tropical soils were proposed by Van den Berg et al. (1997), Tomasella et al. (2000Tomasella et al. ( , 2003)), Hodnett and Tomasella (2002), Oliveira et al. (2002) and Benites et al. (2007).They developed PTFs based on multiple linear regression equations and artificial neural networks to determine the parameters specified in the Brooks-Corey (1964) andvan Genuchten (1980) equations, or to estimate parameters related to water retention and hydraulic conductivity.This study aims to obtain unsaturated K-values from Federal District's shallow aquifers in order to contribute to the sustainable management of the water resources of this region.

STUDY AREA
The Federal District of Brazil is located in the Brazilian central plateau and covers an area of 5.789 km².It includes three of the most important Latin American basins (Paraná, São Francisco and Araguaia-Tocantins) (Fig. 1).Streams are composed by small rivers.Aquifer recharge occurs dominantly in areas covered by natural vegetation, in thick soils and gentle slopes.
The potential of shallow aquifers depends among other features on relief, parental rock, climate and biological actions.Campos and Freitas-Silva (1998) classified them as porous aquifers if they occur in soils, weathered rocks and alluvium materials.Soils and saprolites correspond to unconsolidated materials with large horizontal and vertical extensions, ranging from a few meters to up to 80 m in thickness (average: 15 -25 m).They also point out that these systems are isotropic and heterogeneous.
The Brazilian Corporation of Agricultural Research (Embrapa 1978) produced a soil distribution map of the area.Different soil types show distinct conditions of infiltration and recharge because infiltration is controlled by thickness, textural and morphological structures, conductivity, vegetation, humidity, and slope.Textural and morphological structures are the main parameters to determine the dynamic processes in unsaturated zones of this environment.
The most representative soils of the region are Red Oxisols, Red-Yellow Oxisols, and Inceptisols (Embrapa 2006), which cover about 85% of the study area.Alfisols, Mollisols, Histosols and Entisols are also found in the study site.Oxisols present the best conditions for infiltration and recharge, which are frequently associated to savannas and occurring in flat to gentle topography.They occupy around 55% of the total area, show average thickness of 25 m, and present sandy to clayey texture and a highly stable microstructure.Inceptisols have a smaller infiltration capacity.Inceptisols cover about 30% of the area and occur in hilly topography.They are typically thin, show sandy to loamy textures and blocky to granular structures.Grasslands are the most common vegetation cover.Other soil classes occur in small areas and are grouped in: ■ Alfisols and Mollisols; ■ Fluvic Entisols; ■ Histosols; and ■ Sandy Entisols.
The first group shows medium K-values with clayey texture, occurs in gentle reliefs and is strongly structured.The second group is represented by thin soils near streams and in gentle reliefs, and shows often loamy to clayey textures.Histosols have the lowest K-values and are related to wetlands.Sandy Entisols are composed essentially by quartz, derived from the weathering of quartzite.
The mineralogy of the soil layers is mainly controlled by the lithology of parent rocks of the Paranoá, Canastra, Araxá and Bambuí Groups (Campos & Freitas-Silva 1998, Faria 1995).The most important minerals are iron oxides and hydroxides, gibbsite, kaolinite and quartz.The rocks of the Paranoá Group constitute 65% of the Federal District and are widespread in the western and northwestern portions of the territory.The R3 unit (sandy metarrithmite) and the Q3 unit (medium quartzite) are responsible for keeping the high plateaus, with deep weathered mantles, sandy texture and high permeability.The Clay Metarrithmite unit and the Psamo-Pellitic Carbonate unit also produce Oxisols and Inceptisols, depending on the intensity of weathering and parent material.Soils derived from carbonates and pellitic rocks are thick, while quartzites are more resistant to hydrolysis.Similarly, soil texture is controlled by the parent materials: carbonates and pellitic rocks result in clayey texture; quartzites often produce sandy texture.
The Canastra Group occupies about 15% of the area and occurs mainly in valleys of Bartolomeu River.Chlorite phyllites of this Group are commonly covered by a thick weathered soil layer, but it can also occur in hilly regions, associated with thin Inceptisols.The Araxá Group is limited to the southwestern part of the Federal District, mainly along valleys.Main rock types are schist and quartzite lens covered by a thin layer of soil in the majority of the area.The Bambuí Group is distributed throughout in the eastern portion of the Federal District, along the valleys of Preto River as well as in restricted areas in the central and northern portions (valleys of Maranhão River).The main rock types are siltstone, shale and arkoses.Outcrops of this Group can be found along streams and roads.The majority of the area is covered by thick Red Oxisols or by Inceptisols, in the case of a hilly relief.Climate is also important for the Federal District's water resources availability.Besides the relatively high amount of annual rainfall (approx.1,500 mm), there is a marked seasonality over the Federal District.The peak of the rainy summer comprises more than 50% of the rainfall from December to March.Usually there is no rainfall during the peak of the dry season (June to August).The first rainfall usually occurs in late September.This seasonal pattern is also followed by other parameters such as relative humidity (which can go down to approx.5% in the hottest hours of August) and temperature (approx.17°C in colder months; > 30°C in hotter months).Because of this seasonal climate variation, shallow aquifers play a key role in keeping streams perennial.
According to Coimbra (1987), the actual evapotranspiration is around 900 mm per year and 12% of total precipitation infiltrates in the soil, reaching the saturated zone of the aquifers.The presence of different classes of savanna vegetation is associated with different geomorphological compartments and includes grasslands, shrublands and forestlands.Shrublands are commonly developed in thick Oxisols, grasslands are widespread in Inceptisols and forestlands are found mainly along streams.Forestlands are important to facilitate infiltration of rainfall and for aquifer recharge due to bioturbation by roots and protection against soil erosion due to heavy rainfall.
Based on the types of weathering mantles, thickness of saturated soil layer and hydraulic conductivity, Campos and Freitas-Silva (1998) classified the Federal District's shallow aquifers into four systems: P1, P2, P3 and P4.Gonçalves et al. (2007) proposed a different approach for mapping these aquifers, which was based on new infiltration data and on curve-number theory (SCS 1954).P1, P2 and P3 systems are characterized by high thickness and high, moderate and low hydraulic conductivities, respectively.The P4 system presents low thickness and low hydraulic conductivity.Shallow wells operated in these aquifers show low average yields (< 800 L.h -1 ).Theses aquifers feature a local groundwater flow regime and the discharge is associated with contact and slope of springs.
These shallow aquifers are important not only to act as complementary water supply of relatively isolated rural communities, but also to keep streams, rivers, springs and lakes perennial during the dry season.The unsaturated zone of these aquifers functions as an important protection factor to avoid anthropogenic contamination.As the saturated zone stays incontact with fractured aquifers, the porous aquifers, mainly the P1 and P2 systems, play an important role in their recharge.

METHODS
As mentioned before, hydraulic conductivity is a key parameter to characterize porous aquifers.Although many authors have studied it (e.g., Campos & Freitas-Silva 1998, Souza 2001, Lousada & Campos 2005, Almeida et al. 2006, Fiori et al. 2010), none have considered the influence of soil structure in the porous aquifers.According to Williams et al. (1992) and Danalatos et al. (1994), soil structure in tropical countries is essential to understanding the high K-values in soil profiles with high clay contents.
To estimate K-values from porous domains of the Federal District's soils, soil structure was taken into consideration.Infiltration tests were conducted in different soil classes using the open-end-hole tests as well as the permeability tests proposed by Heitfeld (1979).Soil texture was determined in the laboratoryon soil samples.K-values were also estimated by applying pedotransfer functions available in the Rosetta software package (Schaap et al. 2001).

Open-end-hole method
Determinations of hydraulic conductivity by in situ soil permeability (infiltration) tests were used to compare field test data with the results of pedotransfer functions.Shallow boreholes were drilledto conduct unsteady-state permeability tests in unsaturated zones of soils.The evaluation of permeability tests or infiltration tests is based on Darcy's equation.Field results were evaluated using equations as proposed by Heitfeld (1979) and Scheytt and Hengelhaupt (2001).Soil samples were taken from each borehole to determine chemical composition and particle size distribution.

Location And Soil Types
Eight sites located in the surroundings of the urban areas of the Federal District, representing the main soil types found in the study area.Boreholes were drilled to obtain undisturbed, 1-meter depth soil samples (Fig. 3).Additional boreholes near sampling points were drilled in six sites at 30 -50 cm depth and installed for in situ permeability tests in the unsaturated zone.
Due to the thickness of Inceptisols and related saprolithes, which are normally less than 1 m thick, infiltration tests could not be realized.Very high bulk density Alfisols are difficult to befound in the Federal District.Field investigations led to the assumption that the permeability of Alfisols depends more on structural elements like faults and cracks associated with neotectonic events than on soil porosity.Textural analyses were conducted for three depths of each borehole: 0 -20 cm; 20 -50 cm; and 50 -100 cm.The particle size distribution cannot be used for permeability estimations with conventional methods because tropical soils have high clay and silt contents.Common methods of permeability estimations of fine sand to gravel textured materials, as described by Beyer (1964), are therefore not applicable.The particle size distribution was used in this study as input parameter for pedotransfer functions.The rate of infiltration through a soil tends to the partially saturated permeability of the soil.The soil does not reach full saturation because of the presence of a small volume of air in the porosity, for instance, associated with injected water.This may reduce the K-values measured in the field.For soils with low permeability (e.g., loamy and clayey soils), the effect of suction at the wetting front can also affect the test results, especially in unsaturated zones.Therefore, saturation of soils prior to the test itself is recommended (ISO 2012).The test was conducted with a limited saturation phase of 10 -20 minutes, which surely did not avoid the effect of suction.

Permeability tests
For the evaluation of field test with variable pressure height, the following equation for the open-end-hole test was used: where r = radius of the borehole; h 1 , h 2 = height of water column; and t = time.Factor C depends on the screen length and borehole diameter.The equation for the open end test with nominal screen length (L) = 0 can be calculated by: The infiltration test after Heitfeld (1979) and Earth Manual (1951, after Scheytt andHengelhaupt, 2001) for unsaturated soils can be evaluated by the following equation:
Comparing the open-end-hole test with the setup of permeability test after Heitfeld (1979), additional space at the bottom of the cased borehole is used.The space is filled with gravel and works as drainage space instead of an installed borehole with casing and screens.

Pedotransfer Functions (Rosetta Software)
PTFs are recognized as an efficient method to predict soil hydraulic properties such as water retention and hydraulic conductivity.Reviews regarding the development and the use of PTFs can be found in Rawls et al. (1991), Wösten (1997), Pachepsky and Rawls (1999), and Wösten et al. (2001).There are few references on results of research with PTFs that had used database from tropical countries.Cresswell et al. (1999) reviewed the application of PTFs to predict hydraulic properties of soils from Australia.Also in Australia, Minasny and McBratney (2002) applied the neural network method to predict the parameters of water retention of the van Genuchten (1980) equation.McKenzie and Jacquier (1997) presented a regression tree that predicts saturated hydraulic conductivity from soil morphology classes, including texture, structure and porosity.Minasny et al. (1999) predicted the parameters of the van Genuchten equation from basic soil properties and broad classes of soil texture.
In Brazil, estimations of hydraulic conductivity of kaolinitic Yellow Oxisols in the Amazonian region, based on water retention curves, were obtained by Tomasella and Hoodnett (1997).Tomasella et al. (2000) used soil texture, organic carbon, equivalent humidity and soil density data to establish relations between the van Genuchten equation and water retention parameters.Gaiser et al. (2000) studied the influence of clay minerals in soil water retention from semiarid areas of western Africa and northwestern Brazil and showed that PTFs of soils with low activity clay minerals differ from soils with high activity clay minerals due to silt and organic carbon content.Hoodnett and Tomasella (2002) developed new PTFs for tropical soils related to water retention parameters used in the van Genuchten equation.
The van Genuchten equation was successfully applied to several different soils [e.g.soils reported by Villagra et al. (1988) and Wösten and van Genuchten (1988)].Van den Berg et al. (1997) found that this equation describes adequately the water retention curves of soils with low activity clay minerals, commonly found in Brazilian soils.The parameters of the van Genuchten equation adapted to Brazilian soils showed high accuracy in the estimation of hydraulic conductivity, as observed by Hoodnett and Tomasella (2002).
Most of these PTFs were compiled by Schaap et al. (2001) and integrated into the Rosetta software package, which estimates the parameters of water retention and the hydraulic conductivity according to models proposed by van Genuchten (1980) and Mualem (1976).Rosetta offers five PTFs that predict hydraulic properties of the soil from limited sets of input data.Hierarchic models use the following input data: soil textural classes; percentages of sand, silt and clay, density and a point of water retention in 330 cm (33 kPa) and the same parameters and water retention points in 330 and 15,000 cm (33 and 1500 kPa).The first PTF is based on average hydraulic parameters for varying types of soil texture.Other four models are based on analyses of neural network and usually provide more accurate forecasts, depending on the quality of the input data.Besides the hierarchical approach, the model provides prediction of hydraulic conductivity based on van Genuchten (1980) and Mualem (1976).
In this study, 25 soil profiles from the Brazilian Soil Survey (Soil Map of the Federal District -Embrapa 1978) were used to estimate unsaturated K-values by applying PTFs available in the Rosetta software.Six types of more representative soils found in the study area (Red Oxisols, Red-Yellow Oxisols, Alfisols, Inceptisols and Histosols) were described according to the textural and structural variations of each diagnostic horizon.In order to estimate the hydraulic conductivity more appropriately and to adapt the PTFs for soils of the Federal District, values of soil density (Dd) and gravimetric water contents (W)were obtained.PTFs developed by Benites et al. (2007), shown in equation ( 4) and Tomasella et al. (2003), shown in equation ( 5 The PTF used to estimate the water retention parameters from van Genuchten's equations in the range between -33 kPa (W33) and -1500 kPa (W1500) is as follow (Tomasella et al. 2003 (5) where CS = coarse sand (%); FS = fine sand (%); SI = silt (%); CL = clay (%); Db = bulk density (g.cm -3 ); Me = moisture equivalent (%); W 10kPa , W 33kPa , W 1500kPa gravimetric water content at 10, 33, and 1500 kPa, respectively (g.g -1 ); X 14 -X 17 and Z 9 -Z 11 and Z 13 , auxiliary variables.
Reference PTFs were chosen according to the following criteria: ■ involving only textural properties; ■ being developed from a vast amount of samples of Brazilian soils; and ■ having been extensively evaluated in terms of its applicability.

RESULTS AND DISCUSSIONS
To reach the objectives of this study, the main classes of soils from the Federal District were characterized in terms of textural and structural features observed in field and described by Embrapa (1978) (Tab. 1, Fig. 6).The diagnostic horizons (approx.30 cm depth) were considered as the basis for each profile classification.
The results demonstrated a significant increase of clayey texture at a depth of between 50 and 100 cm in Oxisols.Profiles showed that more than 50% of Oxisols have clayey texture.Such characteristics reveal strong development of finer granular structure in Oxisols, which reflects the high macroporosity of these soils.In average, these soils present 30% of sand, 10% of silt and 50% of clay.Alfisols also indicated the increase of clayey texture but in a lesser extent in relation to Oxisols.These soils present about 30% of sand, silt and clay for all horizons and are also structured in blocks, with strong aggregation.The clayey texture is more apparent in the first 20 cm of Histosols and Inceptisols.These soils have slight decrease of clay content between 20 and 50 cm.These soils did not occur at depths higher than 50 cm.The structures of these soils are prismatic, in blocks and, in some cases, massive.

Open-End-Hole and Permeability Test after Heitfeld (1979)
Results for open-end-hole tests and permeability tests after Heitfeld (1979) showed good correlation.The highest difference between the two methods was found in the site ID 6 (14%).Results for P1, P2, and P4 hydrogeological systems allowed estimating permeability values for different soils types, including Histosols and Alfisols.The hydraulic conductivity of these soils ranged from 10 -8 to 10 -6 m.s -1 (Tab.2).
Because of the amount of field tests, comparison between the results of this study and those obtained by previous investigations in the Federal District and in Goiás State was necessary.Souza and Campos (2001), Reatto et al. (2007) andFiori et al. (2010) reported K-values of soils based on in situ infiltration tests and laboratory tests.Souza and Campos (2001) directly referred the results to the hydrogeologic units P1 to P4. Mean values of these studies of comparable depths and soils were taken to visualize differences in results obtained by these studies and current study (Fig. 7).There was a good approximation with the results obtained by Reatto et al. (2007).K-values obtained by Souza and Campos (2001) were around one decimal power higher.K-values from Fiori et al. (2010) matched at good to very good approximation when compared with the soils characteristics observed in situ.

K-values Obtained by the Rosetta Method
K-values from Embrapa's data were obtained after adapting PTFs for the Federal District's soils by using appropriated values of soil bulk density (Db) and gravimetric water content in the Rosetta software (Tab.3).Db of Oxisols ranged from 1.44 to 1.92 g.cm -3 .Db from Alfisols ranged from 1.62 to 1.79 g.cm -3 .Finally, Db from Inceptisols and Histosols ranged from 1.60 to 1.65 g.cm -3 .There was an increase of Db values with increasing depth in most of soil profiles which can be explained by the increase of clay content and consequent decrease of total porosity and hydraulic conductivity.Gravimetric water contents at -33 kPa (W33) and -1500 kPa (W1500) from Oxisols ranged from 0.180 to 0.993 g g -1 and 0.048 to 0.834 g.g -1 , respectively.For Alfisols, gravimetric water contents ranged from 0.049 to 0.514 g.g -1 and from 0.004 to 0.196 g.g -1 for the W33 and W1500, respectively.Gravimetric water contents from Inceptisols and Histosols ranged from 0.272 to 0.466 g.g -1 and from 0.079 to 0.161 g.g -1 , respectively.In general, these results showed that the percentage of variance gradually decreased from 73 to 26% in Oxisols, from 75 to 25% in Alfisols, and from 75 to 24% in Inceptisols and Histosols.It is related to the decrease in the proportion of porosity that retains water when the water potential decreases.The same evidence was noticed for varying soil types considered by Bruand et al. (1988), Dexter (2004), andReatto et al. (2007).
K-values calculated for Oxisols ranged from 3.44 * 10 -6 to 1.71 * 10 -7 m.s -1 in surface zone and from 2.08 * 10 -6 m.s -1 to 2.53 * 10 -7 m.s -1 in deeper zones.For Alfisols, they ranged from 1.45 * 10 -6 to 5.73 * 10 -8 m.s -1 in surface zone and 2.88 * 10 -5 to 4.25 * 10 -8 m.s -1 in deep zone.Inceptisols and Histosols ranged from 1.28 * 10 -6 to 8.96 * 10 -8 m.s -1 in surface zone and from 6.33 * 10 -7 to 9.15 * 10 -8 m.s -1 in deeper zones (Tab.4).Again, these results showed a decreasing trend of K-values with depth due to the decrease of porosity and the increase of clay contents in Oxisols.This was also observed in estimations based on tests after Heitfeld.Furthermore, these soils commonly occur in high plateaus (chapadas).When found in natural conditions, they impose a great control of conditions of infiltration, regardless of the parent material.High correlations were also observed in aquifers belonging to P2 system.Hydraulic conductivity from Alfisols showed a similar range of values which was observed in tests after Heitfeld (1979) and by the Rosetta software.The average value estimated by Rosetta was 3.57 * 10 -6 m.s -1 , while the value observed on the tests after Heitfeld was 1.87 * 10 -6 m.s -1 .The behavior of this soil class was depth-dependent.There were both increase and decrease of hydraulic conductivity in depth, but Alfisols showed the same behavior.It followed the expected pattern: clay content increased with depth and this led to a decrease in hydraulic conductivity.
As noticed by previous studies (Souza & Campos 2001, Cadamuro et al. 2002, Fiori et al. 2010) and also by this study, Inceptisols presented the least favorableconditions of infiltration and drainage.These soils are poorly developed, derived from low-degree metasedimentary rocks, and occur in hilly regions (intense runoff and limited infiltration).Based on infiltration tests and on comparison and addition of data derived by previous studies, Tab. 4 shows K-values for different aquifer systems.These are considered the typical values for the different shallow aquifers and the best fit between the parameters and aquifer types.

CONCLUSIONS
Good correlation between observed and estimated data was found when field permeability tests and PTF methods were used.Values of the P1, P2, P3 and P4 hydrogeological systems varied from 10 -8 ms -1 to 10 -6 ms -1 .Highest K-values were found in Oxisols, while lowest rates were found in Inceptisols.A decreasing conductivity trend was noticed as soil depth increased, mainly due to the presence of materials with loamy texture.
Contrary to the previous studies, results obtained from the Rosetta software indicated that porous aquifers of the Federal District are anisotropic, once the K-values changing in space (vertically and also in different points for the same soil type), i.e. vertical and horizontal conductivity are commonly different.This evidence is related to variations of K-values within a single class of soil.Heterogeneity of porous aquifers is conditioned by the texture, structure and macroporosity of soils.On the other hand, anisotropic behavior of soil layers can influence the behavior of porous aquifers significantly.Depending on parent material, structural elements (for instance, faults and fractures), and geomorphological characteristics, different vertical water flows during the recharge can induce variable rates of inward flow.
The thickness of porous aquifers in the Federal District ranged from few centimeters to 80 m, with predominance (> 60%) of thickness between 15 -25 m.P1 and P2 systems presented thickness higher than 20 m.The thickness of P3 system was lower than 10 m and, for the P4 system, lower than 1 m.In this system, the lack of saturation zone is common in saprolites, especially when developed on clayey rocks.Inceptisols usually have low K-values because their unsaturated zone lays directly on the rocky substratum.

Figure 1 .
Figure 1.Location map of the study area (Federal District of Brazil).
Geomorphology controls local and regional hydrogeologic flow systems and influences the recharge and discharge and overall flow patterns.Five different compartments were recognized by Novaes Pinto (1994): High Plateau (elevation > 1,100 m; drainage density: low; topography: flat to gentle; major soil type: Oxisols); Intermediate Plateau (elevation: 850 -1,100 m; topography: flat to gentle; major soil type: Oxisols); Incased Valley (topography: flat to gentle; drainage density: high; major soil type: Inceptisols); Soft Border (transition from High Plateau to Intermediate Plateau); and Strong Border (abrupt transition from High Plateau to Incased Valley) (Fig. 2).

Figure 2 .
Figure 2. Geology groups, soils and elevation of the Federal District.
results of the infiltration tests represent the hydraulic characteristics of the sediments surrounding the boreholes.In other words, K-values for permeability are essentially point-based, in which accuracy is shown only for small areas surrounding the boreholes.Three types of permeability tests in open boreholes can be conducted: ■ constant flow rate test with change in hydraulic head in the borehole; ■ constant head test with stabilized infiltration rate; and ■ variable head test with instant change in the hydraulic test due to quick injection.The change of hydraulic head time is measured.Stable and variable head tests were chosen as the preferable infiltration test methods.The water level above the infiltration point corresponds to the filter screen or to the open end of boreholes and varies over time.

Figure 3 .Figure 5 .
Figure 3. Porous aquifer systems and the locations of soil samples in the study area.

Table 1 .
Textural and structural characterization of major soil types from Federal District.

Table 4 .
Porous aquifer systems, soil types and K-values of Federal District soils.