Heat Flow Data in an Area of the Eastern Southern Basin and Range in Arizona Contribute to an Analysis of Neogene Lithosphere Thinning Greater than 100 Km

Obtaining data-driven estimates of lithosphere thickness ~25Ma and Neogene lithosphere thinning in different geological provinces of the western United States may provide valuable boundary conditions for geologic models of Neogene tectonics. Geothermal information providing subsurface temperature estimates relevant to a past geologic time can be considered with present heat flow data to estimate the related past lithosphere thickness as well as lithosphere thinning over the corresponding geologic interval. For example, mineralogy and heat flow studies were used to estimate heat flow ~25Ma for the Four Corners Area of the Colorado Plateau [1, 2] and to suggest Neogene lithosphere thinning of ~100 km for the area. A heat flow-temperature-depth expression, the same as in this study, was used to match upper mantle depth-temperature estimates at ~25Ma consistent with the mineralogy data relevant at ~70 and~140 km depth [1, 2]. From this analysis, a heat flow estimate of ~47mWm for 25Ma was possible, and therefore, an estimate of lithosphere thickness could be made [1]. Heat flow analyses also provided estimates of present lithosphere thickness from the present heat flow data for the Colorado Plateau, which was consistent with seismic studies [1]. Lithosphere thinning as proposed for the Colorado Plateau results from a mantle delaminating process or upper mantle layer separation [1, 3]. Interestingly, Axen et al. [4] present a very different mechanism for Laramide tectonics in which flatslab subduction plows off 20-50 km of lower continental lithosphere across large regions of the western United States. In the present study, I employ thermochronologic derived near-surface geothermal gradients relative to the GeoScienceWorld Lithosphere Volume 2020, Article ID 8872618, 9 pages https://doi.org/10.2113/2020/8872618


Introduction
Obtaining data-driven estimates of lithosphere thickness 25 Ma and Neogene lithosphere thinning in different geological provinces of the western United States may provide valuable boundary conditions for geologic models of Neogene tectonics. Geothermal information providing subsurface temperature estimates relevant to a past geologic time can be considered with present heat flow data to estimate the related past lithosphere thickness as well as lithosphere thinning over the corresponding geologic interval. For example, mineralogy and heat flow studies were used to estimate heat flow~25 Ma for the Four Corners Area of the Colorado Plateau [1,2] and to suggest Neogene lithosphere thinning of 100 km for the area. A heat flow-temperature-depth expression, the same as in this study, was used to match upper man-tle depth-temperature estimates at~25 Ma consistent with the mineralogy data relevant at~70 and~140 km depth [1,2]. From this analysis, a heat flow estimate of 47 mW m -2 for 25 Ma was possible, and therefore, an estimate of lithosphere thickness could be made [1]. Heat flow analyses also provided estimates of present lithosphere thickness from the present heat flow data for the Colorado Plateau, which was consistent with seismic studies [1]. Lithosphere thinning as proposed for the Colorado Plateau results from a mantle delaminating process or upper mantle layer separation [1,3]. Interestingly, Axen et al. [4] present a very different mechanism for Laramide tectonics in which flatslab subduction plows off 20-50 km of lower continental lithosphere across large regions of the western United States.
In the present study, I employ thermochronologic derived near-surface geothermal gradients relative to the beginning of extension at a location in the southern Basin and Range~87 km north of Tucson, Arizona (the Grayback normal fault block in the Tortilla Mountains; [5,6]; Figure 1). Extension is thought to have begun~25 Ma or 25-20 Ma [5,6]. This location is quite unique in that more than a dozen present-day heat flow measurements are also available in the area. Differing from the study in the Four Corners area, near-surface paleogeothermal gradients combined with appropriate thermal conductivity data (from [7,8]) allow preextension heat flow to be estimated at the study location in the southern Basin and Range. Using the derived heat flow estimate for~25 Ma, the cotemporal lithosphere thickness is calculated. From present-day heat flow data over the study area, the present lithosphere thickness is also calculated, which agrees with seismic studies of lithosphere depth. These calculations show that a significant increase in heat flow over the past~25 Myr is likely coupled with a major decrease in lithosphere thickness.
The present study area in the southern Basin and Range has experienced considerable Neogene extension and crustal thinning whereas the Four Corners area of the Colorado Plateau has not experienced Neogene extension [9,10]. In this study, it will be shown that these two very different geologic regions have comparable pre-Neogene (25 Ma) heat flow and rather similar pre-Neogene (25 Ma) lithosphere thickness, as well as a large amount of Neogene lithosphere thinning. Lithosphere separation is proposed as the Neogene lithosphere thinning process in the Four Corners region and may have occurred in the southern Basin and Range as well. However, additional advective mechanisms such as a primary thermal plume, convection cells, and magma migration along conduits are likely present for southern Basin and Range development [11,12]. Crustal and mantle lithosphere extensions are considered to show the necessity of upper mantle asthenosphere and lithosphere advection in the southern Basin and Range. The present study will develop numerical estimates for southern Basin and Range Neogene lithosphere thinning in Arizona.

