benchmark_area_intensity function

Benchmark area intensity against regional data