add script to aggregate data and get number of records of each taxon by year