Panama City Fisheries Resources Office: FY 2008 Annual Report

HIGHLIGHTS FOR FY 2008 1. Completed the first of a two-year Gulf sturgeon population study on the Choctawhatchee River, Florida. The sub adult and adult Gulf sturgeon population was estimated at 2,800 fish. 2. Gulf sturgeon eggs were collected at three hard bottom sites in the Apalachicola River, Florida; two sites were previously confirmed spawning areas and one was a newly confirmed spawning area. 3. Documented 55 potential environmental threats to Gulf sturgeon spawning habitat in the Pea River, Florida and Alabama. 4. Assigned the Eglin AFB Road-Stream Crossing Working Group to guide the closure, repair and maintenance of roads and road stream crossings that impact threatened and endangered species. 5. Conducted 81 assessments of fish and stream invertebrates on and in watersheds surrounding Eglin AFB. 6. Provided technical support for the 5-year status review and reclassification proposed rule for the Okaloosa darter. 7. Initiated an intensive population genetic analysis of the Okaloosa darter throughout its range. Tissues from over 200 Okaloosa darters were collected and analyzed. 8. Established a GIS database to serve as a host for data from any sites sampled for mussels in Northeast Gulf of Mexico drainages. 9. Conducted habitat surveys at 115 locations in the Apalachicola River to assess the effects of drought-related mussel mortality and strandings, evaluate habitat conditions, and assess population demography. 10. A land use/aerial imagery threats assessment data analysis was completed for the Chipola River. A total of 266 impoundments/borrow pits and 471 unpaved road crossings were identified among the threats. 11. Okaloosa darters marked with elastomeric dyes were monitored in Mill Creek, Eglin AFB, to determine movement and habitat use following completion of a fish passage project. 3 12. Partners for Fish and Wildlife funded a streambank and riparian restoration project on Econfina Creek consisting of 3,900 feet of streambank fencing to exclude cattle access. One acre of riparian floodplain was planted with native trees. 13. We provided design and on-the-ground assistance for restoring surface hydrology at St. Vincent NWR. The project restored approximately 1.5 miles of tidal stream and 100 acres of wetlands. 14. A study was completed on 11 coastal streams to document large wood debris relationships with fluvial geomorphic characteristics. 15. We developed a Population Viability Analysis model for the fat threeridge mussel to determine current and future risk of extinction. 17. A Gulf Sturgeon Friends Group, “Gulf Sturgeon Preservation Society” was organized in FY 08. 18. Multiple outreach projects were completed to detail aquatic resource conservation needs and opportunities, including National Fishing Week, Earth Day, several festivals and school outreach.

Limited resources? Let the habitat priority planner help!

In this time of scarce resources, coastal resource managers must find ways to prioritize conservation, land use, and restoration efforts. The Habitat Priority Planner (HPP) is a free geospatial tool created by the National Oceanic and Atmospheric Administration’s Coastal Services Center that has received wide praise for its ease of use and broad applicability to conservation strategic planning, restoration, climate change scenarios, and other natural resource management actions. Not a geographic information system (GIS) user? Don’t worry―this tool was designed to be used in a team setting. One intermediate-level GIS user can push the buttons to show quick results while a roomful of resource managers and stakeholders provide input criteria that determine the results. The Habitat Priority Planner is a toolbar for ESRI’s ArcGIS platform that is composed of three modules: Habitat Classification, Habitat Analysis, and Data Explorer. The tool calculates basic ecological statistics that are used to examine how habitats function within a landscape. The tool pre‐packages several common landscape metrics into a user‐friendly interface for intermediate GIS users. In addition, HPP allows the user to build queries interactively using a graphical interface for demonstrating criteria selections quickly in a visual manner that is useful in stakeholder interactions. Tool advocates and users include land trusts, conservation alliances, nonprofit organizations, and select National Estuarine Research Reserves and refuges of the U.S. Fish and Wildlife Service. Participants in this session will learn the basic requirements for HPP use and the multiple ways the HPP has been applied to geographies nationwide. (PDF contains 5 pages)

Predicting the effects of climate change on sea turtle nesting habitat in Florida

