Chlorophyll-protein complexes are ideal model systems for protein energy landscape research. Here pigments, used in optical spectroscopy experiments as sensitive probes to local dynamics, are built into protein by Nature (in a large variety of local environments; without extraneous chemical manipulations or genetic engineering). Distributions of the tunneling parameter, λ, and/or protein energy landscape barrier heights, V, have been determined for (the lowest energy state of) the CP43 core antenna complex of photosystem II. We demonstrate that spectral hole burning (SHB) and hole recovery (HR) measurements are capable of delivering important information on protein energy landscape properties and spectral diffusion mechanism details. In particular, we show that tunneling rather than barrier hopping is responsible for both persistent SHB and subsequent HR at 5-12 K, which allows us to estimate the md(2) parameter of the tunneling entities as ~1.0 × 10(-46) kg·m(2). The subdistributions of λ actually contributing to the nonsaturated spectral holes (and affecting their recovery) differ from the respective full true distributions. In the case of the full λ-distribution being uniform (or the barrier height distribution ~1/√V, a model which has been widely employed in theories of amorphous solids at low temperatures and in HR analysis), the difference is qualitative, with λ subdistributions probed in the HR experiments being highly asymmetrical, and barrier V subdistributions deviating significantly from ~1/√V. Thus, the distribution of λ for the protein energy landscape tier directly probed by SHB is likely Gaussian and not uniform. Additionally, a Gaussian distribution of barriers, with parameters incompatible with those of the landscape tier directly probed by SHB, contributes to the thermocycling results.