**1. Introduction**

Previous studies have shown fault slip analysis at the outcrop scale provides a means to deduce the orientation of the principal stress fields and their evolution through successive tectonic events [1–7]. Additional information obtained from other structures, such as joints [8] tension gashes, and stylolites [9], is also important but will not be presented here. In this paper, we define a fault as simply a parting in rock with no claim whether it formed as a Mode 1 (opening), Mode 2 (shearing), or Mode 3 (tearing) [10]. If a fault shows shear o ffset (Mode 2). The input for fault slip analysis is field data collected on the surfaces of individual faults which includes the orientation of the, slip direction, and sense of slip. The latter two are determined by one or more of the following displacement indicators visible on the fault surface: slickensides, asperity ploughing, slickolite spikes, crescent marks, the growth of mineral patches on the lee side of hills on a rough fault surface, mineral fibers and steps, and Reidel shears [11].

The basic assumptions behind fault slip analysis are that: (1). conjugate fault sets result from a single brittle deformation event, and (2). slip on a fracture surface occurs in the direction of maximum resolved shear stress. The first step in the analysis consists of reconstructing the "reduced stress tensor". The reduced stress tensor di ffers from the actual stress tensor only in that the absolute magnitudes of the principal stresses: σl (maximum compressional stress), σ2 (intermediate stress), and σ3 (minimum stress) are not determined, only their relative magnitudes. However, the relative magnitude, order, and orientation of the three principal stresses are the same as for the actual stress tensor and enable one to define the directions of compression and extension which prevailed during tectonic events. Knowing the stress state, one determines the shear stress and hence the slip orientation expected on any plane. The first attempt at formulating and solving the inverse problem was [12]. Numerical

methods have since been developed for reconstructing paleo-stress orientations from fault slip data. In the general case illustrated in this paper, any planar discontinuity in a rock may be activated as a fault. The discontinuity may be either a pre-existing fracture activated or reactivated (inherited fault) by the tectonic stress. The basic properties of the reduced stress tensor and its determination is summarized below. The indicators of the direction and sense of shear on a discontinuity reactivated in shear is to collect and analyze fault slip data. The method of direct inversion used in this paper can be found in [1–6,13,14] and in [15].

### **2. Geologic and Tectonic History of the Southwestern White Mountain Region**

The bedrock geology of the study site and Hubbard Brook valley has been mapped (sheet 1) at a scale of 1:12,000 [16] and at a scale of 1:10,000 [17]. The study site is also included in earlier bedrock maps [18,19] and at a scale of 1:250,000 [20]. The bedrock underlying the study site (Figure 1) consists of metamorphic rocks intruded by igneous rocks and belongs to the Central Maine Trough [20]. The metamorphic rocks of the Rangeley and Perry formations were deposited as sandy, clay-rich marine sediments [21] on a continental shelf, rise, and abyssal plain of the Rheic Ocean [22] over a time spanning the Silurian (approximately 443–428 Ma). These sediments were then buried and multiply folded by at least two deformation episodes [17] in the Acadian orogeny (early Devonian, 410–390 Ma), as the Rheic Ocean closed and Avalon collided with eastern North America. At this time, the rocks were metamorphosed to the lower sillimanite grade (approximately 600 ◦C and 4 kb pressure, equivalent to a burial depth of approximately 15 km), which resulted in local melting (migmatization). At approximately 410 Ma, these rocks were at or near the conditions of maximum pressure and temperature and were intruded locally by the southern portion of a large pluton of the Kinsman granodiorite of the New Hampshire plutonic series. The Kinsman granodiorite underlies most of the western half of the Hubbard Brook Valley immediately to the west of the study site [16] see Figure 1. Near the bottom of Figure 1 is a low angle thrust called the Thornton Fault on [20]. This fault thrusts older Silurian rocks over younger Devonian rocks. The fault is cuto ff by and therefore must be older than the intrusion of the Kinsman Granodiorite, but younger than the Littleton Formation. This fault extends below the study site and below the Rangeley Formation at the study site. This fault may have formed during Tectonic Event 1 in Tables 1 and 2. During the late Devonian (370–365 Ma) the metamorphic rocks and the Kinsman intrusion were multiply intruded by small tabular dikes and small discordant bodies of Concord granite, also of the New Hampshire plutonic series. The Concord granite is shown on the map at a scale of 1:200 [16] and in well logs [23], but is not shown in Figure 1 which is modified from the map of [20] (scale of 1:500,000).

The Alleghenian orogeny (325–260 Ma) created the Appalachian Mountains principally by collision with North Africa. While it may not have resulted in large scale deformations at the study site, it was strong enough to create or reactivate fractures.

From early to late Jurassic (194 to 155 Ma) the are immediately to the east and north of the study site was a region of extensive granitic intrusion expressed by the huge batholiths and ring dikes of the White Mountain plutonic/volcanic series [24]. During that time or possibly later (130 to 100 Ma), the metamorphic and igneous intrusive rocks at the study site were intruded by tabular diabase dikes, emplaced as part of continental rifting associated with the opening of the present Atlantic Ocean basin.

From the time of the Acadian Orogeny to the present, erosion and uplift have brought the bedrock from a depth of approximately 15 km to at or near the Earth's surface. The last episode of deformation was the loading and unloading of the bedrock by the advance and retreat of multiple glacial ice sheets over this region in the past 100,000 years [25]. Reconstructions of the thickness of the Laurentide ice sheet yield a glacial ice loading and unloading of three kilometers for New England [26].
