Perez Model

​

Summary

The Perez transposition model estimates sky diffuse irradiance on a tilted surface by decomposing DHI into three components: isotropic background, circumsolar brightening near the solar disk, and horizon brightening near the horizon band. It uses empirical coefficients derived from extensive sky radiance measurements. Sky clearness and brightness indices determine which coefficient set is applied. PlantPredict supports multiple Perez coefficient sets derived from different locations and time periods.

​

Inputs

---

​

Outputs

---

​

Detailed Description

​

Sky Clearness and Brightness Indices

The Perez model uses two atmospheric parameters to characterize sky conditions: Sky Clearness Index ($\varepsilon$): Represents the clarity of the atmosphere, accounting for the ratio of total to diffuse horizontal irradiance and solar zenith angle: $$\varepsilon = \frac{\frac{DHI + DNI}{DHI} + \kappa \theta_z^3}{1 + \kappa \theta_z^3}$$ where $\kappa = 5.535 \times 10^{-6}$ is an empirical constant that corrects the clearness index for zenith angle dependence. Sky Brightness Index ($\Delta$): Represents the amount of diffuse irradiance relative to extraterrestrial irradiance, normalized by air mass: $$\Delta = \frac{DHI \cdot AM}{DNI_{extra}}$$ where $AM$ is the (non pressure-corrected) air mass calculated using the configured air mass model.

​

Perez Coefficient Lookup

The clearness index $\varepsilon$ is binned into 8 categories corresponding to sky conditions from overcast to clear. Each bin has associated empirical coefficients used to compute brightness coefficients $F_1$ and $F_2$: For each bin, coefficients $f_{11}, f_{12}, f_{13}, f_{21}, f_{22}, f_{23}$—from which the Brightness Coefficients $F_1$ and $F_2$ are computed—are defined. PlantPredict supports multiple coefficient sets, each tailored to specific climatic conditions and derived from empirical measurements: Composite Coefficient Sets:

• PlantPredict: Default coefficient set, matches the All Sites Composite 1990 set to the third digit
• All Sites Composite 1990: Comprehensive set derived from data across multiple locations using 1990 methodology.
• All Sites Composite 1988: Earlier composite set based on 1988 methodology from various sites.
• Sandia Composite 1988: Developed using data from Sandia National Laboratories.
• USA Composite 1988: Based on data collected from various U.S. locations.

Location-Specific Coefficient Sets (1988):

• France 1988: Derived from data collected in France.
• Phoenix 1988: Derived from data collected in Phoenix, Arizona.
• El Monte 1988: Derived from data collected in El Monte, California.
• Osage 1988: Derived from data collected in Osage, Iowa.
• Albuquerque 1988: Derived from data collected in Albuquerque, New Mexico.
• Cape Canaveral 1988: Derived from data collected in Cape Canaveral, Florida.
• Albany 1988: Derived from data collected in Albany, New York.

The selection of coefficient set can influence model accuracy depending on site-specific climatic conditions. Composite sets provide general applicability, while location-specific sets may offer improved accuracy for sites with similar climate characteristics. The complete coefficient values for all sets can be found in the pvlib documentation.

​

Brightness Coefficients

The brightness coefficients $F_1$ and $F_2$ are computed using the binned coefficients: $$F_1 = f_{11} + f_{12} \Delta + f_{13} \cdot \frac{\pi}{180} \theta_z$$ $$F_2 = f_{21} + f_{22} \Delta + f_{23} \cdot \frac{\pi}{180} \theta_z$$ A quality control constraint is applied: $F_1 \geq 0$. Physically, $F_1$ quantifies circumsolar brightening (which cannot be negative), while $F_2$ quantifies horizon brightening or darkening (can be negative under clear sky conditions).

​

Diffuse Irradiance Components

DHI is then decomposed into three components (isotropic sky, circumsolar, and horizon) using $F_1$ and $F_2$. Isotropic Component: Uniform diffuse irradiance from the sky dome, reduced by the circumsolar fraction $F_1$: $$G_{isotropic} = DHI \cdot (1 - F_1) \cdot \frac{1 + \cos(\beta)}{2}$$ The term $(1 + \cos(\beta))/2$ is the view factor from the tilted surface to the sky dome. Circumsolar Component: Directional diffuse irradiance concentrated near the solar disk: $$G_{circumsolar} = DHI \cdot F_1 \cdot \frac{\cos(\theta)}{\cos(\theta_z)}$$ The ratio $\cos(\theta)/\cos(\theta_z)$ is the projection ratio that converts horizontal circumsolar irradiance to the tilted plane, identical to the Hay-Davies model. When the sun is behind the module ($\cos(\theta) < 0$), the circumsolar component is set to zero. To avoid numerical instabilities near the horizon, a minimum threshold is applied: $\cos(\theta_z) \geq 0.0872$ (corresponding to $\theta_z \leq 85°$). Horizon Component: Diffuse irradiance from the horizon band: $$G_{horizon} = DHI \cdot F_2 \cdot \sin(\beta)$$ The horizon component can be negative when $F_2 < 0$, which can occur under clear sky conditions when the horizon band is darker than the average sky dome (horizon darkening effect).

​

Circumsolar Allocation

The circumsolar component can be allocated to beam or diffuse POA irradiance, similar to the Hay-Davies model: Direct Allocation: $$G_{sky} = G_{isotropic} + G_{horizon}$$ $$G_{beam} = DNI \cdot \cos(\theta) + G_{circumsolar}$$ Diffuse Allocation: $$G_{sky} = G_{isotropic} + G_{circumsolar} + G_{horizon}$$ $$G_{beam} = DNI \cdot \cos(\theta)$$ Circumsolar irradiance originates from near the solar disk and is blocked by obstructions the same way direct beam is. Allocating it to beam ensures that row-to-row shading calculations apply appropriate losses to circumsolar. Allocating to diffuse treats circumsolar as unaffected by direct shading, which may be appropriate for unshaded systems or when shading is negligible.

​

Ground-Reflected Component

The ground-reflected component accounts for irradiance reflected from the ground onto the tilted surface: $$G_{ground} = \rho \cdot GHI \cdot \frac{1 - \cos(\beta)}{2}$$ The term $(1 - \cos(\beta))/2$ is the view factor from the tilted surface to the ground.

​

Quality Control

Physical constraints are enforced by clamping values:

• If $F_1 < 0$ → $F_1 = 0$
• If $G_{isotropic} < 0$ → $G_{isotropic} = 0$
• If $G_{circumsolar} < 0$ → $G_{circumsolar} = 0$
• If $G_{sky} > 800$ W/m² → $G_{sky} = 800$ W/m²

The horizon component $G_{horizon}$ is not clamped and may be negative.

---

​

References

• Perez, R., Seals, R., Ineichen, P., Stewart, R., & Menicucci, D. (1987). A new simplified version of the Perez diffuse irradiance model for tilted surfaces. Solar Energy, 39(3), 221–231.
• Perez, R., Ineichen, P., Seals, R., Michalsky, J., & Stewart, R. (1990). Modeling daylight availability and irradiance components from direct and global irradiance. Solar Energy, 44(5), 271–289.
• pvlib python. Irradiance module source code. https://pvlib-python.readthedocs.io/en/latest/_modules/pvlib/irradiance.html