Resin Production in Natural and Artificial Red-Cockaded Woodpecker Cavity Trees

Resin flow was measured in red-cockaded woodpecker (Picoides borealis Vieillot) clusters in longleaf pine (Pinus palustris Mill.) in the southern region of the Angelina National Forest, Texas. Resin flow (ml) at 1.4 m height over 24 hrs was measured from one 2.5 cm punch through the phloem between 0700 and 1000 hrs from March 1999 to September 2000, for a total of 9 measurements per tree. Resin was sampled in naturally active cavity trees, artificial (insert) active, natural inactive, artificial inactive and control pines (84 sample trees). Resin flow pattern was significantly different during the year, but not significantly different in the cavity tree type resin flow. Cavity trees in the 90th percentile (>33.0 ml resin in 24 hrs) were defined as super resin producing. High average resin flows in August 1999 and September 2000 indicate when to sample resin for potential cavity trees. Regression equations were produced to estimate future resin production.


Introduction
The federally listed endangered red-cockaded woodpecker (RCW hereafter) (Picoides borealis Vieillot) (U. S. Department of Interior, 1985) is endemic to pine forests in the southeastern United States from Virginia to Florida and west to Texas (Jackson, 1971, Hooper et al., 1980Lennartz et al., 1983;James, 1995). The RCW is a cooperative breeder where young males remain in their natal groups as helpers (Ligon, 1970;U. S. Department of Interior, 1985;Lennartz et al., 1987;Walters et al., 1988;Walters, 1990;Conner and Rudolph, 1995). Redcockaded woodpeckers excavate roost and nest cavities into the heartwood of living pines, generally 85 -350 years old . Upon cavity completion, RCWs excavate shallow holes in the xylem, termed resin wells, and peck them daily to stimulate copious resin flow (Ligon, 1970;Dennis, 1971), to protect against predators and some nest competitors (Jackson, 1974;Jackson, 1978;Rudolph et al., 1990). Bark beetles are less successful in invading and killing southern pines with high resin flow (Lorio et al., 1990).
Artificial cavity construction (Copeyon, 1990;Copeyon et al., 1991;Taylor and Hooper, 1991) is an effective method to induce red-cockaded woodpeckers to occupy previously unoccupied habitat (Hooper et al., 1990;Watson et al., 1995). A tree characteristic not directly considered in the tree selection process for artificial cavities is the pine's resin yielding capabilities. When selecting a potential artificial cavity tree, tree diameter at cavity height (38 cm minimum), trees with larger crowns and stand location is important. Resin yield can be sampled following the methods of Lorio and Sommers (1986) and Lorio (1990). Resin production from preformed resin ducts is a primary defense against bark beetle attack (Hodges et al., 1979;Paine et al., 1985;Cook and Hain, 1986;Mitchell et al., 1991;Bowman and Huh, 1995;Coulson et al., 1995;Ross et al., 1995;Ross et al., 1997;Conner et al., 1998). Resin sampling to estimate resin from preformed resin ducts follows the methods of Lorio (1990) for measurement from 0700 -1000 hrs. Ross et al. (1997) examined resin production in RCW longleaf pine cavity and non-cavity trees on the Angelina National forest. Bowman and Huh (1995) examined resin flow in Florida in two vegetation communities. Losses to hurricanes and bark beetles can cause loss of RCW cavity trees (Hooper et al., 1990;Conner et al., 1998), so both methods of evaluation for resin production for potential RCW trees and current resin status of existing cavity can assist in long term management of the species.
The objectives of this project were to evaluate resin production in longleaf pine, Pinus palustris, in redcockaded woodpecker clusters in the Angelina National Forest, Texas. Comparisons in resin production were made between natural active and inactive, and artificial natural and inactive cavity trees. Similar sizes of longleaf pines were sampled for resin as control trees.

