Geotools indexes shapefiles in memory when they are being written. So during our analysis, when we have large datasets in memory, and try to write them to disk, the geotools indexer adds lots of extra usage to the heap, and we run out of heap space.

The temporary fix is to turn off indexing while writing analysis results, and have them be indexed when they are loaded.