add negative binomial for abundance simulation