Module 4

Wood Formation & Cell Wall Mechanics

Xylogenesis, cellulose/lignin biosynthesis, microfibril angle, viscoelastic mechanics, and reaction wood

4.1 Cambial Activity & Xylogenesis

Wood is produced by the vascular cambium, a lateral meristem one to a few cells thick wrapped around the trunk. Cambial initials divide periclinally to produce xylem mother cells (centripetally) and phloem mother cells (centrifugally), then undergo further differentiation.

Auxin Gradient Control

Indole-3-acetic acid (IAA, auxin) is the primary cambial activator. It is synthesised in young leaves and transported polarly downward via the PIN-family efflux carriers (PIN1, PIN2) and AUX1 influx carriers. The auxin gradient across the cambial zone controls cell differentiation:

  • High IAA (\(> 150\) ng g\(^{-1}\) FW): cambial cell division, early xylem expansion
  • Medium IAA (~80–150 ng g\(^{-1}\)): secondary wall deposition begins
  • Low IAA (\(< 40\) ng g\(^{-1}\)): programmed cell death (autolysis), mature wood cell

The local auxin concentration follows a reaction-diffusion model with active transport. The steady-state profile in a tissue of width \(L\) with degradation rate \(k\) and active transport velocity \(v\):

\[ D\frac{d^2[IAA]}{dx^2} - v\frac{d[IAA]}{dx} - k[IAA] = 0 \]

with solution \([IAA](x) = A e^{\lambda_1 x} + B e^{\lambda_2 x}\),\(\lambda_{1,2} = (v \pm \sqrt{v^2 + 4Dk})/(2D)\).

Stages of Xylogenesis

1. Cell Division
Cambial initial divides; new xylem mother cell formed; 1–3 h per cell cycle in active season
2. Cell Expansion
Primary wall loosening by expansins (pH 4.5 activation); Lockhart growth equation; radial expansion 5–20× diameter increase
3. Secondary Wall
Cessation of expansion; cellulose synthase (CESA) complexes deposit S1, S2, S3 layers; lignin deposition begins simultaneously
4. Lignification
Monolignol monomers transported, oxidised by peroxidases/laccases; middle lamella lignified first, then S1, S2
5. Programmed Cell Death
Vacuolar processing enzymes (VPE) trigger autolysis; nucleus, organelles, tonoplast degrade; dead cell remains as conduit
6. Maturation
Pit aspiration (conifers) or vessel perforation (angiosperms); final cell dimensions fixed

4.2 Cell Wall Structure & Cellulose

CESA Complex & Microfibril Angle

Cellulose is synthesised by CESA (Cellulose Synthase A)complexes, visible as hexameric “rosettes” of 18–36 catalytic subunits in freeze-fracture electron microscopy. Each complex extrudes one glucan chain; the 18–36 chains bundle into a elementary microfibril of 3–5 nm diameter. The complex moves through the plasma membrane, driven by the force of microfibril crystallisation, its direction guided by cortical microtubules.

The Microfibril Angle (MFA) is the angle between cellulose microfibrils and the cell's longitudinal axis in the S2 layer. It profoundly affects wood mechanics:

\[ E_{cell} \approx E_{cellulose} \cos^4\theta_{MFA} + E_{matrix} \sin^4\theta_{MFA} \]

Small MFA (\(\theta < 10^\circ\)) gives high stiffness and low shrinkage (mature wood). Large MFA (\(\theta > 30^\circ\)) gives flexibility and high longitudinal shrinkage (juvenile wood, tension wood G-layer). In the S2 layer of normal softwood, MFA \(\approx 5\text{--}20^\circ\).

Cell Wall Layer Architecture

Cell LumenS3MFA ~70-90 degThin (~0.1 um)S2MFA 5-30 degThick (~2-5 um)80% of wallS1MFA 50-70 degThin (~0.2 um)P + MLPrimary wall+ Middle lamella~70-80% ligninWood Cell Wall Layers (cross-section)

4.3 Lignin Biosynthesis

Lignin is a complex aromatic polymer providing compression strength and waterproofing. It is synthesised from three monolignol precursors (H, G, S) via the general phenylpropanoid pathway:

PhenylalaninePALtrans-CinnamateC4Hp-Coumarate4CLp-Coumaryl-CoA
HCT/CSTCaffeate → Ferulate (G-unit precursor)→ CCoAOMT → CCR → CAD

The three monolignol types and their structural characteristics:

H-lignin
p-Hydroxyphenyl
Methoxy groups: 0
Mainly in: Gymnosperms (minor)
G-lignin
Guaiacyl
Methoxy groups: 1
Mainly in: Gymnosperms (dominant)
S-lignin
Syringyl
Methoxy groups: 2
Mainly in: Angiosperms (dominant)

