Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

We use the Szekeres inhomogeneous relativistic models in order to fit supernova combined data sets. We show that with a choice of the spatial curvature function that is guided by current observations, the models fit the supernova data almost as well as the LCDM model without requiring a dark energy component. The Szekeres models were originally derived as an exact solution to Einstein's equations with a general metric that has no symmetries and are regarded as good candidates to model the true lumpy universe that we observe. The null geodesics in these models are not radial. The best fit model found is also consistent with the requirement of spatial flatness at CMB scales. The first results presented here seem to encourage further investigations of apparent acceleration using various inhomogeneous models and other constraints from CMB and large structure need to be explored next.