Study Area
The 61,988 ha of the Angelina National Forest (31˚15'N, 94˚15'W), is divided by Sam Rayburn Reservoir into approximately equal areas forming the southern and northern sections. Located 45 km east of Lufkin, Texas, portions of the forest are in Angelina, Jasper, San Augustine and Nacogdoches counties. The study was conducted on the southern portion dominated by longleaf pine (Pinus palustris) on xeric soils of the Tehran (loamy, siliceous, thermic, Grossarenic Palendult) and Letney (loamy, siliceous, thermic, Arenic Palendult) soil series (USDA Natural Resource Conservation Service, 2010). Low water holding capacity makes regeneration difficult. The study took place in six national forest compartments: 78-1, 78-2, 88-1, 88-2, 92-1 and 92-2 located within the USDA Forest Service RCW Habitat Management Areas (HMAs). Permission to access these HMAs for study was given by the US Department of the Interior, Fish and Wildlife Service, Federal Fish and Wildlife Permit, permitted under Dr. Richard N. Conner.

Status of Cavity Trees Examined
Longleaf pines in the RCW clusters used for resin sampling were categorized according to their current activity status and method of cavity construction. The naturally excavated active cavity trees (NA) were those excavated and currently occupied by a RCW (n NA = 13). The artificial active cavity trees (AA) contained an artificial insert and was currently occupied by a RCW (n AA = 15). The artificial inactive cavity trees (AI) were pines with an insert that was newly constructed and either had been occupied or abandoned by a RCW (n AI = 16). The natural inactive cavity tree (NI) was a naturally excavated tree not currently used by a RCW (n NI = 16). The control tree (C) (n c = 24) was one that had external structural characteristics similar to a cavity tree including minimum age of 60 years; height, live crown ratio (LCR) and diameter at breast height (Conner and O'Halloran, 1987) but no evidence of RCW activity, defined as resin well pecking stimulating the flow of resin. Resin flow from 83 longleaf pine (one was lost to mortality) were evaluated. There were a minimum of two control trees sampled per cluster.

Resin Flow
Measurements on pines in the study included diameter at breast height (1.4 m); total height (m); height (m) to the lowest live limb; and bole length (m). Live crown ratio was computed as the percentage of the total height of the pine covered with living branches. Tree age was determined from past data (R. N. Conner, data on file). Rainfall for one month prior to each resin sampling event was obtained from the U. S. Department of Agriculture, Forest Service, National Forests and Grasslands in Texas, Lufkin, Texas.

Analysis
A repeated measures-split plot analysis of variance (ANOVA) was designed and then analyzed using SAS Version 8 software to compare means among sample trees. The ANOVA was used to test the null hypothesis of no differences among cavity-tree types with respect to resin flow, time and its effects on tree resin abilities, and the interaction of time and tree type with respect to resin flow with an alpha of 0.05 for testing. A repeated measures design was used as there are "essentially split plot designs allowing for correlation within each random effect" (Yandell, 1997).
Preliminary checks for normality and homogeneity of variances were performed prior to the ANOVA. The statistical package, JMP, was used to create distributions of resin per group (tree category × time period). Homogeneity was checked by graphing the average resin amount per group by the variance of the groups. Median resin yield instead of the mean was also graphed against time to locate outliers that hamper statistical analysis. The outliers were found above the 90 th percentile of the sample population and were considered important to the study due to high resin yields. Therefore, those trees were classified as "super resin producers" because they yielded a consistently higher amount of resin than other trees involved in the study.
A variance-stabilizing transformation was used as a step to correct heteroscedastcity of the variances. These transformations can aid in achieving the normality requirement in ANOVA. A square root transformation was used in the study because the plot of the sample averages against the sample variances resulted in a 45˚ linear trend. This transformed data was used in the resin differences. Regression analysis was used to determine differences between slopes of each of the cavity tree types sampled for resin. The regression equations also allow for general predictions for that year's resin between sample dates for the same study trees.

Resin Production
The distribution analysis showed a heavy skew with many outliers. The asymmetry indicated that the sample mean and the sample median did not coincide; therefore, departing from normality and suggesting the need to transform the data prior to analysis.
The statistical analysis showed no tree type effect on resin ( Table 1) (F = 1.58, df = 4, P = 0.1875), or any effect due to an interaction of tree type and time (F = 0.71, df = 32, P = 0.8849). Although not significantly different, the artificial cavities produced resin yields with the highest mean (AA) and the lowest mean (AI) ( Table  1). The control trees used in the study has a mean only slightly higher than AI cavities and slightly lower that NI, both of the inactive categories (Table 1). Table 2 reveals that the AA and C tree categories had more trees in the 90 th percentile and high resin producers in the month of August 1999. The other tree categories, NA, NI and AI each only had one tree that produced resin over 33.02 ml ( Table 2). Conversely, both active cavity tree categories, NA and AA, had the smallest number of cavity trees below the 33.02 threshold ( Table 2). Therefore, there are few trees present in the lower percentiles of resin production ( Table 3). There were more AI trees in the lower percentiles of resin production, while the NI was located more in the middle and the C category in the mid to upper percentiles ( Table 3). Table 1. Descriptive statistics for twenty-four hour resin yield (ml) sampled on red-cockaded woodpecker, longleaf pine cavity trees (natural active (NA), artificial active (AA), natural inactive (NI), artificial inactive (AI), and control (C)) sampled in the southern region of the Angelina National Forest, Texas, 1999-2000. There was no significant difference among the means (F = 1.58, df = 4, P = 0.1875).  The factor of time, however, did effect resin yield (F = 47.46, df = 8, P = 0.0001) through the sampling periods. Results indicate resin yield peaks in August and September (Table 4). It then tapers off during the late fall and through the early spring months. In both late summer measurements, the total rainfall is highest two months prior but very low one month prior to the sample date. It is during the beginning of the growing season, around April, that resin yield begins to increase again. The rainfall measurements between March and May show different trends for the two years involved. The amount of resin produced in the first year, 1999, started lower in March but peaked higher in the late summer that the year start of 2000 and last measured in September 2000 ( Table 4).
A test of homogeneity of slopes using regression analysis was performed on a yearly basis. The first year, March through August 1999, showed no difference in tree types (F = 1.924, df = 4, 0.10 < P < 0.25), thus all year 1 regression slopes can be called statistically equal. The second year, November 1999 through September 2000, however, showed a tree type effect (P < 0.0003) on the resin yield of at least one tree type. The differences were found between the active trees and the inactive (including control) trees ( Table 5). Those that had the same activity status showed no difference in resin yields. Table 6 is similar to a multiple comparisons for fixed   Table 6. Regression equations formulated with the 24-hour data (March 1999-August 1999) from the red-cockaded woodpecker resin-cavity (natural active (NA), artificial active (AA), natural inactive (NI), artificial inactive (AI), and control (C)) study in the southern region of the Angelina National Forest, Texas, 1999-2000. The slopes (F = 1.924, df = 4, 0.10 < P < 0.25 from these equations can be used to estimate the gain in resin per unit of time (days). effects, but it is a coefficient of determination used for random effects (Dowdy and Wearden, 1991), such as was time in this study. Any differences in tree type could only be determined in year two (November 1999 through September 2000) but testing in the two years together reduces the effect made by year 2000 tree sample differences.
General predictions for resin yield can be made for these study trees using the regression equations in Table 5 and Table 6. The values of R 2 exceed 0.8 and was as high as 0.9 (Table 5 and Table 6). These represent a strong correlation between the actual resin values and what the linear regression model predicted. Further general inferences can be made about the covariance structure for resin yields over time.

Super Rein Producer Analysis
Six trees were identified as consistent outliers and were labeled super resin producers. The six trees included one from each tree type category (NA, AA, NI, AI) and two from the control (C) ( Table 2, Table 7 and Table 8).
Removing the super producers did little to change the results (F = 1.76, df = 4, P = 0.1414) for tree type effect. The quantiletest performed on the allocated super resin producers had a grand overall median of 33.02 in the 90 th percentile in August 1999 ( Table 2). The special application of the chi-square test for homogeneity found strong evidence (P = 0.726) for accepting the null hypothesis that the August 1999 90 th percentile yields of resin were equal across all tree types ( Table 3). The median test resulted in no substantial evidence (P = 0.298) to reject the null hypothesis.
The majority of forest stand characteristics around the trees sampled differed significantly among treatments ( Table 9). The diameters (F = 2.71, P = 0.0359), live crown ratio (F = 2.66, P = 0.0386), and ages (F = 7.15, P = 0.0045) of the sample trees varied among the treatments but the total height (F = 0.45, P = 7722) did not. The differences found among the tree diameters and live crown ratio existed between the naturally excavated cavity trees and the inactive naturally excavated cavity trees. The large tree spacing in the cluster areas allows most trees to obtain higher dominance class as measured by live crown ratio.   Table 9. Diameter at breast height, total height, live crown ratio of longleaf sample trees by cavity type (natural active (NA), artificial active (AA), natural inactive (NI), artificial inactive (AI), and control (C)) study in the southern region of the Angelina National Forest, Texas, 1999-2000.

Discussion
The seasonal pattern of resin production observed coincides with that previously reported (Lorio, 1985;Ross et al., 1993Ross et al., , 1997. With the growth-differentiation balance at a low in the winter, there is little oleoresin flow from wounds while photosynthates accumulate in storage. The spring brings the start of new growth but little differentiation occurs because few vertical resin ducts are formed. Resin yield from wounds slowly begins during this time. The formation of summer latewood and the onset of water deficits results in differentiation that yield high resin flow. The shorter day lengths and cooler temperatures in the fall months show a halting of radial growth and a slowing of resin production and yields (Lorio, 1985). The seasonal dynamics in resin yield supports the statistical design output that there was no difference in the tree types or their interaction with time but that there was a difference in the effect of the season on resin.
With the seasonal differences (Table 4), biologists may track a tree at any time of the year to gauge its resin abilities before the decision is made for selection of potential trees for insert installation. Super rein producers occurred over the five types of sample trees indicating sampling from a homogeneous population. Depending on the limiting environmental factors, growth and differentiation processes become competitive sinks for photosynthates (Lorio, 1985). Wider spacing of trees can yield larger crowns, and faster growth until senescence is approached with tree age, and increased synthesis and flow of resin that produces conditions that favor the RCW but hinder the southern pine beetle, Dendroctonus frontalis Zimmermann, especially in longleaf pine.
Extreme moisture stress can affect the duration and quantity of resin flow decreasing the host tree defenses against bark beetles and also decreasing the success of the RCW with less resin for protection from predation by rat snakes. Adequate resin flow is affected by moisture stress reducing the flow and duration of resin (Lorio and Hodges, 1977). Droughts during the two years of this study and may be a contributing reason for less resin in 2000 ( Table 4, Table 8 and Table 9).
Each tree type showed no significant statistical variance in resin flow compared to other tree types even during periods of increased resin yield. Trees selected for artificial cavity installation for use by the RCW should have a minimum threshold for resin production if these cavities are to be protected from rat snakes by an adequate rein barrier. The AI category had the lowest mean ( Table 1, Table 8 and Table 9) for most sampling periods, lower than the C category. Conversely, the AA, artificial active trees had the highest mean ( Table 1, Table 8 and Table 9) and produce more resin in the upper percentiles ( Table 2 and Table 3). The active categories (NA, AA) separate from the inactive categories in the second year's spring resin yields (Table 5) with the actives holding the minimum 11 ml amount, the highest mean in the peak resin yielding month ( Table 1). Since the AA tree type mean was higher that the NA, it can be assumed that some trees chosen were very high resin producers, beneficial to the red-cockaded woodpecker. Increased observations are needed to see if high resin yields can negatively affect the RCW. The NI had a higher mean that the AI ( Table 1, Table 8 and Table 9) indicating inactive natural trees can potentially produce more resin that an inactive artificial tree. In the AI category, there was one tree above the 90 th percentile, but 46 % were in the 25 th percentile or below (Table 3).
There were droughts during the two years the study took place. The cumulative effect of two years of moisture stress related to drought may be the cause of the second year of resin measurement ( Table 8 and Table 9). The same increasing trend of resin production exists from late spring to late summer during 2000, but the extent of the increase in 2000 does not reach the peaks observed in 1999 ( Table 4). The start of the year 2000 resin production was higher than that of 1999. The overall trend of resin production needs to be sampled over a longer period of time to determine resin amounts over time. With guidance from this study and Ross et al. (1997) in the same area, sampling can indicate the best resin producing longleaf pine to select for artificial cavity insertion. Strom et al. (2002) investigating loblolly pine, (Pinus taeda), promoted identifying, producing and planting high resin yielding loblolly pine as a silvicultural approach to reduce damage from the southern pine beetle.
Although longleaf pine is not the primary host of the southern pine beetle, the amount and seasonal production of resin is a deterrent against the bark beetle. Reduced resin production in RCW cavity trees may make them vulnerable to bark beetle attack (Conner et al., 2001). The presence of resin volatiles around active cavity trees may be a factor in bark beetle induced morality in active cavity trees Conner et al., 1995;Hain et al., 2011). Active RCW cavity trees were more susceptible to bark beetle induced mortality than inactive cavity trees in loblolly, shortleaf (Pinus echinata) and longleaf pine Conner and Rudolph, 1995;Conner et al., 1998;Hain et al., 2011). RCW breeding males selected nest cavities in trees that were heavier resin producers than cavity trees used for roosting by other group members (Conner et al., 1998). Newer cavity trees were selected for the nest cavity in longleaf pine, but the RCW did not switch to the newest artificial cavity indicating sufficient resin was available (Conner et al., 1998). Results of our study can provide times to sample and indications of higher producing resin trees. Thresholds of resin production ( Table 3) indicate higher producing trees are available. For example, if 22 ml of resin in 24 hrs in deeded adequate, then this indicates the top 25 percent of the measured trees would be potential candidates for RCW cavi-ties. Mean resin production is this study was 16.4 ml in August 1999 and 13.3 ml in September 2000 (Table 4) and can be used as a guideline for artificial cavities. The months of July, August and September are used as reference months as they are typically the peak resin yielding months (Lorio, 1985;Bowman and Huh, 1995;Ross et al., 1997). Bowman and Huh (1995) found RCWs excavate cavities in pines with maximum crown-bole ratios. Use of the regression models gives guidance to resin production ( Table 5 and Table 6) over time and these can be used depending on resin measurements. RCW management is aimed at mimicking natural processes that shape the habitat for the bird.

Conclusion
Although no significant differences in resin yield were detected among cavity tree types, some pines produced more resin and were termed super producers. Pines chosen for cavity insert installation should be those pines that produce yields above the overall mean, but closer to the optimum mean of the AA and NA categories ( Table 1). Any resin production by a pine that is closer to 11 ml in the optimum resin producing months of July, August and September would have less chance for placing the RCW at risk for snake predation, or the host tree for bark beetle attack. It is not known if the RCW chooses super resin producers for nesting or roosting as super resin producers were found in each cavity tree type indicating not all super resin producers were colonized by RCW.
Mean resin production in this study was 16.34 ml (standard deviation of 9.53) in August 1999, and 13.33 ml (standard deviation of 4.99) in September 2000. These mean values could be used as a guideline to select the appropriate artificial cavity tree in the southern Angelina National Forest. The mean, however, will fluctuate according to moisture stress. The months of July, August and September can be used as reference months for cavity selection based on resin production. Regression equations were developed to predict resin production later in the season from early season sampling.