Lidar terrain analysis · Middleborough, Massachusetts

150 wooded acres.
A withdrawn warehouse.
What does the land say?

In 2023 a developer proposed 673,000 square feet of warehouses on forest near the Route 44 / I-495 interchange. The proposal was withdrawn in 2024. This planning-level study uses public LiDAR to ask what terrain, mapped constraints, drainage, and conceptual grading could have contributed before a final design.

Site location in Massachusetts

Scope: public-data demonstration for planning and comparison. Mapped wetlands and terrain signatures are not field delineations.

246 acres of this study area sit inside wetland and stream setbacks — over half the land here already carries regulatory constraints.
84,437 yd³
LiDAR-based screening estimate for the conceptual 16-acre pad
789,453 yd³
apparent earthwork if graded to the treetop surface
2.0 cm
how closely my model matches the state's published DEM
49.5M
laser points processed for this site
Start here

Drag the line

Left of the line is the treetop surface — what camera-only mapping is most likely to recover over closed forest. Right of it is the bare ground, recovered from laser returns that made it through gaps in the canopy. Same flight, same week, same points. The yellow rectangle is the building pad analyzed below.

Bare-earth lidar hillshade
Canopy surface hillshade
TREETOP SURFACE (CAMERA-STYLE)
BARE EARTH (LIDAR)
USGS/MassGIS lidar, flown spring 2021 · 1.5 × 1.2 km area · 0.5 m resolution
What it costs to get this wrong

One pad, three surfaces

A 16.3-acre pad sized for the proposed warehouses, sited by an optimization pass over the constraint screens (95.9% clear of wetlands and steep ground — my first hand-drawn attempt turned out to be half swamp, which is exactly what the screens are for). It's 69% forest with canopy to 23 m. Grading it to a level 12.82 m requires:

Conceptual level-pad earthwork sensitivity by input surface
Surface used for the estimateCutFillTotal moved@ $20 / yd³
Lidar bare earth42,27342,16484,437 yd³$1.69M
Treetop surface (canopy counted as ground)775,62613,827789,453 yd³$15.79M
10 m national DEM (all you had before 2021)52,17734,46186,638 yd³$1.73M
Canopy error hiding in the treetop surface705,016 yd³≈ $14.10M

The old 10-meter data is a sneakier failure than the treetop surface. Its total looks fine — but it calls the site 17,700 yd³ cut-heavy when the site actually balances, mislabels cut vs fill on 7% of the pad, and is locally off by up to 2.9 m. A grading plan based on it would include a haul-off that doesn't exist. (Also a data-diligence note: the "current" national DEM has quietly been rebuilt from the 2021 lidar — to compare against what was actually available before, you have to pull the archived 2019 version — worth checking where any comparison surface actually came from.)

A uniform ±10 cm vertical shift changes the modeled volume by ±8,628 yd³ ≈ ±$173k. This is a data-accuracy sensitivity—not a full construction contingency. Design, soils, water, rock, material handling, and field conditions can create larger changes.

Cut and fill map of the pad
Where the dirt moves: cut in red, fill in blue, against the balanced design grade.
Site analysis from the same data

More answers from the same flight

The earthwork number is the headline, but the same dataset answers a bunch of other early-stage questions:

Wetland-risk discrepancy map
Where does the map disagree with the terrain? Cells whose wetness signature matches the mapped wetlands, but that the state map calls dry: 112 acres here — including 10 of the pad's 16.3. Sandy ground drains, so this is a screening flag, not a delineation; it shows where a wetland scientist should look first.
Buildable land screen
Where does the terrain screen pass? Of 441 study-area acres, 175 remain after excluding the stated slopes and mapped wetland/stream buffers. Zoning, access, utilities, habitat, soils, ownership, and field delineation are not included.
Access road grade profile
Can you get a road in? Yes — and that's worth knowing too. The route from Route 44 to the pad's far corner runs 333 m, climbs 3.2 m, and never exceeds 5.1% grade. That route has no terrain-grade flag at this screening level; final access geometry and cost still require design.
Canopy height model
What does clearing cost? Tree heights for every cell. The pad carries 11.2 acres of forest — roughly $34–67k to clear at typical rates. Stumpage revenue can offset part of that — species- and market-dependent (UMass/DCR Southern New England Stumpage Price Survey); not netted here.
Stormwater flow accumulation
Where does water go? Runoff concentration paths on the real terrain. On a parcel this wet, drainage will be a big part of permitting, and this is a first look at where the water wants to go.
Viewshed before and after clearing
Who will see it? From the pad today: nothing — 0.0% of the surrounding area is visible through the forest. After clearing: 18.2%, reaching the Route 44 frontage and the homes south of the parcel. A retained tree buffer on those edges would soften most of it — easier to plan for now than to argue about later.
Why this parcel