Data Presentation and Site Relevance
Near-surface geothermal gradients, determined by thermochronologic methods, relating to pre-Neogene (just before crustal extension began~25 Ma) in the study area ( Figure 1) are rather low (17 ± 5°C km −1 , [5]; 14 ± 7°C km −1 , [6]). The two values average 15:5 ± 2:1 ± 1:5°C km −1 (average ± standard deviation ± standard error of the mean; more study values may provide a more certain mean and reduce the mean uncertainty). Although the study location has undergone a high degree of extension, the location is not within a metamorphic core complex [6,13]. The present crustal thickness at the study area is~28.5 km, almost the same as the characteristic thin crust of 28:2 ± 0:5 km for the southern Basin and Range in Arizona [10]. The location is amagmatic,  Dots represent heat flow sites considered in study area, X represents radiogenic sites considered in study, and triangles represent sites where temperature logs were made in boreholes > 500 m in plutonic rocks [7,8]. Squares in study area indicate additional considered data [13]. 2 Lithosphere as opposed to two other paleogeotemperature gradient study sites (with higher calculated temperature gradients) located in the lower Colorado River area in a region of major extension [5,14]. It is suggested therefore that the study locale is a rather representative location of Neogene geologic development occurring in many areas of the southern Basin and Range of Arizona. Sixteen present-day heat flow measurements at sites within~50 km of the study locale yield a mean value of 83:1 ± 14:4 ± 3:7 mW m −2 [7,8,13]. Including values within the 99.9% confidence range, the mean of thirteen heat flow values is 84:5 ± 7:7 ± 2:2 mW m −2 ( Figure 1). This value happens to be very similar to the average of seven heat flow values taken over much of the southern Basin and Range in Arizona in plutonic rocks to depths greater than 500 m (84:5 ± 6:1 ± 2:5 mW m −2 ; [7,8]; Figure 1; deeper tests in plutonic rocks typically allow longer conductive temperature gradient intervals to be observed). Neogene erosion of~1 km is suggested for the study locale but is also questioned [15]. The present heat flow value for the study area is therefore reduced by about 0-3% (calculated after [16]), giving a heat flow of~83 mW m -2 . This value is quite similar to the average of a large number of values across the southern Basin and Range (82 mW m −2 ; [13]; 157 values). Therefore, the study location may also be rather representative of geothermal evolution over the past 25 Myr in many areas of the southern Basin and Range in Arizona.
Although the mean thermal conductivity for the thirteen heat flow sites in various rock types for the study area is 3:26 ± 0:30 ± 0:10 W m −1°C−1 (giving a present-day heat flow of 83 mW m -2 ), the thermochronologic study site is located specifically in plutonic-granitic rocks [6]. As such, it is suggested that the mean thermal conductivity value for the seven deep tests greater than 500 m in plutonic rocks across the southern Basin and Range in Arizona best estimates the thermal conductivity at the thermochronologic study location (discussed above). The mean thermal conductivity value for the seven plutonic rock heat flow sites is 3:04 ± 0:44 ± 0:17 W m −1°C−1 . It follows the pre-Neogene heat flow is~47 ± 5 mW m -2 , using the statistical calculation for (15:5 ± 1:5°C k m −1 × 3:04 ± 0:17 W m −1°C−1 ). Therefore, it is proposed that the near-surface heat flow at the study site in the southern Basin and Range has changed from~47 to~83 mW m -2 over the past~25 Myr; the change in the heat flow at the Four Corners area over the same time interval is~47 tõ 64 mW m -2 [1].