Rising global temperatures threaten the survival of many plant and animal species. Having already risen at an unprecedented rate in the past century, temperatures are predicted to rise between 0.3 and 7.5C in North America over the next 100 years (Hawkes et al. 2007). Studies have documented the effects of climate warming on phenology (timing of seasonal activities), with observations of early arrival at breeding grounds, earlier ends to the reproductive season, and delayed autumnal migrations (Pike et al. 2006). In addition, for species not suited to the physiological demands of cold winter temperatures, increasing temperatures could shift tolerable habitats to higher latitudes (Hawkes et al. 2007). More directly, climate warming will impact thermally sensitive species like sea turtles, who exhibit temperature-dependent sexual determination. Temperatures in the middle third of the incubation period determine the sex of sea turtle offspring, with higher temperatures resulting in a greater abundance of female offspring. Consequently, increasing temperatures from climate warming would drastically change the offspring sex ratio (Hawkes et al. 2007). Of the seven extant species of sea turtles, three (leatherback, Kemp’s ridley, and hawksbill) are critically endangered, two (olive ridley and green) are endangered, and one (loggerhead) is threatened. Considering the predicted scenarios of climate warming and the already tenuous status of sea turtle populations, it is essential that efforts are made to understand how increasing temperatures may affect sea turtle populations and how these species might adapt in the face of such changes. In this analysis, I seek to identify the impact of changing climate conditions over the next 50 years on the availability of sea turtle nesting habitat in Florida given predicted changes in temperature and precipitation. I predict that future conditions in Florida will be less suitable for sea turtle nesting during the historic nesting season. This may imply that sea turtles will nest at a different time of year, in more northern latitudes, to a lesser extent, or possibly not at all. It seems likely that changes in temperature and precipitation patterns will alter the distribution of sea turtle nesting locations worldwide, provided that beaches where the conditions are suitable for nesting still exist. Hijmans and Graham (2006) evaluate a range of climate envelope models in terms of their ability to predict species distributions under climate change scenarios. Their results suggested that the choice of species distribution model is dependent on the specifics of each individual study. Fuller et al. (2008) used a maximum entropy approach to model the potential distribution of 11 species in the Arctic Coastal Plain of Alaska under a series of projected climate scenarios. Recently, Pike (in press) developed Maxent models to investigate the impacts of climate change on green sea turtle nest distribution and timing. In each of these studies, a set of environmental predictor variables (including climate variables), for which ‘current’ conditions are available and ‘future’ conditions have been projected, is used in conjunction with species occurrence data to map potential species distribution under the projected conditions. In this study, I will take a similar approach in mapping the potential sea turtle nesting habitat in Florida by developing a Maxent model based on environmental and climate data and projecting the model for future climate data. (PDF contains 5 pages)

An analysis of Nigerian freshwater fishes: those under threat and conservation options

The study assessed qualitatively the threat status of all nigerian freshwater fishes using such criteria as rarity, size at maturity, mode of reproduction, human population density, habitat degradation, pollution and range of each species among others. The biology of 48% (129n) of nigerian freshwater species is not well known. Of the 266 known freshwater fishes, 47 species represented 17% are critically endangered, 15 (5%) are endangered , 8(3%), are vulnerable while 23(8%) are near threatened. The paper suggests increased basic knowledge of threatened species and conservation policy along three lines public awareness, legislation and creation of national parks, aquaria and reserves as measures needed to ensure the conservation of the fishes

I. Impact spallation experiments: fracture patterns and spall velocities. II. Craters in carbonate rocks: an electron paramagnetic resonance analysis of shock damage

This work is divided into two independent papers.

PAPER 1.

Spall velocities were measured for nine experimental impacts into San Marcos gabbro targets. Impact velocities ranged from 1 to 6.5 km/sec. Projectiles were iron, aluminum, lead, and basalt of varying sizes. The projectile masses ranged from a 4 g lead bullet to a 0.04 g aluminum sphere. The velocities of fragments were measured from high-speed films taken of the events. The maximum spall velocity observed was 30 m/sec, or 0.56 percent of the 5.4 km/sec impact velocity. The measured velocities were compared to the spall velocities predicted by the spallation model of Melosh (1984). The compatibility between the spallation model for large planetary impacts and the results of these small scale experiments are considered in detail.

The targets were also bisected to observe the pattern of internal fractures. A series of fractures were observed, whose location coincided with the boundary between rock subjected to the peak shock compression and a theoretical "near surface zone" predicted by the spallation model. Thus, between this boundary and the free surface, the target material should receive reduced levels of compressive stress as compared to the more highly shocked region below.

PAPER 2.

