Module 8: Conservation — DDT, Lead, and Recovery

The mid-twentieth-century collapse of bald-eagle, peregrine and osprey populations across North America and Europe provided the single most important public demonstration of ecological toxicology. Dichloro-diphenyl-trichloroethane (DDT), via its metabolite p,p’-DDE, was identified by Ratcliffe (1967) as the causal agent of eggshell thinning through inhibition of uterine Ca-ATPase transport. The subsequent ban (USA 1972, most of Europe 1970s–1980s) yielded remarkable recoveries — bald eagles delisted in 2007, peregrines recovered via captive breeding (Cade 1982). The current-generation threats are lead ammunition (inhibiting ALA-D through Pb²⁺ displacement of Zn²⁺) and, for Old-World Gyps vultures, the veterinary NSAID diclofenac. We develop the quantitative toxicology, the demographic consequences, and the policy responses.

1. The DDT Era: Biomagnification and Eggshell Thinning

DDT was first synthesised by Othmar Zeidler in 1874; its insecticidal properties were discovered by Paul Hermann Müller in 1939 (Nobel Prize 1948). Deployment began 1942 (malaria and typhus control during WWII), and agricultural application exploded postwar. Production peaked at ~8 × 10⁴ tonnes/year by 1963. Rachel Carson’s Silent Spring (1962) catalysed public and regulatory concern; the US Environmental Protection Agency banned agricultural use in 1972.

Bioaccumulation and biomagnification

DDT is lipid-soluble (log \( K_{\text{ow}} \) ≈ 6.9) and persistent (environmental half-life ~2–15 years depending on substrate). At each trophic step, lipid-normalised concentration amplifies by a biomagnification factor (BMF):

\[ C_{n+1} \;=\; \text{BMF}_{n \to n+1} \cdot C_n, \qquad C_\text{top} \;\approx\; C_\text{water} \cdot \prod_{i} \text{BMF}_i \]

For a typical four-step food-web (water → plankton → fish → piscivorous bird), the total amplification reaches 10⁵–10⁷. A water concentration of 10 ng/L thus yields 1–100 ppm in eagle adipose tissue.

Mechanism: Ca-ATPase inhibition

DDE binds the oviductal Ca-ATPase with \( K_i \approx 3 \) ppm wet-weight adipose equivalent, reducing transported calcium to the shell-gland lumen in proportion. Because shell thickness is determined by the time-integral of Ca-flux over the ~19-hour calcification window, DDE produces proportional thinning:

\[ \frac{L_\text{shell}(\text{DDE})}{L_\text{shell}(0)} \;\approx\; \frac{1}{1 + [\text{DDE}]/K_i} \]

Ratcliffe (1967, 1970), measuring museum egg collections, documented a stepped decline of 15–25 % in British peregrine, sparrowhawk and golden-eagle shell thickness from 1946 onwards. Similar declines were reported for North American bald eagles, osprey, and peregrines. Thin shells crack under incubating-parent weight, and reproductive success collapses.

Bald eagle crash and recovery

Pre-DDT bald-eagle populations in the contiguous USA are conservatively estimated at 100 000 breeding pairs (pre-1940). By 1963 the population had crashed to a documented 417 pairs (USFWS 1983). Following the 1972 ban and decades of protection under the Endangered Species Act, population recovery proceeded: 500 pairs by 1970, 2 500 by 1980, 12 000 by 2000, and >70 000 by 2020. Bald eagle was delisted from the US Endangered Species list in June 2007 — one of the most successful ESA recoveries on record.

DDT biomagnification through a lake food web

2. Peregrine Recovery and the Model of Captive Breeding

Peregrine falcon (Falco peregrinus) eastern-North-American populations were extirpated by 1964 (Hickey 1969). Tom Cade and the Peregrine Fund initiated a captive-breeding and reintroduction programme in 1970, operating from the Cornell hack-box facility. By 1999 the peregrine had been removed from the US Endangered Species list. The programme pioneered techniques widely used subsequently:

Hacking: Young are released from elevated boxes, fed in absentia through a pipe, and achieve wild flight in situ — avoiding human-imprinting.
Cross-fostering: Eggs from captive pairs are swapped into active wild nests of conspecifics or closely related species.
Double-clutching: Induced egg- replacement elicits a second clutch from a captive pair, doubling output.
Puppet-feeding: Life-sized parental-model puppets feed chicks to avoid mis-imprinting on human handlers, critical for later-used condor and Spanish imperial eagle programmes.

By 2020, US peregrine population exceeded 3 000 breeding pairs and the species is commonly encountered in major cities (urban ledge-nesting).

3. Lead Ammunition: ALA-D Inhibition and the Condor Case

