Skip to content

Research

My ongoing work on solar radiation data in West Africa: how far a ground calibration stays valid away from the station that produced it, and what a satellite product erases where its grid cell is too coarse.

01In progress

Spatial calibration of solar radiation

The transferability of a ground calibration is governed by climatic similarity between sites, not by geographic distance.

Progress. Analysis complete over the 31 stations and 930 pairs, confirmed on two independent satellite products (CAMS then SARAH-3), then mapped across West Africa. Article write-up in progress.

At a glance

31
ESMAP/WAPP stations
930
source-target pairs
2
independent products
+0.32
residual vs Δclim (GHI)
01

The problem

Satellite products (GHI, DNI) are locally biased; they are corrected by calibrating against ground measurements. But measurement stations are scarce in West Africa, and sizing is almost always done far from any station.

The question: is a calibration established at one station transferable elsewhere, and by what criterion? Common intuition assumes geographic proximity governs transferability. This work tests the competing hypothesis: it is the climatic similarity between sites, written Δclim\Delta_{\mathrm{clim}}, that governs it: two sites with close climates but geographically distant calibrate one another better than two neighbours with different regimes.

02

Hypotheses

H1

The transferability of a calibration is governed by climatic similarity Δclim\Delta_{\mathrm{clim}} between sites, not by their geographic distance.

H2

The transfer residual grows with Δclim\Delta_{\mathrm{clim}}, and this law is stable from one independent satellite product to another.

H3

Accounting for aerosol (AOD) enriches the law and reduces the transfer residual.

03

Method & maths

A network of 31 stations, i.e. 930 ordered source–target pairs. For each pair, the product is calibrated on the source (bias correction by quantile matching), the calibration is applied to the target, and the residual is measured. The residual is then compared against two competing distances: geographic distance and climatic distance Δclim\Delta_{\mathrm{clim}}.

Tools: a relationship between residual and Δclim\Delta_{\mathrm{clim}}, a variogram pitting climate against distance, station-by-station cross-validation, an aerosol-aware enriched version, then replication on a second independent product (SARAH-3 after CAMS) and a map of West Africa at 1° resolution. The whole thing follows a deterministic pipeline, traced to the formulas and independently checked.

Fig. 1The ESMAP/WAPP network: 31 stations across West Africa (top) and the distribution of the 930 inter-station distances, from 100 to 2,800 km (bottom).
04

Results & figures

Transferability is indeed governed by Δclim\Delta_{\mathrm{clim}} and not by distance: the climatic variogram structures the residual where geographic distance does not. The result is validated on 31 stations and 930 pairs, with two independent satellite products (CAMS, SARAH-3), then mapped across West Africa, applied to the ESMAP/WAPP network, the Guinea case, with a field anchor at Kankan (ESMAP/WAPP campaign).

The resulting method, CLIM3, calibrates each site on its three best climatic analogues, weighted by the inverse of Δclim\Delta_{\mathrm{clim}}. On Guinea, it brings the expected residual down to 2.4% for GHI and 6.9% for DNI. A few stations well spread across climate zones are then enough to cover a territory, provided the least-instrumented areas are validated locally.

Fig. 2Variogram. The residual semivariance rises clearly with climatic distance Δclim (r = +0.53 for GHI, +0.62 for DNI) yet stays flat with geographic distance (r ≈ +0.05). Climate, not distance, structures the transfer.
Fig. 3Point by point over the 930 pairs: the transfer penalty correlates with Δclim (r = +0.38 GHI, +0.26 DNI) and never with distance (r ≈ +0.05). The law holds across the whole cloud.
Fig. 4Expected residual mapped across West Africa (1° grid), for GHI and DNI: the method extends to where no station exists.
Fig. 5Zoom on Guinea: application to the use case, with the Kankan field anchor (ESMAP/WAPP campaign) and the Tarambaly site.
05

Open questions

  1. 01DNI stays the least well-constrained component: its climatic trend is there, but still depends on the satellite product. It needs evidence as strong as GHI's.
  2. 02The result rests on a two-year window, a single network and a single region: beyond that, generality remains to be shown.
  3. 03The map is validated station by station, with no ground measurement inside a cell: a short campaign in a target zone would confirm it.
  4. 04The aerosol optical depth comes from a coarse reanalysis; a finer product would sharpen Δclim\Delta_{\mathrm{clim}}.