Carbonate samples from the nuclear explosion crater, OAK, and a terrestrial impact crater, Meteor Crater, were analyzed for shock damage using electron para- magnetic resonance, EPR. The first series of samples for OAK Crater were obtained from six boreholes within the crater, and the second series were ejecta samples recovered from the crater floor. The degree of shock damage in the carbonate material was assessed by comparing the sample spectra to spectra of Solenhofen limestone, which had been shocked to known pressures.

The results of the OAK borehole analysis have identified a thin zone of highly shocked carbonate material underneath the crater floor. This zone has a maximum depth of approximately 200 ft below sea floor at the ground zero borehole and decreases in depth towards the crater rim. A layer of highly shocked material is also found on the surface in the vicinity of the reference bolehole, located outside the crater. This material could represent a fallout layer. The ejecta samples have experienced a range of shock pressures.

It was also demonstrated that the EPR technique is feasible for the study of terrestrial impact craters formed in carbonate bedrock. The results for the Meteor Crater analysis suggest a slight degree of shock damage present in the β member of the Kaibab Formation exposed in the crater walls.

A new ground motion intensity measure, peak filtered acceleration (PFA), to estimate collapse vulnerability of buildings in earthquakes

In this thesis, we develop an efficient collapse prediction model, the PFA (Peak Filtered Acceleration) model, for buildings subjected to different types of ground motions.

For the structural system, the PFA model covers modern steel and reinforced concrete moment-resisting frame buildings (potentially reinforced concrete shear wall buildings). For ground motions, the PFA model covers ramp-pulse-like ground motions, long-period ground motions, and short-period ground motions.

To predict whether a building will collapse in response to a given ground motion, we first extract long-period components from the ground motion using a Butterworth low-pass filter with suggested order and cutoff frequency. The order depends on the type of ground motion, and the cutoff frequency depends on the building’s natural frequency and ductility. We then compare the filtered acceleration time history with the capacity of the building. The capacity of the building is a constant for 2-dimentional buildings and a limit domain for 3-dimentional buildings. If the filtered acceleration exceeds the building’s capacity, the building is predicted to collapse. Otherwise, it is expected to survive the ground motion.

The parameters used in PFA model, which include fundamental period, global ductility and lateral capacity, can be obtained either from numerical analysis or interpolation based on the reference building system proposed in this thesis.

The PFA collapse prediction model greatly reduces computational complexity while archiving good accuracy. It is verified by FEM simulations of 13 frame building models and 150 ground motion records.

Based on the developed collapse prediction model, we propose to use PFA (Peak Filtered Acceleration) as a new ground motion intensity measure for collapse prediction. We compare PFA with traditional intensity measures PGA, PGV, PGD, and Sa in collapse prediction and find that PFA has the best performance among all the intensity measures.

We also provide a close form in term of a vector intensity measure (PGV, PGD) of the PFA collapse prediction model for practical collapse risk assessment.

Three-dimensional nonlinear inelastic analysis of steel moment-frame buildings damaged by earthquake excitations

The Northridge earthquake of January 17, 1994, highlighted the two previously known problems of premature fracturing of connections and the damaging capabilities of near-source ground motion pulses. Large ground motions had not been experienced in a city with tall steel moment-frame buildings before. Some steel buildings exhibited fracture of welded connections or other types of structural degradation.

A sophisticated three-dimensional nonlinear inelastic program is developed that can accurately model many nonlinear properties commonly ignored or approximated in other programs. The program can assess and predict severely inelastic response of steel buildings due to strong ground motions, including collapse.

Three-dimensional fiber and segment discretization of elements is presented in this work. This element and its two-dimensional counterpart are capable of modeling various geometric and material nonlinearities such as moment amplification, spread of plasticity and connection fracture. In addition to introducing a three-dimensional element discretization, this work presents three-dimensional constraints that limit the number of equations required to solve various three-dimensional problems consisting of intersecting planar frames.

Two buildings damaged in the Northridge earthquake are investigated to verify the ability of the program to match the level of response and the extent and location of damage measured. The program is used to predict response of larger near-source ground motions using the properties determined from the matched response.

A third building is studied to assess three-dimensional effects on a realistic irregular building in the inelastic range of response considering earthquake directivity. Damage levels are observed to be significantly affected by directivity and torsional response.