Temperature-Depth Calculations
The expression used to arrive at first-order approximations to lithosphere temperatures is the steady-state onedimensional heat flow equation (from equation 4, [17]): where K is thermal conductivity, TðzÞ is temperature at depth z, q is near-surface heat flow, A o is radiogenic heat production at the top of a granitic crustal radiogenic layer where radioactive heat production theoretically decreases exponentially with depth, and D equals 10 km (D is the depth above which all but 1/e th of radioactivity occurs). The value presently calculated for A o in the southern Basin and Range of Arizona is 2:6 ± 1:9 ± 0:8 μW m −3 (average for six granite and granodiorite sample values from boreholes ≥ 200 m depth, average depth 482 m) or 2:3 ± 0:8 ± 0:4 μW m −3 for four of the sites if values within the 99% confidence level are considered [7,8]. Sass et al. [13] arrive at a somewhat lower value for radiogenic heat production from crystalline rock over the southern Basin and Range (2.1 μW m -3 ).
Using Equation (1) implies that lithosphere temperatures have been relatively constant for an extended geologic time, which may be hypothesized as a limiting approximation for the pre-Neogene (~25 Ma) initial time condition but will be further discussed for the present. As well as model considerations, the thermal conductivity and radiogenic distributions are always somewhat speculative. Sass et al. [13] discuss nonlinear variability of radiogenic heat production A o vs. heat flow in the Basin and Range resulting from extension and magmatism. Many models of crustal radiogenic distribution (~500, [18]), as well as limited data, contribute to calculated lithosphere temperature uncertainty. The lithosphere thermal conductivity values in the present study are estimated to 200 km depth from average rock type values weighted by the appropriate interval length proportional to the total depth considered (rock type given and extrapolated from [10]; thermal conductivity values from [19]). Interestingly, temperature and pressure effects on mafic lithosphere rock thermal conductivity appear small compared to variation within rock type at room temperature [20].

Lithosphere Temperatures and the Lab
(Lithosphere-Asthenosphere Boundary) Figure 2 shows the calculated lithosphere temperatures (using Equation (1)) for heat flows of 47 and 83 mW m -2 .
The temperatures related to 47 mW m -2 estimate the initial temperature distribution proposed for the present study site at about 25 Ma or just before the start of extension. The LAB is suggested to occur at 1300°C or 1573 K [18,21]. From Figure 2, the LAB appears to occur at~184 km depth corresponding to near-surface heat flow of 47 mW m -2 . As mentioned above, the calculated temperature distribution is sensitive to the near-surface radiogenic heat concentration (A o in Equation (1)). For example, a valve of 2.6 μW m -3 , instead of 2.3 μW m -3 , would indicate 1300°C at about 210 km depth for a heat flow of 47 mW m -2 , instead of 184 km. In both cases, the preextension lithosphere is suggested to be quite deep, generally consistent with (and possibly the same as) estimates of LAB depths for early Proterozoic areas (200-220 km, [18]; the study site is located in a Paleo-Proterozoic area, [22]). If A o equals 2.1 μW m -3 [13] instead of 2.3 μW m -3 as used above, the lithosphere thickness using Equation (1) would still be quite deep for a heat flow of 47 mW m -2 , 170 km. Figure 2 also shows the calculated present-day temperature-depth profile derived from Equation (1) (1)). This estimate of the LAB is in quite good agreement with seismic study estimates of the LAB for the study location (~70-72 km; [23,24]). Agreement of the LAB depth between several seismic studies and the steady-state conduction temperature model suggests that advection of heat is moving the temperature distribution much faster than conduction alone. Lachenbruch and Sass [17] also present a steady-state thermal model of lithosphere extension by mantle stretching and crustal intrusion; the predicted LAB is the same as above. Agreement of the LAB depth between the two thermal models likely results from different thermal conductivity and A o estimates. From these estimates of heat flow~25 Ma and present day, considerable thinning of the lithosphere in the southern Basin and Range of Arizona over the past 25 Myr is suggested, >~100 km (~114 km = 184 -70 km). This is somewhat larger than predicted for the Four Corners region in the Colorado Plateau (~100 km; [1]).
In an area of about 50 km in radius centered at 111°W and 32°N, just south of the study location, there are 19 heat flow values [7,13]. Considering the data within a 99.9% confidence level and applying a small Neogene erosion correction as discussed above, the average heat flow is 87 mW m -2 . The resulting LAB calculated as above is 65 km. From seismic studies, the LAB is also~65 km [23,24]. Although we have no preextension thermal gradient estimate at this second location, the elevated presentday heat flow data correlating with a relatively thin lithosphere supports the notion that thinning of the lithosphere relates to a thermal source which may be the fundamental cause for extensional tectonic development of the southern Basin and Range.

