CyberFaCES

This module walks participants through how to calculate several heat metrics and how to calculate changes in labor productivity from changes in these heat metrics. You will run four notebooks.

• Notebook Summaries

• **01_WBGT_fundamentals.ipynb **— This is a teaching notebook that walks through the physics of Wet Bulb Globe Temperature (WBGT) calculation from the ground up, following Kong & Huber (2022). It derives relative humidity from ERA5 temperature/dew point, computes solar geometry, and shows the source code for the natural wet-bulb and black-globe energy-balance solvers (using wbgt_physics.py). It then runs a single synthetic ERA5-like case end-to-end to compute explicit WBGT, compares it against simplified WBGT and ESI approximations, and demonstrates ISO-7243 labor-productivity limits with a sensitivity analysis.

• 02_ERA5_Africa_WBGT.ipynb — This notebook scales the WBGT calculation from a single point to a gridded domain, reading one week of hourly ERA5 reanalysis data over Africa (with synthetic, local-file, or direct AWS S3 download modes). It standardizes variable names/units, derives humidity, wind, and solar terms across the grid, then computes both outdoor and indoor WBGT fields using the same physics as Notebook 1. Results are visualized as time series and daytime-mean maps, and a compact baseline NetCDF product (wbgt_outputs/era5_africa_wbgt_baseline_week.nc) is saved for downstream notebooks, with an optional step to export per-country crops.

• 03_country_labor_productivity.ipynb — This notebook takes the gridded WBGT baseline from Notebook 2, masks it to a single country (Ghana by default) using Natural Earth boundaries via regionmask, and computes ISO-style labor productivity at four metabolic work rates. It carefully applies the nonlinear productivity function at each daytime timestep before averaging, then produces area-weighted country summaries comparing baseline productivity to a teaching-only "+2°C" future scenario. Outputs include productivity bar charts, a spatial loss map, and a country-level CSV/NetCDF export to country_outputs/, along with notes on how to extend the approach toward the full Saeed et al. (2022) methodology.

• 04_country_capacity_loss_for_GTAP.ipynb — This notebook implements the final step of the Saeed et al. (2022) heat-stress-to-labor-capacity chain, computing the grid-cell labor-capacity loss quantity (z^{i,c}{r,g}) that feeds Equation 4 of their GTAP economic shock. It combines baseline outdoor/indoor WBGT with a bias-corrected "+3°C" CMIP6-derived anomaly, applies the ISO 7243 capacity function per work intensity and indoor/outdoor condition, and computes the relative capacity loss between baseline and future climates. The notebook optionally crops results to one country and saves a gridded hand-off NetCDF product, documenting the remaining steps (BLS work-share weighting, GDP adjustment, spatial aggregation) needed to reach the actual GTAP shock.