Several strong recorded ground motions clearly exceed code-based levels. Properly designed buildings can have drifts exceeding code specified levels due to these ground motions. The strongest ground motions caused collapse if fracture was included in the model. Near-source ground displacement pulses can cause columns to yield prior to weaker-designed beams. Damage in tall buildings correlates better with peak-to-peak displacements than with peak-to-peak accelerations.

Dynamic response of tall buildings shows that higher mode response can cause more damage than first mode response. Leaking of energy between modes in conjunction with damage can cause torsional behavior that is not anticipated.

Various response parameters are used for all three buildings to determine what correlations can be made for inelastic building response. Damage levels can be dramatically different based on the inelastic model used. Damage does not correlate well with several common response parameters.

Realistic modeling of material properties and structural behavior is of great value for understanding the performance of tall buildings due to earthquake excitations.

Habitat use by 0+ cyprinid fish in the River Great Ouse, East Anglia

This study was designed to examine the habitat use of several species of 0+ cyprinid in the regulated River Great Ouse and to determine the reasons for specific habitat use. In general, all fish species were found associated with the marginal zone, with little diel variation. Use of shallow habitats in the presence of macrophytes correlated well with the distribution of zooplankton in the river channel, the preferred food source of 0+ cyprinids. During the early to late larval phase, all species fed upon rotifers and diatoms. Cladocera, particularly Alona spp. and Chydorus spp., and early instar larvae of Chironomidae, then became prevalent in the diet along with small numbers of Copepoda. Models were developed to determine habitat availability over a range of discharges, using the physical habitat simulation (PHABSIM) component of the Instream Flow Incremental Methodology (IFIM). The results of this analysis revealed that habitat suitable for 0+ fishes comprised a relatively small percentage of the main channel and generally decreased with discharge.

Produção in vitro e criopreservação de raízes de Cleome rosea Vahl (Capparaceae)

Cleome rosea é uma espécie nativa, de porte herbáceo, ocorrente em restingas brasileiras. Estudos recentes têm revelado o potencial medicinal da espécie para importantes propriedades farmacológicas, como por exemplo, as atividades anti-inflamatória, antigenotóxica, antiviral e antibacteriana. Porém, nos últimos anos, C. rosea não tem sido encontrada em várias regiões de seu ambiente natural, devido, principalmente, às ações antrópicas. Dessa forma, torna-se relevante o desenvolvimento de métodos de conservação que permitam o estudo e exploração das propriedades medicinais da espécie. O cultivo in vitro de raízes representa uma forma eficiente para produção de biomassa, devido ao rápido crescimento, produção estável de metabólitos, além de representar uma potencial fonte de explantes para a propagação em massa de diferentes espécies. O presente trabalho teve como objetivo a produção in vitro de culturas de raízes de C. rosea, associada à criopreservação, como forma de manutenção em longo prazo das culturas, monitorada através da análise de estabilidade genética. As culturas estabelecidas a partir de explantes radiculares de plantas propagadas in vitro de C. rosea demonstraram excelente capacidade de multiplicação de raízes em meio de cultura suplementado com o fitorregulador ANA, com manutenção dessa capacidade ao longo de sucessivas subculturas. Associado a esses resultados, o estabelecimento de protocolos de criopreservação pelo método de vitrificação resultou em elevados valores de frequência de recuperação do material após congelamento em nitrogênio líquido com as soluções de vitrificação PVS2 e PVS3. Os estudos de monitoramento da estabilidade genética, pela técnica de marcadores moleculares RAPD, revelaram a presença de polimorfismos significativos em uma das três culturas iniciadas a partir de raízes de C. rosea criopreservadas. Esses resultados demonstram as possibilidades de produção de raízes de C. rosea e conservação em longo prazo através da criopreservação, iniciando estudos inéditos para a espécie.

Lock Acquisition and Sensitivity Analysis of Advanced LIGO Interferometers

Laser interferometer gravitational wave observatory (LIGO) consists of two complex large-scale laser interferometers designed for direct detection of gravitational waves from distant astrophysical sources in the frequency range 10Hz - 5kHz. Direct detection of space-time ripples will support Einstein's general theory of relativity and provide invaluable information and new insight into physics of the Universe.

Initial phase of LIGO started in 2002, and since then data was collected during six science runs. Instrument sensitivity was improving from run to run due to the effort of commissioning team. Initial LIGO has reached designed sensitivity during the last science run, which ended in October 2010.