Scavenging raptors ingest lead shot fragments or bullet jacketing when feeding on carrion from hunter-killed ungulates (deer, wild pigs) left in the field as “gut piles”. A single 0.30-06 bullet can leave 15–50 fragments per carcass after expansion. Bald eagles, golden eagles, turkey vultures, and California condors are all affected; condors are uniquely vulnerable because of their obligate scavenging and ~9 kg body mass (small lead doses relative to body).

Mechanism: ALA-D inhibition

δ-aminolevulinic acid dehydratase (ALA-D, EC 4.2.1.24) catalyses the condensation of two molecules of δ-ALA to porphobilinogen in heme biosynthesis. The enzyme is a zinc metalloprotein; Pb²⁺ competitively displaces Zn²⁺ at the catalytic site with half-maximal inhibition near 0.3 μM Pb in erythrocyte assays. Pain (2019) reviewed the avian toxicology literature: clinical signs appear at blood-Pb > 60 μg/dL, lethal blood-Pb > 100 μg/dL. Subclinical effects — impaired cognition, reduced hunting success — are detectable below 20 μg/dL.

\[ \frac{v(\text{Pb})}{v_0} \;=\; \frac{1}{1 + ([\text{Pb}]/K_i)^n}, \qquad K_i \approx 0.3\ \mu\text{M},\ n \approx 1.8 \]

The downstream consequences are accumulation of ALA (neurotoxic), decreased heme synthesis (anaemia), and eventual multi-organ failure.

Blood-Pb pharmacokinetics

A lead shot fragment lodged in the proventriculus dissolves at the gastric pH (~2–3) at a rate of ~15 mg/day. Blood-Pb half-life is ~13 days; bone-Pb half-life is 5–20 years, providing a long-term depot that re- enters circulation during calcium stress (laying, incubation). Finkelstein et al. (2012) published the definitive California condor blood-Pb survey: more than half of wild-sampled condors exceeded 20 μg/dL, and a third exceeded 60 μg/dL during deer-hunting seasons.

California condor reintroduction

The last wild California condor was captured in 1987; total population stood at 27 individuals, all in captivity. The San Diego Zoo and Los Angeles Zoo implemented captive breeding, with first releases in 1992. By 2020, ~500 condors existed, with roughly half wild-flying across California, Arizona, Utah, and Baja California. However, the population is not self-sustaining without continued management: annual lead exposure rates remain high, chelation therapy is required for many wild birds, and captive releases continue to offset ongoing mortality.

California AB 711 (2013) and federal context

California Assembly Bill 711 (signed 2013, fully phased in 2019) prohibits lead-based ammunition for taking any wildlife in California. US federal lands adopted variable restrictions. The effect is real but partial: compliance is imperfect, and imported ammunition continues to spread fragments. Sustainable condor populations in the long term require broad regional shift to copper- based ammunition.

4. Diclofenac and the Collapse of South Asian Vultures

Three South Asian Gyps species (G. bengalensis, G. indicus, and G. tenuirostris) collapsed by >95 % between 1993 and 2007. Oaks et al. (2004) identified the veterinary non-steroidal anti-inflammatory drug diclofenac as the causal agent: vultures feeding on cattle carcasses treated shortly before death received a lethal renal dose (visceral gout from uric-acid deposition on kidney, liver, spleen). The LD₅₀ is ~0.1–0.2 mg/kg, far below typical carcass residues.

India, Pakistan, and Nepal banned veterinary diclofenac in 2006 and substituted meloxicam, which is vulture-safe. Recovery has been slow; populations remain >95 % below historic baseline. The incident is the textbook demonstration of how a veterinary pharmaceutical can devastate a wild apex scavenger guild and cascade through ecosystem function: carcass-disposal collapse has been linked to the rise of feral-dog populations and increases in human rabies incidence across the subcontinent (Markandya et al. 2008).

Eagle parallels

While diclofenac has not been directly implicated in eagle collapses at the same scale, veterinary NSAIDs remain a concern for facultatively scavenging eagles (Aquila chrysaetos, A. heliaca, Haliaeetus albicilla). The European white-tailed eagle population has shown scattered diclofenac- attributable mortality following Spanish approval of veterinary diclofenac in 2013, later restricted.

5. Collision and Electrocution

Modern-era anthropogenic mortality sources for eagles:

Wind-turbine strikes

At California’s Altamont Pass Wind Resource Area (5 400 turbines at peak), golden-eagle mortality was ~60–100 adults/year (Smallwood & Thelander 2008). Subsequent repowering with fewer, larger turbines reduced fatality rates per-megawatt but did not eliminate mortality. Rotor-tip velocities exceed 80 m/s, making tip strikes effectively unavoidable once an eagle is committed to the swept disk. Because golden-eagle populations are buffered by floater recruits (Module 7), sustained mortality at Altamont was masked for years before the pool depleted and territory occupancy fell.