02In progress

Pixel degeneracy of solar products

Naming and quantifying the case where a satellite product assigns the same value to genuinely distinct sites, and measuring where spatial information is erased.

Progress. Diagnosis established on NASA POWER (34 communes of Guinea); extension to finer products and article write-up to come.

At a glance

34
communes of Guinea
22
distinct values
20
co-celled communes
504
identical months (42 yrs)
01

The problem

A satellite or reanalysis product delivers radiation on a grid. Two genuinely distinct sites that fall into the same cell receive the same value, to the byte. The founding note documented this on a single pair: Kindia and Mamou, 92 km apart and 378 m in altitude, to which NASA POWER assigns strictly identical series over decades.

This work generalises the observation from a pair to a territory, and quantifies it: how much spatial information a product erases over Guinea, and where. The phenomenon is named here pixel degeneracy: a product's inability to tell apart genuinely distinct sites because they share a cell.

To be distinguished from the spatial representativeness of a ground measurement (Hakuba 2013, Schwarz 2018) and from resource interpolation: this is not a point-to-area error, it is an intrinsic property of the product, measured against a reference of real sites.

02

Hypotheses

H1

At NASA POWER's resolution (1°), a substantial share of the 34 communes share their cell and receive byte-identical series.

H2

Degeneracy decreases with cell size: strong for NASA POWER (1°), lower for ERA5 (0.25°), negligible for fine products (CAMS, SARAH-3, ERA5-Land).

H3

The erased separations are not harmless: co-celled communes differ in altitude and regime, hence in genuinely distinct resource.

H4

Where the coarse product is degenerate, fine products separate the sites but disagree on the sign of the difference: the real gap exists without being resolved by satellite alone.

03

Method & maths

For a set of N=34N = 34 sites, a product PP tiles space into cells and assigns each commune to the cell containing it. Identity property, verified rather than assumed: the product returns an identical series to all sites in one cell.

Degeneracy indicators are defined from the number of distinct values DPD_P (occupied cells):

The rate ρP=0\rho_P = 0 if the product resolves every site, and ρP1\rho_P \to 1 if it collapses them all: it is the fraction of spatial information lost. We also measure the maximal collapse MPM_P and, for each shared cell, the erased separation (geographic distance and altitude range).

Deterministic pipeline: raw data read from disk, outputs dated and traceable down to the computation.

Indicators

ρP=1DPN,δP=1N{s:ss,  cP(s)=cP(s)}.\rho_P = 1 - \frac{D_P}{N}, \qquad \delta_P = \frac{1}{N}\,\bigl\lvert\,\{\, s : \exists\, s' \neq s,\; c_P(s') = c_P(s)\,\}\,\bigr\rvert.
04

Results

On NASA POWER and the 34 communes, only 22 distinct values are returned:

A third of the spatial information is erased: 20 of 34 communes share their cell, only 14 are resolved alone. Within each shared cell, the monthly GHI series are strictly identical over 504 months (42 years), without exception.

The note's Kindia-Mamou case is in fact a triplet with Dalaba, erasing 371 m of altitude range over 97 km; the coastal cell collapses the capital and three neighbouring prefectures into a single value. A telling irony: NASA POWER knows different altitudes for these sites (289 m, 660 m) yet gives them the same radiation. Hypotheses H1H_1 and H3H_3 are established.

Key result

D=22,ρ=122340.35,δ=20340.59,M=4.D = 22, \qquad \rho = 1 - \tfrac{22}{34} \approx 0.35, \qquad \delta = \tfrac{20}{34} \approx 0.59, \qquad M = 4.
05

Open questions

  1. 01Extend the diagnosis to ERA5-Land, CAMS then SARAH-3 and show that degeneracy ρP\rho_P decreases as the cell tightens (H2H_2).
  2. 02On the pairs NASA POWER conflates, measure whether fine products agree on the sign of the inter-site difference (H4H_4).
  3. 03Define an effective resolution: the cell size below which degeneracy vanishes on this reference.
  4. 04Add ERA5 at 0.25° as an intermediate resolution point between 1° and ~5 km.