Geothermal Considerations of Lithosphere Thinning
Why the steady-state conduction thermal model (Equation (1)) appears to present accurate estimates of the present LAB depth is rather enigmatic (of course more study sites are needed). Lachenbruch et al. [25] suggest that various advection phenomena during southern Basin and Range extension followed by conductive cooling may be consistent with present heat flow in the southern Basin and Range. The 83 mW m -2 temperature profile represents steady-state conduction ( Figure 2). Advection processes in the asthenosphere and mantle lithosphere that would bring temperatures of 1300°C from~184 km depth (25 Ma LAB estimate) tõ 70 km depth (present LAB estimate, Figure 2) beginning 25 Ma must be accompanied by additional advection processes in the mantle lithosphere which bring warmer temperatures closer to the earth's surface than 70 km. This is because the conduction time constant for 70 km is~39 Myr (the conduction time constant is the time it takes for the near-surface temperatures to reach within 1/e th of steady state assuming a step continuous temperature increase). If conduction were the sole means of heat transport from 70 km to the surface, while 1300°C was initiated and maintained at 70 km depth

Crustal and Mantle Lithosphere Thinning Calculation
Pre-Neogene (25 Ma) crustal thickness estimates depend on Neogene extension estimates and present crustal thickness. The characteristic crustal thickness for the southern Basin and Range of Arizona is 28.2 ± 0.5 km, at the study site I esti-mate~28.5 km [10]. Extension over the entire southern Basin and Range from~25 Ma (El Paso TX to the North American Plate margin) is~57% (change in width/initial width; [9]); the ratio of the present width to the initial width is then 1.57. If the cross section area of the crust along cross section remains constant during extension ( [10], suggest from seismic data a simple extension of the upper and lower crust),  57%; [9,27]). The preextension crustal thickness in these areas for the present model (~74 km) would be quite large even compared with much of the present southern Rocky Mountains' crustal thickness (48 -50 km; [27]). This suggests that additional processes contributed to crustal thickness and crustal extension in these regions of the northern Basin and Range, e.g., possibly extensive crustal under plating and sill and dike formation.
The above calculations suggest much greater thinning in the mantle lithosphere than in the crust for the southern Basin and Range. If a given cross section area of the mantle lithosphere remains the same during the thinning process, then as the mantle lithosphere thins, it should extend much more than the crust extends; alternatively, additional vertical advection in the mantle lithosphere accompanies extension and contributes to its thinning.