Power-line electrocution

Electrocution on medium-voltage distribution lines (most hazardous: 12–33 kV, where conductor spacing is shorter than typical eagle wingspan of 2 m) is a major cause of mortality for perching raptors. Mitigation by insulation of live parts on poles (Avian Power Line Interaction Committee 2006) reduces mortality by >90 % on retrofitted infrastructure.

Secondary rodenticide poisoning

Second-generation anticoagulant rodenticides (SGARs: brodifacoum, bromadiolone, difenacoum) persist in the liver of poisoned rodents for weeks; predators and scavengers accumulate sublethal and occasionally lethal doses. Golden eagles, barn owls, and great-horned owls are all documented affected. California and the EU have tightened SGAR restrictions through the 2010s.

6. Reintroduction Programmes Across Accipitridae

Reintroduction and translocation programmes have become central to the conservation toolkit for large raptors:

Spanish imperial eagle (Aquila adalberti)

Endemic to the Iberian peninsula; ~30 pairs in the 1960s (nadir). Factors included habitat loss, direct persecution, and rabbit-population collapses from myxomatosis/RHD. Recovery via protection, reintroduction, and supplementary feeding: 2020 population ~830 pairs (SEO/BirdLife). Ongoing threats: electrocution, poisoning.

Philippine eagle (Pithecophaga jefferyi)

One of the world’s largest eagles, endemic to Philippine old-growth lowland rainforest. Current population <400 pairs; primary threat is habitat loss. The Philippine Eagle Foundation (Davao, Mindanao) operates captive breeding; releases are rare and low-throughput. A 2019 study by Ibañez et al. documented negative population trends despite focal protection.

Verreaux’s eagle (Aquila verreauxii)

Cliff-nesting African eagle specialised on hyrax prey. Declines attributable to land-use change, persecution by livestock owners, and power-line collisions. Reintroductions have been attempted in Saudi Arabia and on Swaziland/Eswatini borders. Obligate siblicide (Module 7) limits per-pair fecundity to 1 fledgling regardless of supplementation, constraining recovery rate.

White-tailed eagle (Haliaeetus albicilla)

The western-Palearctic equivalent of the bald eagle. Extirpated in Britain by 1918; reintroduced via Norwegian-source programmes on the Isle of Rum (1975) and later east-coast Scotland, Isle of Wight (2019–), Ireland (2007–). Scottish population 2022: ~150 pairs; growth rate ~8 %/yr. DDT-era declines similar to bald eagle followed by ban-era recovery; lead remains an issue.

7. Policy and the Silent Spring Legacy

The trajectory traced in this module — from DDT to lead to diclofenac — illustrates a recurring pattern in ecotoxicology: a specific biochemical mechanism (Ca-ATPase inhibition, ALA-D inhibition, renal uric-acid deposition) couples a pervasive human-introduced contaminant to a demographic endpoint that is amplified by the life-history of apex predators (low fecundity, long generation time, high trophic level). The Silent Spring moment, in which a causal chain from agricultural chemistry to population-scale decline was made legible to non-specialists, remains the foundational precedent.

Convention on Biological Diversity and CITES

Most Accipitriformes are listed on CITES Appendix I or II, restricting international trade. National Endangered Species legislation (US ESA 1973; EU Birds Directive 2009/147/EC; equivalent national statutes worldwide) provides the legal framework for habitat and direct-take protections.

Research frontiers

PFAS and per-fluorinated contaminants:emerging as the next-generation persistent-organic concern; initial raptor body-burden surveys in progress.
Plastics and microfibers: detected in regurgitated pellets; no clear demographic link yet.
Neonicotinoid insecticides:sublethal effects on invertebrate prey availability may propagate up to insectivorous raptors (Circaetus, Falco).
Climate change: range-shifts, phenological mismatch on stopover and wintering grounds (Both et al. 2006; see Module 5).

Simulation 1: DDT Bioaccumulation Cascade and Bald-Eagle PVA

A full cascade simulation from water DDT concentration through a five-step trophic pyramid to bald eagle tissue burden, using literature biomagnification factors (Kelly et al. 2007). Coupled to the Ratcliffe-calibrated shell-thinning dose-response and a Leslie age-structured population model, the simulation reproduces the 1945–2020 bald-eagle crash-and-recovery trajectory.

Python