In parallel with commissioning and data analysis with the initial detector, LIGO group worked on research and development of the next generation detectors. Major instrument upgrade from initial to advanced LIGO started in 2010 and lasted till 2014.

This thesis describes results of commissioning work done at LIGO Livingston site from 2013 until 2015 in parallel with and after the installation of the instrument. This thesis also discusses new techniques and tools developed at the 40m prototype including adaptive filtering, estimation of quantization noise in digital filters and design of isolation kits for ground seismometers.

The first part of this thesis is devoted to the description of methods for bringing interferometer to the linear regime when collection of data becomes possible. States of longitudinal and angular controls of interferometer degrees of freedom during lock acquisition process and in low noise configuration are discussed in details.

Once interferometer is locked and transitioned to low noise regime, instrument produces astrophysics data that should be calibrated to units of meters or strain. The second part of this thesis describes online calibration technique set up in both observatories to monitor the quality of the collected data in real time. Sensitivity analysis was done to understand and eliminate noise sources of the instrument.

Coupling of noise sources to gravitational wave channel can be reduced if robust feedforward and optimal feedback control loops are implemented. The last part of this thesis describes static and adaptive feedforward noise cancellation techniques applied to Advanced LIGO interferometers and tested at the 40m prototype. Applications of optimal time domain feedback control techniques and estimators to aLIGO control loops are also discussed.

Commissioning work is still ongoing at the sites. First science run of advanced LIGO is planned for September 2015 and will last for 3-4 months. This run will be followed by a set of small instrument upgrades that will be installed on a time scale of few months. Second science run will start in spring 2016 and last for about 6 months. Since current sensitivity of advanced LIGO is already more than factor of 3 higher compared to initial detectors and keeps improving on a monthly basis, upcoming science runs have a good chance for the first direct detection of gravitational waves.

Studies of response to earthquake ground motion

A study is made of the accuracy of electronic digital computer calculations of ground displacement and response spectra from strong-motion earthquake accelerograms. This involves an investigation of methods of the preparatory reduction of accelerograms into a form useful for the digital computation and of the accuracy of subsequent digital calculations. Various checks are made for both the ground displacement and response spectra results, and it is concluded that the main errors are those involved in digitizing the original record. Differences resulting from various investigators digitizing the same experimental record may become as large as 100% of the maximum computed ground displacements. The spread of the results of ground displacement calculations is greater than that of the response spectra calculations. Standardized methods of adjustment and calculation are recommended, to minimize such errors.

Studies are made of the spread of response spectral values about their mean. The distribution is investigated experimentally by Monte Carlo techniques using an electric analog system with white noise excitation, and histograms are presented indicating the dependence of the distribution on the damping and period of the structure. Approximate distributions are obtained analytically by confirming and extending existing results with accurate digital computer calculations. A comparison of the experimental and analytical approaches indicates good agreement for low damping values where the approximations are valid. A family of distribution curves to be used in conjunction with existing average spectra is presented. The combination of analog and digital computations used with Monte Carlo techniques is a promising approach to the statistical problems of earthquake engineering.

Methods of analysis of very small earthquake ground motion records obtained simultaneously at different sites are discussed. The advantages of Fourier spectrum analysis for certain types of studies and methods of calculation of Fourier spectra are presented. The digitizing and analysis of several earthquake records is described and checks are made of the dependence of results on digitizing procedure, earthquake duration and integration step length. Possible dangers of a direct ratio comparison of Fourier spectra curves are pointed out and the necessity for some type of smoothing procedure before comparison is established. A standard method of analysis for the study of comparative ground motion at different sites is recommended.

I. Nuclear magnetic resonance analysis of ferrocenylcarbonium ion. II. Ferrocene catalyzed photochemistry

I. Nuclear magnetic resonance spectra of appropriately substituted ferrocenylcarbonium ions reveal the α-protons of the substituted ring to be more shielded than β-protons. The observation is discussed in terms of various models proposed for the ferrocenylcarbonium ion and is found to support a model in which the iron is bonded to all six carbona of the substituted ring.

II. Ferrocene catalyzes the photoisomerization of the piperylenes and the photodimerization of isoprene. Our results suggest a mechanism in which a complex of ferrocene and diene is excited to its second singlet state which dissociates to a triplet-state ferrocene molecule and a triplet-state diene molecule. The triplet-state diene, then, proceeds to isomerize or attack ground-state diene to form dimers.

Investigation of strong earthquake ground motion

