1. Introduction, previous research, and study purpose
The alluvial fan system located beneath and around Ciudad Juárez, Chihuahua, northern México (Figure 1) developed during the late Tertiary to Quaternary throughout the southwestern USA and northern México region. Much research has been published (in the USA) about the history of the Río Grande, with key papers published by Hawley et al. (1976) and Mack et al.(2006). However, there is little knowledge of the southern continuation of the system into northern México where the Río Grande is called the Río Bravo. This paper presents the first attempt to link the chronological processes of the alluvial system in the Ciudad Juárez area to USA sites, using Optically Stimulated Luminescence (OSL) dating, and aims to advance understanding the history of the alluvial fan depositional system of northern México.
There is a globally recognized need to understand the behavior of regional scale alluvial fan systems due to flood risk and sediment transport associated with climate change. Furthermore, understanding the regional effects of climate on alluvial fan deposits is a key research frontier (Owen et al., 2014). However, little previous research has been conducted on Mexican regional alluvial fans. Sparsely inhabited areas of the El Fresnal Playa within the Los Muertos Basin, about 100 km west-southwest of our study area in Chihuahua, have well-exposed sections that allowed Ortega-Ramírez et al. (2004) to identify three generations of apparent late Quaternary alluvial fans in the hanging wall and footwall of a major fault. In contrast, the high population density has limited the geologic outcrops in Ciudad Juárez. Although small road cuts are abundant across the fans and expose undisturbed deposit outcrops, these outcrops cannot be easily correlated to build up a clear stratigraphic context for the fan history without geochronological age control. Studies using modern geochronological tools, such as those presented here, are critical for reaching a better understanding of alluvial fan dynamics (Owen et al. , 2014).
One such geochronological tool is OSL dating (Aitken, 1998; Rhodes, 2011). OSL dating involves measuring electrons that have been displaced by ionizing radiation and in the process became trapped in crystal lattice defects within a mineral (Huntley et al., 1985). The concentration of these electrons increases with time while buried and decreases while in sunlight (Gray and Mahan, 2015). This property has allowed grains of quartz and feldspar to be used as geo-chronometers for sedimentary units in a variety of geomorphic environments (Huntley and Johnson, 1976; Wintle and Huntley, 1979; Huntley et al. , 1985; Rhodes, 2011). Quartz OSL is particularity useful in that the populations of electrons can be easily removed by sunlight such that there is a high probability of obtaining an accurate age. This is advantageous for alluvial fans where deposition can occur during flash floods, which have limited light exposure for sedimentary grains. Additionally, robust statistical techniques such as the Minimum-Age-Model (MAM) can be used to further control for the uneven light exposure of grains (Galbraith et al., 1999).
In an effort to establish these stratigraphic correlations and obtain age control, we present three new OSL dates on sediment in the alluvial fan system. We use these data to constrain the ages of formation of the fans and discuss the implications of these dates in relation to regional fan processes and late Quaternary climate. This study will also help improve understanding of the recent geological history of the region in order to place the current problems of flooding in Ciudad Juárez into context. This is part of long-range planning for mitigation of the frequent floods affecting Ciudad Juárez (Esparza et al., 2007).
2. Río Bravo and fans: features and geological setting
The pre-Neogene drainage network of the study area was derived during the compressive Laramide orogeny during Paleocene to Eocene time. After that, deformation of Cretaceous rocks created a synclinorium structure and four thrust sheets (oriented NW) were emplaced in the area and formed the rock massif now represented as the Juárez Mountains. Then, during Eocene to Oligocene time, intrusive and extrusive andesite rocks were deposited at the surface. Therefore, many reverse, normal and oblique faults typical of the Basin and Range physiographic province were created at this time, in probable response to magmatic events and these faults extend across Chihuahua State, Ciudad Juárez, El Paso (Texas) and New México (García, 1970; Wacker, 1972; Nodeland, 1977; Drewes and Dryer, 1993; Carciumaru and Ortega, 2008). During this time, the ancestral Chihuahua trough was destroyed and its deposits were transported, filling the northern Hueco and Mesilla bolsons, which form typical Basin and Range intermontane basins. Upon this antecedent surface, alluvial fans have been developed, generating the topography upon which Ciudad Juárez is constructed.
The geomorphic evolution is spatially and temporally associated with both, tectonic and climatic controls. Down-faulted areas produced basins, which accumulated the sedimentary Lower, Middle and Upper Santa Fe Formations in the Hueco and Mesilla bolsons (Hawley and Kottlowsky, 1969; Hawley et al. , 1976). Furthermore, recent tectonic and/or climatic events associated with the arrival of major ephemeral distributaries of the Río Bravo during Pleistocene to Holocene time likely produced river terraces. Finally, ages of ash derived from Lava Creek B (from the Yellowstone volcanic system), preserved in the highest layer 75 m above the Río Bravo floodplain, in the area of the University of Texas at El Paso (near Mesa Street), demonstrate incision of the Río Bravo channel around 660 ka before present along the Paso del Norte between the Franklin Mountains and the Sierra de Juárez uplift (Hawley and Kottlowsky, 1969; Hawley et al. , 1976).
3. Río Bravo system evolution during early Pleistocene to Holocene time:
The Río Bravo system and surrounding mountain streams were developed during Pliocene-Pleistocene to Holocene time (Mack and Seager, 1990, 1995; Mack et al. , 2006) leaving a varied regional landscape which can be broadly divided into three major areas (Figure 1): A) The northwest sector composed of the Anapra alluvial fan, where the fan processes are active transporting sediment via ephemeral streams. There is little human modification in this area, with only one major protective structure (the Benito Juárez dyke) to prevent flooding in the lower reaches of the fan. Otherwise, there is no interruption of natural processes. Sector B) the central sector comprised of: West Snake Stream; East Snake Stream (VAF on Figure 1); the alluvial fan called Colorado Alluvial fan (CAF) and the piedmont mountain front fan system (PCAF; Figure 1); finally, sector C) The southeast sector, which integrates the Palo Chino (PCH) and Jarudo fan system (PCHJAF; Figure 1). Sectors B) and C) are areas of major construction of buildings and roads.
In the streams, numerous minor dykes trap downstream movement of floodwater and alluvial sediment; and many streams are lined and shaped with concrete and act as both roadways and flood spillways. Therefore in these sectors, the natural fan processes are severely interrupted, and fan action does not operate in its normal mode. However, despite recent human interruptions of fan processes, there is evidence, produced by new dates, that the fan system was also interrupted naturally by dissection and lateral erosion in the Late Pleistocene. The lower parts of the fans were subsequently covered by Río Bravo sediments, which were later eroded to form terraces, broadly correlating with documented Pleistocene to Holocene climatic processes of New México and El Paso, Texas (Mack and Leeder, 1998). It is worth noting that Mack and Leeder (1998) commented that the Río Bravo is a special case, because channel migration and avulsion were frequent, and these processes might impact the degree of intercalation between the river and alluvial terrace sediments.
Figure 1 shows a Digital Elevation Model (DEM) of the study area with the distribution of alluvial fans derived from the Juárez Mountains. OSL ages were obtained from the CAF (Figure 2). At least five shifts of the Río Bravo tributaries occurred during Pleistocene to Holocene time, providing these possible scenarios: i) the Cabeza de Vaca Lake of Plio-Pleistocene time (> 2.5 Ma) (Strain, 1966; Hawley and Kottlowsky, 1969; Vanderhill, 1986; Gustavson, 1991); ii) the ancestral Río Bravo path changed again during the Middle to Late Pleistocene (0.66 Ma to 0.16 Ma); iii) the older Río Bravo tributaries built a dissected terrace during Late Pleistocene time (0.16 to 0.10 Ma) (Connell et al., 1998); iv) another terrace system was formed during Late Pleistocene time from 78 to 28 ka (Connell et al. , 1998) and again from 15 to 22 ka (Hawley and Kottlowsky 1969; Hawley et al. , 1976; Allen and Anderson, 1993). During this period of time, younger Río Bravo tributaries were active again. Finally, v) from 15 ka to the present, the last shift of the Río Bravo has been in progress. Connell and Love (2001) described a series of late Quaternary sands and gravels deposited above the Santa Fe Group (making them Post-Santa Fe Group sediments). These form a series of river terraces adjacent to the modern Río Grande in New México. Of interest to our study, the upper terraces (informal depositional units Qra and Qay of Connell and Love (2001)) are approximately age-equivalent to the dates that we have obtained from the alluvial fans above the Ciudad Juárez (see Figures 3 and 6 of Connell and Love, 2001) noting that these authors did not provide numerical ages).
4. Methods
Fan distribution was studied and mapped using a combination of field observations and GIS-based topographic mapping, with the caveat that comprehensive fieldwork in the Ciudad Juárez area is precluded by the large-scale ground coverage of buildings and roads. Mapping was used to determine the most appropriate locations for OSL sediment sampling (Personius and Mahan, 2000, 2003; Nelson et al., 2015). OSL was applied to samples collected from three specific sites mostly composed of fine aeolian quartz sands (see Figures 3A-C). The samples were collected by the first author using the procedure recommended in Cole et al. (2007). Care was taken to keep the samples light-tight during collection and shipment to the OSL laboratory. The samples were processed using protocols set forth in Gray et al. (2015). We used the 250 - 180 µm fraction because this one yielded the largest amount of quartz particles useable of the sample. However, even using this largest fraction of quartz, we were limited in the number of aliquots available for analysis due to small quartz yields in the lab. Samples were treated to standard HCl, H2O2, and HF treatments, magnetic and density separation, and dose rates were measured using gamma spectrometry at the U.S. Geological Survey TRIGA facility (Gray et al. , 2015) cosmic ray contribution to dose rate for luminescence was achieved using the method of Prescott and Hutton (1994).
As noted above, a key problem in fan analysis is that the system has been dissected and extensively covered by the modern development of Ciudad Juárez. Determining detailed geometric relationships between the various parts of the fans is effectively impossible. However, it is expected that the uppermost exposed deposits at any site will provide the best opportunity to assemble age data for more recent fan evolution. Consistency of the material was achieved by sampling of aeolian sands at all the three sites. In each outcrop location, the local stratigraphy shows that the sands sampled for OSL dating are overlain by conglomerates (see Figures 3A-C) so episodes of aeolian deposition were followed by renewed fan building activity. The three sample sites help to constrain the overall general history of the alluvial fans but are specifically tied to the CAF, which was selected in order to understand the most recent history of the fan system. Therefore, aeolian layers should be the best candidates for a complete bleaching during the depositional processes and would function as proxies for arid or windy period. Samples were collected at the following sites on the CAF (see Figure 2):
P1 collected at 13R: 354,694.63 m E; UTM 3,512,048.16 m W; elevation 1252 masl.
P3 collected at 13R: 355,593.59 m E; UTM 3,512,861.18 m W; elevation 1201 masl.
P2 collected at 13R: 355,730.21 m E; UTM 3,512,678.11 m W; elevation 1241 masl.
5. Results
The OSL data (Table 1) and results (Figure 4) show the ages distribution for the three sites on the CAF. The most notable result from the OSL analyses is that, there are two sediment populations in the two younger samples. Furthermore, the oldest sediment populations can be observed in each of the succeeding samples. Sample P1 (S1) shows a population at 31 ka and another population at 45 ka P1 (S2) while sample P2 shows a population at 41 ka (S1) and a second, older population at 69 ka (S2). This 69 ka population overlaps with the population measured from the oldest sample P3 (S1)) at 74 ka (see Table 1).
aField moisture, with figures indicating the complete sample saturation %. Ages calculated using approximately 15 % of total measured saturation.
bAnalyses obtained using laboratory gamma spectrometry (high-resolution Ge detector).
cIncludes cosmic dose rate of 0.23 Gy/ka for P1, 0.13 Gy/ka for P2, and 0.15 Gy/ka for P3 (using methods of Prescott and Hutton, 1994; depth and attentuation)
dNumber of replicated equivalent dose (De) estimates used to calculate the age. Figures in parentheses indicate total number of measurements included in calculating the represented De and age using the central age model (CAM; bottom age) and the minimum age model (MAM; top age). Preferred ages in bold font.
eDefined as "overdispersion" of the DE values. Obtained by taking standard deviation over the average. Values > 30 % reflect poorly bleached or mixed sediments.
fAge for 250 - 180 microns quartz sand. Exponential fit used on equivalent dose, errors to one sigma.
The results coincide with the period of time covering Oxygen Isotope Stages (OIS) 3 to 5 as follows: samples P1 was measured at 30 and 45 ka (late OIS Stage 3), samples P2 was measured at 41 and 69 ka (around the boundary between OIS 4 and 3), and sample P3 was measured at 74 ka (very latest OIS 5a). These ages represent deposition during both glacial and interglacial phases.
The sample with the tightest equivalent dose distribution was P3 (Table 1). However, P3 is also the oldest sample and came very close to having all quartz OSL trap sites completely filled (usually quartz saturates around 150 - 200 Grays in traditional OSL measurements; Mahan et al., 2015). The two younger samples were not saturated for two reasons: firstly, because the growth curves (see Figures 4A, 4B, 5A and 5B and table 5 of supplementary material) did not exhibit saturation behavior and secondly because the dose rate was moderate enough (2.80 to 3.26 Gy/ka) to allow natural accumulation to occur for some length of time (Table 1). The quartz in sample P3 was approaching to saturation at 208 Grays, so this age could be even older, but is certainly not younger than 74 ka (Figure 6B of supplementary material). P1 and P2 have two equivalent dose populations with correspondingly higher dispersions, close to or over 30 % to 40 % (P1 = 29 % and P2 = 39 %; (see Table 1).
It is not entirely unexpected to have several populations of grains within poorly sorted alluvial fans sediment (Mahan et al., 2007, 2015; Miller et al., 2010) due to loose sand and gravel packing, fluvial reworking; partial bleaching, due to rapid depositional processes, or reactivation of unconsolidated or uncemented sediment layers. However, the aeolian sediment within the alluvial fan material is most appropriate for OSL dating because the nature of these deposits ensures that maximum exposure to sunlight before deposition occurred (Bateman et al., 2012) and because targeting finer-grained layers within the fan also increases the chances that a majority of the sediments were zeroed or bleached before final deposition.
6. Discussion
Relict topography has influenced the alluvial and fluvial history of the study area, producing a complex sedimentary system comprised of river sediments from the Río Bravo in the northeast, and the alluvial fan material derived from the Juárez Mountains in the southwest. The timing and extent of the inter-fingering of the Río Bravo terrace sediments with the Juárez Mountains alluvial material is currently unknown and represents an area for future exploration; we will concentrate on evidence provided by the OSL ages and what they indicate about alluvial fan incision and aggradation. Corroborating evidence for substantial, long-term erosion of the fans, using the OSL dates presented in this paper, is shown in Figure 6, where it can be seen that the fans are incised to a substantial degree, so that the city is built on an eroded fan surface. Furthermore, it is possible that increased Holocene aridity, developed after and throughout the last glacial phase, may have led to substantial degradation of the fans.
6.1. Ciudad Juárez alluvial fan age
The younger of the two ages of P1 and P2 are more reliable than the older ages of these two samples, and the younger dates are obtained using a (MAM) (Galbraith et al., 1999). The MAM thus marks three episodes of aeolian deposition (at approximately 30 ka, 40 ka and 75 ka) contained within the uppermost (youngest) fan structure. At least three explanations for the deposition and the preservation of the sands are possible: A) sediment became available early in the two last glacial maxima, when there was more erosion of the highlands due to increased moisture, so that when arid periods impacted the basin there was a ready supply of sand to be deposited on the fans (Ellwein et al., 2011, 2015); B) dust or sediment became available because of dune activation or reactivation of the easily deflated playa, river terrace, and lower alluvial fan sediments (Muhs, 2013). It is possible that these periods contain episodes of low effective moisture and high dust accumulation due to drying winds as the glacial phase develops, in which case the formation of aeolian sands for P1 and P2 match global climate changes. However, this is a broad interpretation, and it is more likely that there were C) simply more localized controls operating in this region, such as vegetation and wind shifts (Cole et al., 2007; Hall et al., 2008).
Vegetative cover change was highlighted by Harvey et al. (1999) where attention was drawn to the differences in erosion characteristics of late Quaternary fans in California, which are influenced by the degree to which they are covered by vegetation. In the Ciudad Juárez area, throughout the fan system sparse vegetation is obvious, and it has proved difficult to find outcrops or cores containing pollen, seeds, or other physical evidence that the region was previously well vegetated. Ceballos et al. (2010, Figure 6) provided maps of Pleistocene vegetation throughout México, and identified xeric scrub as the dominant vegetation in the Ciudad Juárez study area, indicating poorly developed land cover only minimally sufficient to support a land mammal migration corridor between northern central México and southern central USA. Thus, erosion may be sustained in the highlands of the Juárez Mountains leading to aggrading fan development in the late Quaternary, with subsequent erosion in more recent times in the lower Juárez Mountains. Furthermore, there are studies of global teleconnections that suggest monsoon rains during the late Quaternary glaciation and possibly earlier (Coats et al., 2008; Miller et al., 2010; Muhs, 2013) along the Pacific coasts of North America, which would be compatible with shifts between aeolian and conglomeratic deposits in the Ciudad Juárez fans.
Unfortunately, because of the urban sprawl of buildings and roads, the history of the fans cannot be fully determined, but one key point emerges: the OSL dates reveal that the sand has older ages than expected, and suggest that large parts of the fan system became inactive, with the youngest preserved activity being around 30 ka. The amount of erosion on the fans cannot be totally known, and so it remains possible that the fans could have aggraded in the Holocene, with subsequent erosion. To avoid the problem of younger aggradation and subsequent removal, however, the samples were collected from the highest preserved points where it is not likely to have been easily eroded (Figures 1, 2 and 6, for visual inspection of the geomorphic relationship between the fans and the surrounding hills).
In more recent times, the alluvial system was partly reactivated in the Anapra fan. The Anapra fan is active in the north part of the study area and this observation suggests shifting of fan aggradation, while abandoning the older parts (where the city is built). It is reasonable to presume, from the geometry and activity of the Anapra fan (Figure1) and the continued flooding of the city after heavy rainfall events that the hydrological system has been intermittently to continually active through the Pleistocene to modern time, but has shifted position, abandoning the older fan system by starving it of sediment. The older (OSL dated) fans were likely degraded by more recent erosive activity but this cannot be satisfactorily investigated because of the paucity of outcrops in the city. Although making overall interpretations from only three dates is tenuous, these dates do provide valuable preliminary constraints on the history of the area.
6.2. Comparison with other alluvial fans and climate
Comparison between the three OSL ages and the pattern of the global oxygen isotope curve (Figure 5) shows that there is not a clear relationship between glacial/interglacial cycles and broad-scale regional sedimentary process in the northern México area. The proximity of the fans to the Juárez Mountains, and the repeated local stratigraphy in the three sample sites, suggests instead that localized climatic changes, perhaps influenced by the steep topography of the mountains and the prevailing winds (Muhs, 2013), were a major control on fan formation. It is possible that cool-season Pacific frontal storms caused river flow, ephemeral lakes, and fan incision, whereas periods of intense warm-season storms caused hillslope erosion and alluvial fan aggradation, as suggested by Miller et al. (2010) and Muhs (2013). Alternatively, the regional topography from the Juárez Mountains, across the Hueco Basin to the Franklin Mountains to the north, created a geographically constrained climatic environment that was locally more influential than the global glacial/interglacial cycles. Support for this suggestion comes from recent work on the Río Bravo terraces in the Albuquerque area of New México (Cole et al., 2007), which shows that there is no clear relationship between terrace age and climate in the Albuquerque area: dates on terraces from a variety of dating methods show scatter throughout OIS-6 to OIS-1, with no firm pattern. Furthermore, compilations of alluvial fan cosmogenic geochronology ages demonstrated no clear pattern with climate (Owen et al., 2014), further challenging the idea that alluvial fans are a climate indicator. Hall and Goble (2012, their Figure 2) described Late Pleistocene paleosols in New México, and measured two OSL dates (81 ± 6 ka and 91 ± 7 ka) for the lower aeolian sand deposited after the Mescalero paleosol, but before the Berino paleosol (this paleosol is dated as late OIS Stage 3 to early OIS Stage 2). Thus the lower aeolian sand represents a dry period between two wetter periods as inferred for the paleosols. This aeolian sand is older than the Ciudad Juárez P3 sample at 74.3 ka but the two do overlap within error at 75 - 79 ka. The aeolian sand from which P3 was collected is overlain by coarse-grained conglomerates that were presumably deposited by high or rapid water flow indicating a humid phase after aridity. Thus we make a tentative correlation between the oldest age of the three dates on Juárez fans on one hand, and depositional events in New México on the other, reflecting a regional climate feature.
McDonald et al. (2003) described the results of climate shifts between latest Pleistocene to Holocene in the Mojave Desert, California, highlighting the increase in wetness that spurred enhanced erosion and transport of sediment, leading to fan development. They also noted that aeolian and alluvial sediments in the region were deposited at approximately the same time, indicating widespread climate effects on erosion and sedimentation processes. Dorn (2009, Figure 24.1) drew attention to larger-scale climatic events in relation to the late Quaternary. In particular Dansgaard-Oeschger Cycles 11 and 12 from the Greenland ice cores that reflect moisture changes in the northern hemisphere are close in time to the Ciudad Juárez P2 sample. In contrast to interpretations of aridity discussed above for the 74 ka date in this paper, there is evidence in Greenland of temperature rise (implying increase moisture) in the study of Dansgaard-Oeschger cycles by Landais et al.(2004, see their Figure 3) at the same time, suggesting therefore that the Ciudad Juárez fans do respond to regional climate controls.
It remains unclear what the influence of short-term climate events, such as ENSO (El Niño - Southern Oscillation) episodes, is on the deposition of sediments in the Ciudad Juárez alluvial fans. Magaña and Conde (2000) demonstrated a tentative link between ENSO events and winter flooding in southern USA and northwestern México. On the other hand, Pavia et al.(2006) showed finer resolution of identification of wet phases involving both ENSO and Pacific Decadal Oscillations. Thus the Ciudad Juárez region may have been influenced more by Pacific than by Atlantic Ocean processes, and is an avenue for further investigation, although whether it will be possible to extend such interpretations back into the Pleistocene of the study area or not, remains an open question.
The implications of these results for the alluvial-fluvial system studied here in the Ciudad Juárez area provide an avenue for future work. Thus the OSL ages reported here may be indicators of climatic fluctuations in northernmost México through the late Quaternary that relate to broader scale climate changes, but so far our dated sediments do not allow a firm correlation with wider scale events, and change into the Holocene. Thus we regard our results as preliminary, but showing great promise, and emphasize the need for further dating work.
7. Conclusions
This study has provided the following outcomes:
The first numerical ages of quaternary alluvial fans in northern Mexico as measured by OSL in three samples are: a) sample P1 contains a minimum age (S1) of 31 ka and a second population (S2) at 45 ka; b) sample P2, a minimum age (S1) 41 ka and second population (S2) at 69 ka; and c) sample P3 has one dominant population (S1) at 74 ka. We consider the younger of the two ages for P1 and P2 to be the more reliable. These ages are within the range of OIS 3 to 5. They demonstrate that fan development in the area, now occupied by the city, either terminated in the Late Pleistocene, after which there was apparently an extended period of aridity, or was not preserved when the Río Bravo incised the distal fan layers.
The oldest age (74 ka) corresponds to the developing glaciation of OIS Stage 4, but the younger two ages relate broadly to the development of the last glacial. These results imply that aridity was locally to regionally controlled, with limited evidence of a global overprint. However, this interpretation is tentative and further research is necessary. Alluvial fan margins inter-finger with fluvial terrace sediments derived from the Río Bravo indicating an additional component of fan dissection by the Río Bravo lateral erosion, presumed to be active during earlier times than our OSL ages, but these are not yet dated.
The age of the fans are broadly correlated with sand sheets and alluvial fans of other desert areas of the USA, including the Mojave Desert, and there are preliminary indications of relationships between the Ciudad Juárez fans and climate change in the region.
The fan system has been potentially reactivated to form the Anapra fan in the northern part of the study area, but detail is lacking because incision in the distal fans has been affected by the Río Bravo eroding the stratigraphy as the river changed course.
These preliminary dates provide an avenue of investigation for further dating to ascertain the controls on the alluvial fan and fluvial system.