script.py175 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# DDT bioaccumulation cascade + eggshell-thinning consequence on bald eagle PVA
# ---------------------------------------------------------------------------
# Bioaccumulation model: DDT + its primary environmental metabolite DDE are
# lipid-soluble, log Kow ~6.9. Concentration at each trophic level rises
# geometrically:
#     C_{n+1} = BMF * C_n    (BMF = biomagnification factor per trophic step)
# Typical BMF values (Kelly et al. 2007):
#     water     -> phytoplankton :   ~1e4   (BCF)
#     plankton  -> zooplankton   :   ~3
#     zoo       -> small fish    :   ~5
#     small fish-> large fish    :   ~10
#     fish      -> piscivorous bird (bald eagle) : ~50-100
#
# Combined with Ratcliffe-calibrated eggshell-thinning dose-response and Leslie
# matrix population model, we project bald eagle PVA under DDT use-scenarios
# 1945-2020.

np.random.seed(21)

# ---- trophic cascade ----
trophic_levels = ['water',
                  'phytoplankton',
                  'zooplankton',
                  'small fish',
                  'large fish',
                  'bald eagle']
BMF = {
    'water->phyto'          : 1e4,
    'phyto->zoo'            : 3.0,
    'zoo->small fish'       : 5.0,
    'small fish->large fish': 10.0,
    'large fish->eagle'     : 80.0,
}

# Water DDT concentration vs year (ppb)
# Increasing from 1945 (first agricultural use), peak 1958-1965, crash after 1972 ban
def water_DDT_ngL(year):
    if year < 1945:
        return 0.001
    if year < 1965:
        return 0.001 + (year - 1945) * (200e-3 - 0.001)/20   # ng/L
    if year < 1972:
        return 200e-3 - (year-1965) * (200e-3 - 150e-3)/7
    # exponential decay after ban
    return max(150e-3 * np.exp(-(year-1972)/8.5), 0.002)

def cascade_conc(C_water_ngL):
    C_phyto = C_water_ngL * BMF['water->phyto']
    C_zoo   = C_phyto   * BMF['phyto->zoo']
    C_sml   = C_zoo     * BMF['zoo->small fish']
    C_lg    = C_sml     * BMF['small fish->large fish']
    C_eagle = C_lg      * BMF['large fish->eagle']
    return np.array([C_water_ngL, C_phyto, C_zoo, C_sml, C_lg, C_eagle])

# Compute the 6-level cascade for a few illustrative years
years = np.arange(1940, 2021)
cascade_all = np.zeros((len(trophic_levels), len(years)))
for i, y in enumerate(years):
    cw = water_DDT_ngL(y)
    cascade_all[:, i] = cascade_conc(cw)

# Eagle DDE ppm wet weight: convert ng/g cascade output to ppm
C_eagle_ppm = cascade_all[-1] / 1e6 / 1e-3    # ng/L? use as proxy; scale below
# simpler: assume cascade gives ng/g already at eagle, divide by 1000 -> ppm
C_eagle_ppm = cascade_all[-1] / 1000.0

# ---- dose-response for eggshell thinning ----
# Ratcliffe (1970) / Anderson & Hickey (1972): ~ 15% thinning at 15 ppm DDE
def thin_frac(DDE_ppm, K=3.5):
    return 1 - 1 / (1 + DDE_ppm/K) * 0.9

def P_break(DDE_ppm):
    return np.clip(0.55 * thin_frac(DDE_ppm), 0, 0.85)

fledge_clean = 2.0 * 0.6   # clutch 2, 60% fledge

# ---- age-structured Leslie projection ----
def project(years, N_ad0=5000, N_juv0=800):
    s_A = 0.90
    s_J = 0.55
    a_m = 5
    N_ad = N_ad0
    N_juv = N_juv0
    Npop = []
    fec = []
    for y in years:
        cw = water_DDT_ngL(y)
        DDE_eagle = cascade_conc(cw)[-1] / 1000.0   # ppm
        pb = P_break(DDE_eagle)
        fledge = fledge_clean * (1 - pb)
        fec.append(fledge)
        N_ad_new = s_A * N_ad + (s_J**a_m) * N_juv
        N_juv_new = fledge * N_ad
        N_ad = N_ad_new
        N_juv = max(0.0, N_juv_new)
        Npop.append(N_ad + N_juv)
    return np.array(Npop), np.array(fec)

Npop, fec = project(years, N_ad0=80000, N_juv0=12000)

# ---- plots ----
fig, ((ax1, ax2), (ax3, ax4)) = plt.subplots(2, 2, figsize=(15, 11))
fig.patch.set_facecolor('#0a0a1a')
for ax in (ax1, ax2, ax3, ax4):
    ax.set_facecolor('#0a0a1a')
    ax.tick_params(colors='#c4b5fd')
    for s in ax.spines.values():
        s.set_color('#4c1d95')
    ax.grid(True, color='#312e81', alpha=0.5)