The pattern of energy release during the Imperial Valley, California, earthquake of 1940 is studied by analysing the El Centro strong motion seismograph record and records from the Tinemaha seismograph station, 546 km from the epicenter. The earthquake was a multiple event sequence with at least 4 events recorded at El Centro in the first 25 seconds, followed by 9 events recorded in the next 5 minutes. Clear P, S and surface waves were observed on the strong motion record. Although the main part of the earthquake energy was released during the first 15 seconds, some of the later events were as large as M = 5.8 and thus are important for earthquake engineering studies. The moment calculated using Fourier analysis of surface waves agrees with the moment estimated from field measurements of fault offset after the earthquake. The earthquake engineering significance of the complex pattern of energy release is discussed. It is concluded that a cumulative increase in amplitudes of building vibration resulting from the present sequence of shocks would be significant only for structures with relatively long natural period of vibration. However, progressive weakening effects may also lead to greater damage for multiple event earthquakes.

The model with surface Love waves propagating through a single layer as a surface wave guide is studied. It is expected that the derived properties for this simple model illustrate well several phenomena associated with strong earthquake ground motion. First, it is shown that a surface layer, or several layers, will cause the main part of the high frequency energy, radiated from the nearby earthquake, to be confined to the layer as a wave guide. The existence of the surface layer will thus increase the rate of the energy transfer into the man-made structures on or near the surface of the layer. Secondly, the surface amplitude of the guided SH waves will decrease if the energy of the wave is essentially confined to the layer and if the wave propagates towards an increasing layer thickness. It is also shown that the constructive interference of SH waves will cause the zeroes and the peaks in the Fourier amplitude spectrum of the surface ground motion to be continuously displaced towards the longer periods as the distance from the source of the energy release increases.

Influence of local geology on earthquake ground motion

As a simplified approach for estimating theoretically the influence of local subsoils upon the ground motion during an earthquake, the problem of an idealized layered system subjected to vertically incident plane body waves was studied. Both the technique of steady-state analysis and the technique of transient analysis have been used to analyze the problem.

In the steady-state analysis, a recursion formula has been derived for obtaining the response of a layered system to sinusoidally steady-state input. Several conclusions are drawn concerning the nature of the amplification spectrum of a nonviscous layered system having its layer stiffnesses increasing with depth. Numerical examples are given to demonstrate the effect of layer parameters on the amplification spectrum of a layered system.

In the transient analysis, two modified shear beam models have been established for obtaining approximately the response of a layered system to earthquake-like excitation. The method of continuous modal analysis was adopted for approximate analysis of the models, with energy dissipation in the layers, if any, taken into account. Numerical examples are given to demonstrate the accuracy of the models and the effect of a layered system in modifying the input motion.

Conditions are established, under which the theory is applicable to predict the influence of local subsoils on the ground motion during an earthquake. To demonstrate the applicability of the models to actual cases, three examples of actually recorded earthquake events are examined. It is concluded that significant modification of the incoming seismic waves, as predicted by the theory, is likely to occur in well defined soft subsoils during an earthquake, provided that certain conditions concerning the nature of the incoming seismic waves are satisfied.

Landscape Genetics for the Empirical Assessment of Resistance Surfaces: The European Pine Marten (Martes martes) as a Target-Species of a Regional Ecological Network

Coherent ecological networks (EN) composed of core areas linked by ecological corridors are being developed worldwide with the goal of promoting landscape connectivity and biodiversity conservation. However, empirical assessment of the performance of EN designs is critical to evaluate the utility of these networks to mitigate effects of habitat loss and fragmentation. Landscape genetics provides a particularly valuable framework to address the question of functional connectivity by providing a direct means to investigate the effects of landscape structure on gene flow. The goals of this study are (1) to evaluate the landscape features that drive gene flow of an EN target species (European pine marten), and (2) evaluate the optimality of a regional EN design in providing connectivity for this species within the Basque Country (North Spain). Using partial Mantel tests in a reciprocal causal modeling framework we competed 59 alternative models, including isolation by distance and the regional EN. Our analysis indicated that the regional EN was among the most supported resistance models for the pine marten, but was not the best supported model. Gene flow of pine marten in northern Spain is facilitated by natural vegetation, and is resisted by anthropogenic landcover types and roads. Our results suggest that the regional EN design being implemented in the Basque Country will effectively facilitate gene flow of forest dwelling species at regional scale.