Heat Input to the New Lithosphere and Crust as Heat Flow Increases from 47 to 83 mW M -2
One may estimate the immense heat input to the lithosphere at 70 km accompanying the LAB change from 184 to 70 km depth and the heat flow increase from 47 to 83 mW m -2 by subtracting the area under the 47 mW m -2 temperaturedepth curve, from the surface to 70 km, from the corresponding area under the 83 mW m -2 curve ( Figure 2); i.e., ΔQ ðchange in heatÞ = ðΔ integrals of Tdz from 0 to 70 kmÞ × c × ρ, c is specific heat and ρ is density. Bashir et al. [10] suggest Vs/Vp and R data for the southern Basin and Range in Arizona, referenced to a typical cratonic region, indicate simple extension of both the upper and lower crusts (R is amplitude of P to S converted wave). Therefore, as a first-order approximation allow the pre-Neogene (~25 Ma) upper and lower crust thicknesses to have the same ratio as today (upper crust/lower crust~1.5, but does vary along profile, [10] To better illustrate the magnitude of this amount of heat, compare with a surface basalt flow of height "z", the associated heat is Q = ρ × c × T × z + L × ρ × z. Choosing density and specific heat as above, temperature T = 1250°C, and latent heat of solidification L = 4 × 10 5 J kg −1 [29], the amount of heat per horizontal square meter in the basalt flow is~4:62 × 10 9 J m −3 × z. Equating to the amount of heat input at 70 km when heat flow is increased from 47 to 83 mW m -2 , the estimate z is~2:43 × 10 4 m, or 24.3 km, a remarkable height. Proceeding as above, consider the heat input into the new crust at 28.5 km as heat flow changes from 47 to 83 mW m -2 . The calculated heat input at 28.5 km is~16:6 × 10 12 J m −2 , equivalent to~3.6 km of basalt flow. As heat flow changes from 47 to 83 mW m -2 , heat input into the mantle lithosphere would be~9:54 × 10 13 J m −2 (11:2 × 10 13 -1:66 × 10 13 ), or about 5.75 times the heat input into the crust.

Discussion
Heat flow and seismic data integrated with thermochronologic paleogeotemperature gradient data provide an estimate of considerable lithosphere thinning (~114 km or~62%) over the past~25 Myr at the study location in the eastern southern Basin and Range of Arizona. The difference between crustal thinning and mantle lithosphere thinning proposed from the above calculations, as well as the conduction time constants, suggests considerable vertical upward advective heat transfer to thin the mantle lithosphere. The study location is an amagmatic area with considerable extension, but not in a metamorphic core complex, has a heat flow almost the same as the southern Basin and Range average heat flow (83 vs. 82 mW m -2 ), and a crustal thickness essentially the same as that characterizing the southern Basin and Range in Arizona (28:2 ± 0:5 vs.~28.5 km; [5,10,13]). The change from subduction to divergence between the North American and Pacific plate boundary is followed by the beginning of extension across the southern Basin and Range and at the study area as well [5,6,9,30]. For these reasons, I suggest the study area as a potential location for first-order representation of geothermal and extensional processes in many areas over the southern Basin and Range of Arizona. Figure 3 shows depth estimates from several studies for the Moho and LAB along 33°N across southern Arizona. The LAB shallows dramatically west of 114°W approaching the Salton Sea. East of 114°W, the LAB depth shallows gradually to 110°W and then deepens somewhat to 109°W. The depth of 70 km at 111°W (the study location) appears to be an approximate average depth east of 114°W, although the LAB depths do vary along the profile east of 114°W about 5 Lithosphere 10 to 16 km depending on the study. It is noted that the depth estimates are ± several km while interpolating the colorcoded maps for LAB depth. Further study may better correlate heat flow with the LAB; this has been done in the present study for another site near 111°W and 32°N (above). The crustal depth estimates shallow approaching 115°W and are within several km of one another eastward to 111°W at which point depths appear to diverge going eastward. The study by Bashir et al. [10] shows only a couple km difference in crustal depths east of 114°W along 33°N; the map by Gilbert [27] shows rather uniform crustal thickness of about 28 km across the entire southern Basin and Range of Arizona.
To further examine lithosphere thinning in the southern Basin and Range, consider a model where mantle lithosphere extension across the southern Basin and Range is limited by the western boundary of the North American plate and the eastern boundary of the Pacific Plate and is therefore the same as the crustal extension. One may calculate potential thinning of mantle lithosphere due to extension or spreading for a straightforward geometric model versus separate vertically upward advection not explicitly a component of lithosphere extension (extension or spreading and thinning the mantle lithosphere transports heat both horizontally and vertically, see [17], their Figure 9-8). As before, allow the model cross sectional area along the cross section to remain constant during mantle lithosphere extension and the LAB temperature to also remain constant. In this model, the upward advection associated with mantle lithosphere spreading results from the LAB rising with mantle lithosphere thinning inversely proportional to extension. As mentioned before, overall east-west extension across the southern Basin and Range during the past~30 to 25 Myr is~57% (El Paso, TX to the North American plate margin; from [9]). Of course this is an average, but from Figure 3 and previous discussions, it would appear the crustal and LAB depths, as well as extension and heat flow, at the study site are reasonable average approximations across the southern Basin and Range in Arizona. If extension was 57%, the lithosphere thickness after extension may be calculated as the initial depth divided by final depth z: (184/z = 1:57) giving the LAB depth z as 117 km. This represents a first-order estimate, as per the present model, of the effect of vertical heat transport associated with mantle lithosphere extension and thinning. As such, in the present model, additional processes advecting heat vertically upward and involved with lithosphere thinning are necessary to account for the LAB at 1300°C to move an additional~47 km upward from~117 km to 70 km ( Figure 2). Some of the fundamental advection mechanisms suggested are a rising primary mantle plume inducing convection cells, delaminating parts of the mantle with associated upward advection, and magma movement along conduits (a few references are, respectively, [3,11,12]).
The same basic geometric model can be applied to the crust and upper mantle lithosphere. If crustal extension and thinning occurred in such a manner as to preserve Moho temperatures, then preextension temperatures at~44.5 km (429°C for 47 mW m -2 ) would be brought to~28.5 km. However, the present heat flow suggests a temperature of~656°C at 28.5 km; this requires additional heat transfer involving advection in the thinning mantle lithosphere (Figure 2).
Seismic studies of rifted regions in southern California suggest that mechanisms of lithosphere deformation and strain accommodation are unclear; however, a lack of systematic offset between the lower lithosphere and crustal surface deformation agrees with the notion of lithosphere symmetric extension [26]. Colocation of crustal surface rifting evidence and a flat and shallow LAB suggests that deep lithosphere extension and crustal extension are related [26].
Increased elevation and crustal thickness as well as lithosphere thinning can provide deviatoric stresses promoting southern Basin and Range extension; lithosphere thinning with asthenosphere replacement can also increase elevation [30]. Lithosphere thinning warms the remaining lithosphere which decreases its viscosity allowing for easier deformation. Thinning of lithosphere thickness by over 100 km at the study location should provide considerable deviatoric stresses and initially preextension elevation increase [30]. The change from a convergent to divergent plate boundary just beginning~30 Ma [9] appears to have allowed deviatoric tension to proceed with extension.
Yuen and Fleitout [12] suggest that a primary thermal plume causes strong small-scale convection that thins the lithosphere and reduces viscosity. Their model predicts observed rapid rates of uplift and lithosphere thinning. Lachenbruch et al. [25] suggest that the history of extension in the southern Basin and Range (occurring mainly between 28 and 16 Ma; [30]) indicates advection rapidly changing  The above discussions indicate mantle lithosphere thinning results because of both mantle lithosphere extension and additional vertical heat advection; but mantle lithosphere extension is in turn dependent on lithosphere thinning. A thick crust and a thinning lithosphere result in a potential energy configuration most likely to develop deviatoric stresses [30]. Both the crust and the lithosphere in the southern Basin and Range of Arizona were much thicker~30-25 Ma (~44.5 and~184 km, respectively, as calculated before) when the Pacific and North American plates just began to diverge. Because a thick crust is conducive for deviatoric stress, it would seem reasonable that lithosphere thinning began with thermal erosion near the LAB caused by a large primary thermal source. The seismic data indicating warmer conditions to depths of at least 200 km under almost all of the southern Basin and Range suggest a very large thermal anomaly at depth ultimately associated with a primary thermal upwelling or plume, subsequent convection to contribute to lithosphere thinning as well as magma movement along conduits is also likely occur [11,12,32].
With the thick crust and divergent plate motion, the thinning mantle lithosphere may promote sufficient deviatoric stress to initiate extension. One may suggest that a feedback mechanism could have been operating in the southern Basin and Range. As the mantle lithosphere continues to thin, more extensional stress is generated. However, as the crust extends and thins, extensional stress is decreased. A decrease in thermal activity associated with mantle lithosphere thinning may have brought about a balance between decreasing stresses generated by crustal thinning and increasing stresses generated by mantle lithosphere thinning, arriving at a cessation in southern Basin and Range extension.
Processes that may produce the present lithosphere depth as well as the estimated Moho temperature (Equation (1)) would establish the present LAB by asthenosphere and mantle lithosphere advection and lithosphere extension; upward advection in the thinned mantle lithosphere would establish a new elevated shallower Moho temperature. The crustal temperatures could respond conductively rather quickly to the new Moho temperature (relatively short time constant as discussed above). The high thermal gradient just below the Moho caused by advection may supply heat conductively to maintain Moho temperature after advection and extension stop (after [25]). Crustal intrusion and sill formation have not presently been included. Temperatures shown in Figure 2 between the present Moho and LAB are probably not linear as shown but are instead convex upward as per temperature profiles showing upward fluid movement across a layer. Long-term thermal relaxation should approach temperatures shown in Figure 2 as long as the LAB remains at 70 km.

Conclusion
The increase in heat flow sometime during the past~25 Myr as measured at the study site in the southern Basin and Range of Arizona identifies thinning of the lithosphere caused by thermal phenomena as a probable fundamental driving force of southern Basin and Range extension. The straightforward model presented above suggests lithosphere thinning may be accomplished by both extension and additional thermal advection, while extension is allowed by the change to divergent plate motion.

Data Availability
The data sources are all given in the references cited in the text and presented in the reference list. Several of the primary references are (1)

Conflicts of Interest
There is no conflict of interest with respect to the present manuscript.