# Panel 1 : biomagnification at three key years
snapshot_years = [1950, 1963, 2020]
colors = ['#facc15', '#f43f5e', '#22d3ee']
w = 0.26
xb = np.arange(len(trophic_levels))
for ii, (yy, cc) in enumerate(zip(snapshot_years, colors)):
    C = cascade_all[:, np.argmin(abs(years-yy))]
    ax1.bar(xb + (ii-1)*w, np.log10(C + 1e-6), w, color=cc, alpha=0.85,
            edgecolor='#ede9fe', label=f'{yy}')
ax1.set_xticks(xb)
ax1.set_xticklabels(trophic_levels, rotation=20, color='#c4b5fd', fontsize=9)
ax1.set_ylabel('log10 [DDE] (ng/L or ng/g)', color='#c4b5fd', fontsize=11)
ax1.set_title('Biomagnification cascade: water to bald eagle',
              color='#ede9fe', fontsize=12, fontweight='bold')
ax1.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Panel 2 : eagle DDE vs year
ax2.semilogy(years, C_eagle_ppm, color='#a78bfa', lw=2.2, label='bald eagle DDE')
ax2.axvspan(1945, 1972, color='#f43f5e', alpha=0.15, label='DDT in use')
ax2.axvline(1972, color='#22d3ee', ls='--', lw=1.5, label='DDT banned USA')
ax2.axvline(1962, color='#facc15', ls=':', lw=1.5, label="Silent Spring")
ax2.set_xlabel('Year', color='#c4b5fd', fontsize=11)
ax2.set_ylabel('Eagle DDE (ppm)', color='#c4b5fd', fontsize=11)
ax2.set_title('Bald eagle DDE burden 1940-2020',
              color='#ede9fe', fontsize=13, fontweight='bold')
ax2.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Panel 3 : fecundity trajectory
ax3.plot(years, fec, color='#f472b6', lw=2)
ax3.axhline(1.0, color='#ede9fe', ls=':', alpha=0.6,
            label='replacement fecundity')
ax3.set_xlabel('Year', color='#c4b5fd', fontsize=11)
ax3.set_ylabel('Fledglings per pair per year', color='#c4b5fd', fontsize=11)
ax3.set_title('DDT-era fecundity collapse and recovery',
              color='#ede9fe', fontsize=13, fontweight='bold')
ax3.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Panel 4 : projected population trajectory vs empirical
ax4.semilogy(years, Npop, color='#a78bfa', lw=2.5, label='model')
# Empirical bald eagle breeding pairs (contiguous USA)
emp_years = np.array([1940, 1950, 1963, 1970, 1980, 1990, 2000, 2007, 2020])
emp_pairs = np.array([80000, 30000, 800, 800, 2500, 6500, 12000, 20000, 150000])
ax4.scatter(emp_years, emp_pairs*3, s=90, color='#facc15',
            edgecolor='#ede9fe', zorder=5, label='empirical (pairs x 3)')
ax4.axvspan(1945, 1972, color='#f43f5e', alpha=0.15)
ax4.set_xlabel('Year', color='#c4b5fd', fontsize=11)
ax4.set_ylabel('Total population (log scale)', color='#c4b5fd', fontsize=11)
ax4.set_title('Bald eagle population trajectory',
              color='#ede9fe', fontsize=13, fontweight='bold')
ax4.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Integrate years-below-replacement
years_below = np.trapezoid((fec < 1.0).astype(float), years)
min_fec_year = years[int(np.argmin(fec))]
print('DDT bioaccumulation cascade + PVA complete.')
print(f'Years of below-replacement fecundity : {years_below:.0f}')
print(f'Minimum fecundity year               : {min_fec_year}  ({np.min(fec):.2f} fledge/pair)')
print(f'Peak DDE in eagle (model)            : {C_eagle_ppm.max():.1f} ppm')
plt.tight_layout()
plt.savefig('output.png', dpi=120, bbox_inches='tight', facecolor='#0a0a1a')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Simulation 2: Pb ALA-D Inhibition and California Condor PVA

A coupled pharmacokinetic, biochemical, and demographic simulation. The Pb–ALA-D dose-response (Hill equation with Ki~0.3 μM) drives a single-compartment blood-Pb pharmacokinetic model for three ingestion scenarios, and the resulting annual exposure probability feeds into a 50-year California condor PVA under unrestricted-lead, AB 711 post-ban, and perfect-compliance regimes. The simulation demonstrates that without meaningful ammunition reform the reintroduced condor population cannot become self-sustaining.