The analysis that wasn't in the hearing room

When the warehouse proposal was argued in 2023–24, most of the technical information in the room came from studies the developer commissioned. Everything on this page comes from free public data and could have been on the table for either side: how much land passes stated terrain screens, where water concentrates, what the neighbors might see, and how conceptual grading compares. The parcel is still forest, still zoned, still on the market in practice — which makes it a live case study rather than a postmortem.

Aerial imagery of the parcel today
The same ground today: still forest, with logistics buildings and trailer yards on both sides. Basemap: Esri World Imagery; provider credits are listed in methodology and sources.
The bottom line

What this means for the deal

Planning-level terrain subtotal for the conceptual pad
Screening itemEstimate
Earthwork, graded to balance$1.69M ± $173k
Clearing, 11.2 forested acres$34–67k
Access terrain flagscreened route is 333 m and peaks at 5.1%; final geometry not designed
Terrain-screening subtotal: modeled earthwork + clearing only≈ $1.74M

Not included: final grading design, topsoil, shrink/swell, unsuitable soils, rock, haul/disposal, retaining walls, utilities, stormwater construction, erosion control, mobilization, permitting, survey, or engineering.

The modeled terrain subtotal is lower than in the Hopkinton case; the open questions here are regulatory and field-based. Diligence should cover:

  1. Wetland delineation before the offer is priced — the pad clears the state's wetland map, but 10 of its 16.3 acres share the terrain-wetness signature of the mapped wetlands. If the field delineation moves the line, the pad moves — worth making the offer contingent on the delineation.
  2. Check the vintage of any elevation data in the pro forma — the pre-2021 public DEM shows a 17,700 yd³ haul-off that isn't really there, and older underwriting may still be carrying it.
  3. Vernal pools — depression screen found no candidates in or near the pad; 5 state-mapped potential pools sit elsewhere in the study area. "Screened, clear" at the pad, field confirmation at permitting.
Methodology

The numbers are checkable

  1. Data — USGS 3DEP lidar (2021 Central-Eastern Massachusetts, block 2, published accuracy 10 cm RMSE). 49.5 million points for this site. Wetlands and streams: MassGIS DEP layers with 100 ft / 200 ft buffers.
  2. Processing — same PDAL pipeline as Project 01: ground/vegetation separation, 0.5 m ground model, treetop model, canopy heights. Same scripts, different coordinates — which is most of what I wanted to prove.
  3. Processing consistency — my ground model compared with the state-published DEM derived from the same flight: 1 mm average difference and 2.0 cm RMSE over open ground. That checks processing agreement, not independent field accuracy. The pre-2021 baseline was pulled from the archived 2019 national DEM after finding the "current" one had been rebuilt from this same lidar.
  4. Volumes — cell-by-cell raster math, cross-checked in QGIS: agreement to five significant figures. Cost basis is an illustrative $20/yd³ screening assumption, with a $10–$40 sensitivity range; a real project would require current contractor or estimator input.

Public LiDAR includes published accuracy specifications, but these case-study values remain planning-level. Survey-grade or construction-facing work requires licensed control and project-specific professional oversight. Screens run for every project: terrain-wetness discrepancy, closed-depression / vernal-pool candidates, surface roughness. Also available: intensity analysis, contours to CAD, per-tree inventory. Scripts and pipeline files available on request. See the shared method, sources, assumptions, and limitations.