We have measured the muon flux at the underground Gran Sasso National Laboratory (3800 m w.e.) to be (3.41 \pm 0.01) \times 10-4m-2s-1 using four years of Borexino data. A modulation of this signal is observed with a period of (366\pm3) days and a relative amplitude of (1.29 \pm 0.07)%. The measured phase is (179 \pm 6) days, corresponding to a maximum on the 28th of June. Using the most complete atmospheric data models available, muon rate fluctuations are shown to be positively correlated with atmospheric temperature, with an effective coefficient {\alpha}T = 0.93 \pm 0.04. This result represents the most precise study of the muon flux modulation for this site and is in good agreement with expectations.

Borexino is a large-volume liquid scintillator detector installed in the underground halls of the Laboratori Nazionali del Gran Sasso in Italy. After several years of construction, data taking started in May 2007. The Borexino phase I ended after about three years of data taking. Borexino provided the first real time measurement of the $^{7}$Be solar neutrino interaction rate with accuracy better than 5% and confirmed the absence of its day-night asymmetry with 1.4% precision. This latter Borexino results alone rejects the LOW region of solar neutrino oscillation parameters at more than 8.5 $\sigma$ C.L. Combined with the other solar neutrino data, Borexino measurements isolate the MSW-LMA solution of neutrino oscillations without assuming CPT invariance in the neutrino sector. Borexino has also directly observed solar neutrinos in the 1.0-1.5 MeV energy range, leading to the first direct evidence of the $pep$ solar neutrino signal and the strongest constraint of the CNO solar neutrino flux up to date. Borexino provided the measurement of the solar $^{8}$B neutrino rate with 3 MeV energy threshold.

Borexino, a liquid scintillator detector at LNGS, is designed for the detection of neutrinos and antineutrinos from the Sun, supernovae, nuclear reactors, and the Earth. The feeble nature of these signals requires a strong suppression of backgrounds below a few MeV. Very low intrinsic radiogenic contamination of all detector components needs to be accompanied by the efficient identification of muons and of muon-induced backgrounds. Muons produce unstable nuclei by spallation processes along their trajectory through the detector whose decays can mimic the expected signals; for isotopes with half-lives longer than a few seconds, the dead time induced by a muon-related veto becomes unacceptably long, unless its application can be restricted to a sub-volume along the muon track. Consequently, not only the identification of muons with very high efficiency but also a precise reconstruction of their tracks is of primary importance for the physics program of the experiment. The Borexino inner detector is surrounded by an outer water-Cherenkov detector that plays a fundamental role in accomplishing this task. The detector design principles and their implementation are described. The strategies adopted to identify muons are reviewed and their efficiency is evaluated. The overall muon veto efficiency is found to be 99.992% or better. Ad-hoc track reconstruction algorithms developed are presented. Their performance is tested against muon events of known direction such as those from the CNGS neutrino beam, test tracks available from a dedicated External Muon Tracker and cosmic muons whose angular distribution reflects the local overburden profile. The achieved angular resolution is 3-5 deg and the lateral resolution is 35-50 cm, depending on the impact parameter of the crossing muon. The methods implemented to efficiently tag cosmogenic neutrons are also presented.