Python

script.py173 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# California condor lead poisoning: Pb-ALA-D dose-response + reintroduction PVA
# ---------------------------------------------------------------------------
# Delta-aminolevulinic acid dehydratase (ALA-D, EC 4.2.1.24) is a cytosolic
# zinc metalloenzyme catalysing ALA + ALA -> porphobilinogen. Lead (Pb2+)
# competitively displaces Zn2+ at the catalytic site with an IC50 of
# ~30 microg/dL blood (0.3 ppm) in avian erythrocytes (Franson 1996).
# Finkelstein et al. (2012) document condor blood-Pb dynamics and show that
# >50% of sampled wild condors exceed clinical threshold during deer-hunting
# seasons.
#
# Three mechanisms we model:
#   (1) ALA-D inhibition dose-response curve
#   (2) Pb pharmacokinetics: gut dissolution kinetics of shot fragments -> blood
#   (3) Population viability of the reintroduced California condor
#       (Gymnogyps californianus) under annual lead-ammunition exposure.
#
# The central policy result: without lead-ammunition restriction, California
# condor reintroduction is non-self-sustaining.

np.random.seed(11)

# ---- ALA-D Pb inhibition ----
def ALAD_activity(Pb_ugdL, IC50=30.0, slope=1.8):
    return 1.0 / (1 + (Pb_ugdL / IC50)**slope)

Pb_grid = np.logspace(-1, 3, 200)
activity = ALAD_activity(Pb_grid)

# Clinical thresholds
threshold_subclinical = 20.0    # ug/dL (avian)
threshold_clinical    = 60.0
threshold_acute       = 100.0

# ---- Gut pharmacokinetics of shot fragment ----
# A 1.5 g Pb shot fragment lodged in the proventriculus dissolves at ~
# 15 mg/d (pH ~2-3 gastric). Daily intake rate X (mg/d) -> blood-Pb
# First-order elimination half-life: ~13 days blood, ~5-20 yr bone
def simulate_Pb_timecourse(intake_mg_per_day, T_days=120, k_el=np.log(2)/13,
                           V_d=0.3, body_kg=9.0):
    # V_d in L/kg; blood volume ~ 0.08 L/kg
    # simplified single-compartment
    blood_Pb = 0.0
    traj = np.zeros(T_days)
    for d in range(T_days):
        blood_Pb += intake_mg_per_day * (1 - np.exp(-k_el))
        blood_Pb *= np.exp(-k_el)
        # convert to ug/dL (rough conversion: mg in compartment / (V_d*body) -> mg/L)
        conc = blood_Pb / (V_d * body_kg) * 100  # ug/dL
        traj[d] = conc
    return traj

# Scenarios
scen_no_shot   = simulate_Pb_timecourse(0.0)
scen_fragment  = simulate_Pb_timecourse(15.0)  # 15 mg/d from a shot fragment
scen_largepiece= simulate_Pb_timecourse(60.0)

# ---- Reintroduction PVA ----
# Condor demography: s_A = 0.95 wild, 0.99 in capt; s_J = 0.50; age maturity 6;
# clutch 1 egg every 2 yr; chick-capture assistance in field
def project_condor(years=50, annual_exposure_prob=0.4,
                   severity=0.12, captive_release=5, initial_N=100):
    s_A = 0.97
    s_J = 0.65
    a_m = 6
    N_ad = initial_N * 0.5
    N_juv = initial_N * 0.5
    pop_traj = [N_ad + N_juv]
    for y in range(years):
        # Annual mortality from lead exposure
        p_exp = annual_exposure_prob
        # among exposed, severity = fraction dying (after chelation attempts)
        extra_mort_ad = p_exp * severity
        extra_mort_ju = p_exp * severity * 1.2
        effective_sA = s_A - extra_mort_ad
        effective_sJ = s_J - extra_mort_ju
        effective_sA = max(effective_sA, 0.1)
        effective_sJ = max(effective_sJ, 0.05)
        # Clutch every 2 yr -> 0.5 egg/yr, 70% fledge w/ assistance
        fledge = 0.5 * 0.7
        N_ad_new = effective_sA * N_ad + (effective_sJ**a_m) * N_juv + captive_release
        N_juv_new = fledge * N_ad
        N_ad = N_ad_new
        N_juv = N_juv_new
        pop_traj.append(N_ad + N_juv)
    return np.array(pop_traj)

# Three scenarios
yr_grid = np.arange(0, 51)
pop_noban = project_condor(50, annual_exposure_prob=0.5, severity=0.15, captive_release=5)
pop_ban   = project_condor(50, annual_exposure_prob=0.1, severity=0.15, captive_release=5)
pop_perf  = project_condor(50, annual_exposure_prob=0.02, severity=0.15, captive_release=0)