Monolignols are exported to the cell wall where they are oxidised by class III peroxidases (H\(_2\)O\(_2\)-dependent) and laccases (O\(_2\)-dependent) to resonance-stabilised radicals, which couple randomly to form the branched polymer network. Key inter-unit bonds: \(\beta\)-O-4 aryl ether (48–60%),\(\beta\)-\(\beta\) (9–12%), \(\beta\)-5 phenylcoumaran (6–12%).

S/G ratio varies: conifers (0/80:20 G:H), hardwoods (40–70% S). Higher S content correlates with lower cross-linking frequency (S-units cannot form 5-5 bonds) and higher pulp yield for paper making. Average molecular mass of lignin in wood: \(\bar{M}_n \approx 5\)\(20\) kDa.

Lignin Monomer Structures

H unitp-OHCH=CH-CH2OHOHG unitp-OHm-OMeCH=CH-CH2OHOHOCH3S unitp-OH2x OMeCH=CH-CH2OHOHOCH3OCH3Lignin Monolignol Units (H, G, S)Each bears a propenyl side chain (shown above ring)

4.4 Viscoelastic & Composite Mechanics

Lockhart Growth Equation (revisited)

During expansion (Stage 2 of xylogenesis), the growing cell obeys the Lockhart equation where the wall extensibility \(\phi\) is set by expansin activity and wall pH:

\[ \frac{1}{V}\frac{dV}{dt} = \phi(\Psi_p - Y), \quad \phi = \phi_0 \cdot 10^{(pH_{opt} - pH)/n} \]

Expansin activity is maximal at pH 4.5 (acid growth hypothesis: H\(^+\)-ATPase acidifies the wall). The wall yield threshold \(Y \approx 0.2\text{--}0.5\) MPa for expanding wood cells.

Viscoelastic Creep

Wood is viscoelastic: under constant stress, it deforms over time (creep). The creep compliance\(J(t)\) (strain per unit stress) is described by the standard linear solid (Kelvin–Voigt model with spring in series):

\[ J(t) = J_0 + J_1\!\left(1 - e^{-t/\tau}\right) + \frac{t}{\eta_v} \]

where \(J_0 = 1/E_0\) is the instantaneous elastic compliance, \(J_1\) is the delayed elastic compliance (recoverable), \(\tau = \eta_v J_1\) is the retardation time (minutes to days for wood), and the last term represents viscous flow (permanent). At room temperature and 12% MC,\(\tau \approx 10\text{--}100\) min for small clear wood specimens.

Composite Stiffness: Voigt & Reuss Models

Wood as a fibre composite: cellulose fibrils (\(E_c \approx 135\) GPa, volume fraction \(V_c \approx 0.4\)) embedded in a matrix of hemicellulose + lignin (\(E_m \approx 2\text{--}4\) GPa):

\[ E_{Voigt} = E_c V_c + E_m V_m \quad \text{(upper bound, isostrain)} \]
\[ \frac{1}{E_{Reuss}} = \frac{V_c}{E_c} + \frac{V_m}{E_m} \quad \text{(lower bound, isostress)} \]

Real wood lies between these bounds. The Halpin–Tsai equation for the longitudinal modulus accounts for fibre aspect ratio \(\xi = 2l/d\):

\[ E_L = E_m \frac{1 + \xi\eta V_c}{1 - \eta V_c}, \qquad \eta = \frac{E_c/E_m - 1}{E_c/E_m + \xi} \]

4.5 Reaction Wood

Trees respond to mechanical displacement (leaning, gravity sensing by statoliths in endodermis) by producing anatomically modified wood that generates restoring forces:

Tension Wood (Angiosperms)

  • Forms on the upper side of leaning stem
  • Contains a G-layer (gelatinous layer) of almost pure cellulose (\(> 95\%\))
  • G-layer has very small MFA (\(< 5^\circ\)) and high crystallinity
  • Generates large tensile maturation stress (\(\sigma_T \approx\) 10–70 MPa) on setting
  • Mechanism: cellulose crystallisation shrinkage pulls fibre to smaller length

Compression Wood (Gymnosperms)

  • Forms on the lower side of leaning stem
  • Tracheids: rounded, intercellular spaces, thick S2, helical checks
  • High lignin content (\(> 30\%\) vs. normal 25%)
  • Large MFA (\(> 30^\circ\)) in S2 layer
  • Generates compressive maturation stress (\(\sigma_C \approx\) 2–20 MPa)
  • Mechanism: lignin swelling during polymerisation pushes fibre longitudinally

Together, tension wood pulling from above and compression wood pushing from below straighten the leaning stem in a months-long process. The balance determines righting efficiency and contributes to stem form optimisation.

4.6 Python: Wood Mechanics & MFA Effects

We model wood cell longitudinal modulus as a function of microfibril angle and lignin content, then generate stress–strain curves for normal wood, tension wood, and compression wood using a viscoelastic composite model.

Python
script.py106 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server