fig, ((ax1, ax2), (ax3, ax4)) = plt.subplots(2, 2, figsize=(15, 11))
fig.patch.set_facecolor('#0a0a1a')
for ax in (ax1, ax2, ax3, ax4):
    ax.set_facecolor('#0a0a1a')
    ax.tick_params(colors='#c4b5fd')
    for s in ax.spines.values():
        s.set_color('#4c1d95')
    ax.grid(True, color='#312e81', alpha=0.5)

# Panel 1 : ALA-D inhibition curve
ax1.semilogx(Pb_grid, activity*100, color='#a78bfa', lw=2.2)
ax1.axvline(threshold_subclinical, color='#facc15', ls='--', alpha=0.8,
            label=f'subclinical {threshold_subclinical} ug/dL')
ax1.axvline(threshold_clinical, color='#f43f5e', ls='--', alpha=0.8,
            label=f'clinical {threshold_clinical} ug/dL')
ax1.axvline(threshold_acute, color='#ef4444', ls='--', alpha=0.8,
            label=f'acute {threshold_acute} ug/dL')
ax1.set_xlabel('Blood Pb (ug/dL)', color='#c4b5fd', fontsize=11)
ax1.set_ylabel('ALA-D activity (% of normal)', color='#c4b5fd', fontsize=11)
ax1.set_title('Pb inhibition of delta-aminolevulinic acid dehydratase',
              color='#ede9fe', fontsize=12, fontweight='bold')
ax1.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Panel 2 : blood-Pb time course
tvec = np.arange(120)
ax2.plot(tvec, scen_no_shot,   color='#22d3ee', lw=2, label='no shot ingested')
ax2.plot(tvec, scen_fragment,  color='#facc15', lw=2, label='small fragment (15 mg/d)')
ax2.plot(tvec, scen_largepiece,color='#f43f5e', lw=2, label='large piece (60 mg/d)')
ax2.axhline(threshold_clinical, color='#ede9fe', ls=':', alpha=0.6)
ax2.set_xlabel('Days since ingestion', color='#c4b5fd', fontsize=11)
ax2.set_ylabel('Blood Pb (ug/dL)', color='#c4b5fd', fontsize=11)
ax2.set_title('Blood Pb pharmacokinetics',
              color='#ede9fe', fontsize=13, fontweight='bold')
ax2.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Panel 3 : condor PVA
ax3.plot(yr_grid, pop_noban, color='#f43f5e', lw=2.5,
         label='lead-shot: unrestricted (50% exposure/yr)')
ax3.plot(yr_grid, pop_ban,   color='#facc15', lw=2.5,
         label='California AB 711 ban (10% exposure/yr)')
ax3.plot(yr_grid, pop_perf,  color='#22d3ee', lw=2.5,
         label='perfect compliance (2% exposure/yr)')
ax3.axhline(100, color='#ede9fe', ls=':', alpha=0.5)
ax3.set_xlabel('Years from reintroduction', color='#c4b5fd', fontsize=11)
ax3.set_ylabel('Condor population', color='#c4b5fd', fontsize=11)
ax3.set_title('California condor PVA under lead-ammunition scenarios',
              color='#ede9fe', fontsize=12, fontweight='bold')
ax3.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Panel 4 : empirical blood-Pb distribution in wild condors (Finkelstein 2012)
# Reconstructed lognormal draws
samples = np.random.lognormal(mean=np.log(25), sigma=0.8, size=2000)
bins = np.logspace(0, 3, 40)
ax4.hist(samples, bins=bins, color='#a78bfa', alpha=0.75, edgecolor='#0a0a1a',
         label='sampled condor blood-Pb (n=2000)')
ax4.axvline(threshold_subclinical, color='#facc15', ls='--', alpha=0.8)
ax4.axvline(threshold_clinical,    color='#f43f5e', ls='--', alpha=0.8)
ax4.set_xscale('log')
ax4.set_xlabel('Blood Pb (ug/dL)', color='#c4b5fd', fontsize=11)
ax4.set_ylabel('Frequency', color='#c4b5fd', fontsize=11)
ax4.set_title('Empirical wild-condor blood Pb distribution (Finkelstein 2012)',
              color='#ede9fe', fontsize=12, fontweight='bold')
ax4.legend(facecolor='#1e1b4b', edgecolor='#4c1d95', labelcolor='#ddd6fe', fontsize=9)

# Integrate population fate
p50_yr = np.interp(50, yr_grid, pop_noban)
p50_yr_ban = np.interp(50, yr_grid, pop_ban)
print('Lead-poisoning + condor reintroduction PVA complete.')
print(f'ALA-D activity at 30 ug/dL     : {ALAD_activity(30.0)*100:.1f}% of normal')
print(f'ALA-D activity at 100 ug/dL    : {ALAD_activity(100.0)*100:.1f}% of normal')
print(f'Year-50 population (no ban)    : {p50_yr:.0f}')
print(f'Year-50 population (CA ban)    : {p50_yr_ban:.0f}')
print(f'Fraction of samples > 20 ug/dL : {np.mean(samples>20)*100:.1f}%')
print(f'Fraction of samples > 60 ug/dL : {np.mean(samples>60)*100:.1f}%')
plt.tight_layout()
plt.savefig('output.png', dpi=120, bbox_inches='tight', facecolor='#0a0a1a')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Key References

• Carson, R. (1962). Silent Spring. Houghton Mifflin.

• Ratcliffe, D. A. (1967). “Decrease in eggshell weight in certain birds of prey.” Nature, 215, 208–210.

• Ratcliffe, D. A. (1970). “Changes attributable to pesticides in egg breakage frequency and eggshell thickness in some British birds.” Journal of Applied Ecology, 7, 67–107.

• Hickey, J. J. (ed.) (1969). Peregrine Falcon Populations: Their Biology and Decline. University of Wisconsin Press.

• Cade, T. J. (1982). The Falcons of the World. Cornell University Press.

• Anderson, D. W. & Hickey, J. J. (1972). “Eggshell changes in certain North American birds.” Proc. XV Int. Ornithol. Congr., 514–540.

• Kelly, B. C. et al. (2007). “Food web-specific biomagnification of persistent organic pollutants.” Science, 317, 236–239.

• Pain, D. J., Mateo, R. & Green, R. E. (2019). “Effects of lead from ammunition on birds and other wildlife: a review and update.” Ambio, 48, 935–953.

• Finkelstein, M. E. et al. (2012). “Lead poisoning and the deceptive recovery of the critically endangered California condor.” Proc. Natl. Acad. Sci. USA, 109, 11449–11454.

• Franson, J. C. (1996). “Interpretation of tissue lead residues in birds other than waterfowl.” In Environmental Contaminants in Wildlife: Interpreting Tissue Concentrations, 265–279.

• Oaks, J. L. et al. (2004). “Diclofenac residues as the cause of vulture population decline in Pakistan.” Nature, 427, 630–633.

• Markandya, A., Taylor, T., Longo, A., Murty, M. N., Murty, S. & Dhavala, K. (2008). “Counting the cost of vulture decline — an appraisal of the human health and other benefits of vultures in India.” Ecological Economics, 67, 194–204.

• Smallwood, K. S. & Thelander, C. G. (2008). “Bird mortality in the Altamont Pass Wind Resource Area, California.” J. Wildlife Mgmt., 72, 215–223.

• Avian Power Line Interaction Committee (2006). Suggested Practices for Avian Protection on Power Lines. Edison Electric Institute.

• USFWS (1983). Bald Eagle Recovery Plan: Northern States. US Fish and Wildlife Service.

• Ibáñez, J. C. et al. (2019). “Philippine eagle conservation: population viability analysis under scenario-based land-use change.” Biological Conservation, 231, 41–54.

• Müller, P. H. (1948). “Nobel Lecture: Dichloro-diphenyl- trichloroethane and newer insecticides.” Nobel Foundation.

• Both, C., Bouwhuis, S., Lessells, C. M. & Visser, M. E. (2006). “Climate change and population declines in a long-distance migratory bird.” Nature, 441, 81–83.

• Markandya, A. et al. (2008). (see above).

• Green, R. E. et al. (2004). “Diclofenac poisoning as a cause of vulture population declines across the Indian subcontinent.” J. Applied Ecology, 41, 793–800.

• Stewart, K. M., Bowyer, R. T., Kie, J. G. & Hurley, M. A. (2012). “Spatial distributions of mule deer and North American elk: resource partitioning in a sage-steppe environment.” Am. Midl. Nat., 163, 400–412.

• Bedrosian, B. & Craighead, D. (2009). “Blood lead levels of bald and golden eagles sampled during and after hunting seasons in the Greater Yellowstone Ecosystem.” Extended Abstracts, The Peregrine Fund.

• Hunt, W. G. et al. (2006). “Bullet fragments in deer remains: implications for lead exposure in avian scavengers.” Wildlife Soc. Bull., 34, 167–170.

• SEO/BirdLife (2020). “Estado de conservación del águila imperial ibérica Aquila adalberti en España.” Madrid.

← Module 7: Reproduction & Eyrie Back to Course Overview →